2000坐标系转经纬度工具_用编程赋能工作系列——百度VS高德经纬度互转

上一篇刚学会了如何获取地址所在经纬度以及经纬度对应地址，于是信心满满的准备在老板面前露一手，准备花点心思做个漂亮的可视化地图放在报告亮眼的位置。

当地图跑出来的那一刻顿时傻眼了，卧槽这些点定位的位置明显不对呀，酒店直定位到湖里了，这让老板看到你说后果严重不。

其实是因为百度地图的经纬度和高德地图经纬度所使用的坐标系编码不同，所以如果你制图的软件或平台是基于高德地图服务的话，那么就需要使用对应的高德标准的经纬度来进行绘图，同理如果你用的服务是基于百度地图的，必须使用百度标准的经纬度坐标，如果拿到的是高德坐标则必须经过算法转为百度标准。

今天这篇内容就给大家分享如何使用R和Python对百度坐标系和高德坐标系的经纬度进行互转，解决地图绘制最后一道坎儿。

内容同样是两个模块四小节，使用R语言和Python分别进行百度经纬度转高德经纬度以及高德标准转百度，这样我们在经纬度获取和转化处理上的技能掌握就比较系统了，不会再受制于工具和平台服务的标准差异而苦恼。

一、R语言方案

1）百度经纬度转腾讯&高德

library("leaflet")
library('baidumap')
library('ggplot2')
library('ggmap')


dt1 <- read.table(pipe("pbpaste"), sep="t", header=T,stringsAsFactors =FALSE)
#百度转腾讯&高德
bMapTransQQMap <- function (lng,lat) {
  x_pi = 3.14159265358979324 * 3000.0 /180.0
  x = lng - 0.0065
  y = lat - 0.006
  z = sqrt(x^2 + y^2) - 0.00002 * sin(y * x_pi)
  theta = atan2(y,x) - 0.000003 * cos(x * x_pi)
  lngs = z * cos(theta)
  lats = z * sin(theta)
  return(data.frame(lng,lat,lngs,lats)) 
}

result <- bMapTransQQMap(dt1$lon,dt1$lat)

使用高德地图服务来呈现百度坐标系下的经纬度，可以看到偏移量非常明显。

leaflet(result) %>