利用PCI进行遥感异常提取流程.docVIP

遥感异常提取工作手册 第一节、利用PCI进行遥感异常提取流程 遥感异常提取的过程可以用下面的流程图来描述： 去干扰异常主分量阈值化技术流程图 注：* OHA即HA,为羟基异常的简称; FCA即FA,为铁染异常的简称. **阈值化,在张玉君教授以前的著作中称为门限化。 （本流程以14231景ETM数据为例进行处理） 一、了解该景的基本信息 拿到该数据后，首先查看该景数据的基本信息。在写字板中打开后缀为“HRF-FST”的数据，出现下面内容： 记录下划线表示的信息：获取日期（ACQUISITION DATE）、图像的像素数（PIXELS PER LINE、LINES PER BAND、图像的位置坐标（UL）、太阳高度角（SUN ELEVATION）、太阳方位角（SUN AZINUTH ANGLE）。 在写字板中打开后缀为“MLT-FST”的数据，查出： 记录下各波段获取时的灵敏度“L”或“H”。 二、计算准归一化因子 经试验推算，卫星在近日点和远日点时所获取数据相差达6.7%，因此日地纠正不可忽略。该景位于东经80-90度之间，90/15=6，即在第六时区。卫星过境约为上午十时，实为格林尼治10-6=4时，查日地距离表，可知获取该景数据时的日地距离d=1U。 日地校正因子公式：d2/sinG。其中：d=1见附表1），=41.20，G=0.2（G为数，也可以取0.25，视图像而定），代入后得校正因子为7.6054726。 三、建立数据库 （一）建立“.pix”数据库 利用“HRF.FST”引导文件导入1，-6（实际为1，-5波段，7波段）。点击“”，弹出如下 “ImageWorks Configuration”对话框，设置“Image Size”“Visible Window”“Image Planes”参数，为了保证显示窗口放大而且图像的行列不失调，显示窗口（Visible Window）的X=1024*6967/5965，其中6967=PIXELS PER LINE，5965=LINES PER BAND，另外“Image Planes”选择“8”，即控制面板中出现8个电子枪。 点击“Accept”。在“File”下拉菜单中选择“Utility…”, 选择“Tools”下拉菜单中的“Subset…”，弹出下面的对话框，输入原文件（_HRF.FST）的路径名称，确定输出的文件路径名称，选择需要添加的通道，如下所示： 输入原文件（_HRF.FST）的路径名称，确定输出的文件路径名称，在“Selectable Layers”框下将需要的波段从源文件层添加到目标层，点击“Subset”，如上图框所示。 （二）“.pix”数据库中增加25个通道 点击“”，在“View”下拉菜单下，选择Aplications/Database Utilities，选择“Tasks”下的“PCIMOD”，出现右下边的对话框，填写“FILE”，“PCIOP”（点击后，出现右上图框，选择“ADD”），“PCIVAL”（数字表示8U、8S、16U、32R四种类型的通道数目）。如下所示填好参数，点击“Status”“Run”。 四、全图统计 在“View”下拉菜单下，选择“Histograms…” 弹出“Multi Histogram Panel”控制面板，双击直方图可以弹出如下统计图，从图中读取并记录该通道像素的最大最小值，其它的通道的最值依次读取。（注：选中“From File”，即在文件中而非视图中进行统计） 表1 各波段的灰度最值统计表 波段（Band） DN Min DN Max 1 29 108 2 19 114 3 16 147 4 8 111 5 4 143 7 3 124 五、去干扰 1、去除边框 ETM+各通道组成的图像边框不整齐，因此需要进行去边框处理。对整幅图像来说，最靠左边数据是Band1采集的，最右边数据Band5采集的，所以可以利用Band1*Band5处理。 点击右上图的“”，出现右下图的“Xpace”工具栏，选择好命令后，出现右下图的对话框，设置好参数，点击“Status”“Run”。经过运算，去边框后图像存在第8通道中。其中：“OPER”表示操作的运算；“CNST”表示运算系数；“DBIC”表示参与运算的通道；“DBOC”表示输出的通道；“ZERODIV”表示被0除的数输出值；“RVAL”填写输入范围与输出范围（之间用逗号间隔）。 2、去除白泥地 首先看白泥地在各波段的波谱特征，目视解译及现场考察白泥地的分布，找出典型白泥地的位置进行其波谱特征分析。如下： 在图像上选择含白泥地的某一剖面如上图所示，在“View”下拉菜单中，点击“Profile Windows…