今天我们来学习一下关于获取栅格数据的一些信息,首先我在这里呢是遍历某个工作空间下面所有数据,然后读取这些栅格信息,最后将这些信息写入到txt文件中。其中有关arcpy获取栅格信息的标记可以参考esri的相关文档。http://resources.arcgis.com/zh-cn/help/main/10.1/index.html#//0017000000m7000000,下面我们来看一下具体的实现代码。
#coding=utf-8
import arcpy
outputTxtFile = r'D:\ArcpyProject\GetRasterPropertiesWriteTXT\result.txt'
arcpy.env.workspace = r'D:\TIF\dem\data_84'
rasters = arcpy.ListRasters("*")
print rasters
for raster in rasters:
minResult = arcpy.GetRasterProperties_management(raster, "MINIMUM")
maxResult = arcpy.GetRasterProperties_management(raster, "MAXIMUM")
meanResult = arcpy.GetRasterProperties_management(raster, "MEAN")
img_SR = arcpy.Describe(raster).spatialReference.name
minRes = minResult.getOutput(0)
maxRes = maxResult.getOutput(0)
meanRes = meanResult.getOutput(0)
# minRes=str(minResult)
# maxRes = str(maxResult)
# meanRes = str(meanResult)
f = open(outputTxtFile, 'a+')#打开txt文件,并追加信息
f.write(raster + "," + minRes + "," + meanRes + "," + maxRes + "\n")
print "finish"
最后将信息写入到txt文件中。
当然在写这个程序的时候遇到如下的错误。
错误一:
这个提示是因为我的txt文件路径下没有创建好。
错误二:
这个错误有点让人一下子找不到解决问题,根据提示,显示统计没有是无效的。这时我突然想起用arcgis desktop打开来看一下文件信息。
如上图所示,显示在统计一栏中都是空白的,目前没有找到什么的办法统计出红色框里面的数据,于是我重新导出了一份,便有了下面的信息。
更多内容,请微信扫二维码关注公众号,或者加入arcpy开发qq学习群:487352121