weixin_42187487的博客

使用该功能来恢复已经被存为文本文件的显示组。参见第四章的“恢复保存的显示组”以获得详细介 绍。 （5） 关闭所有文件 在 ENVI中完成对文件的操作后，通过选择 File > Close All Files来关闭它们。 注意：ENVI 将关闭所有磁盘文件并删除所有内存中的数据项。在选择该按钮前，确认你已保存了内存中 所有重要的数据项。 将出现一个要求你确认的警告消息。点击“OK”。 2.17 ENVI队列管理 使用 ENVI Queue Manager 选项可以执行被排序的操作程序（参见第 18 页的“将未执行功能列入队 列”）。 注意：一个排列程序的结果不能用作另一个程序的输入。 先要确定运行队列功能所需的所有文件都已经被打开并列在可用波段列表中。在 ENVI主菜单中，选 择 File >ENVI Queue Manager，将出现 ENVI Queue Manager对话框，其中包含所有队列的操作程序。 要浏览一个队列功能的输入文件名，点击该功能名称。关于空间、波谱维数和所用子集的信息将显示 在“Procedure Information”区域中。要删除一个队列的操作程序，点击该功能名称，然后点击“Remove Selected”按钮 图 2-20：ENVI Queue Manager对话框 在“Queued Procedures List”中，点击功能名称。点击“Execute Selected”按钮。被执行功能的名称

（4） 跨平台文件的可移植性 当使用 ENVI命名文件时，要考虑的一个附加因素就是跨平台的可移植性。在 UNIX系统上的文件名 后附加一个 .hdr。在Windows系统上，用.hdr 代替当前的扩展名。这对于保持图像文件和它们相应的（.hdr） 头文件之间的关系是特别重要的 。例如，一个名为 image_1.img 的 UNIX图像文件将 image_1.img.hdr作 为它相应的头文件。一个名为 image_1.img的 Windows图像文件将 image_1.hdr作为它相应的头文件。在 Windows下，两幅图像 image_1.img和 image_1.dat将有同样的头文件名 image_1.hdr，若这两幅图像大小 和特征不同，那么 ENVI将会出现问题。在 UNIX下，同样的两个文件将有相应的头文件名 image_1.img.hdr 和 image_1.dat.hdr，因此不会发生混乱。若图像的跨平台可移植性存在问题，那么最简单的解决方案就是 不用扩展名来命名图像文件。这样，图像文件 image_1 在 UNIX 和 Windows 下有同样的头文件名 （image_1.hdr）。此外，将图像从 UNIX系统移到Windows系统前，可以按Windows约定对图像和头文件 进行重命名。 1.11 ENVI支持的输入文件格式 ENVI 支持的输入文件格式如下所示。有关这些格式的更多信息在第二章的“打开图像文件”，“打开 外部文件”，“打开矢量文件”中可以找到。ENVI支持多种输入文件格式： 二进制文件： ● BSQ（波段顺序格式）

13.8 纹理滤波器 ENVI包括几个纹理滤波器，它们允许从 SA 或其它数据类型中抽取纹理信息。这些滤波器是根据数 据范围、RMS、数据的一阶矩和二阶矩进行滤波的。 要选择一个滤波器，点击Radar > Texture > 所需滤波器。这些滤波器也可以从ENVI主菜单中的Filters 菜单里调用，详细介绍，请参阅第515页的“纹理滤波器”。 13.9 合成彩色图像 使用Synthetic Color Image选项可以将一幅灰阶图像转换成一幅彩色合成图像。该转换通常用于在保 留有用细节的情况下，增强雷达数据中大比例尺细微特征的显示。 要生成一幅彩色合成图像，选择Radar > Synthetic Color Image。 注意：详细介绍，请参阅第515页的“合成彩色图像”。 13.10 极化工具 ENVI提供一套极化雷达分析工具，用于对NASA/JPL、AIRSAR、TOPSAR和SIR-C类型的数据进行 分析。 由于这些极化雷达数据被存储为一种压缩的Stokes格式或散射矩阵格式（scattering matrix format），所 以不能直接对它们进行浏览，而需要专用工具。 ENVI提供的极化工具允许进行如下操作，包括：浏览头文件、解压和合成图像；多视(multilook) SIR-C 压缩数据；计算相位和消隐脉冲高度图像；运行AIRSAR散射分类技术；以及提取极化信号。 ‧ JPL AIRSAR数据合成 在AIRSAR图像用于ENVI标准处理程序之前，必须先被合成。 使用Synthesize AIRSAR Data选项可以从压缩的Stokes文件中合成标准的和特定的发射和接收极化图 像以及总功率图像。输入的AIRSAR数据必须处于JPL stokes矩阵格式（标准文件扩展名为 .stk）。 JPL AIRSAR数据通常使用压缩的Stokes矩阵格式，作为三个独立的文件（P、L、C-波段）在9-轨、 8mm或4mm 磁带上传输。 详细介绍，请参阅JPL AIRSAR文件描述。 van Zyl, J.J., H.A. Zebker, and C. Elachi: Imaging radar polarization signatures: theory and observation. Radio Science 22(4):529-543, 1987.

（2） 构建新的决策树 在ENVI主菜单中，选择Classification >Decision Tree >Build New Decision Tree。将出现ENVI Decision Tree窗口，其中包含一个决策树节点和两个类别（分支）。 点击“Node”按钮，并在Edit Decision Tree Properties对话框中输入节点名和表达式（参见第315页 的“决策树变量和表达式”）。点击“OK”，将出现Variable/File Pairings对话框。下图显示了Edit Decision Tree Properties对话框和Variable/File Pairings对话框。 图 6-11：Edit Decision Tree Properties对话框（左）和 Variable/File Pairings对话框（右） 点击变量名，选择与之对应的输入文件或波段。要添加新的节点，在“Class”按钮上点击鼠标右键， 选择“Add Children”。点击新的节点并在Edit Decision Tree Properties对话框中输入一个表达式，然后点击 “OK”。如果需要，重复步骤4和5，添加足够多的节点。决策树构建完毕后，选择Options > Execute。 （3） 运行决策树 要交互式地运行和编辑决策树：在ENVI Decision Tree窗口中，选择Options > Execute。当出现Decision Tree Execution Parameters对话框时，若需要，选取一个空间子集，并输入一个输出文件名。如果数据经过 地理坐标定位，选择哪幅图像将作为基图像，键入输出像元尺寸，并在下拉菜单中选择重采样方法。 在决策树运行时，每个节点都将变为绿色。分类结果出现在一个新的显示窗口中。系统将为每个节点 保存一个临时文件，避免在决策树每次运行时都要重复计算结果。

（9） IAR反射率定标 IAR (Internal Average Relative) Reflectance定标工具用于将图像标准化成一景平均波谱。这对于在没有 地面测量的鲜为人知的区域，将高光谱数据减小到“相对反射率”非常有效。该工具特别适用于没有植被 的干旱区域。从整个景区中计算出平均波谱，并将它作为参照波谱，然后被分成图像中每一个像元的波谱。 选择Basic Tools >Preprocessing > Calibration Utilities > IAR Reflectance。当出现Calibration Input File 对话框时，选择输入文件，若需要，构建任意子集。点击“OK”。当出现IARR Calibration Parameters对话 框时，选择输出到“File”或“Memory”。点击“OK”。如果不存在统计文件，当统计信息被计算时，将 出现一个运行状态窗口。定标完成后，将出现另外一个状态窗口。定标完毕后，定标后的数据将被加载到 可用波段列表中。 （10） 经验线性定标(Empirical Line Calibration) Empirical Line 定标技术被用来将波谱数据与所选区域的反射波谱相匹配。线性回归被用在每一个波 段中使DN值和反射率相等，这等价于消除了太阳幅亮度和大气程辐射。下面的等式介绍了经验线性定标 中的偏移是如何计算的。 反射率（区域波谱）= 增益 * 幅照度（输入数据）+ 偏移 ENVI经验线性定标至少需要一个区域、波谱库或其它参照波谱，这些可以来自波谱剖面或波谱图、 波谱库、感兴趣区、统计表和ASCII文件。输入的波谱将自动被重采样，以与选择的数据波长相匹配。如 果使用了多个波谱，则对每个波段的回归将通过穿越所有波谱的回归行来计算。如果只用到一个波谱，则 假定回归行只穿越初始区域（0反射率等于0 DN）。可以用已经存在的系数对数据集进行定标。 ‧ 计算系数并定标 在应用经验行定标时，可以选择一个暗区和一个亮区应用（假定这些区域中的参照波谱是可以获得 的），这提供了一个更加精确的线性回归。使用尽可能多的数据/区域波谱也可以提高定标精度，至少需要 一组数据/区域波谱。 注意：要使用来自感兴趣区的波谱，在运行本功能前先定义一个感兴趣区。 选择Basic Tools >Preprocessing > Calibration Utilities > Empirical Line > Compute Factors and Calibrate。 当出现Empirical Line Input File对话框时，选择输入文件。若需要，构建任意子集。点击“OK”，将出现 Empirical Line Spectra对话框。

