使用pyshp包進(jìn)行shapefile文件修改的例子-創(chuàng)新互聯(lián)

由于最近在處理shp文件,想要跳出arcpy的限制,所以打算學(xué)習(xí)一下pyshp包的使用方法。在使用《Python地理空間分析指南(第2版)》的時(shí)候發(fā)現(xiàn)書(shū)中部分代碼由于版本更新,無(wú)法運(yùn)行。開(kāi)貼記錄踩過(guò)的雷。

創(chuàng)新互聯(lián)公司是一家專(zhuān)注于做網(wǎng)站、成都網(wǎng)站建設(shè)與策劃設(shè)計(jì),瀘水網(wǎng)站建設(shè)哪家好?創(chuàng)新互聯(lián)公司做網(wǎng)站,專(zhuān)注于網(wǎng)站建設(shè)10余年,網(wǎng)設(shè)計(jì)領(lǐng)域的專(zhuān)業(yè)建站公司;建站業(yè)務(wù)涵蓋:瀘水等地區(qū)。瀘水做網(wǎng)站價(jià)格咨詢(xún):028-86922220

這個(gè)問(wèn)題是出現(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)
#版本修改,路徑和類(lèi)型都要在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()

標(biāo)題名稱(chēng):使用pyshp包進(jìn)行shapefile文件修改的例子-創(chuàng)新互聯(lián)
標(biāo)題URL:http://muchs.cn/article14/dhgige.html

成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián),為您提供網(wǎng)頁(yè)設(shè)計(jì)公司、全網(wǎng)營(yíng)銷(xiāo)推廣服務(wù)器托管、面包屑導(dǎo)航做網(wǎng)站、App開(kāi)發(fā)

廣告

聲明:本網(wǎng)站發(fā)布的內(nèi)容(圖片、視頻和文字)以用戶(hù)投稿、用戶(hù)轉(zhuǎn)載內(nèi)容為主,如果涉及侵權(quán)請(qǐng)盡快告知,我們將會(huì)在第一時(shí)間刪除。文章觀(guān)點(diǎn)不代表本網(wǎng)站立場(chǎng),如需處理請(qǐng)聯(lián)系客服。電話(huà):028-86922220;郵箱:631063699@qq.com。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載,或轉(zhuǎn)載時(shí)需注明來(lái)源: 創(chuàng)新互聯(lián)

外貿(mào)網(wǎng)站制作