1. 为什么你需要掌握栅格数据提取?
如果你正在处理遥感影像、环境监测数据、气象数据或者任何一张由无数小方格(像元)组成的“图片”,那么你迟早会遇到一个核心需求:把图上每个小格子的数值,以及它在地球上的精确位置,给批量“抠”出来,整理成一张规规矩矩的表格。这个表格,大概率就是Excel。
我遇到过太多这样的场景了。比如,生态学的朋友要分析一整片森林的NDVI(植被指数)变化,他需要知道每个像元点的指数值和经纬度;搞环境评估的同事,手里有一张土壤重金属污染浓度分布图,他得把每个位置的浓度值提取出来,才能做后续的空间统计分析。这时候,如果对着屏幕一个个点去手动记录,那简直是灾难。数据量动辄几十万、上百万个点,这活儿根本没法干。
所以,我们今天要聊的,就是如何用ArcGIS这套几乎人手一套的工具,不写一行代码,自动化、批量化地完成“栅格数据像元值与坐标值提取,并导出到Excel”的全过程。这听起来像是个高阶操作,但我跟你保证,它的核心步骤就三步,比想象中简单得多。关键在于理解每一步工具在“干什么”,以及如何避开那些我踩过的“坑”。掌握了这个流程,你就能把任何栅格图层,无论是TIFF、IMG还是GRID格式,都变成一份结构清晰、随时可用的表格数据,为你的数据分析、报告撰写或者模型输入扫清障碍。
2. 核心武器库:认识你将用到的ArcGIS工具
在开始动手之前,我们得先熟悉一下“兵器”。整个流程主要依赖ArcGIS工具箱里的几个核心工具,它们就像一条流水线上的不同工位,各司其职。
第一个工位:栅格转点(Raster to Point)。这是整个流程的起点,也是最关键的一步。它的作用,顾名思义,就是把一张完整的、由像元矩阵构成的栅格图,打散成一个个独立的点要素。你可以想象成,把一张由乐高积木紧密拼成的底板,拆分成一块块独立的积木块。工具会为每一个非NoData(空值)的像元生成一个点,这个点就位于该像元几何中心的位置。在输出的点要素属性表里,默认会有一个叫做GRID_CODE的字段,这个字段里存放的,就是该像元原始的像元值(DN值)——可能是温度、高程、浓度,或者任何你的栅格所代表的数值。
第二个工位:添加XY坐标字段。栅格转点后,我们得到了带着值的点,但这些点还没有明确的坐标信息。虽然它们在地图上有位置,但属性表里并没有数字化的X和Y字段。这一步,我们需要手动在属性表里创建两个新字段,通常命名为“X”和“Y”,然后利用ArcGIS的“计算几何”功能,批量将每个点的坐标值(可以是经纬度,也可以是投影坐标系下的X,Y米制单位)计算出来并填入。这是将空间位置信息表格化的关键。
第三个工位:表转Excel(Table to Excel)。当点要素的属性表里同时拥有了GRID_CODE(值)和“X”、“Y”(坐标)字段后,这张属性表本身就是一张完美的数据表了。最后一步,就是使用“表转Excel”工具,将这张表从地理数据库的格式,导出为标准、通用的.xlsx或.xls文件,方便在Excel、SPSS、R或Python中进行下一步的数据处理和分析。
理解了这个工具链,我们就知道,整个过程其实是一个“空间数据属性化、表格化”的思路。接下来,我们进入实战环节,我会用一个真实的DEM(数字高程模型)数据作为例子,带你一步步走通,并分享每个环节的细节和注意事项。
3. 第一步:从栅格到点——生成携带像元值的矢量点
好,现在我们打开ArcGIS Pro或者ArcMap,加载你需要处理的栅格数据。我这里用的是一个dem.tif文件,它表示一块区域的地形高程。
3万+




被折叠的 条评论
为什么被折叠?