（9） 散点图管理 使用散点图菜单栏上的File下拉菜单，可以进行二维散点图的输入和输出，选项包括输出散点图、导 入保存的感兴趣区，以及退出散点图功能。 ‧ 启动新的散点图 要启动使用不同波段的新散点图：从下列选项中选择： 在散点图窗口，选择File >New Scatter Plot。在主图像窗口，选择Tools >2D Scatter Plot当出现Scatter Plot Band Choice对话框时，选择X和Y波段。点击“OK”。除了最初的散点图之外，将出现一个新的散 点图。 ‧ 散点图的输出 可以将散点图存为脚本（PostScript）文件或图像文件（RGB 二进制文件）：BMP、HDF、JPEG、PICT、 SRF、TIFF或XWD文件。要存储一个散点图，选择File >Save Plot As >Postscript或Image File。 详细介绍请参阅第225页的“显示输出选项”。 ‧ 将 ROIs导入到散点图 要把原先保存的感兴趣区导入到散点图：选择File > Import ROIs。当出现Import ROIs to Scatter Plot 对话框时，从现有的感兴趣区列表中选择。被选择的感兴趣区将被绘制在图像窗口中，相应的点突出显示 在散点图中。在同一散点图区域内的其它像元将突出显示在图像窗口中。 4.34 3-D曲面浏览 通过使用ENVI的3-D SurfaceView TM 功能，可以对3-D数据进行曲面飞行。

4.10 交互式密度分割 使用 Density Slice功能可以为灰阶图像中突出显示的区域选择数据范围和颜色。用于密度分割着色的 数据范围可以来自显示图像或其他相同尺寸的图像。 4.11 网格线 Grid Lines功能允许在一幅图像上覆盖一个或多个网格。网格可以基于像元或地图坐标和/或基于经纬 度（对于地理坐标定位的图像）。每个图像显示可以都有它自己的一套网格，显示在主图像、滚动和缩放 窗口。 提示：使用 ENVI快速制图功能可以在地理坐标定位的图像上快速叠加网格线、标题、磁偏角图、指北针 和边界。详细介绍请参阅第 159页的“快速制图”。 图 4-18：基于像元的网格图像

（7） 关闭所有显示窗口 要关闭所有显示窗口，从显示窗口的菜单或 ENVI主菜单中，选择Window >Close All Display Windows。 （8） 关闭所有图表窗口 要关闭所有图表窗口，选择Window >Close All Plot Windows。

（2） 选择注记类型 Annotation对话框中的 Object 下拉菜单选项允许选择多种注记类型，包括文本、符号、形状、地图信 息和彩色阶跃斜面（color ramp）。每种类型都包含各自的注记选项，更多细节描述如下。 • 文本注记 在 Annotation对话框中，选择 Object > Text。通过在窗口底部的文本框内键入，来输入文本。一次可 以输入多行文本。通过选择以下特定选项来更改文本特征。 注意：在对话框中的文本特征不会改变，要显示改变，在图像窗口上点击即可。注意在对显示的文本特征 满意之前不要确认注记的放置位置。 要更改文本颜色，点击“Color”按钮菜单，选择所需颜色。要把文本放置在实心的彩色背景上，使用 “Back”按钮菜单选择一种颜色；要把文本直接放置在图像上而不使用背景（系统默认），选择“Off”。 要调整文本宽度，在“Thick”文本框中输入一个新宽度，或使用增减箭头来更改“Thick”参数框内的值。 要更改文本字体，点击“Font”按钮菜单，然后选择一种字体。要更改字体的大小（按像元），在“Size” 参数框中更改数值，或使用增减箭头按钮。要更改屏幕上的文本方向（书写的角度），在“Orien”参数框 中输入角度的度数值（相对于水平线（0 度）逆时针方向），或使用增减箭头按钮。要居中或者左对齐、 右对齐文本，在“Align”按钮菜单下选择适当的选项。 按照第 120页“注记定位”中所描述的方法来放置和确认注记。 注意（关于字体）：在文本注记窗口的“Font”按钮菜单中包含多种字体类型，也包含特定的 ENVI字体。 使用“ENVI Symbols”字体类型来放置制图符号和 ENVI、IDL的标识语。 在创建注记时，ENVI 提供了很强的灵活性，允许使用基于文本或符号的“TrueType”字体。如果下 载了该字体，应该安装它，以便在 ENVI中使用。 为 ENVI安装 TrueType字体 IDL和 ENVI通过源文件（ttfont.map文件，位于 IDL安装目录下的 resource/fonts/tt子目录下）提供了 TrueType字体的接口。ttfont.map文件包含了一系列可用的字体文件。下图是一个字体源文件的示例。

3.10 ENVI矢量窗口 为了显示矢量数据及合成简单的矢量地图，ENVI提供了一个独立于（stand-alone）操作系统的 GIS 图 表窗口。ENVI 也在标准 ENVI 显示上提供矢量覆盖，包括在所有窗口内（包括缩放窗口）覆盖图的真矢 量化（true vectorization）。矢量数据的精确度被保留以避免 “Pixellation”。矢量可以从多种输入文件中读 取，也可以在一个矢量窗口或在图像上创建和绘制新的矢量层。当交互地跟踪矢量时，能够显示纬度/经度 和地图坐标信息。当光标跟踪每个矢量时，一个 Vector Information窗口可以实时地显示属性信息。你可以 直接查询矢量的 GIS 属性信息，生成具有已选定属性信息的新层。由 ENVI 内部的 .evf 格式，创建 ArcView Shape文件以及相应的 .dbf 属性文件和索引，或 DXF文件。用 ENVI强大的图像处理能力生成 的新矢量层，以及 ENVI中矢量层的变化，很容易导出工业标准的 GIS格式。 ENVI的矢量窗口提供了浏览下列数据的一种方式，包括：USGS Digital Line Graphs（DLG）、空间数 据转换标准格式（SDTS）的 USGS DLGs、DXF文件、ARC/INFO 交换格式文件及 ArcView Shape文件等。 你可以浏览、编辑和查询与 ArcView Shape文件相关的属性，也可以创建你自己的矢量文件和属性。ENVI 的矢量窗口包括：一个矢量显示窗口或图像窗口，及一个#N Vector Parameters对话框，该对话框控制矢量 的显示以及与矢量属性的交互。

2.16 处理文件和脚本文件 （1） 更改输出目录 要更改输出的缺省目录：选择 File > Change Output Directory。当出现“Change Output Directory”对话 框，键入你想让 ENVI输出文件到某一目录的完整路径。点击“OK”来把当前输出目录更改为新的路径。 （2） 将Session存为脚本文件（Save Session to Script） 该选项允许你把当前打开的图像文件、图像显示和它们的波段保存到 ENVI 的启动脚本文件。要执行 这个启动文件，请参阅本页的“Execute Startup Script”。 选择 File > Save Session to Script。当出现 Output ENVI Script Filename对话框时，输入启动文件名（通 常扩展名为.ini）。 （3） 执行启动脚本文件（Execute Startup Script） ENVI的 Execute Startup Script 允许你当启动 ENVI 时，打开图像文件、把波段导入显示、打开矢量 文件，以及打开感兴趣区文件。如果一个文件名在 envi.cfg 文件中被设置，当 ENVI 启动时该脚本被自 动地执行。通过使用该功能，它也可以在任何时候被执行。附录 A中的 “ENVI Startup Script commands” 提供了 ENVI 启动脚本命令的说明。

