处理步骤

  • 1. 目标:
    • (1)获取最新中国县级行政区划矢量数据
    • (2)学习如何使用R语言将.json数据转为shp文件
  • 2. 步骤:
    • 将json数据处理为shapefile
  • 3. 结果
    • 数据分享

数据来源:高德开放平台

1. 目标:

(1)获取最新中国县级行政区划矢量数据

(2)学习如何使用R语言将.json数据转为shp文件

2. 步骤:

将json数据处理为shapefile

library(pacman)
p_load(sf,tmap,maptools,rgdal)
china<-st_read("https://geo.datav.aliyun.com/areas_v2/bound/100000_full.json")
class(china)
# [1] "sp"
chinaqtm()qtm(china, fill = "name")+tm_text("name")crs(china)
# CRS arguments: +proj=longlat +datum=WGS84 +no_defs
china$name
# [1] "北京市"           "天津市"           "河北省"           "山西省"           "内蒙古自治区"
# [6] "辽宁省"           "吉林省"           "黑龙江省"         "上海市"           "江苏省"
# [11] "浙江省"           "安徽省"           "福建省"           "江西省"           "山东省"
# [16] "河南省"           "湖北省"           "湖南省"           "广东省"           "广西壮族自治区"
# [21] "海南省"           "重庆市"           "四川省"           "贵州省"           "云南省"
# [26] "西藏自治区"       "陕西省"           "甘肃省"           "青海省"           "宁夏回族自治区"
# [31] "新疆维吾尔自治区" "台湾省"           "香港特别行政区"   "澳门特别行政区"   ""


导出为shapefile

st_write(china, "G:/Rdata/China/GeoAltas/China.shp", layer_options = "ENCODING=UTF-8",delete_layer = TRUE)
**province <- paste0("https://geo.datav.aliyun.com/areas_v2/bound/",china$adcode,"_full.json")
province1 <- province[1:34]
province1 count=0
for (i in province1){x <- st_read(i)count <- count + 1st_write(x,dsn = "G:/Rdata/China/GeoAltas/1",layer = china$name[count],layer_options = "ENCODING=UTF-8",delete_layer = TRUE, driver = "ESRI Shapefile")
}county <- list.files("G:/Rdata/China/GeoAltas/1/省市",full.names = TRUE,pattern = ".shp$")library(plyr)
data <- NULL
for(aFile in county){data <- rbind.fill(data, st_read(aFile))
}
data <- st_as_sf(data)
class(data)
# [1] "sf"         "data.frame"
qtm(data,fill = "name")
data#https://stackoverflow.com/questions/64569432/st-union-sf-polygons-with-data-table
data=data[1:nrow(data),]
st_write(data,"G:/Rdata/China/GeoAltas/city.shp",layer_options = "ENCODING=UTF-8",delete_layer = TRUE)library("dplyr")
data <- data%>%dplyr::filter(chldrnN >0 )
qtm(data)city_county <- paste0("https://geo.datav.aliyun.com/areas_v2/bound/",data$adcode,"_full.json")
head(city_county)
# [1] "https://geo.datav.aliyun.com/areas_v2/bound/340100_full.json"
# [2] "https://geo.datav.aliyun.com/areas_v2/bound/340200_full.json"
# [3] "https://geo.datav.aliyun.com/areas_v2/bound/340300_full.json"
# [4] "https://geo.datav.aliyun.com/areas_v2/bound/340400_full.json"
# [5] "https://geo.datav.aliyun.com/areas_v2/bound/340500_full.json"
# [6] "https://geo.datav.aliyun.com/areas_v2/bound/340600_full.json"count1 = 0
for(i in city_county){x <- st_read(i)count1 <- count1 + 1st_write(x,dsn = "G:/Rdata/China/GeoAltas/1/市县",layer = data$name[count1],layer_options = "ENCODING=UTF-8", delete_layer = TRUE,driver = "ESRI Shapefile")
}**

3. 结果

数据分享

https://pan.baidu.com/s/1yhLJodusCwx9GTK858Uxuw
提取码:1111

欢迎关注个人公众号GeoSuper

