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.0x = lng - 0.0065y = lat - 0.006z = 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) %>% setView(116.2938,40.00939, zoom = 13) %>%addTiles('http://webrd02.is.autonavi.com/appmaptile?lang=zh_cn&size=1&scale=1&style=8&x={x}&y={y}&z={z}',tileOptions(tileSize=256, minZoom=9, maxZoom=17),group="高德地图") %>%addCircleMarkers(lng = ~  lng,lat = ~lat,stroke = FALSE,fillOpacity = 1,radius =8)

leaflet(result) %>% setView(116.2938,40.00939, zoom = 13) %>%addTiles('http://webrd02.is.autonavi.com/appmaptile?lang=zh_cn&size=1&scale=1&style=8&x={x}&y={y}&z={z}',tileOptions(tileSize=256, minZoom=9, maxZoom=17),group="高德地图") %>%addCircleMarkers(lng = ~  lngs,lat = ~lats,stroke = FALSE,fillOpacity = 1,radius =8)

# 腾讯&高德地图转百度
qqMapTransBMap <- function(lngs,lats){x_pi = 3.14159265358979324 * 3000.0 / 180.0x = lngs y = lats z = sqrt(x^2 + y^2) + 0.00002 * sin(y * x_pi)theta = atan2(y,x) + 0.000003 * cos(x * x_pi)lng_b = z * cos(theta) + 0.0065lat_b = z * sin(theta) + 0.006  return(data.frame(lngs,lats,lng_b,lat_b))
}
result_b <- qqMapTransBMap(result$lngs,result$lats)

options(baidumap.key = '***************')
# 这里启用的baidumap服务，后台调用百度地图api，
# 需要使用百度api开发秘钥，需自己申请
bjMap <- getBaiduMap(location = c(116.284028,40.001732), width = 800, height =800, zoom = 16, scale = 2 , messaging = TRUE)
ggmap(bjMap) + geom_point(data = result_b, aes(x =  lngs, y = lats),shape = 21 , col = 'white',fill = 'red',size = 5) +theme_nothing()

ggmap(bjMap) + geom_point(data = result_b, aes(x =  lng_b, y =  lat_b),shape = 21 , col = 'white',fill = 'red',size = 5) +theme_nothing()

import pandas as pd
import numpy as np
import folium
from folium import plugins
df=pd.read_clipboard()
def bMapTransQQMap(lng,lat):x_pi = 3.14159265358979324 * 3000.0 /180.0x = np.array(lng) - 0.0065y = np.array(lat) - 0.006z = np.sqrt(np.power(x,2) + np.power(y,2)) - 0.00002 * np.sin(y * x_pi)theta = np.arctan2(y,x) - 0.000003 * np.cos(x * x_pi)lngs = z * np.cos(theta)lats = z * np.sin(theta)return(pd.DataFrame({'lng':lng,'lat':lat,'lngs':lngs,'lats':lats}))
result = bMapTransQQMap(df['lon'],df['lat'])

lng = np.array(result["lng"],dtype=float)
lat = np.array(result["lat"],dtype=float)
data1 = [(lat[i],lng[i]) for i in range(len(result))]map_osm = folium.Map(location=[39.996710,116.281012],zoom_start=115,tiles='http://webrd02.is.autonavi.com/appmaptile?lang=zh_cn&size=1&scale=1&style=8&x={x}&y={y}&z={z}',attr="&copy; <a href="http://ditu.amap.com/">高德地图</a>")
marker_cluster = plugins.MarkerCluster().add_to(map_osm)
for i in data1:folium.Marker(i).add_to(marker_cluster)
display(map_osm)

lng = np.array(result["lngs"],dtype=float)
lat = np.array(result["lats"],dtype=float)
data2 = [(lat[i],lng[i]) for i in range(len(result))]map_osm = folium.Map(location=[39.996710,116.281012],zoom_start=115,tiles='http://webrd02.is.autonavi.com/appmaptile?lang=zh_cn&size=1&scale=1&style=8&x={x}&y={y}&z={z}',attr="&copy; <a href='http://ditu.amap.com/'>高德地图</a>")
marker_cluster = plugins.MarkerCluster().add_to(map_osm)
for i in data2:folium.Marker(i).add_to(marker_cluster)
display(map_osm)

#腾讯&高德地图转百度
def qqMapTransBMap(lngs,lats):x_pi = 3.14159265358979324 * 3000.0 / 180.0x = np.array(lngs)y = np.array(lats )z = np.sqrt(np.power(x,2) + np.power(y,2)) + 0.00002 * np.sin(y * x_pi)theta = np.arctan2(y,x) + 0.000003 * np.cos(x * x_pi)lng_b = z * np.cos(theta) + 0.0065lat_b = z * np.sin(theta) + 0.006  return(pd.DataFrame({'lngs':lngs,'lats':lats,'lng_b':lng_b,'lat_b':lat_b}))result2 = qqMapTransBMap(result['lngs'],result['lats'])