（1） 在UNIX下启动ENVI 要启动 ENVI，在 UNIX命令行中键入“envi”；要启动 ENVI RT，在 UNIX命令行中键入“envi_rt”。 （2） 从Macintosh系统启动ENVI 指定使用 OroborOSX安装器。 双击 OroborOSX图标。OroborOSX发出 Xdarwin并在一个 OS X窗口显示一个 UNIX X-Windows命令 行。 操作下列步骤之一：要启动 ENVI，在 OroborOSX 命令提示行中键入“envi”；要启动 ENVI RT，在 OroborOSX命令提示行中键入“envi_rt”。 （3） 从Windows系统启动ENVI 选择 Start >Programs >ENVI 4.2 >ENVI。程序成功加载和运行 ENVI以后，ENVI主菜单将出现。 （4） ENVI启动的命令行选择 你也可以改变 ENVI和 IDL的启动方法，通过设置一个相应的命令行开关来启动 ENVI。 （5） ENVI启动的UNIX开关 可以在 UNIX平台使用如下命令来启动 ENVI： 以 32-bit模式启动 IDL，否则对于那些支持 64-bit的平台，将以默认的 64-bit模式 启动 IDL。如果你没有安装 64-bit译本，32-bit译本将被自动打开。如果你没有安装 32-bit 译本，那么这个“flag”将无法工作。例如：%envi –32。 1.5 ENVI图形用户界面（GUI）

（2） 打开SPOT文件 根据如下步骤来打开 old SPOT文件、GeoSPOT 文件、ACRES SPOT、SPOT植被文件和 SPOT DIMAP 文件。 • 打开 SPOT文件 要直接从光驱或磁盘读入 1A、2A和 1B的 SPOT数据：选择 File > Open External File >Spot > SPOT。 当出现文件选择对话框时，选择要读取的图像文件。ENVI 自动地提取头文件信息，并把图像波段输入到 可用波段列表。 • 打开 GeoSPOT - ArcView Raster（.bil）文件 SPOT 提供的数据大多为 GeoSPOT格式。ArcView 栅格图像文件有一个相似的格式规范。在从 SPOT 图像得到的文档中有 GeoSPOT 格式的详细描述。虽然 GeoSPOT 格式提供很多种栅格和矢量数据，但是 目前 ENVI 只支持 GeoSPOT 栅格图像。这些图像的文件扩展名为 .bil，并且有一个相应的扩展名为 .hdr 的头文件。ENVI 将从头文件中读取 UTM 和 State Plane projection georeferencing 信息。 选择 File > Open External File > Spot > GeoSPOT。当出现文件选择对话框时，选择要读取的文件。 ENVI 自动地提取头文件信息，包括 UTM 和 State Plane projection georeferencing 信息，并把图像波段输 入到可用波段列表。 • 打开 ACRES SPOT文件 要打开澳大利亚遥感中心(ACRES) CCRS和 SPIM格式的 Spot文件：选择 File > Open External File > Spot > ACRES SPOT。当出现文件选择对话框时，选择 Imag_xx .dat文件。ENVI 自动地提取头文件信息，

图5．6基于单幅图像建模的流程图 奉文的重建对象是剑桥大学的一栋建筑物，无法通过测量它的实际尺寸验证本文 方法的可行性。为此对实验室中一个简单的长方体茶盒进行重建，图像实际像素大小为 640x481，如图5 7所示，选择相应的世界坐标系后对相关数据进行具体分析。 曲原困 M设定世界坐标系的位置 圈5 7原阁lea．b佩 对图像tea．box进行特征检测，如图5 8所示，根据设定的世界坐标系，分别在三 个主方向选取平行线的投影线段确定灭点，选取的线段坐标见表5 6。 a)canny算子检测边缘 W概率Hough变换检测直线 圉5．8 tea-box特征检测结果

