利用pyshp包给shapefile文件添加字段的实例

yipeiwu_com6年前Python基础

在已有的shapefile文件的基础上增加字段:

 # -*- coding:gb2312 -*-
import shapefile

r=shapefile.Reader(r"C:\Users\skfzh\Documents\python地理空间分析指南\第五章\5.5.4\output\NYC_MUSEUMS_UTM")
w=shapefile.Writer(r"C:\Users\skfzh\Documents\python地理空间分析指南\第五章\5.5.4\output\NYC_MUSEUMS_UTM_update",shapeType=r.shapeType)
w.fields=list(r.fields)
#新增加两个字段
w.field("LAT","F",8,5)
w.field("LON","F",8,5)

#将另外一个文件中的坐标点的信息存入新增加的两个字段
geo=shapefile.Reader(r"C:\Users\skfzh\Documents\python地理空间分析指南\第五章\5.5.4\NYC_MUSEUMS_GEO\NYC_MUSEUMS_GEO")
i=0  
for rec in r.iterShapeRecords():
 ls=rec.record
 lon,lat=geo.shape(i).points[0]  #可以再此处增加判断语句,限制i的大小。防止增加的字段数目不匹配问题
 # print(geo.shape(i).points)
 ls.extend([lat,lon])
 # print(ls)
 w.record(*ls)
 w.shape(rec.shape)
 i+=1

w.close()

以上这篇利用pyshp包给shapefile文件添加字段的实例就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持【听图阁-专注于Python设计】。

相关文章

Python 硬币兑换问题

硬币兑换问题: 给定总金额为A的一张纸币,现要兑换成面额分别为a1,a2,....,an的硬币,且希望所得到的硬币个数最少。 # 动态规划思想 dp方程式如下 # dp[0] =...

详细解读Python的web.py框架下的application.py模块

本文主要分析的是web.py库的application.py这个模块中的代码。总的来说,这个模块主要实现了WSGI兼容的接口,以便应用程序能够被WSGI应用服务器调用。WSGI是Web...

Python实现基于PIL和tesseract的验证码识别功能示例

Python实现基于PIL和tesseract的验证码识别功能示例

本文实例讲述了Python实现基于PIL和tesseract的验证码识别功能。分享给大家供大家参考,具体如下: 之前搞这个搞了一段时间,后面遇到了点小麻烦,导致识别率太低了,最多也就百分...

跟老齐学Python之有点简约的元组

关于元组,上一讲中涉及到了这个名词。本讲完整地讲述它。 先看一个例子: >>>#变量引用str >>> s = "abc" >>>...

python list转矩阵的实例讲解

如下所示: <pre name="code" class="python">#list转矩阵,矩阵列合并 x = [[1.2,2.2,1.4],[1.3,2.4,2.1]...