某位仁兄问我一个问题:
有一幅栅格数据影像,想要提取每一个小栅格的经纬度和属性值,该怎么操作?
对于新手来说,最简单的就是用ArcGIS进行提取操作,本着学习无界,一起进步的思想,写出来分享一下。
当然了,除了用GIS操作以外,还可以用python,matlab等编程实现这个目标,感兴趣的可以自己探究,其实就是对数据矩阵的处理,当然这里面包含空间位置,需要一些地理知识。
长话短说,今天重点讲提取,掐头去尾,重点说属性值的提取部分,总体步骤如下,图片看不清楚可保存放大看:
1.本例数据来源:地理空间数据云,四川省某地区DEM高程图。
所以,这个栅格数据影像里面的value值就是高程值。
掩膜提取自己感兴趣的区域,如图:
2.提取完毕后,可以使用识别工具查看栅格属性信息,如下图所示,Pixel value为1482,就表明这个点的高程值为1482。Location是经纬度信息。我们的目的就是提取Location(经纬度)和Pixel value(高程值).某仁兄要提取LAI,其实也是一样的道理,就是DEM变成了LAI或者NDVI等等…
3.打开arctoolbox下的转换工具下的栅格到点工具,如下图所示:
输入栅格为上一步提取的感兴趣区域,字段(field)选择value(其实就是属性值比如高程,叶面积指数什么的),输出随意,然后运行,本次得到点文件convert_dem。
4.放大影像可以看到,得到的点文件位于每一个栅格的中心位置,如下图所示。
5.打开arctoolbox下的提取工具下的提取值到点工具,如下图所示:
输入点文件为convert_dem,输入栅格为感兴趣区域栅格,输出随意,本次输出为Point_value。
6.打开point_value的属性数据表,如下图,可以看见高程值已经包含进去,也就是红框部分。
7.新建两个字段分别存放经纬度,字段选择DOUBLE类型。然后右键点击字段,选择计算几何工具,如下图所示,计算经纬度。
提醒:X Coor是经度,Y Coor是纬度。
8.计算结果如下图所示,至此,每一个像素的属性值(这里是DEM高程),经纬度,都已经存在属性数据表中了。
9.导出属性数据表.dbf文件,用Excel愉快的处理数据吧。
如有任何错误或疑问,可联系邮箱:chinesevoice@163.com(本人邮箱)
公众号:科研萌新
,