1．2．2国内外研究现状 近些年，单幅图像建模的研究越来越受到关注，由于一幅图像无法完全重构出对 应的三维模型，针对这一点，众多研究者利用图像中己知的几何信息对目标物体进行三 维重建【121。 在国外，Caprile和Torre首先提出了利用灭点进行摄像机标定的思想【I31。在给出 并证明了灭点属性的基础上，根据灭点与摄像机内参数的关系通过单幅图像确定摄像机 的内参数，再由根据特定模板计算摄像机的外参数。Guillou等人提出了一种利用灭点 实现单幅图像中三维物体重建的方法【14】，假设相机在投影面上的投影中心即为图像的中 心，而不是利用灭点的属性进行计算，且在对相机内参数的估算时也没有充分利用灭点 的属性，因而计算过程稍显复杂【15】。Horry等人在1997年提出的基于单幅图像的场景 漫游方法(Tour into the picture)D6】。Shum等人利用几何约束条件重建三维模型【17】。 Debevec等将目标物体的几何形状作为相应参数，实现结构化场景的交互式建模【l 81。 ZHANG[19】提出利用一些用户指定的约束，如目标物体表面的位置、轮廓、皱痕等生成 满足这些约束的三维表面。 在国内，对利用单幅图像进建模工作的研究也很多，张祖勋等通过分析灭点几何， 2

表 4.29  带符号饱和运算的示例运算结果  输入(R0)  输出(R1)  Q标志位  0x2000(8192)  0x7FF(2047)  1  0x537(1335)  0x537(1335)  无变化  0x7FF(2047)  0x7FF(2047)  无变化  0  0  无变化  0xFFFFE000(‐8192)  0xFFFFF800(‐2048)  1  0xFFFFFB32(‐1230)  0xFFFFFB32(‐1230)  无变化  如果需要把 32 位整数饱和到无符号的 12 位整数（0‐4095），则可以如下使用 USAT指令    USAT{.W} R1, #12, R0 该指令的执行情况如图 4.4演示  图 4.4  无符号饱和运算  表 4.30  无符号饱和运算的示例运算结  输入(R0)  输出(R1)  Q标志位  0x2000(8192)  0xFFF(4095)  1  0xFFF(4095)  0xFFF(4095)  无变化  0x1000(4096)  0xFFF(4095)  1  0x800(2048)  0x800(2048)  无变化  0  0  无变化  0x80000000(‐2G)  0  1  0xFFFFFB32(‐1230)  0  1  CM3 中的一些有用的新指令 75   

图9.2 异常咬尾示意图 为进一步帮助读者理解，译者从另外文献上截取并改编下图：  图9.2B 异常咬尾与常规处理的比较（以ARM7TDMI为例）  晚到（的高优先级）异常   CM3的中断处理还有另一个机制，它强调了优先级的作用，这就是“晚到的异常处理”。 当CM3对某异常的响应序列还处在早期：入栈的阶段，尚未执行其服务例程时，如果此时收 到了高优先级异常的请求，则本次入栈就成了为高优先级中断所做的了——入栈后，将执行 高优先级异常的服务例程。可见，它虽然来晚了，却还是因优先级高而受到偏袒，低优先级 的异常为它“火中取栗”。    比如，若在响应某低优先级异常#1的早期，检测到了高优先级异常#2，则只要#2没有太 晚，就能以“晚到中断”的方式处理——在入栈完毕后执行ISR #2，如图9.3所示。如果异常 #2来得太晚，以至于已经执行了ISR #1的指令，则按普通的抢占处理，这会需要更多的处理 器时间和额外32字节的堆栈空间。    在ISR #2执行完毕后，则以刚刚讲过的“咬尾中断”方式，来启动ISR #1的执行。  138   

2 分分分类类类型型型数数数据据据聚聚聚类类类 2.1 相相相关关关定定定义义义 设属性集A = {A1, A2, · · · , Am},属性值域集合 为D = {D1, D2, · · · , Dm},则属性Ai的值域为Di = {vi,1, vi,2, · · · , vi,pi}, i = {1, 2, · · · ,m}. 其中: vi,r为 属性Ai的第 r个值,即属性Ai的一个域值, r = {1, 2, · · · , pi}, pi为属性Ai的不同取值数; Di中取值的次 序只是为了计算方便,没有其他含义. 设数据集X = {x1, x2, · · · , xn},其中第 i个样本 的m个属性值表示为xi = [xi1, xi2, · · · , xim]T, xij ∈ Dj , j = {1, 2, · · · ,m}. 定义 1 设一对象x, 若其在属性Ai取 vi,r的情 况下,属性Aj取 vj,θ(vj,θ表示属性Aj的第 θ个取值), 则称域值 vi,r和 vj,θ在对象x中共生, 此时将形如 (vi,r, vj,θ)的形式称为值对. 定义 2 设值对 (vi,r, vj,θ)在数据集中出现的次 数为 tmvj,θvi,r , θ ∈ {1, 2, · · · , nAj},简记为 tm. 设Ai的两个取值为 vi,rx, vi,ry, 且有 vi,rx 6=vi,ry, rx, ry ∈ {1, 2, · · · , pi},统计出 vi,rx和 vi,ry与Dj中各 元素构成的值对出现的次数. 设 vi,rx和 vi,ry与Dj中 各元素均能构成值对, 且各 tm均一致. 假设现在增 加了Ak, 分别统计出 vi,rx和 vi,ry与Dj中各元素构 成的值对的 tm.一方面,若 vi,rx和 vi,ry均能与 vk,α(α = {1, 2, · · · , pk})构成值对且两者的 tm相差较小, 则 vi,rx和 vi,ry (在Ai, Aj和Ak条件下)的差异性较小, 反之较大; 另一方面, 若 vi,rx与Dk中的部分值构成 值对, 则 vi,ry与Dk中剩余的元素构成值对, 分析结 果同上. 定义 3 假如属性Ai的第 r个值为 vi,r, 计算其 与剩余m − 1个属性的各个值共同出现的次数,以此 构造一个行向量为 V Aj∈A,j 6=i vi,r =(tm v1,1 i,r , tm v1,2 i,r , · · · , tm v1,p1 i,r , · · · , tmvi−1,1i,r , · · · , tm vi−1,p(i−1) i,r , · · · , tmvi+1,1i,r , · · · , tm vi+1,p(i+1) i,r , · · · , tmvm,1i,r , · · · , tm vm,pm i,r ), (1) 将该行向量简记为VIJ. 定义 4 将同一个属性的两个不同值之间的距 离简记为值距离. 本文采用余弦相似性来度量两个行向量, 即属 性Ai的第 p个值与第 q个值之间的值距离为 σAi(vi,p, vi,q) = 1− V Aj∈A,j 6=i vi,p V Aj∈A,j 6=i vi,q ‖V Aj∈A,j 6=ivi,p ‖‖V Aj∈A,j 6=ivi,q ‖ , (2) 其中 ‖V Aj∈A,j 6=ivi,r ‖ = √√√√ m∑ j=1,j 6=i pj∑ θ=1 (tmvj,θvi,r ) 2. 若对象x和 y在属性Ai(i = 1, 2, · · · ,m)下的取 值分别为 vi,p和 vi,q,则x和 y的距离为 d(x, y) = m∑ i=1 σAi(vi,p, vi,q). (3) 定义 5 对于每个类,以该类中所拥有的各个属 性的各个属性值在该类中所占的比例来代表该类[9], 记第 i类的类中心为 ci = 1 |Ci| ((n1,1,c, n1,2,c, · · · , n1,p1,c), (n2,1,c, n2,2,c, · · · , n2,p2,c), · · · , (nm,1,c, nm,2,c, · · · , nm,pm,c)), (4) 其中nu,v,c表示在类Ci中第u个属性的第 v个取值 所占的个数. 那么对象x与类中心 ci的距离为 d(x, ci) = 1 |Ci| (n1,1,cσA1(v1,1, xj) + n1,2,cσA1(v1,2, xj)+ · · ·+ n1,p1,cσA1(v1,p1, xj) + · · ·+ ni−1,1,cσAi−1(vi−1,1, xj) + · · ·+ ni−1,p(i−1),cσAi−1(vi−1,p(i−1), xj) + · · ·+ ni+1,1,cσAi+1(vi+1,1, xj) + · · ·+ ni+1,p(i+1),cσAi+1(vi+1,p(i+1), xj) + · · ·+ (nm,1,cσAm(v1,1, xj) + · · ·+ (nm,p,mcσAm(vm,pm, xj)). (5) 2.2 算算算法法法描描描述述述 下面给出有效的分类型数据聚类 (ECDC)算法 的主要步骤: 输入: 数据集X = {x1, x2, · · · , xm},聚类数 k; 输出:已聚好的类C1, C2, · · · , Ck. Step 1: 扫描数据集,指定聚类数 k; Step 2: 计算每个属性值对应的VIJ; Step 3: 根据式 (2)计算属性值间的距离; Step 4: 为对象随机分配 1, 2, · · · , k类标号; Step 5: 按照定义 5计算类中心; Step 6: 按照式 (5)重新计算每个对象到各个类中 心距离, 将其分配到离其最近的类, 待分配完后重新 计算类中心; Step 7: 若此时有大于或等于一个类的对象成员 为空,则返回Step 4; 否则,重复Step 6,直到类成员不 再变化或者达到指定的迭代次数为止,算法结束. 3 算算算例例例分分分析析析 文献 [9]已表明,对于如表 1所示的数据集,将其 聚为两个类最佳;而从表 1可以看到,对于类 c1,属性 Genre和Actor都只有 1个值, 只有属性Director有 2 个值, 类 c2中每个属性都有 2个值, 这便使得 c2的紧

专用实践 1.2 确定哪些培训需要是组织的职责 确定哪些培训需要是组织的职责，哪些是各个项目与支持组的职责。 关于项目与支持组特定的培训计划的更多信息，参见项目策划过程域。 组织的培训除了满足战略培训需要之外，还要满足跨项目和支持组的共同培训需要。项目与支持组对标识和 处理他们的特定培训需要负主要责任。组织的培训工作人员只负责处理跨项目与支持组的共同培训需要（例如， 有关多个项目的共同工作环境的培训）。然而在某些情况下，组织的培训工作人员也可以解决项目与支持组的附 加培训需要，但必须在培训资源可用且组织的培训优先的条件下，并按照与项目与支持组协商的结果进行。 典型工作产品： a) 项目与支持组共同的培训需要 b) 培训承诺 子实践： a) 分析由各项目与支持组所标识的培训需要。 分析项目与支持组的培训需要的目的是，标识出能在组织范围内 有效处理的共同培训需要。这些需要 分析活动，用于预计在项目与支持组级首先察觉到的未来的培训需要。 b) 与各项目和支持组协商如何满足他们的特定培训需要。 组织的培训工作人员能提供哪些支持，取决于培训资源是否可用，以及组织培训的优先级。 适合由项目与支持组实施的培训之例，如： 项目的应用领域或服务领域方面的培训 项目或支持组所使用的独特工具与方法方面的培训 在安全性、安全保密和人员因素方面的培训 c) 将向项目和支持组提供培训支持的承诺文档化。 专用实践 1.3 制定组织的战术培训计划 建立和维护组织的战术培训计划。 组织的战术培训计划是关于组织负责提供培训的计划，并且这种培训对个人有效履行其角色是必要的。此计

专用实践 1.2 使用组织过程资产策划项目活动 使用组织的过程资产与测量仓库，来估计和策划项目活动。 关于组织的过程资产与组织测量仓库的更多信息，参见组织过程定义过程域。 典型工作产品：

