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

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

让客户满意是我们工作的目标，不断超越客户的期望值来自于我们对这个行业的热爱。我们立志把好的技术通过有效、简单的方式提供给客户，将通过不懈努力成为客户在信息化领域值得信任、有价值的长期合作伙伴，公司提供的服务项目有：国际域名空间、虚拟主机、营销软件、网站建设、柳州网站维护、网站推广。

#scrip

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

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

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

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

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数据，正在处理第"+str(shps.index(sh)+1)+"个："+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,"")#注释则裁剪的栅格不删除，不注释删掉。

    # except:

    print(outname+'cannot delete')

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

print "Finish!"

python如何替换shp字段内容

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

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

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又有啥区别呢？还是第一种吧

网页题目：python的shp函数 python shape0
链接分享：http://scyanting.com/article/dosjsed.html