使用pyshp包進(jìn)行shapefile文件修改的例子
由于最近在處理shp文件,想要跳出arcpy的限制,所以打算學(xué)習(xí)一下pyshp包的使用方法。在使用《Python地理空間分析指南(第2版)》的時(shí)候發(fā)現(xiàn)書中部分代碼由于版本更新,無法運(yùn)行。開貼記錄踩過的雷。
這個(gè)問題是出現(xiàn)在5.5.4的shapefile文件修改中:
# -*- coding:gb2312 -*-
import shapefile
import utm
file_path=r"C:\Users\skfzh\Documents\python地理空間分析指南\第五章\5.5.4\NYC_MUSEUMS_GEO\NYC_MUSEUMS_GEO.shp"
r=shapefile.Reader(file_path)
# print(list(r.fields))
# print(r.shapeTypeName)
#版本修改,路徑和類型都要在writer里面定義
w=shapefile.Writer(r"C:\Users\skfzh\Documents\python地理空間分析指南\第五章\5.5.4\output\NYC_MUSEUMS_UTM",\
shapeType=r.shapeType)
w.fields=list(r.fields[1:])
for rec in r.iterShapeRecords():#新版本已經(jīng)刪除了w.records
# print(*rec.record)
w.record(*rec.record)
for sha in r.iterShapes():
# print(sha.points[0])
lon,lat=sha.points[0]
y,x,zone,band=utm.from_latlon(lat,lon)
w.point(x,y)
w.close()
以上這篇使用pyshp包進(jìn)行shapefile文件修改的例子就是小編分享給大家的全部內(nèi)容了,希望能給大家一個(gè)參考,也希望大家多多支持腳本之家。
相關(guān)文章
python中p-value的實(shí)現(xiàn)方式
今天小編就為大家分享一篇python中p-value的實(shí)現(xiàn)方式,具有很好的參考價(jià)值,希望對大家有所幫助。一起跟隨小編過來看看吧2019-12-12
Python3中對json格式數(shù)據(jù)的分析處理
這篇文章主要介紹了Python3中對json格式數(shù)據(jù)的分析處理,本文給大家介紹的非常詳細(xì),對大家的學(xué)習(xí)或工作具有一定的參考借鑒價(jià)值,需要的朋友可以參考下2021-01-01
深入理解Python虛擬機(jī)中常見魔術(shù)方法的使用
本文主要給大家介紹在 python 當(dāng)中與數(shù)學(xué)計(jì)算相關(guān)的一些常見的魔術(shù)方法,是在很多科學(xué)計(jì)算的包當(dāng)中都使用到的魔術(shù)方法,感興趣的小伙伴可以了解一下2023-05-05
python用pickle模塊實(shí)現(xiàn)“增刪改查”的簡易功能
本篇文章主要介紹了python用pickle模塊實(shí)現(xiàn)“增刪改查”的簡易功能,具有一定的參考價(jià)值,感興趣的小伙伴們可以參考一下。2017-06-06
Python實(shí)現(xiàn)監(jiān)視程序的內(nèi)存使用情況
我們使用Python和它的數(shù)據(jù)處理庫套件進(jìn)行大量數(shù)據(jù)處理時(shí)候,可能使用了大量的計(jì)算資源,那么如何監(jiān)視程序的內(nèi)存使用情況就顯得尤為重要,下面我們就來了解一下具體實(shí)現(xiàn)方法吧2023-12-12