13.2 程序调用 到现在为止，我们介绍的都是如何在 Tcl解释器内部进行编程，其实 Tcl是一种通用的粘连脚 本语言，能将其他的程序和软件包按照求组装在一起。为了完成这种功能，Tcl提供了一些调用 其他程序的方法。主要有两中方法在 Tcl中调用一个程序：  open打开一个连接到文件描述符的进程管道来运行其他程序  exec 将一个程序作为子进程运行 这些命令以前是为 UNIX系统设计的，Tcl7.5以后的Windows版本也可以实现他们。为了能 保证操作正确性，建议本章内 exec命令和 open命令的程序在 UNIX系统内完成。 13.2.1 用 open命令打开一个进程管道 此处的 open调用和打开一个文件的命令相同，也返回一个文件描述符。如果文件名变量的 第一个字符是“ |“ 的话，open命令则将此变量剩余部分当作是一个程序名，并将运行 exec 来处理，而输入、输出被重定向到了文件描述符。管道则可被子进程打开用以只读、只写或者 读写。 进 程 A 进 程 B 写 读 读 写

第2章 多矿物模型 传统的测井资料分析方法与处理程序，基本上是按固定的解释模型，采用有限的几种测井 资料分析计算储层参数进行地层评价及油气分析的。多矿物分析模块采用体积理论模型，最优 化数学计算方法和面向对象的程序设计，充分利用现有的测井资料，采用不同的解释模型和解 释方法，对不同地质剖面所测的测井资料进行有效评价，对计算结果提供自我检验手段。使用 多矿物模型分析模块处理测井资料可提供以下成果： 岩性体积剖面； 孔隙度、含水饱和度和渗透率等多种储层参数； 地层多种物理特征参数（如视骨架时差，视颗粒密度等）； 连续的测井相剖面。 多矿物模型分析模块属于多步处理程序，主要包括模型分析、矿物计算、曲线合成、岩相 划分和岩相修正五个分析程序，这些程序关系紧密，一个程序的输出可作为另一个程序的输入， 其关系图如下： 这是一个多步处理的程序，根据地质需要分步进行。 2.1 建立测井响应参数数据库 测井解释响应参数数据表中包括各种矿物的各种测井响应值，用于产生线性方程组。处理 测井资料直接从数据表中调用。数据表中的数据主要来源于岩心分析、区域地质资料、文献和 过去的资料。测井解释响应参数数据库可由用户不断扩充和完善。 界面菜单工具栏下设有“模型参数库编辑”功能键。

第4章 地层倾角 模块采用可视相关对比法对地层倾角测井资料进行数字处理。倾角可视化地质应用工具包 括加速度校正、倾角可视相关对比、倾角合成、沉积构造分析、沉积相分析五个模块，能够处 理各种 4臂、6臂和 8臂地层倾角资料。 4.1 方法原理 随着油田勘探开发形势的不断发展，地层倾角测井资料的应用越来越广泛，它主要应用于 地质构造解释和沉积环境研究。Forward.net平台下的倾角可视化地质应用工具包括加速度校正、 倾角可视相关对比、倾角合成、沉积构造分析、沉积相分析五个模块，能够处理各种 4 臂、6 臂和 8臂地层倾角资料，其处理流程为： 1、数据加载 通过本平台提供的工具，将地层倾角原始数据加载到平台，形成相应的WIS数据格式文件。 （详见数据管理部分） 2、.深度校正 由于地层倾角测井的结果要同其他常规测井资料的地质解释结合使用，而地层倾角测井和 常规测井通常是两次测井，故要求两者深度一致。用倾角对比曲线或地层倾角带测的自然伽马 曲线与常规曲线对比，通过“曲线校深”模块进行深度校正。（详见曲线校正部分） 3、处理 用可视相关对比法对地层倾角测井资料进行数字处理。 4.1.1 原始测井曲线 倾角可视化地质应用模块可以处理国外斯仑贝谢公司（4臂、8臂）和阿特拉斯、哈里伯顿 公司（4臂、6臂）以及国内各公司的地层倾角测井资料。以下以 4臂倾角资料为例，进行具体 的说明。 前阿特拉斯的 CLS（3700）倾角仪和斯仑贝谢的 CSU测井系列 HDT的探测器都有一个四 臂极板系统，其外表相同而结构有差别，四个极板顺时针编号 1，2，3，4，两两相隔 90°。四 臂极板测井仪每一个极板上有一个微聚焦电极系，测四条电导率曲线，两条井径和其它曲线。 其原始曲线有： FC1 第一个极板测的电导率曲线。 FC2 第二个极板测的电导率曲线。 FC3 第三个极板测的电导率曲线。 FC4 第四个极板测的电导率曲线。 PIAZ I号极板方位角。

（2）电磁制动的特性 电机型号 项目 HF75B, HF105B HF54B, HF104B HF154B HF204B, HF354B 形式 （注 1） 弹簧关闭无励磁动作电磁制动 [保持・紧急控制用] 额定电压 DC24V 20℃时的额定电流 (A) 0.38 0.8 1.4 容量 (W) 9 19 34 静摩擦扭矩 (N･m) 2.4 8.3 43.1 惯量力矩（注 2） (kg･cm2) 0.2 2.2 9.7 释放延迟时间（注 3） (s) 0.03 0.04 0.1 制动延迟时间(DC断) （注 3） (s) 0.03 0.03 0.03 每次制动 (J) 64 400 4,500 允许制动 工作量 每小时 (J) 640 4,000 4,5000 电机轴的制动松动 (度) 0.1 ～ 0.9 0.2 ～ 0.6 0.2 ～ 0.6 次数 (回) 20,000 20,000 20,000 制动器寿命 （注 4） 每次制动量 (J) 32 200 1,000 (注 1) 没有手动释放装置。如需要对机械进行轴芯取出操作时，应另备 DC24V电源以电气方式释放制动。 (注 2) 对无制动的伺服电机累加的数值。 (注 3) 20℃条件下初始接通时的吸合间隙值。 (注 4) 制动间隙会因制动时制动闸片的磨损而增大，但间隙不是可调节的。因此以需要调节前的期限作为制动寿命。 (注 5) 有电磁制动的伺服电机在轴端会发生漏磁。 (注 6) 在低速区域运转时，有时会发出闸片的松动声，此现象并非功能上出现问题。

4.6 常见问题 指针可以对内存进行读写操作，对指针错误的理解和应用可能造成内存数据的错误。本 节将讲解指针操作时常见的错误。用指针读取数据时，不正确的访问可能读取不到需要的数 据。当用指针进行写操作时，错误的操作可以破坏内存中其他程序的数据。 4.6.1 错误的写操作 在使用指针时，需要牢记，不可以随意向不确定的内存单元写入数据。这种写操作可能 破坏内存中其他的数据，引起其他程序的错误。正确的方法是，一个指针先取得正确的内存 地址，然后访问一个确定的内存单元。例如下面的程序就是错误地向一个未知的内存单元写 入数据。 #include <stdio.h> /*包含一个头文件。*/ main(void) /*主函数。*/ { int i; /*定义一个整型变量。*/ int *p; /*定义一个指向整型变量的指针。*/ i=10;

图 9.69 创建图层动画对话框 图 9.70 导入路径生成动画工具 图 9.71 沿路径飞行动画对话框

