python的shp函数 python中的shape函数

使用shp数据批量裁剪栅格数据并统计均值

# -*- coding: UTF-8 -*-

#scrip

# 使用shp数据批量裁剪栅格数据并统计均值

print"使用shp数据批量裁剪栅格数据并统计结果均值"

###########修改部分###########

ras_file= r"D:\python\nc\GLEAM\Transpiration_TIF_OUT"  #待裁剪python的shp函数的原始栅格数据存储文件夹路径

suffix= 'tif'  #待裁剪的原始栅格数据后缀

bvalue=-999.0  #待裁剪的原始栅格数据背景值

shp_file= r"D:\python\shp"  #裁剪所需shp模板数据存储文件夹路径

ras_file_cut= r"D:\python\nc\GLEAM\Transpiration_YR_SWAT_mask"  #裁剪后栅格数据结果存储文件夹路径

txtname=r"D:\python\nc\GLEAM\Transpiration_YR_SWAT_mask"    #输出统计文本路径

#############################

#计算部分

import fenqutongji_arcpy

fenqutongji_arcpy.env.workspace=shp_file

shps=fenqutongji_arcpy.ListFeatureClasses()

fenqutongji_arcpy.env.workspace=ras_file

ras=fenqutongji_arcpy.ListRasters('*', suffix)

print "共有"+str(len(ras))+"个栅格数据"

#

print "Processing......"

for shin shps:

    shtmp=sh.encode('cp936')

shpfile=shp_file+"\\"+shtmp

print "共有"+str(len(shps))+"个shp数据python的shp函数，正在处理第"+str(shps.index(sh)+1)+"个python的shp函数："+shtmp

result=[]

for rsin ras:

        rstmp=rs.encode('cp936')

outname=ras_file_cut+"\\"+rstmp[0:len(rstmp)-4]+shtmp[0:len(shtmp)-4]+".tif"

        #arcpy.Clip_management(rstmp,"#",outname,shpfile,"#","ClippingGeometry")

        fenqutongji_arcpy.Clip_management(rstmp, "#", outname, shpfile, str(bvalue), "ClippingGeometry")#忽略无效值

        stats=fenqutongji_arcpy.GetRasterProperties_management(outname, "MEAN")

result.append(rstmp+'  '+str(stats)+"\n")

# try:

        #    arcpy.Delete_management(outname,"")#注释则裁剪的栅格不删除python的shp函数，不注释删掉。

        # except:

        print(outname+'cannot delete')

file(txtname+"\\"+shtmp[0:len(shtmp)-4]+".txt",'w').writelines(result)

print "Finish!"

相关问答

Q1: python如何替换shp字段内容

python实现字符串替换时，可利用replace函数来实现，

具体代码为：stringold.replace（strfrom,strto），其中stringold就是需要更改的字符串，strfrom是需要替换的子字符串，strto是需要转换成的子字符串。Python是一种跨平台的计算机程序设计语言，也是一种面向对象的动态类型语言，最初被设计用于编写自动化脚本。随着版本的不断更新和语言新功能的添加，越来越多被用于独立的.大型项目的开发。Python语言具有简洁性.易读性以及可扩展性，在国外用Python做科学计算的研究机构日益增多，一些知名大学已经采用Python来教授程序设计课程。

Q2: python函数调用

inname = r"C:\Python27\esri.shp"

outname = "outname.cst"

# 在此处调用该函数。函数体定义必须放在调用以前。可以通过import

read_ESRT_……(file = inname, fileOut = outname)

# 这两个参数只是字符串而已，指明你的文件路径。注意在python中，若有 \ 号，则最好使用 \\ 双斜杠，或者如上例，加上前缀 r

Q3: python转义双引号

为啥一定都用双引号呢？

gp.CalculateField_management(a, "Type", '"林地"') #外面是单引号，里面可以直接使用双引号，不必转义

如果使用变量代替可以这样：

a = “林地”

b = '"%s"' % a

gp.CalculateField_management(a, "Type", b)

也可以这样：

a = “林地”

b = a

gp.CalculateField_management(a, "Type", '"%s"' % b)

这与直接用a又有啥区别呢？还是第一种吧

Q4: 如何使用python对shp文件进行切片

取一个list的部分元素是非常常见的操作。比如，一个list如下：

L = ['Adam', 'Lisa', 'Bart', 'Paul']

取前3个元素，应该怎么做？

对这种经常取指定索引范围的操作，用循环十分繁琐，因此，Python提供了切片（Slice）操作符，能大大简化这种操作。

对应上面的问题，取前3个元素，用一行代码就可以完成切片：

L[0:3]

['Adam', 'Lisa', 'Bart']

L[0:3]表示，从索引0开始取，直到索引3为止，但不包括索引3。即索引0，1，2，正好是3个元素。

只用一个 : ，表示从头到尾：

L[:]

['Adam', 'Lisa', 'Bart', 'Paul']

因此，L[:]实际上复制出了一个新list。

切片操作还可以指定第三个参数：

L[::2]

['Adam', 'Bart']

第三个参数表示每N个取一个，上面的 L[::2] 会每两个元素取出一个来，也就是隔一个取一个。

把list换成tuple，切片操作完全相同，只是切片的结果也变成了tuple。

python的shp函数的介绍就聊到这里吧，感谢你花时间阅读本站内容，更多关于python中的shape函数、python的shp函数的信息别忘了在本站进行查找喔。