(二十三)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文件中。

(二十三)arcpy开发&利用GetRasterProperties_management获取栅格数据相关信息_第1张图片

当然在写这个程序的时候遇到如下的错误。

错误一:

(二十三)arcpy开发&利用GetRasterProperties_management获取栅格数据相关信息_第2张图片

这个提示是因为我的txt文件路径下没有创建好。

错误二:

(二十三)arcpy开发&利用GetRasterProperties_management获取栅格数据相关信息_第3张图片

这个错误有点让人一下子找不到解决问题,根据提示,显示统计没有是无效的。这时我突然想起用arcgis desktop打开来看一下文件信息。

(二十三)arcpy开发&利用GetRasterProperties_management获取栅格数据相关信息_第4张图片

如上图所示,显示在统计一栏中都是空白的,目前没有找到什么的办法统计出红色框里面的数据,于是我重新导出了一份,便有了下面的信息。

(二十三)arcpy开发&利用GetRasterProperties_management获取栅格数据相关信息_第5张图片


                          更多内容,请微信扫二维码关注公众号,或者加入arcpy开发qq学习群:487352121

                                                                  

                                                                  

                                        

你可能感兴趣的:(arcpy,python,arcgis,插件开发)