（5） 运行和使用模型 经过模型验证并对模型结果满意后，即可双击模型图标 ，实现基于 DEM 的水系 自动提取。 12.3脚本文件 12.3.1简介 空间数据的处理常常是复杂而费时的，而且其过程具有很强的重复性，有必要进行自 动化的处理。任何可以支持 COM 的脚本语言都可以执行 ArcGIS9 的地理处理工具，如 Python、Jscript和 VBScript等。这些脚本语言都是公开的，而且易学易用。 脚本可以通过一个工具或多个工具实现一个简单或者复杂的处理，也可以通过循环操 作对输入数据进行批处理。因为数据不是特定的，所以脚本可以重用。脚本的高效性还体 现在可以独立于 ArcGIS9桌面程序执行。 熟悉 ArcInfo Workstation的 AML用户转而使用一种新的脚本语言是非常容易的。在 ArcGIS9 中 AML 不仅可以执行 Arc 命令，而且任何 AML 程序都可以作为脚本添加到 ArcToolbox中。对不熟悉脚本语言的用户，模型生成器是构建脚本的方便工具，只要先构 建一个模型再输出成脚本即可。但是模型不可以独立于 ArcToolbox运行，脚本却可以脱离 ArcGIS的环境独立运行。而这个脚本可以是任何支持 COM的脚本语言。 12.3.2 脚本编写基础 通过脚本调用 ArcToolbox进行空间处理，实际上是通过 GPDispatch组件类来实现的。 n GPDispatch组件类的作用就是调用和执行各种地理处理工具。它实现了 IGPComHelper 接口和 IGPDispatch接口。关于这两个接口的具体使用可以参考 ArcGIS的在线帮助。 n 在 VBS 脚本中通过 WScript.CreateObject("esriGeoprocessing.GPDispatch.1");语句可以 创建一个 GPDispatch对象。其中 esriGeoprocessing指的是 GPDispatch所在的类型库。 后面的“1”表示的是版本号。GPDispatch 最常用的方法是 CheckOutExtension()检查 扩展许可、AddToolbox（）添加工具、RemoveToolbox()移除工具、Execute()执行工具。 n 其中 Execute()用法较为特殊下面为它的函数原型： Set variable = object.Execute (ToolName, params, trackcancel ) ToolName为工具名称，params为工具的一系列参数参数，trackcancel用于跟踪 cancel 键是否被按下。

7.7 EMI问题 DDR 由于其速度快，访问频繁，所以在许多设计中需要考虑其对外的干扰性，在设计时需要注意一下几点 原理有性能指标要求的，易受干扰的电路模块和信号，如模拟信号，射频信号，时钟信号等，防止 DDR 对其干扰，影响指标。 � DDR 的电源和不要与其他易受干扰的电源模块使用同一电源，如必须使用同一电源，要注意使用电 感、磁珠或电容进行滤波隔离处理。 � 在时钟及 DQS 信号线上，预留一些可以增加的串联电阻和并联电容的位置，在 EMI 超出标准时，在 信号完整性允许的范围内增大串联电阻或对地电容，使其信号上升延变缓，减少对外的辐射。 � 进行屏蔽处理，使用金属外壳的屏蔽结构，屏蔽对外辐射。 � 注意保持地的完整性。 7.8 测试方法 注意示波器的探头和示波器本身的带宽能够满足测试要求。 测试点的选择要注意选到尽量靠近信号的接受端。 由于 DDR 信令比较复杂，因此为了能快速测试、调试和解决信号上的问题，我们希望能简单地分离 读/写比特。此时，最常用的是通过眼图分析来帮助检查 DDR 信号是否满足电压、定时和抖动方面的要求。 触发模式的设置有几种，首先可以利用前导宽度触发器分离读/写信号。根据 JEDEC 规范，读前导的 宽度为 0.9 到 1.1 个时钟周期，而写前导的宽度规定为大于 0.35 个时钟周期，没有上限。第二种触发方 式是利用更大的信号幅度触发方法分离读/写信号。通常，读/写信号的信号幅度是不同的，因此我们可以 通过在更大的信号幅度上触发示波器来实现两者的分离。 测试中要注意信号的幅度，时钟的频率，差分时钟的交叉点，上升沿是否单调，过冲等。 时序中最重要，最需要注意的就是建立时间和保持时间。 1. 写建立分析： DQ vs. DQS 2. 写保持分析： DQ vs. DQS 3. 读建立分析： DQ vs. DQS 4. 读保持分析： DQ vs. DQS

10.15 ASCII坐标转换 ASCII Coordinate Conversions选项用于将一种投影和参数的ASCII坐标转换为另外一种。输入的ASCII 文件必须包含X、Y坐标栏，以及可选的Z 值（高程）。如果在输入的坐标中不包含高程值，则假定高程 为0。输出文件将包含处于新选投影中的 X、Y（Z）栏。可以将坐标输出到ENVI的地面控制点文件中。 选择Map > ASCII Coordinate Conversions。选择输入的ASCII坐标文件。 要向“Selected Input Files”列表中添加文件，点击“Input Additional File”。要从列表中删除文件，在 “Selected Input Files”列表中选择一个文件名，然后点击“Delete”。 键入一个输出文件名，或用“Choose”按钮选择一个输出文件名。点击“OK”继续。当出现 ASCII

13.11 TOPSAR工具 使用TOPSAR Tools选项可以读取综合的TOPSAR数据、浏览头文件，以及把综合的 TOPSAR文件 转换为实际单位。综合的TOPSAR数据通常包括一个数字高程模型 (DEM)、C-波段VV 数据、入射角图 像、相关性图像，以及L-和P-波段极化AIRSAR数据。其它分析工具的详细介绍，请参阅第515页的“极 化工具”。 关于综合TOPSAR数据的详细介绍，请参阅以下参考文献： AIRSAR Integrated Processor Documentation, Data Formats, Version 0.01, May 3, 1995, JPL Publication, Jet Propulsion Laboratory, Pasadena, CA. （22） 打开TOPSAR文件 要读取原始的TOPSAR格式数据文件 (Cvv、入射角、相关性图像、或DEM)：选择Radar > TOPSAR Tools > Open TOPSAR File。 当出现文件选择对话框时，选择需要读取的文件。ENVI自动提取所需的头文件信息，并把图像波段 放置在可用波段列表中。 注意：要读取所有的TOPSAR文件，并自动地把它们转换为实际单位，请参阅下一节的内容。要合成AIRSAR 图像，请参阅第515页的“JPL AIRSAR数据合成”。 （23） TOPSAR数据转换 使用Convert TOPSAR Data选项可以将原始TOPSAR格式数据文件读取和转换为实际单位(浮点数据)。 Cvv文件被转换为sigma zero。入射角文件被转换为0-180度角。相关性文件被转换为0-1值。DEM被转 换为米。 选择Radar > TOPSAR Tools > Convert TOPSAR Data。当出现文件选择对话框时，选择一个数据文件。 当出现TOPSAR Conversion Parameters对话框时，通过点击文件名旁的复选框，选择要转换的文件。 若需要，点击“Spatial Subset”按钮，然后使用标准的ENVI构造子集方法选取任意空间子集。选择 输出到“File”或“Memory”。 对于每个输入文件，ENVI都将所有的转换文件放置在包含一个波段的单个输出文件中。 （24） 替换DEM坏值 使用Replace Bad Values或DEM Replace Bad Values选项，利用表面拟合技术计算所得到的数值，来 填充数字高程模型(DEMs) 中的坏数值。

12.7 不规则点栅格化 使用Rasterize Point Data选项可以将不规则栅格数据内插为一幅栅格图像。ENVI 的栅格化功能是使 用平面点的Delaunay三角测量来实现的。在不规则栅格数据点被三角测量后，它们将被内插到一个规则的 网格上，可以选用线性内插或连续的五次多项式内插。对三角形区域外的点可以选用外插方法。栅格点读 自 ASCII 文件，支持不同的输入和输出投影。 详细介绍，请参阅第489页的“不规则点栅格化”。 12.8 矢量地形图转化为栅格 DEMs 矢量和栅格格式的地形数据通常都是可用的。不同的格式有不同的优点，ENVI为从栅格数字高程模 型数据中提取附加信息提供了一些工具。例如：栅格DEMs可用于生成附加地形数据集，其中包括：坡度、 坡向、阴影地貌（参见第491页的“地形建模”）；或用于生成地形分类图像（参见第493页的“提取地形 特征”）。当把图像拖放到地形高程数据上，用于3-D曲面浏览时，也需要栅格DEM。 Convert Contours to DEM选项利用矢量记录结点为DEM进行栅格的三角测量。ENVI建议设定的输出 像元尺寸要适当大于矢量结点的采样距离。如果输出像元尺寸与结点采样距离差距很大，矢量等高线将自 动变的非常密集，从而减少内插次数。除非相应的高程属性值处于用户指定的有效高程范围之外，否则包 含在文件中的所有矢量记录都会被用于转化中。 ‧ 输入所需值 从EVF文件输入。矢量数据必须是ENVI矢量格式数据（.evf），如果矢量地形包含在外部格式文件中 （例如：ArcView Shape文件），把它转换成ENVI矢量格式也很简单——只需在ENVI中读取数据即可， 输入程序将提供数据到ENVI矢量格式的转换（参见第42页的“打开矢量文件”）。 所需高程属性。EVF文件必须包含用于指定每个矢量等高线的高程的属性。 ‧ 将等高线转换成 DEM 在ENVI主菜单中，选择Topographic >Convert Contours to DEM，或Vector > Convert Contours to DEM。 在文件选择对话框中，选择包含地形等高线的EVF文件和高程属性。将出现Convert Elevation Contours to Raster DEM对话框。在“Elevation Attribute Column”标签旁的下拉菜单中，选择包含等高线高程的属性。 如果需要，输入有效高程范围（与高程属性单位相同），高程处于有效范围外的矢量在构建DEM时将 被忽略。键入输出像元尺寸和数据类型。 注意：像元尺寸在处理开始之前可以更改（参见下面的第 9步）。建议设定的输出像元尺寸要适当大于矢

