水滴石穿，绳锯木断

在使用arcpython时，经常要用到坐标系的信息，直接获取比较麻烦。但是，如果手头已经有了一个你想要的坐标系文件，可以通过读取它的坐标系文件，直接用就可以了。下面就是使用arcgis的重投影的函数，就用到了你想要投影到的坐标系，比如wgs1984#获取特定坐标系的方法#方法1：直接读取已有的文件的坐标系信息dataset = r'D:\aa\weihe\RiverBorder.shp'spatial_ref = arcpy.Describe(dataset).spatialReference.

说明：现实中可能获取了很多存在excel中的点坐标信息，当需要获取shape文件时，常规的方法是通过arcgis中add XY,生成event图层，然后导出为shape的点文件，如果需要再生成line或者polygon则需要其他的工具。比较繁琐，基于此，本代码实现批处理的有excel转为Polygon文件。# -*- coding: utf-8 -*-import arcpy#

1、编译器为了更好的编译，调试等，使用好的IDE很关键，开发python时推荐使用PyCharm软件，下载地址： http://www.xt70.com/soft/show/22323.html安装后就可以系统管理代码，包括按F8逐步调试。

雪水当量数据(EASE-GridData)处理流程1、 数据描述该数据描述的是冰雪覆盖融化后相当于的水量值。是由欧空局(ESA)数据用户要素支持的全球雪盖项目支持，覆盖了北半球雪要素：雪范围（SE）和雪水当量（SWE），时尺度分别为17年和34年。其中SWE是被动微波和地面气象观测融合而成，时间覆盖为1979-2013日尺度。空间分辨率是25km，全球覆盖北纬35-85。       

图像归一化就不多说了，就是(数值-min)/(max-min)，把结果都划归到0-1范围，便于不同变量之间的比较，取消了不同数量差别。第一个方法，需要对栅格数据预先知道取值范围。第二种方法，比较好点，直接读取属性

1、 setnull函数 arcpy.CheckOutExtension("Spatial") outSetNull = SetNull(inputRaster,inputRaster,"VALUE < 0") outSetNull.save(outputRaster)注意：一定要检查许可否则会报错；另外，其中的whereclaus...

1、获取一个文件的子串如：  str1="myfirst.tif",想获取不带扩展名的文件名最简单方式为str2=str1[0：-4]2、获取指定文件夹下的所有指定类型文件arcpy.env.workspace=inputDirfileList=arcpy.ListFiles("*.tif")即可获取所有tif文件的链表3、获取指定文件路径下的所有文件名称d

对于Python不能直接调用arcgis中的栅格计算器工具。但是可以通过使用Raster函数实现。如下面代码实现了每年12个月做一个累加，然后输出一个年值的功能：for year in range(1960,1965):    print year    temp=0     #必须提前定义要不报错    for month in range(1,13):        p

生成字符串变量str='python String function' 字符串长度获取：len(str)例：print '%s length=%d' % (str,len(str))连接字符串sStr1 = 'strcat'sStr2 = 'append'sStr1 += sStr2print sStr1复制字符串#strcpy(sStr1,sS

1、  获取指定位置的点值：OutputFile = open(statisticResultTXT, 'w')  cellvalue=arcpy.GetCellValue_management(inputfilepath+filenname+".tif",staionXY, "1")OutputFile.write(stationID+"_"+filenname+""+str(cel

首先通过批量裁切为arcgis的extractbymask工具，代码如下： 其次，利用arcgis的Python脚本的GetRasterProperties_management方法，代码如下：