![](https://img-blog.csdnimg.cn/20201014180756927.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
Python
忠言睿长
从事3S技术与水文气象研究
展开
-
ArcPython常用函数操作记录
在使用arcpython时,经常要用到坐标系的信息,直接获取比较麻烦。但是,如果手头已经有了一个你想要的坐标系文件,可以通过读取它的坐标系文件,直接用就可以了。下面就是使用arcgis的重投影的函数,就用到了你想要投影到的坐标系,比如wgs1984#获取特定坐标系的方法#方法1:直接读取已有的文件的坐标系信息dataset = r'D:\aa\weihe\RiverBorder.shp'spatial_ref = arcpy.Describe(dataset).spatialReference.原创 2021-01-22 19:56:36 · 640 阅读 · 0 评论 -
excel中的点文件直接生成polygon文件
说明:现实中可能获取了很多存在excel中的点坐标信息,当需要获取shape文件时,常规的方法是通过arcgis中add XY,生成event图层,然后导出为shape的点文件,如果需要再生成line或者polygon则需要其他的工具。比较繁琐,基于此,本代码实现批处理的有excel转为Polygon文件。# -*- coding: utf-8 -*-import arcpy#原创 2017-10-09 17:10:17 · 1709 阅读 · 0 评论 -
Python学习
1、编译器为了更好的编译,调试等,使用好的IDE很关键,开发python时推荐使用PyCharm软件,下载地址: http://www.xt70.com/soft/show/22323.html安装后就可以系统管理代码,包括按F8逐步调试。原创 2017-03-13 16:32:04 · 493 阅读 · 0 评论 -
雪水当量数据(EASE-Grid Data)处理方法
雪水当量数据(EASE-GridData)处理流程1、 数据描述该数据描述的是冰雪覆盖融化后相当于的水量值。是由欧空局(ESA)数据用户要素支持的全球雪盖项目支持,覆盖了北半球雪要素:雪范围(SE)和雪水当量(SWE),时尺度分别为17年和34年。其中SWE是被动微波和地面气象观测融合而成,时间覆盖为1979-2013日尺度。空间分辨率是25km,全球覆盖北纬35-85。原创 2016-09-16 20:41:10 · 10567 阅读 · 7 评论 -
arcpython设置处理边界extent
原创 2016-08-02 16:13:23 · 1798 阅读 · 0 评论 -
Arcgis+Python实现对栅格归一化处理
图像归一化就不多说了,就是(数值-min)/(max-min),把结果都划归到0-1范围,便于不同变量之间的比较,取消了不同数量差别。第一个方法,需要对栅格数据预先知道取值范围。第二种方法,比较好点,直接读取属性原创 2016-01-15 10:25:53 · 12268 阅读 · 12 评论 -
arcPython细节汇总-setll函数和地图计算Raster Calculator函数
1、 setnull函数 arcpy.CheckOutExtension("Spatial") outSetNull = SetNull(inputRaster,inputRaster,"VALUE < 0") outSetNull.save(outputRaster)注意:一定要检查许可否则会报错;另外,其中的whereclaus...原创 2015-05-30 06:56:51 · 3752 阅读 · 0 评论 -
Python常用常见知识汇总
1、获取一个文件的子串如: str1="myfirst.tif",想获取不带扩展名的文件名最简单方式为str2=str1[0:-4]2、获取指定文件夹下的所有指定类型文件arcpy.env.workspace=inputDirfileList=arcpy.ListFiles("*.tif")即可获取所有tif文件的链表3、获取指定文件路径下的所有文件名称d原创 2015-06-16 23:30:07 · 710 阅读 · 0 评论 -
python实现栅格计算器功能
对于Python不能直接调用arcgis中的栅格计算器工具。但是可以通过使用Raster函数实现。如下面代码实现了每年12个月做一个累加,然后输出一个年值的功能:for year in range(1960,1965): print year temp=0 #必须提前定义要不报错 for month in range(1,13): p原创 2015-05-21 22:20:39 · 6048 阅读 · 1 评论 -
python字符串函数
生成字符串变量str='python String function' 字符串长度获取:len(str)例:print '%s length=%d' % (str,len(str))连接字符串sStr1 = 'strcat'sStr2 = 'append'sStr1 += sStr2print sStr1复制字符串#strcpy(sStr1,sS转载 2015-01-21 17:55:15 · 808 阅读 · 1 评论 -
python获取栅格点和面值
1、 获取指定位置的点值:OutputFile = open(statisticResultTXT, 'w') cellvalue=arcpy.GetCellValue_management(inputfilepath+filenname+".tif",staionXY, "1")OutputFile.write(stationID+"_"+filenname+""+str(cel原创 2015-04-08 08:55:14 · 4569 阅读 · 0 评论 -
arcgis脚本Python批量裁切与统计实现
首先通过批量裁切为arcgis的extractbymask工具,代码如下: 其次,利用arcgis的Python脚本的GetRasterProperties_management方法,代码如下:原创 2015-01-17 11:33:43 · 2597 阅读 · 0 评论