（3） 波谱分析信息解译 波谱分析的输出是所分析的波谱与输入波谱库中每种要素的波谱的匹配得分排序或权重的得分。最高 得分表示匹配最好，波谱相似度的可信度高。相似的要素得分应该相对较高，但是不相关的要素得分应该 较低。 （4） 成功使用波谱分析的提示 波谱分析是建立在波谱匹配技术之上的，对成功运作有具体的要求。下面的条目总结了要成功使用该 工具应考虑的事项。 ‧ 波长范围 许多要素在一个波长范围内是类似的，但是在其它范围内具有较大差异。包含特定吸收特征的波长范 围将产生最好的结果。当显示一个波谱时，波谱分析工具将对显示在图表窗口中的整个范围进行分析。如 果需要使用一个特定范围，在图中用鼠标中键将波谱图放大到所需的波长范围，然后点击Spectral Analyst 窗口中的“Apply”。 （5） 方法的选择 判定要素是否具有吸收特征。如果它们有，最好使用波谱特征拟合方法。此外，波谱角分类和二进制 编码也能得到较好的结果。 ‧ 库中要素缺乏时的情况(与其它要素的类似) 如果感兴趣的要素波谱在波谱库中不存在，它就不能被识别出来，此时，类似的材料也许得分值相对 较高。

（6） 卷积滤波和形态学滤波的应用 可以将滤波结果应用到一个临时显示窗口或当前存在的文件中。选择“Quick Apply”方法可以快速试 验不同的滤波参数和浏览滤波结果；选择“Apply to File”方法可以将滤波结果应用到一个特定的图像文件 中，并进行输出。 ‧ 快速浏览滤波结果 在Convolutions and Morphology Tool对话框中，使用“Quick Apply”按钮可以快速检验不同的滤波结 果。可以根据需要进行更改并重新显示图像。点击“Quick Apply”按钮，将所选滤波应用到一个单独波段。 在Convolution Quick Apply Input Band对话框中，选择输入波段，若需要，选取任意子集。 注意：如果选择一个输入波段，无法选取波谱子集。 点击“OK”。结果文件出现在一个显示窗口中，该文件是一个临时文件，被放置在ENVI默认的临时 文件目录下（参见第82页的“设置ENVI参数”）。当关闭Convolutions and Morphology Tool对话框时，显 示窗口和文件都将被删除。 如果在Convolutions and Morphology Tool对话框中更改了参数设置，重新点击“Quick Apply”按钮。 根据新参数生成的结果将覆盖当前结果。 更改快速应用的输入波段 要更改用于快速应用输出结果的波段：Convolutions and Morphology Tool对话框中，选择 Options >Change Quick-Apply Input Band。文件选择对话框中，选择输入波段，若需要，选取任意子集。 注意：如果选择一个输入波段，无法选取波谱子集。 点击“OK”。 保存快速应用结果 要将快速应用结果保存到一个文件：选择File >Save Quick Result to File。当出现Output Convolution Filename对话框时，键入一个输出文件名。点击“OK”。 ‧ 将滤波应用到文件 在Convolutions and Morphology Tool对话框中，选择“Apply to File”按钮可以将当前滤波应用到一个 文件中，并进行输出。点击“Apply to File”按钮。当出现文件选择对话框时，选择输入文件，若需要，选 取任意子集。点击“OK”。将出现Convolution Parameters对话框。选择输出到“File”或“Memory”。作 为结果的文件将被列在可用波段列表中。 8.4 纹理滤波器 使用 Texture选项可以应用基于概率统计或二阶概率统计的纹理滤波。许多图像包含的区域以亮度变 化为特征，而不仅仅局限于亮度值。纹理是指图像色调作为等级函数在空间上的变化。被定义为纹理清晰

（5） 决策树节点选项 决策树节点和类别节点包含的选项有：剪除、恢复、删除和添加子节点；编辑节点属性；以及将剩余 像元保存到一个掩膜文件中。在节点中点击鼠标右键，在出现的快捷菜单中也可以访问这些选项。 ‧ 剪除/恢复子节点 要暂时剪除一个节点下的子节点并将该节点内的所有类归为一个单独的类别，在节点上点击鼠标右 键，在出现的快捷菜单中选择“Prune Children”。 要恢复剪除的子节点，在节点上点击鼠标右键，在出现的快捷菜单中选择“Restore Pruned Children”。 ‧ 删除子节点 要彻底删除决策树节点下的子节点，并将该节点变为类别节点：在节点上点击鼠标右键，在出现的快 捷菜单中选择“Delete Children”。 ‧ 添加子节点 要添加一个子节点，并将类别节点变为决策节点：在一个类别上点击鼠标右键，在出现的快捷菜单中 选择“Add Children”。 ‧ 编辑决策节点属性 要更改节点名称和表达式，在节点上点击鼠标右键，在出现的快捷菜单中选择“Edit Properties”。 ‧ 编辑类别节点属性 要更改类别名称和颜色，在一个类别上点击鼠标右键，在出现的快捷菜单中选择“Edit Properties”。 ‧ 将剩余像元保存到掩膜文件 要将剩余像元保存到一个掩膜文件中，在节点或类别上点击鼠标右键，在出现的快捷菜单中选择“Save Survivors to Mask”，并键入一个文件名。 （6） 决策树文件管理 使用File菜单中的选项可以进行如下操作，包括：构建新的决策树、将决策树保存到ASCII文件、恢 复保存的决策树以及退出该功能。 选择File >New Tree，构建新的决策树。选择File >Save Tree并键入一个输出文件名，将决策树保存到 ASCII文件。选择File >Restore Tree并选择一个输入文件，恢复先前保存的决策树选择File >Cancel，关闭 决策树对话框。 （7） 编辑现存的决策树 使用Edit Existing Decision Tree选项可以恢复一个现存的决策树，并根据需要对它进行编辑。

