ArcGIS教程：编辑特征

　　摘要

　　通过合并、又一次编号和删除类特征来编辑和更新特征文件。

　　使用方法

　　· 编辑特征工具同意您通过下面全部操作或某一操作来改动现有特征文件：

　　合并一组特征类

　　又一次编号特征类 ID

　　删除不须要的特征

　　· 输入特征文件必须是 ASCII 特征文件。该文件能够是不论什么用于生成包括所需统计信息的文件的“多元分析”工具(如 Iso聚类和创建特征)的输出。

该文件必须至少包括两个类。可通过 .gsg 扩展名来识别此类文件。

　　· 输入特征重映射文件是 ASCII 文件，其每一行有两列值与之相应，以冒号分隔。第一列是原始类 ID 值。

第二列包括用于在特征文件里更新的新类 ID。文件里的全部条目必须基于第一列以升序进行排序。

　　要合并一组类。请为该组各类 ID 的第二个值放置同样的新类 ID。仅仅有须要编辑的类才必须被放入特征重映射文件;不论什么在重映射文件里不存在的类将保持不变。要删除类特征。请使用 -9999 作为该类第二列的值。此外，还可将类 ID 又一次编号为某个不存在于输入特征文件里的值。下面是输入特征重映射文件的演示样例：

　　2 : 3

　　4 : 11

　　5 : -9999

　　9 : 3

　　上例将使用 3 合并类 2 和 类 9。使用 11 合并类 4，并将删除类 5。

　　· 假设输入特征文件里包括类特征名称。而且输入特征重映射文件里的特征要进行合并，则合并值的关联名称将被传输至输出特征文件。

　　· 假设输入为创建自多波段栅格(超过三个波段)的图层，则操作将会考虑与源数据集相关联的全部波段。而不不过由图层载入(符号化)的三个波段。

　　· 您能够使用多种方式将多波段栅格的波段子集指定为工具的输入。

　　·  假设要使用工具对话框，则通过输入栅格波段旁边的button导航至多波段栅格。打开栅格并选择所需波段。

　　·  假设多波段栅格是内容列表中的图层，则能够使用创建栅格图层工具来创建一个新的仅仅包括所需波段的多波段图层。

　　·  您还能够用波段合成来创建一个仅仅包括所需波段的新数据集，并将生成的数据集用作工具的输入。

　　·  在 Python 中，能够在工具參数中以列表形式直接指定所需波段。

　　语法

　　EditSignatures (in_raster_bands, in_signature_file, in_signature_remap_file, out_signature_file, {sample_interval})

　　代码实例

　　EditSignatures 演示样例 1(Python 窗体)

　　本例将基于输入重映射文件编辑特征文件。

　　import arcpy

　　from arcpy import env

　　from arcpy.sa import *

　　env.workspace = "C:/sapyexamples/data"

　　EditSignatures("redl123", "c:/sapyexamples/data/zsamp12.gsg",

　　"c:/sapyexamples/data/zsamp7.rmp",

　　"c:/sapyexamples/output/redlremap.gsg", "")

　　EditSignatures 演示样例 2(独立脚本)

　　本例将基于输入重映射文件编辑特征文件。

　　# Name: EditSignatures_Ex_02.py

　　# Description: Edits and updates a signature file by merging, renumbering,

　　# and deleting class signatures.

　　# Requirements: Spatial Analyst Extension

　　# Import system modules

　　import arcpy

　　from arcpy import env

　　from arcpy.sa import *

　　# Set environment settings

　　env.workspace = "C:/sapyexamples/data"

　　# Set local variables

　　inRaster = "redl123"

　　oldSig = "c:/sapyexamples/data/zsamp12.gsg"

　　sigRemap = "c:/sapyexamples/data/zsamp7.rmp"

　　outNewSig = "c:/sapyexamples/output/redlsig.gsg"

　　interval = ""

　　# Check out the ArcGIS Spatial Analyst extension license

　　arcpy.CheckOutExtension("Spatial")

　　# Execute EditSignatures

　　EditSignatures(inRaster, oldSig, sigRemap, outNewSig, interval)