坐标转换,EPSG:4326转换成高德坐标教程
这里先给大家介绍几个坐标系:
1.wgs84:国际坐标系,为一种大地坐标系,也是目前广泛使用的gps全球卫星定位系统使用的坐标系。
2.gcj02:火星坐标系,是由中国国家测绘局制订的地理信息系统的坐标系统。由wgs84坐标系经加密后的坐标系。
工作中有时需要将shape数据转换成高德坐标数据,由于目前数据大部分都是wgs84国际坐标系,因此需要通过使用高德地图提供的api将84坐标转换成gcj02即火星坐标系。
第一步:准备数据
方法1:
从shape图上随机提取一些点,根据shape图形大小提取,如果图形较大,建议多提取一些。一般在拐点和图形的四周均匀选择。
点提取完成后获取坐标点的经纬度值。可以使用arcgis的计算工具获取x,y值。画红框的地方需要注意,要选择十进制表示的经纬度。
坐标提取完成后新建一个excel表,将x,y值分别复制黏贴在excel表中的前两列。
方法2:
将shape数据存入postgresql,并使用postgis提供的函数随机从shape上提取点,并取出其x,y值。此方法是随机选择,因此选择的点未必均匀分布在图形四周和拐点处。
1 select st_x(geom), st_y(geom) from( 2 select 3 (st_dump (st_generatepoints (nanjing.geom, 20))).geom as geom, 4 md5((random()*random())::text) as id, 5 random()*1000 as val 6 from nanjing where name = '浦口区') 7 k1
第二步:执行下面的python代码调取高德api进行坐标转换
1 import xlrd 2 from xlutils.copy import copy 3 from urllib import request 4 import json 5 6 7 class togd(): 8 def __init__(self, key, coordsys='gps', output='json'): 9 self.key = key # 高德应用的key 10 self.coordsys = coordsys # 原坐标系,默认选择gps,此处实际对应4326的坐标系 11 self.output = output # 设置高德api数据返回类型,可选json和xml 12 self.file_path = input('请输入你的文件路径') # 获取文件路径,即上面新建的excel文件名称 13 14 def split_li(self, locations_li): 15 # 拆分坐标列表 16 li = list() 17 for location in locations_li: 18 location_li = location.split(',') 19 # print(location_li) 20 li.append(location_li) 21 # print(li) 22 return li 23 24 def locations(self): 25 # 读取文件中的数据并返回locations字符串及原坐标列表 26 f = xlrd.open_workbook(self.file_path) 27 index = f.sheet_names()[0] 28 sheet = f.sheet_by_name(index) 29 nrows = sheet.nrows 30 data_str = '' 31 old_locations_li = list() 32 for i in range(nrows): 33 if i == (nrows - 1): 34 data_li = str(sheet.row_values(i)[0]) + ',' + str(sheet.row_values(i)[1]) 35 else: 36 data_li = str(sheet.row_values(i)[0]) + ',' + str(sheet.row_values(i)[1]) + '|' 37 data_str += data_li 38 old_locations_li.append(sheet.row_values(i)) 39 return data_str, old_locations_li 40 41 def make_response(self, locations): 42 # 构造请求url,并获取响应数据,返回拆分后的坐标列表 43 url = 'https://restapi.amap.com/v3/assistant/coordinate/convert?locations=%s&coordsys=%s&output=%s&key=%s' % \ 44 (locations, self.coordsys, self.output, self.key) 45 # 返回值为json 46 resp_json = request.urlopen(url) 47 # json转换成字典 48 resp_dict = json.loads(resp_json.read().decode()) 49 # 提取转换为高德坐标后的坐标 50 locations_str = resp_dict.get('locations') 51 # 字符串分割 52 locations_li = locations_str.split(';') 53 new_locations_li = self.split_li(locations_li) 54 return new_locations_li 55 56 def save(self, old_locations_li, new_locations_li): 57 # 保存成excel或txt文件 58 # new_file = copy(f) 59 # new_sheet = new_file.get_sheet(0) 60 # for m in range(len(locations_li)): 61 # for n in range(len(locations_li[m])): 62 # new_sheet.write(m, n + 2, locations_li[m][n]) 63 # file_path = self.file_path + '2' 64 # new_file.save(file_path) 65 file_path = 'result.txt' 66 with open(file_path, 'w', encoding='utf-8')as f: 67 for m in range(len(new_locations_li)): 68 # for n in range(len(new_locations_li[m])): 69 f.write(str(old_locations_li[m][0])) 70 f.write(' ') 71 f.write(str(old_locations_li[m][1])) 72 f.write(' ') 73 f.write(str(new_locations_li[m][0])) 74 f.write(' ') 75 f.write(str(new_locations_li[m][1])) 76 f.write('\n') 77 78 def run(self): 79 data_str, old_locations_li = self.locations() 80 new_locations_li = self.make_response(data_str) 81 self.save(old_locations_li, new_locations_li) 82 83 84 if __name__ == '__main__': 85 test = togd('154c586add07ef456b90b079935f47a4') # 输入高德应用的key,实例化对象 86 test.run()
代码由我自己编写,水平不高还请理解。
执行上面的代码需要申请高德应用的key,相关申请方法可查看https://lbs.amap.com/dev/key
代码执行完成后会在本地生成一个result.txt的文件,下面需要用这个文件做空间校正。
第三步:空间校正
在arcgis中对需要做坐标转换的数据开启编辑功能,打开空间校正工具栏,输入连接文件。连接文件即第二步得出的result.txt文件。
可以点击view link table查看误差值,如果误差值太大,需要重新做连接文件。
选择需要校正的数据。
点击校正。
校正完成后即可得到与高德地图匹配的shape数据。
上述方法只能用于小范围测试使用,未做精度评估。如需大规模商业化使用,还需再想别的方法。
上一篇: Redmi K20系列4个月全球销量突破300万台:神秘旗舰即将发布
下一篇: 这哥们是穿越了吧