（2） 反向MNF变换 使用Inverse MNF Transforms工具可以将MNF波段变换回原始数据空间。 提示：在运行反向变换之前，应该先对以噪声为主导的图像（高波段号）进行平滑，或在反向变换过程中， 用波谱子集消除噪声波段。平滑或消除噪声波段将减少原始数据空间中的噪声。 选择Transforms > MNF Rotation> Inverse MNF Transform，或Spectral > MNF Rotation> Inverse MNF Transform。当出现标准ENVI文件和子集选择对话框时，用标准ENVI文件选择程序选择以前计算过的 MNF变换图像，若需要，选取任意子集。当出现Enter Forward MNF Stats Filename对话框时，键入正向 MNF统计文件的名称，点击“OK”。当出现Inverse MNF Transform Parameters对话框时，选择输出到“File” 或“Memory”。从“Output Data Type”菜单中，选择所需的数据类型（字节，整型，无符号整型，长整型， 无符号长整型，浮点型，双精度型）。点击“OK”开始处理。 处理完成后，反向MNF变换结果将会被导入到可用波段列表中。 （3） 把正向的MNF变换应用到波谱分析 Apply Forward MNF to Spectra工具可以将端元波谱变换到MNF空间，用于混合调制匹配滤波（参见 第404页的“混合调制匹配滤波”）。该变换要使用混合调制匹配滤波的输入数据文件的MNF统计。使用 一个对话框（类似于端元收集对话框）可以从图表窗口、波谱库、ASCII 文件、感兴趣区和统计文件中收 集将被变换的波谱。 选择Transforms > MNF Rotation> Apply Forward MNF to Spectra。当出现Forward MNF Statistics Filename对话框时，选择为混合调制匹配滤波输入数据文件计算的MNF统计文件名。将出现Forward MNF Convert Spectra对话框。使用Import菜单选项将波谱输入到对话框中，或将波谱拖放到对话框顶部的黑色 拖放小组件中（详细介绍，请参阅第298页“拖放-下拉窗口的使用”、“输入端元波谱”和第301页“端元 Option下拉菜单”）。 注意：在把波谱库中的波谱变换到正向MNF空间，用于混合调制匹配滤波时，来自波谱库的波谱应该被 等比例缩放到使用MNF统计进行正向变换的图像中。 要从对话框中删除波谱，点击波谱名，然后点击“Delete Spectrum”。点击“Apply”。变换后的波谱将 出现在Forward MNF Spectra图表窗口中，并能被输入到Mixture Tuned Matched Filtering Endmember Collection对话框中。 ‧ 变换后的波谱输出 可以打印变换后的波谱图，或将它们保存为ASCII 文件、波谱库、IDL变量、脚本文件或图像文件。 要输出波谱图，选择File >Save Plot As >MNF图表窗口的输出类型。要打印波谱图，选择File >Print。 详细介绍，请参阅第200页的“输出图表数据”。

（4） 浏览图像像元分布 使用Scatter Plot Dancing Pixels TM 功能可以在散点图窗口交互式浏览图像窗口中高亮显示的像元分布。 选择Options >Image :Dance。在主图像窗口内部的任何一处放置指针，当散点图被显示时，点击鼠标 左键，在散点图上用彩色标记像元。这表明了落在与散点图中用到的一样大小的空间方框内的像元的 DN 值分布（见第220页的“设置方框大小”）。由于速度因素，该方框没有被绘制在主图像窗口中。 当在主图像窗口中移动指针时，按住鼠标左键，把这些像元作为跳动像元实时地显示在散点图中。 图 4-62：显示在散点图窗口的“跳动”像元

（3） 多分辨率地理坐标定位后图像的链接 使用 Geographic Link工具来链接包含地理坐标定位数据的显示窗口或矢量窗口。链接以后，当移动指 针时，所有显示的地理坐标定位后的图像和矢量窗口将随指针位置及时更新。该功能不考虑投影、像元尺 寸和每个数据集的旋转因子。 注意：该功能无法进行如下操作，包括：重新定义投影（on-the-fly）、重采样以及动态覆盖。要对数据集 重新定义投影和重采样以达到相同的投影和分辨率，请参阅第 238页的“层的叠加”。 从主图像窗口的菜单栏中选择 Tools >Link >Geographic Link。在 Geographic Link对话框中，通过点击 显示名旁的箭头切换按钮（切换为“On”）来选择要链接的显示。点击“OK”。当在一个经过地理坐标定 位后的图像或矢量窗口移动指针时，其他经过地理坐标定位后的链接图像或矢量窗口中的指针也移动到相 同的地理位置。 注意：要关闭基于地理坐标的链接，选择 Tools >Link >Geographic Link，点击显示名旁的箭头切换按钮， 切换为“Off”。 4.22 交互式剖面和波谱图 ENVI允许抽取水平的（X）、垂直的（Y）、波谱的（对每个像元为 Z ）以及任意的剖面图。 剖面图显示在单独的图表窗口中，X、Y 和 Z剖面图是可以同时激活的。鼠标用来移动一个十字准线 并交互地选择剖面图。剖面图是标准的 ENVI图，有关控件和选项的详细描述，请看第 193页的“交互式

（2） Options下拉菜单 使用 Options下拉菜单可以获取分类的统计信息，以及更改类别的颜色和名称。 • 浏览类别分布 要产生一个记录显示像元数和每类像元的百分比，选择 Options > Class distributions。这个记录被处理 和显示在 Classification Distribution窗口中。在编辑过程中，在类中添加或删除像元后，该记录将被自动更 新。 • 编辑类别名称和颜色 要更改类别颜色和/或名称，选择 Options > Edit class colors/names。将出现 Class Color Map Editing对 话框。详细介绍，请参阅第 206页的“分类彩色制图”。 • 类别合并 要把一个或多个分类合并到一个所选择的基类中：在 Interactive Class Tool对话框中，选择 Options > Merge classes。将出现 Interactive Merge Classes对话框。通过点击所需要的类名，选择基类。通过点击所 需要的类名，选择要合并到基类的类别。点击“OK”。下一次打开 Class Distributions对话框时，它将被更 新。 • 绘制平均波谱曲线 要为活动类绘制平均波谱曲线，选择 Options > Mean for active class。若出现 Input File Associated with Classification Image 对话框，选择计算统计信息所需的输入文件，然后点击“OK”。只有先前没有设置输 入文件，这个对话框才会出现。 要绘制所有分类的平均波谱曲线，选择 Options > Mean for all classes。 ENVI计算统计信息，并在一个 Class Means图表窗口内显示平均波谱曲线。

（3） 保存分类图像 在 Interactive Class Tool对话框中，选择 File > Save Changes to File。 注意：选择“Yes”，将覆盖现有的分类文件。 在询问对话框内，点击“Yes”。当前的文件名将被更改的文件覆盖。 （4） 分类显示输出 要输出被分类图覆盖的图像，选择 File >Save Image As >所需要的输出。要打印被分类图覆盖的图像， 选择 File >Print。 注意：只有那些被启动的类出现在输出中。 分类覆盖图将作为 <Displayed Class Overlay>，列在 Change Graphics Overlay Options 对话框里的 Graphics Overlay Selections列表中。详细介绍，请参阅第 227页的“设置图表覆盖选项”。 （5） 隐藏Interactive Class Tool对话框 要在不擦除分类结果的情况下，隐藏或显示 Interactive Class Tool对话框，参见第 118页的“显示和隐 藏叠加对话框和层”。

（3） 网格Options菜单 Grid Line Parameters对话框中的 Options菜单允许更改三种网格的特征和显示边框。 • 更改边框 要更改显示边框的大小和颜色：选择 Options > Set Display Borders。为左、上、右和底部输入所需要 的边框大小（以像元为单位）。从“Border Color”菜单中，选择所需要的颜色。 • 编辑网格特征 要更改颜色、标注尺寸和字体以及其他网格特征：在 Grid Line Parameters对话框中，选择 Options > Edit Pixel/Map/Geographic Grid Attributes。当出现 Edit Map Attributes对话框时，选择下列步骤，更改标注、线 条、方框和拐角的特征。 要开启和关闭绘制在显示图像上的特定网格要素，切换每个要素名旁的“On/Off”箭头切换按钮。要 在图像外部虚拟边框内绘制网格标注，开启“Labels”。刻度也将沿着图像的边缘，被放置在标注的附近。 要更改颜色，点击“Colors”，然后选择一种颜色。要更改标注的宽度，在“Thick”文本框中输入一个值。 要更改标注的字体，从“Font”菜单中选择。要更改大小，在“Charsize”文本框中输入一个值。要为每种 网格类型，控制标注与图像边缘的距离，使用“Dist”参数，数值越大，离边缘越远。“X-axis Labels”和 “Y-axis Labels”箭头切换按钮，分别为 X和 Y网格轴的标注在水平 H 和垂直（V）显示之间切换。要绘 制 X和 Y网格线，开启“Lines”选项。要绘制围绕图像的一个“Neat Line”方框，开启“Box”选项。要 为每个网格交点绘制刻度线，开启“Corners”。要为网格和方框线设置线条类型（实线、点线、虚线等）， 从“Style”选项中选择。要设置选项的宽度，在“Thick”参数框中输入所需要的值。 当根据需要设置好网格线属性后，点击“OK”。点击“Apply”，在图像上绘制网格线。