Bidirectional UTM_to_WGS84

 1 import utm
 2 import numpy as np
 3 from shapely.wkt import dumps, loads
 4 outfile = "D:/研究/课题资料/申请后提交/验收-准备/coding4-yanshou/data_vec_output/inputTIF.out"
 5 # outfile: LINESTRING (112.87052724867623 32.780393689498865, 112.87060585126358 32.78038236059359, 112.87071589488588 32.7803445975656, 112.87077653116756 32.780337044958074)
 6 
 7 with open(outfile, 'r') as output:
 8     for line_wkt_str in output.readlines():
 9         lr = loads(line_wkt_str)
10         lon_arr, lat_arr = np.array(lr.xy[0]), np.array(lr.xy[1])
11         print(lon_arr), print(lat_arr)
12 
13         # latlon_2_utm
14         easting, northing, zone_number, zone_letter = utm.from_latlon(latitude=lat_arr, longitude=lon_arr)
15         print(easting, northing, zone_number, zone_letter)
16         # utm_2_latlon
17         latitude, longitude = utm.to_latlon(easting, northing, zone_number, zone_letter)
18         print(longitude, latitude)

UTM.zone_number(UTM zone numbers), UTM.zone_letter(UTM zone designators) can be found in url: http://www.jaworski.ca/utmzones.htm

posted @ 2021-09-27 20:24  土博姜山山  阅读(48)  评论(0编辑  收藏  举报