2022/02/draw-china-maps/

[utterances-bot]

R 语言画中国地图 - Xiangyun Huang | 黄湘云

https://xiangyun.rbind.io/2022/02/draw-china-maps/

[chuxinyuan]

http://datav.aliyun.com/tools/atlas/#&lat=30.334491405102664&lng=106.72420352760892&zoom=3.5
这个数据应该是完美的。

[ShiganLuo]

试了一下，确实可以

[XiangyunHuang]

@chuxinyuan 应该是的，在高德地图一节，特意提到了这个地图数据源。

[hongyuanjia]

正愁去哪里下载到绝对正确的中国地图，很有帮助，谢谢！

[ShiganLuo]

太感谢了，帮了我很大忙

[XiangyunHuang]

不客气，我没白写😂

[1junyuan]

厉害的

[qtl0719]

请问怎么得到完整的市县级的地图？笔者描述的使用raster 包下载网站 GADM 上的国家行政边界信息，然而现在无论使用getdata函数或者在GADM网站上下载都无法获取地图了，且缺乏香港澳门台湾的地图，万望哪位老师解答一下，谢谢！

[XiangyunHuang]

你把整篇文章都看完的时候就知道从哪下载了。

[tcgriffith]

最近有个问题是想把一堆地名标在地图上，顺便加一些meta信息，跟着帖子玩了一下，发现leafletCN被cran拿下了，最后还是高德地图+leaflet好使 。

另外高德地图的api识别坐标也挺好用的,问了ds直接拿来用了

library(jsonlite)

# 你的API Key
api_key <- "5a469608137235718fecf7cefcb0c1b3"

# 地址列表
addresses <- c("北京市天安门", "北京师范大学", "北京市通州区梨园镇")


# 初始化结果向量
results <- data.frame(address = character(), lng = numeric(), lat = numeric(), stringsAsFactors = FALSE)

# 循环获取每个地址的坐标
# for (address in addresses) {

for (address in addresses) {
  # 构造请求URL
  url <- paste0("https://restapi.amap.com/v3/geocode/geo?address=", URLencode(address), "&key=", api_key)
  
  # 发送GET请求
  response <- GET(url)
  
  # 解析JSON响应
  data <- fromJSON(content(response, "text"))
  
  # 提取经纬度
  if (data$status == "1" && length(data$geocodes) > 0) {
    location <- strsplit(data$geocodes$location, ",")[[1]]
    lng <- as.numeric(location[1])
    lat <- as.numeric(location[2])
    
    # 将结果添加到数据框
    results <- rbind(results, data.frame(address = address, lng = lng, lat = lat, stringsAsFactors = FALSE))
  } else {
    print(paste("无法获取地址:", address))
  }
  
  # 避免频繁请求，适当延时
  Sys.sleep(1)
}


# 输出结果
print(results)

[XiangyunHuang]

谢谢告知 leafletCN 的情况。我在家乡邵阳的那个例子里就直接使用了高德地图的，没有使用 leafletCN。