快捷導(dǎo)航

python中shapefile庫(kù)讀取shapefile文件信息

更新時(shí)間：2023年06月29日 10:03:35 作者：努力搬磚的giser

本文主要介紹了python中shapefile庫(kù)讀取shapefile文件信息，文中通過(guò)示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值，需要的朋友們下面隨著小編來(lái)一起學(xué)習(xí)學(xué)習(xí)吧

shapefile是GIS中非常重要的一種數(shù)據(jù)類(lèi)型，在ArcGIS中被稱為要素類(lèi)(Feature Class)，主要包括點(diǎn)(point)、線(polyline)和多邊形(polygon)。Python腳本是ArcGIS官方推薦的腳本語(yǔ)言，通過(guò)Python腳本能夠很方便的調(diào)用ArcGIS中的各種工具和函數(shù)批量完成所需操作。本文所介紹的這個(gè)庫(kù)(Python Shapefile Library)是一個(gè)Python庫(kù)，用于在Python腳本中對(duì)ArcGIS中的Shapefile文件(.shp，.shx，.dbf等格式)進(jìn)行讀寫(xiě)操作。

開(kāi)發(fā)準(zhǔn)備

安裝命令：pip install pyshp

使用導(dǎo)入：import shapefile

Shapefile文件的讀操作

通過(guò)創(chuàng)建Reader類(lèi)的對(duì)象進(jìn)行shapefile文件的讀操作。

file = shapefile.Reader('shapefile name')

“幾何數(shù)據(jù)”通過(guò)Reader類(lèi)的shapes( )和shape( )方法來(lái)讀取，二者的區(qū)別在于：shapes()方法不需要指定參數(shù)，其返回值是一個(gè)列表，包含該文件中所有的"幾何數(shù)據(jù)"對(duì)象，而shape( )方法則需要通過(guò)指定參數(shù)返回所需要的"幾何數(shù)據(jù)"對(duì)象。

shape = file.shape(i)#讀取第i+1個(gè)要素，索引序列從0開(kāi)始(參數(shù)是整數(shù)類(lèi)型)
shapes = file.shapes()#讀取全部要素

# -*- coding: utf-8 -*-
import shapefile# 使用pyshp
file = shapefile.Reader("E://arcgisData//行政區(qū)劃數(shù)據(jù)//省界_region.shp")#讀取
#讀取元數(shù)據(jù)
print(str(file.shapeType))  # 輸出shp類(lèi)型
print(file.encoding)# 輸出shp文件編碼
print(file.bbox)  # 輸出shp的文件范圍（外包矩形）
print(file.numRecords)  # 輸出shp文件的要素?cái)?shù)據(jù)
print(file.fields)# 輸出所有字段信息
# print(file.records())  # 輸出所有屬性表

補(bǔ)充 shp類(lèi)型信息

shp屬性表字段信息

字段索引	字段類(lèi)型
C	字符，文字
N	數(shù)字，帶或不帶小數(shù)
F	浮動(dòng)（與“N”相同）
L	邏輯，表示布爾值True / False值
D	日期
M	備忘錄，在GIS中沒(méi)有意義，而是xbase規(guī)范的一部分

展示單個(gè)要素
實(shí)現(xiàn)代碼：

import shapefile  # 使用pyshp
from matplotlib import pyplot as plt
file = shapefile.Reader("E://arcgisData//行政區(qū)劃數(shù)據(jù)//省界_region.shp")
border_shape = file
# 通過(guò)創(chuàng)建reader類(lèi)的對(duì)象進(jìn)行shapefile文件的讀取
# border_points
border = border_shape.shapes()
# .shapes()讀取幾何數(shù)據(jù)信息，存放著該文件中所有對(duì)象的 幾何數(shù)據(jù)
# border是一個(gè)列表
border_points = border[0].points
print(border_points)# 返回第1個(gè)對(duì)象的所有點(diǎn)坐標(biāo)
# border_points = [(x1,y1),(x2,y2),(x3,y3),…]
x, y = zip(*border_points)
# x=(x1,x2,x3,…)
# y=(y1,y2,y3,…)
fig, ax = plt.subplots()  # 生成一張圖和一張子圖
# plt.plot(x,y,'k-') # x橫坐標(biāo) y縱坐標(biāo) ‘k-'線性為黑色
plt.plot(x, y, color='#6666ff', label='fungis')  # x橫坐標(biāo) y縱坐標(biāo) ‘k-'線性為黑色
ax.grid()  # 添加網(wǎng)格線
ax.axis('equal')
plt.show()

展示效果：

shapefile文件的寫(xiě)操作

創(chuàng)建點(diǎn)狀要素類(lèi)

import shapefile  # 使用pyshp
file = shapefile.Writer("E://arcgisData//test//站點(diǎn)_2")#新建數(shù)據(jù)存放位置
#創(chuàng)建兩個(gè)字段
file.field('FIRST_FLD')
file.field('SECOND_FLD','C','40')     #'SECOND_FLD'為字段名稱，C代表數(shù)據(jù)類(lèi)型為字符串，長(zhǎng)度為40
file.point(1,1)
file.record('First','Point')
file.point(300,10)
file.record('Second','Point')
file.point(400,30)
file.record('Third','Point')
file.point(200,200)
file.record('Fourth','Point')
#寫(xiě)入數(shù)據(jù)
file.close()

創(chuàng)建線狀要素類(lèi)

import shapefile  # 使用pyshp
file = shapefile.Writer("E://arcgisData//test//polyline")#新建數(shù)據(jù)存放位置
#創(chuàng)建兩個(gè)字段
file.field('FIRST_FLD')
file.field('SECOND_FLD','C','40')     #'SECOND_FLD'為字段名稱，C代表數(shù)據(jù)類(lèi)型為字符串，長(zhǎng)度為40
file.line([[[1,5],[5,5],[5,1],[3,3],[1,1]]])
file.record('First','polyline')
file.line([[[1,500],[300,30],[1,16]]])
file.record('Second','polyline')
#寫(xiě)入數(shù)據(jù)
file.close()

創(chuàng)建面狀要素類(lèi)

import shapefile  # 使用pyshp
file = shapefile.Writer("E://arcgisData//test//polygon")#新建數(shù)據(jù)存放位置
#創(chuàng)建兩個(gè)字段
file.field('FIRST_FLD')
file.field('SECOND_FLD','C','40')     #'SECOND_FLD'為字段名稱，C代表數(shù)據(jù)類(lèi)型為字符串，長(zhǎng)度為40
file.poly([[[1,5],[5,5],[5,1],[3,3],[1,1]]])
file.record('First','polygon')
file.poly([[[1,500],[300,30],[1,16],[1,500]]])
file.record('Second','polygon')
#寫(xiě)入數(shù)據(jù)
file.close()

到此這篇關(guān)于python中shapefile庫(kù)讀取shapefile文件信息的文章就介紹到這了,更多相關(guān)python讀取shapefile文件信息內(nèi)容請(qǐng)搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持腳本之家！

您可能感興趣的文章: