天地图“带审图号”的行政区划数据

灭火器

看到好几个公众号有发带“2024 年最新审图号”的行政区划数据，是 天地图行政区划可视化 页面提供的 GeoJSON 文件，注册天地图账号后即可下载。我看各公众号提供的都是 shapefile 版的而不是 GeoJSON 版的，因此自己打包一份发在论坛上。



这套数据的特点是：
- 含省、市、县三个级别的数据（国界可以手动用省界合成）
- 由 MultiPolygon、LineString 和 MultiLineString 组成，前者对应行政区划的面元，后两者对应九段线、未定国界线、港澳特区线、海上省界等装饰线，有助于让画出来的地图更接近标准地图网的那种。
- 跟 GCJ-02 的数据对比了一下，坐标应该是 WGS84 的。
- 分辨率略低，没几个字段

地理上的正确性我就不太懂了，可以参考这个公众号链接：首次！！带审图号的"行政区划"矢量数据(GeoJSON)

注意，有审图号的是原链接里那张网络墨卡托投影的地图，而用来画图的数据是不存在审图号的概念的。只不过我们用这套官方的数据来画图能更有底气。

下面附带 Cartopy 读取 GeoJSON 画图的例子

1. import json
   
2. import shapely.geometry as sgeom
   
3. import matplotlib.pyplot as plt
   
4. import cartopy.crs as ccrs
   
5. 
   
6. with open('中国_省.geojson') as f:
   
7.     features = json.load(f)['features']
   
8. geoms = [sgeom.shape(ft['geometry']) for ft in features]
   
9. 
   
10. crs = ccrs.PlateCarree()
    
11. plt.figure(figsize=(10, 10))
    
12. ax = plt.axes(projection=crs)
    
13. ax.set_extent([72, 136, 2, 56], crs)
    
14. ax.stock_img()
    
15. ax.add_geometries(geoms, crs, fc='none', ec='k', lw=0.5)
    
16. ax.add_geometries(geom, crs, fc='r', ec='k', lw=0.5)
    
17. plt.show()

复制代码

天地图行政区划可视化.zip (4.2 MB, 下载次数: 72)

2024-4-30 21:06 上传

点击文件名下载附件

灭火器

口蘑战士 发表于 2024-5-4 15:29
国界可以手动用省界合成
lz可以出个教程不

1. import json
   
2. 
   
3. import shapely.geometry as sgeom
   
4. import matplotlib.pyplot as plt
   
5. import cartopy.crs as ccrs
   
6. from shapely.ops import unary_union
   
7. 
   
8. with open('中国_省.geojson') as f:
   
9.     geoj = json.load(f)
   
10. 
    
11. # 提取省界多边形
    
12. provinces = []
    
13. for feature in geoj['features']:
    
14.     geometry = feature['geometry']
    
15.     if geometry['type'] in ('Polygon', 'MultiPolygon'):
    
16.         provinces.append(sgeom.shape(geometry))
    
17. 
    
18. # 省界合成国界
    
19. country = unary_union(provinces)
    
20. 
    
21. # 构造国界的GeoJSON
    
22. new_geoj = {
    
23.     'type': 'FeatureCollection',
    
24.     'features': [
    
25.         {
    
26.             'type': 'Feature',
    
27.             'geometry': sgeom.mapping(country),
    
28.             'properties': {'name': '国界'}
    
29.         }
    
30.     ]
    
31. }
    
32. 
    
33. with open('中国_国.geojson', 'w') as f:
    
34.     json.dump(new_geoj, f)
    
35. 
    
36. # 画个图看看
    
37. crs = ccrs.PlateCarree()
    
38. ax = plt.axes(projection=crs)
    
39. ax.set_extent((70, 140, 0, 60), crs)
    
40. ax.add_geometries(country, crs, fc='none', ec='k')
    
41. plt.show()

复制代码

用 geopandas 确实会更方便，但我刚把带有 geopandas 的环境搞崩了，所以这里就写个用 shapely 的。

tulalang

顺着lz说一下，geopandas可以直接把geojson转成shp

口蘑战士

国界可以手动用省界合成
lz可以出个教程不

小羊趴在云朵上

有用的分享，学到啦

口蘑战士

灭火器 发表于 2024-5-5 09:41
用 geopandas 确实会更方便，但我刚把带有 geopandas 的环境搞崩了，所以这里就写个用 shapely 的。

大佬！根据这些文件可以产出九段线的geojson么

灭火器

口蘑战士 发表于 2024-5-6 17:37
大佬！根据这些文件可以产出九段线的geojson么

省市县三个文件里都以 MultiLineString 的形式带九段线，不过跟未定国境线、海上省界线等混在一起，名字都被设成了“国境线”。你循环一下把所有线元素提取出来，然后画出来试试，就能从这些线中找出九段线。

1163087557

口蘑战士

灭火器 发表于 2024-5-6 18:42
省市县三个文件里都以 MultiLineString 的形式带九段线，不过跟未定国境线、海上省界线等混在一起，名字 ...

好的好的！