arcgis栅格数据重分类后有八类，进行栅格转面后转出的面要素只有一个分类后获得的grid_code，是什么原因？

Question

想通过计算核密度分析后每个区域的质心坐标，所以要将每个类别形成面。如题，
arcgis的栅格数据进行重分类后，有八类：

但进行栅格转面后转出的面要素只有一个分类“8”获得了grid_code：

请问这是什么原因？查遍全网都没有这种情况

Answer 1

引自免费微信小程序：皆我百晓生

遇到在ArcGIS中进行栅格数据重分类并转换为面要素时，只有一类被成功转换的情况，这可能由几个因素导致。以下是一些可能的原因及解决策略：

1. 栅格数据属性问题

• 确保你的栅格数据在重分类之后确实包含了所有八类值。可以通过查看属性表或使用Identify工具检查不同位置的值。

2. 栅格转面参数设置

• 在执行栅格转面操作时，确保你选择了正确的输出类型和字段映射。特别是grid_code字段应该正确地从输入栅格中的分类值映射到输出的面要素类。

3. 空值处理

• 检查是否有空值或无效值存在。有时，栅格中某些区域可能未被赋值，这可能导致转换时仅识别到部分分类。使用Set Null工具可以将这些值设为一个特定值，确保所有区域都被正确分类。

4. 数据类型不匹配

• 确认栅格数据和输出面要素类的字段类型是否兼容。例如，如果栅格数据中的值是整数类型，而输出面要素类中的grid_code字段是浮点型，可能会导致数据丢失或转换错误。

解决方案步骤：

步骤1: 验证栅格数据

import arcpy
from arcpy.sa import *

# 输入栅格路径
inRaster = "C:/path/to/your/raster"

# 查看栅格统计信息
desc = arcpy.Describe(inRaster)
print("Raster Name:", desc.name)
print("Raster Type:", desc.dataType)

# 使用Identify工具检查栅格值
xy = (100, 200) # 示例坐标
value = arcpy.GetCellValue_management(inRaster, str(xy))
print("Value at", xy, "is", value)

步骤2: 栅格转面

# 输出矢量文件路径
outPoly = "C:/path/to/output/polygon.shp"

# 设置栅格转面参数
arcpy.RasterToPolygon_conversion(inRaster, outPoly, "NO_SIMPLIFY", "VALUE")

步骤3: 检查输出结果

# 检查输出面要素类属性
outPolyDesc = arcpy.Describe(outPoly)
print("Output Polygon Name:", outPolyDesc.name)
print("Output Polygon Type:", outPolyDesc.dataType)

# 打开属性表检查grid_code字段
fields = ["grid_code"]
with arcpy.da.SearchCursor(outPoly, fields) as cursor:
    for row in cursor:
        print("Grid Code:", row[0])

通过上述步骤，你应该能更深入地了解问题所在，并采取相应的措施解决它。如果问题仍然存在，可能需要进一步检查数据源、转换过程中的具体参数设置，或是尝试不同的转换方法。