如何获取最新中国县级行政区划矢量数据相关推荐

  1. 获取最新2020全国行政区划及经纬度

    源码:https://github.com/tyronczt/java-learn/tree/master/Tools/China-Xzqh 最后效果 2020年中华人民共和国县以上行政区划代码 ht ...

  2. python程序获取最新的行政区划名称代码

    一.实现目标           最近由于项目需要,需要获取最新的过去全国县以上行政区划的名称和代码.网上虽然有一些资料,但是不是需要积分就是需要会员,而且担心这些资料不是最新的.因此,想着使用程序从 ...

  3. 最新中国省市区县乡镇 5 级行政区划代码Json数据

    之前我爬了datav的地图选择器中全国省市区的geoJSON数据 为了获取最新省市区的行政区划代码用于级联组件,又爬了国家统计局的内容 点此查看 index_2019_level_3.json:全国省 ...

  4. 实时最新中国省市区县geoJSON格式地图行政边界数据Echarts地图数据(可精确到街道级)

    geojson 数据下载地址:https://hxkj.vip/demo/echartsMap/ 可下载的数据包含省级geojson行政边界数据.市级geojson行政边界数据.区/县级geojson ...

  5. 从百度地图批量获取中国县级以上行政区划边界坐标

    从百度地图批量获取中国县级以上行政区划边界坐标 不知道从什么时候开始,在Google地图或是百度地图中用城市或者省名称搜索的时候,你会发现地图上会勾勒出搜索的行政区的轮廓来,这个功能非常人性化. 比如 ...

  6. 零编程基础获取最新、最权威的县以上行政区划代码信息

    一.行政区划代码 行政区划代码是国家对能够统治的行施管辖区域进行分级分层进行管辖用信息化手段编制的对各层级行政区划编制的替代数码,一般执行两项国家标准: 一是<中华人民共和国行政区划代码> ...

  7. 获取最新的中国IP的脚本,给ROS可以使用的脚本

    获取最新的中国IP的脚本,给ROS可以使用的脚本 wget http://ftp.apnic.net/apnic/stats/apnic/delegated-apnic-latest echo &qu ...

  8. 最新中国省市区县geoJSON格式地图数据Echarts地图数据

    最新全国省市区县geoJSON格式的地图数据,可直接用于echarts地图展示: https://github.com/lyhmyd1211/GeoMapData_CN 全国地图: china.jso ...

  9. 最新中国上市公司市值500强,都分布在哪里?

    本文禁止转载到任何地方! 1. 背景介绍 文中数据,关注公众号:数据分析与统计学之美,后台回复:500,领取数据. 前不久,朋友发给我两张图片,图片中是表格形式的数据,包含了最新中国上市公司市值500 ...

  10. QGIS|三步拥有全国县级行政区划

    "天地图"(http://www.tianditu.gov.cn)是我国自然资源部直属的国家基础地理信息中心建设的国家地理信息公共服务平台,以多种形式面向社会和公众提供地图与地理信 ...

最新文章

  1. Windows上通过VLC播放器搭建rtsp流媒体测试地址操作步骤
  2. python 装饰器 参数-[Python]写个带参数的装饰器
  3. Boost:std ::bind与Boost的_1绑定的测试程序
  4. Spring factoryBeanInstanceCache在哪里进行put?
  5. python实验四_python实验四
  6. Build String(CF-237E)
  7. PTA17、求指定层的元素个数 (10 分)
  8. 交互设计[小插曲]--网站UI配色
  9. Effective_STL 学习笔记(三) 使容器里对象的拷贝操作轻量而正确
  10. 一道笔试题的解法和联想
  11. 谷歌招聘新职员的五大标准
  12. keil5安装_如何安装Keil5
  13. Office基础操作:Word插入visio图片显示不全
  14. CSDN有奖任务答案
  15. 电脑右下角音量键打不开,调不了音量大小的解决方法
  16. java基础回顾之Map中 TreeMap排序原理-二叉树
  17. [Power--IC]电源管理IC-STNS01
  18. Pr:旧版标题设计器
  19. linux实验报告ALU,《linux内核分析》第一次课 实验作业
  20. JavaScript 千分位 货币格式化

热门文章

  1. 【板栗糖GIS】如何将3dmax数据导入到超图软件中
  2. NUC980 DIY项目大挑战 - EtherCAT实现
  3. [趣谈网络协议学习] 17 P2P协议:我下小电影, 99%急死你
  4. 铁路售票系统_高铁铁路运营客票乘务,自动售检票务实训的诞生背景
  5. 三菱plc编程软件gx-developer在win10中无法使用及工程无法建立解决方法
  6. Testbench编写详解
  7. 10年专业导师整理的单片机控制电动机正反转设计类毕业论文文献
  8. 最新卡巴斯基密码管理器:注重便利性和强化密码控制
  9. Delphi2007企业版破解
  10. 返回顶部php代码,页面按需返回顶部代码及注释说明