ArcGIS按照格网提取矢量值

### 如何使用 ArcGIS 按照栅格网格提取对应的矢量属性值 #### 创建与栅格相同大小的网格并赋予唯一标识 为了确保能够准确对应到每一个栅格单元，可以先基于原始栅格创建具有唯一标识符的点集。这一步骤通过将栅格转换成点来实现，并在这些点上附加一个名为`value`的新字段，其值设为各点自身的FID编号[^1]。 ```sql Add Field: value (Type: Long Integer) Calculate Field: value = FID ``` #### 设置处理环境参数 接着，在执行任何地理处理操作前，需配置好ArcGIS中的处理范围(Processing Extent)，使其匹配目标栅格的数据边界；同时启用捕捉(Snap Raster)选项指向同一源栅格文件，从而保证新生成的对象能严格遵循原有栅格的空间布局。 #### 将点数据重构成栅格形式 利用先前准备好的带`value`字段的点层作为输入，采用Spatial Analyst扩展模块下的Point to Raster工具完成向栅格格式的转变过程。这里指定`value`作为分类依据，使得最终输出的每个像元都能继承自唯一的起始位置信息。 ```python import arcpy from arcpy import env from arcpy.sa import * env.workspace = "C:/data" arcpy.CheckOutExtension("spatial") outRaster = PointToRaster_conversion(inPoints="point_layer", valueField="value", cellAssignment="MOST_FREQUENT", priorityField="NONE", cellsize=your_cell_size, outRaster="output_raster") ``` #### 实现栅格至多边形的转换 当获得上述经过特殊编码后的栅格后，下一步就是调用Raster To Polygon功能将其转变为面状要素类。此过程中会自动保留原有的数值型属性，即之前的`value`，进而便于后续与其他空间对象进行关联分析或进一步加工。 ```python inRaster = "path_to_your_prepared_raster" outPolygonFeatures = "polygon_output" field = "" # 如果不需要额外字段可留空 simplifyPolygons = "NO_SIMPLIFY" # 或者"SIMPLIFY"取决于需求 arcpy.RasterToPolygon_conversion(inRaster, outPolygonFeatures, simplifyPolygons, field) ``` #### 提取矢量属性并与栅格关联 对于已经存在的矢量数据（比如由DEM产生的等高线），如果希望获取它们覆盖区域内的平均高度或者其他统计特征，则可以通过Zonal Statistics as Table工具达成目的。它接受两个主要参数——一个是定义分区界限的区层面（Zone Layer），另一个是用来计算统计数据的目标栅格(Raster Dataset)[^2]。 ```python zoneLayer = "zones.shp" # 这里应该是前面得到的多边形 rasterDataset = "elevation.tif" # 原始的高度模型或其他类型的栅格 statisticsTable = "stats.dbf" arcpy.gp.ZonalStatisticsAsTable_sa(zoneLayer, "FID", rasterDataset, statisticsTable, "DATA", "ALL") ``` 以上流程展示了如何借助ArcGIS平台内建的功能和服务，有效地从给定的栅格资源中抽取出结构化的矢量表达及其附带的信息。