python3 将RTK中的点导出后整理为unit为degree的点
"""
Created on Fri Sep 18 23:22:15 2020
@author: Yiting_BAI
"""
import numpy as np
import pandas as pd
import os
input_path = 'E:\\AIR_CAS\\怀来遥感综合试验站实验\\'
output_fn = os. path. join( input_path, 'correctjingweidu.csv' )
x1= pd. read_excel( os. path. join( input_path, 'jingweidu.xls' ) )
contents = pd. DataFrame( x1) . to_numpy( )
lat = [ ]
lon = [ ]
for i in range ( len ( contents[ : , 1 ] ) ) :
temp = contents[ : , 1 ] [ i] . split( '.' )
d = temp[ 0 ]
m = temp[ 1 ]
s = temp[ 2 ]
ls = temp[ 3 ]
du = float ( d) + float ( m) / 60 . + float ( s+ '.' + ls) / 3600 .
lat = np. append( lat, du)
temp1 = contents[ : , 2 ] [ i] . split( '.' )
d = temp1[ 0 ]
m = temp1[ 1 ]
s = temp1[ 2 ]
ls = temp1[ 3 ]
dudu = float ( d) + float ( m) / 60 . + float ( s+ '.' + ls) / 3600 .
lon = np. append( lon, dudu)
ID = contents[ : , 0 ]
elevation = contents[ : , 3 ]
ID_column = pd. Series( ID, name= 'ID' )
Lon_column = pd. Series( lon, name= 'Longitude' )
Lat_column = pd. Series( lat, name= 'Latitude' )
elevation_column = pd. Series( elevation, name= 'elevation' )
predictions = pd. concat( [ ID_column, Lon_column, Lat_column, elevation_column] , axis= 1 )
save = pd. DataFrame( { 'ID' : ID, 'Longitude' : lon, 'Latitude' : lat, 'elevation' : elevation} )
save. to_csv( output_fn, index= False , sep= ',' )