基于python的pyshp库读取.shp数据来获取中国城市边界的经纬度数据，并生成hdf文件

在之前博客中记录了如何获取省界经纬度数据
下面介绍如何获取城市区域的边界经纬度数据，具体思路同上一篇博客
https://blog.csdn.net/weixin_43718675/article/details/93410875

import shapefile
import numpy as np
from mpl_toolkits.basemap import Basemap
import matplotlib
import matplotlib.pyplot as plt
import os
import h5pyfile=shapefile.Reader('gadm36_CHN_shp/gadm36_CHN_2.shp')  #读取城市.shp文件
shapes=file.shapes()   #获取point
records=file.records() #获取省-市名称

打开records文件，如下

pro_city_points=[] #建立城市边界列表
pro_city_names=[] #建立城市名称列表
pro_names=[]  #建立省份列表for i in range(len(shapes)):points=shapes[i].points #h获取经纬度数据    pro_name=file.records()[i][3]  #获取省名称city_name=file.records()[i][6]  #获取市区名称city_ch_name=file.records()[i][8]  #获取市区中文名称pro_city_name=[pro_name,city_name,city_ch_name]lon =[]  lat  =[] #将每个tuple的lon和lat组合起来[lon.append(points[i][0]) for i in range(len(points))]  [lat.append(points[i][1]) for i in range(len(points))]   lon=np.array(lon).reshape(-1,1)lat=np.array(lat).reshape(-1,1) loc=np.concatenate((lon,lat),axis=1)pro_city_points.append(loc)pro_city_names.append(pro_city_name)pro_names.append(pro_name)pro_names=list(set(pro_names))  #使用集合去掉重复的省名称

打开，如下图

将各个城市的边界数据保存为hdf文件，

file_name='province_city_boundary_lon_lat.hdf'  #保存的目标文件名和路径
f=h5py.File(file_name,mode='w')  #创建一个hdf文件for pro_name in pro_names[0:]:#每一个省份创建一个群组 pro = f.create_group(pro_name)for pro_city_name,points in zip(pro_city_names,pro_city_points):#将对应省份的城市的放在一个group下if pro_city_name[0]==pro_name:  city=pro.create_group(pro_city_name[1])city['longitude']=points[:,0]city['latitude']=points[:,1]f.close() #写入完成后，记得关掉。

生成结果如下：

下载使用链接为：
http://bbs.06climate.com/forum.php?mod=viewthread&tid=92137&extra=
如果没有金币下载，可以私聊我，我发给你

基于python的pyshp库读取.shp数据来获取中国城市边界的经纬度数据，并生成hdf文件相关推荐

下载地图包，并基于python的pyshp库读取.shp数据来获取中国省界的经纬度数据
目前画地图的软件都能很方便的调用省界数据,但是如果想要单独对省界做一些个性化设置,如设置宽度.样式.颜色什么的,就需要将省界数据单独拿出来进行设置了. 软件在画地图时,描述的边界都是一些列经纬度格点的 ...
python操作文件的库_Python使用pyshp库读取shapefile信息的方法
通过pyshp库,可以读写Shapefile文件,查询相关信息,github地址为 https://github.com/GeospatialPython/pyshp#reading-shapefil ...
python文件读写用到的库_Python使用pyshp库读取shapefile信息的方法
通过pyshp库,可以读写shapefile文件,查询相关信息,github地址为 import shapefile # 使用pyshp库 file = shapefile.reader(" ...
python中pyshp_Python使用pyshp库读取shapefile信息的方法
通过pyshp库,可以读写Shapefile文件,查询相关信息,github地址为 import shapefile # 使用pyshp库 file = shapefile.Reader(" ...
Python 基于Python从mysql表读取千万数据实践
基于Python 从mysql表读取千万数据实践 by:授客 QQ:1033553122 场景: 有以下两个表,两者都有一个表字段,名为waybill_no,我们需要从tl_waybill_b ...
python人脸识别库_基于Python的face_recognition库实现人脸识别
Python Python开发 Python语言基于Python的face_recognition库实现人脸识别一.face_recognition库简介 face_recognition是Pyt ...
pandas读取csv写入mysql_使用python的pandas库读取csv文件保存至mysql数据库
第一:pandas.read_csv读取本地csv文件为数据框形式 data=pd.read_csv('G:\data_operation\python_book\chapter5\\sales.cs ...
基于python的npcap库与dpkt库实现抓包及存储
基于python的npcap库与dpkt库实现抓包及存储 import pcap import dpkt import socket import sys import getopt import o ...
Python用pydicom库读取dicom文件
Python用pydicom库读取dicom文件并调用方法1:通过group.element读取方法2:通过tag读取方法3:get()方法通用 import pydicompath = 'x ...

基于python的pyshp库读取.shp数据来获取中国城市边界的经纬度数据，并生成hdf文件

基于python的pyshp库读取.shp数据来获取中国城市边界的经纬度数据，并生成hdf文件相关推荐

最新文章

热门文章