第4章空间域处理方法扩展阅读.pdfVIP

空间域处理方法补充案例 （配第4 章）  测试数据 数据目录：附带光盘下的..\chapter04\data\ 文件名 说明 Landsat8_OLI_b1.dat 某地的Landsat 8 OLI 图像某波段，灰度值被拉伸到0 ~ 100 Air_pressure_3d.dat 某地的3 维大气压数据，用于剖面提取，具体参见本章补充材料  IDL 程序 程序目录：附带光盘下的..\chapter04\programs\ 文件名 说明 zhu_local_statistic.sav IDL 代码，用于图像邻域统计，具体参见本章补充材料 4.1 邻域统计 邻域统计常用的指标有多样性（diversity ）、密度（density ）、众数（majority ）、 少数（minority ）、求和（sum ）、均值（mean ）、标准差（standard deviation ）、最 大值（maximum ）、最小值（minimum ）和秩（rank ）等10 个指标（指标含义详 见朱文泉等编著的《遥感数字图像处理——原理与方法》第3 章3.1.1 节中的邻 域统计部分）。目前，ENVI 软件尚无上述几种统计功能，需利用IDL 程序实现， 本书提供封装好的IDL 程序zhu_local_statistic.sav 可实现以上 10 种邻域统计功 能，其基本操作流程如下。 （1）打开图像。在菜单栏中，单击File Open ，选择Landsat8_OLI_b 1.dat 文件，则成功加载图像。 （2 ）将 ENVI 中的图像导出为 IDL 中的变量。在工具箱中，选择 Raster Management IDL Export to IDL Variable ，在Export to IDL Input File 中选择待 处理的波段。单击OK ，弹出Export Variable Name 对话框（图4.1）。 1 图4.1 Export Variable Name 对话框 （3 ）设置变量名称。在Export Variable Name 对话框中，New Variable Name 文本框输入变量名，这里我们将变量命名为DATA 。点击OK 后，我们在IDL 界 面中变量查看器中看到出现一个DATA 的变量（图4.2 中的红框处）。 图4.2 IDL 界面 （4 ）导入邻域统计程序。在IDL 操作界面的IDL 控制台输入下列代码，导 入程序。其中需将 D:\RSDIP\chapter04\documents\zhu_local_statistic.sav 替换为 程序实际所在路径。 2 RESTORE,D:\RSDIP\chapter04\programs\zhu_local_statistic.sav 程序使用及参数说明： 程序示例： result=ZHU_LOCAL_STATISTIC(image, win_size, /diversity) 参数： image: 输入的数字图像，为一个二维数字矩阵。 win_size: 邻域窗口大小，为一个2 元素一维数组，数组元素均为奇数。如[3, 5]， 表示邻域大小为3 列5 行。 关键词： Diversity ：计算邻域的多样性 Density ：计算邻域的密度 Majority ：计算邻域的众数 Minority ：计算邻域的少数 Sum：计算邻域的总和 Mean ：计算邻域的均值 Std：计算邻域的标准差 Maximum ：计算邻域的最大值 Minimum ：计算邻域的最小值 Rank ：计算邻域的秩 （5