python 根据shp文件解析经纬度市县信息

原创已于 2022-06-24 15:19:05 修改 · 1k 阅读

2 ·

CC 4.0 BY-SA版权

文章标签：

#python

于 2022-06-24 14:56:52 首次发布

Python 专栏收录该内容

10 篇文章

订阅专栏

该博客介绍了如何利用Python的Shapely和shapefile库读取和操作Shapefile地理数据。通过示例展示了如何查找指定经纬度坐标点（43.146, 125.247）所在省份、城市和县的名称。通过遍历shapefile的记录和边界，实现了基于空间位置的快速查询。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

from shapely.geometry import Point
from shapely.geometry.polygon import Polygon
from shapely.geometry import shape
import time
import shapefile
try:
    file = shapefile.Reader(os.getcwd()+r"\\shp\HEBEI.shp", encoding='gbk')

except:
    print('No shpfile!')
else:
    attribute_box=[]
    border_shape = file
    border = border_shape.shapes()
    attribute = file.records()
    sf = file.shapeRecords()
    for i in range(len(sf)):
        attribute_box.append(sf[i].shape.bbox)
        # print(attribute[i])
        # print(sf[i].shape.bbox)
    print('Success read shpfile')
    print(attribute_box)
lat = 43.146
lon = 125.247
point = Point(lon, lat)
sheng = '' 
shi = '' 
xian = ''  
for i in range(len(sf)):
    # print(sf[i].shape.bbox)
    [lon_min,lat_min,lon_max,lat_max]=attribute_box[i]
    if lat_min<lat<lat_max and lon_min<lon<lon_max:            ###提高速度
        # if point.within(shape(sf[i].shape)):
        if shape(sf[i].shape).contains(point):
            geo_attr = attribute[i]
            # print(geo_attr)
            xian = geo_attr['XIAN_NAME']
            shi = geo_attr['SHI_NAME]
            sheng = geo_attr['SHENG_NAME']
            print(lg, ly, lc)
            break