python shp文件_python中shapefile学习（pyshp）

python中处理空间地理数据时，往往需要处理shp文件PyShp package可以用于shp文件的读取、写入等操作

https://pypi.org/project/pyshp/

以下用代码形式显示其基本读取操作

import shapefile #导入模块

sf = shapefile.Reader('shapefiles/blockgroups.shp') #读入shp文件

sf.shapeType #查询shp属性，NULL = 0; POINT =1;POLYLINE =3; POLYGON =5;MULTIPOINT =8; POINTZ=11;POLYLINZ=13;POLYGONZ=15;MULTIPOINTZ =18; POINTM=21;POLYLINEM =23;POLYGONM=25;MULTIPOINTM=28;MULTIPATCH=31

sf.shapeTypeName #shp属性的具体名称

sf.bbox #shp的bounds

sf.__geo_interface__ #获取整个shp文件属性(一般将shp文件转换为GeoJSON用到)

sf.shape(index) #利用index获取其中单一shape

shapes = sf.shapes() #读取shape的geometry，输出为一个矩阵

#每个shapes都可以独立查询其属性

shapes[0].shapeType

shapes[0].shapeTypeName

shapes[0].bbox

shapes[0].points #线上的每个点坐标(x,y)

#获得shp文件的fields，每个field含有以下属性

##Field name 描述数据的column index

##Field类型：'C'-字符；'N'-数字；'F'-浮点数；'L'-逻辑数；'D'-时间；'M'-Memo

## Field length 在这个列index下数据长度

## Decimal length 小数点后的有效数位数

fields = sf.fields

# 获取 records

records = sf.records()

# 将records整理为字典型

dct = rec.as_dict()

# 同时读取geometry and records

shapeRecs = sf.shapeRecords()