使用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 實(shí)現(xiàn)調(diào)用子文件下的模塊方法
今天小編就為大家分享一篇python 實(shí)現(xiàn)調(diào)用子文件下的模塊方法,具有很好的參考價(jià)值,希望對(duì)大家有所幫助。一起跟隨小編過來看看吧2018-12-12解決tensorflow由于未初始化變量而導(dǎo)致的錯(cuò)誤問題
今天小編就為大家分享一篇解決tensorflow由于未初始化變量而導(dǎo)致的錯(cuò)誤問題,具有很好的參考價(jià)值,希望對(duì)大家有所幫助。一起跟隨小編過來看看吧2020-01-01Python3實(shí)現(xiàn)統(tǒng)計(jì)單詞表中每個(gè)字母出現(xiàn)頻率的方法示例
這篇文章主要介紹了Python3實(shí)現(xiàn)統(tǒng)計(jì)單詞表中每個(gè)字母出現(xiàn)頻率的方法,涉及Python針對(duì)文件的讀取、遍歷、統(tǒng)計(jì)等相關(guān)操作技巧,需要的朋友可以參考下2019-01-01python函數(shù)缺省值與引用學(xué)習(xí)筆記分享
有關(guān)一個(gè)在函數(shù)參數(shù)設(shè)置缺省值與引用的問題,這個(gè)問題是大多數(shù)Pythoner可能會(huì)忽視的問題,作個(gè)筆記,以備后閱,同時(shí)供需要的朋友參考2013-02-02Python 多線程共享變量的實(shí)現(xiàn)示例
這篇文章主要介紹了Python 多線程共享變量的實(shí)現(xiàn)示例,文中通過示例代碼介紹的非常詳細(xì),對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧2020-04-04python使用正則表達(dá)式替換匹配成功的組并輸出替換的次數(shù)
正則表達(dá)式是一個(gè)特殊的字符序列,它能幫助你方便的檢查一個(gè)字符串是否與某種模式匹配。這篇文章主要介紹了python使用正則表達(dá)式替換匹配成功的組并輸出替換的次數(shù),需要的朋友可以參考下2017-11-11