（二十三）arcpy开发&利用GetRasterProperties_management获取栅格数据相关信息

今天我们来学习一下关于获取栅格数据的一些信息，首先我在这里呢是遍历某个工作空间下面所有数据，然后读取这些栅格信息，最后将这些信息写入到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