# -*- 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实现字符串替换时,可利用replace函数来实现,
具体代码为:stringold.replace(strfrom,strto),其中stringold就是需要更改的字符串,strfrom是需要替换的子字符串,strto是需要转换成的子字符串。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又有啥区别呢?还是第一种吧