遥感数据预处理

本研究以 2008 年 3 月 16 日地震前 IKONOS 遥感影像和 2008 年 9 月 1 日地震后QuickBird 遥感影像为数据源，采用基于多源多时相变化检测技术开展遥感震害信息提取。

基于多源多时相变化检测技术的遥感震害信息提取数据的预处理不同于普通的遥感影像数据的预处理，其对两时相影像质量的要求较高 ( 不管在辐射校正、几何校正还是影像配准等环节中均要求比较高的精度) ，因为这直接关系到后续震害信息提取的准确性。为满足 “快速、高效、准确”的要求，本节提出了基于变化检测技术的多源多时相遥感震害信息提取影像预处理技术流程 ( 图 4 －2) ，通过实验证明取得了良好效果。

( 一) 无控制点共线方程几何校正法

共线方程建立在图像坐标与地面坐标严格的数学变换关系基础上，是对成像空间几何形态的直接描述。该方法的校正过程需要用到数字高程模型，可以在一定程度上修正因地形起伏而引起的投影差和几何变形。当今，所有的卫星遥感数据都附带了卫星传感器的精确位置、高度、速度、太阳高度角和姿态等参数，这些信息一般保存在遥感影像的头文件或 ＲPC 文件里，所以能够十分方便地在没有地面控制点的情况下使用共线方程几何校正法进行较高精度的几何纠正和定位。

高分辨率遥感影像由于幅宽窄、空间分辨率高，受地球切平面、地球曲率等影响相对较小，影像内部几何畸变较小，所以一般在使用了共线方程几何校正法对高分辨率遥感影像进行几何校正后，都能够达到比较好的效果。由于共线方程几何校正法只需要提供卫星传感器飞行的相关参数就能对影像进行纠正，省去了选取控制点的步骤，节省了时间，满足了 “快捷、准确”的要求，所以本研究选择无控制点共线方程几何校正法对影像进行第一次 “粗”校正。

在 ENVI 软件中，分别读取 IKONOS 和 QuickBird 数据的 ＲPC 文件 ( . txt 格式) ，然后在 Georeference 模块中分别对全色和多光谱波段进行无控制点共线方程的几何校正。

图 4 －2 遥感震害信息提取影像预处理技术流程图

( 二) 正射校正

遥感图像成图时，由于受到各种不确定因素的影响，例如传感器的成像方式、外方位元素的变化、地形起伏、地球曲率、大气折射等，图像本身的几何形状与其对应的地物形状往往会不一致，发生几何变形 ( 畸变) 。遥感图像的几何变形是指原始图像上各地物的几何位置、形状、尺寸、方位等特征与在参照系统中的表达要求不一致时产生的变形。为了消除这些因素带来的几何变形，为后续影像配准做好铺垫，还需要利用研究区域的DEM 对影像做数字正射校 正，分别生 成 震 前、震 后 两时相 的 数 字 正 射 影 像 图 ( Digital Orthophoto Map，简称 DOM) 。数字正射纠正的原理就是将中心投影的影像通过数字元纠正形成正射投影的过程 ( 陈文凯，2007) 。

本文的正射校正在 ENVI 软件的 Orthorectification 模块中完成。得到震前、震后两时相的 DOM 后，还应当检查其与 DEM 的匹配情况，与 DEM 同名地物点的点位中误差不能大于表 4 －1 的规定，如果超过规定，需要重新进行正射校正。

表 4 －1 DOM 与 DEM 地物点的点位中误差

( 三) 图像融合

对全色数据与多光谱数据的 DOM 进行融合，形成兼具高分辨率空间信息和多光谱彩色信息的融合影像。融合前须对多光谱数据进行色彩增强处理，拉大不同地类之间的色彩反差，突出其彩色信息同时对影像进行色调调整，提高全色数据的对比度和亮度，增强局部反差，突出纹理细节，降低噪声。融合后须检查影像是否出现重影、模糊等现象，检查影像纹理细节与色彩，判断融合前的处理是否正确，如果存在以上问题，需要返回重新融合。融合后影像如果亮度偏低、灰阶较窄，可采用线性拉伸、亮度对比度等方法进行色调调整，但是应注意尽量保留融合数据的光谱信息和空间信息。

研究为了使融合后的数据仍然保持多光谱特性 ( 四个波段: 红、绿、蓝、近红外) ，便于标准化植被指数 NDVI 的计算，在 EＲDAS 软件下，采用 Subtractive resolution merge 模块( 此种融合方法能够使融合后的数据保持原有的多光谱特性) ，分别对 IKONOS 和 QuickBird数据的全色波段和多光谱波段进行融合，取得了良好的效果 ( 图4 －3、图4 －4) 。

图 4 －3 IKONOS 融合影像 ( 1m)

图 4 －4 QuickBird 融合影像 ( 0. 6m)

( 四) 图像配准

震前、震后 DOM 影像在经过无控制点共线方程几何校正和正射校正之后，已基本实现了叠置，大部分地物都能比较好地重叠在一起，不过也有个别目标存在偏差现象。图 4 －5 中左侧为震后 QuickBird 影像，右侧为 IKONOS 影像，中间黑色划线标注区域的池塘重叠效果存在着偏差，在这种情况下，需要进行影像之间的配准。图像配准也叫影像的精校正，是指消除图像中的几何变形，产生一幅符合某种地图投影或图形表达要求的新图像过程。

图 4 －5 震前、震后 DOM 影像叠置的效果( 黑色划线区域内存在偏差)

本节所指的图像配准是多图像的几何配准，多图像是指同一地区不同时刻的图像 ( 多时相图像) 或者不同传感器获取的多源图像，这里 IKONOS 和 QuickBird 影像就属于多源多时相遥感影像。多图像几何配准就是指将多图像的同名像点通过几何变换实现精确重叠，通常称为相对配准如果将相对配准后得到的多图像归入同一地图坐标系统当中，就叫做绝对配准。

本研究中震前与震后的融合 DOM 影像之间的配准工作在 EＲDAS 软件的 Image Geometric Correction 模块中进行，以震后 QuickBird 影像 DOM 为参照，选择二次多项式校正模型配准震前 IKONOS 影像 DOM，手动选取了 6 个同名控制点建立了多项式模型之后，这时 EＲDAS 软件会根据模型自动找出后面控制点在图像中对应的位置，这时只要在图像窗口中校正其位置即可，节省了时间。同名配准控制点的残差应当满足表 4 － 2 的要求。共选取 20 个地面控制点 ( GCP) ，这 20 个控制点总的均方根误差 ( ＲMSE) 为 ±1. 0773，各个地面控制点坐标值和 ＲMSE 见表4 －3。最后选择最近邻域法 ( Nearest Neighbor) 对影像进行重采样。

表 4 －2 配准控制点残差

表 4 －3 各个地面控制点的坐标和 ＲMSE

续表

DOM 影像几何配准之后，需要对质量进行检查控制。首先校正震前、震后 DOM 影像的同名地物点的点位中误差不能大于表 4 －4 的规定，另外检查两个时相的 DOM 影像与土地利用现状图 ( LUDＲG) 的匹配情况，精度不能大于表 4 － 5 的规定。如果不满足要求，还需要利用土地利用现状图对两时相 DOM 影像进行第二次配准。

表 4 －4 多时相 DOM 同名地物点配准精度

表 4 －5 DOM 相对于土地利用现状图的精度

( 五) 影像辐射增强处理

由于震前、震后 DOM 影像获取的时间不同，地面接收到的太阳辐射度不同，加之高分辨率遥感影像本身像元间光谱的异质性较强，使得两个 DOM 影像在外观上肯定存在一些差异，对震害变化信息检测带来不利影响。为了消除这些不利影响，提高震害信息提取精度，需要事先对震前、震后两时相 DOM 影像做辐射增强处理，主要包括自适应性滤波处理和直方图匹配。

为了控制高分辨率遥感数据中的随机噪声 ( 随机噪声往往影响地类之间的均匀性及边界的稳定性) 和像元间光谱的强异质性，需要对影像进行空间滤波处理。本研究采用 ENVI 软件中的 Frost 自适应性滤波对震前、震后 DOM 影像进行滤波处理，在降低像元光谱异质性，使影像平滑的同时，较好地保持了地类边缘和纹理的清晰。Frost 自适应性滤波是以权重为自适应调节参数的滤波器，对每一个像元都确定一个权重，然后逐个进行滤波。

本研究所说的直方图匹配是指对图像查找表进行数学变换，使一幅多光谱遥感图像所有波段的直方图与另一幅遥感图像所有对应波段相似，其经常用于相邻图像之间的拼接或者多时相遥感图像动态变化信息检测研究的预处理工作，经过直方图匹配可以消除部分由于太阳高度角或者大气辐射造成的多源遥感影像间光谱信息的差异性 ( 党安荣等，2003) 。

本研究在 EＲDAS 软件中以震后 QuickBird 影像 DOM 为标准，对震前 IKONOS 影像DOM 的各个波段完成直方图匹配处理。从上节融合后的结果 ( 图 4 － 3、图 4 － 4) 中可以发现，QuickBird 融合影像在研究区域内存在一片厚厚的云层以及由云层引起的阴影，导致云层和阴影范围内的信息完全丢失，严重影响了数据质量。在这种情况下，使用去除薄云的处理方法无法解决问题，然而也不能使用去除厚云的替补方法，因为后续工作是震害变化信息的提取，如果用其他的影像数据进行替换 ( 例如灾前 IKONOS 数据) ，势必会影响到后面变化信息提取的精度。经过综合考虑，决定在不能修复厚云及其阴影区域的情况下，通过分类单独提出云层和阴影，然后建立掩膜图层，在 QuickBird 影像上剔除上述区域，不参与后续研究。同样在震前 IKONOS 影像上也剔除掉相同的区域。最终经过辐射增强处理以及剔除厚云和阴影部分的前、后两时相 DOM 影像效果见图 4 －6 和图 4 －7。

图 4 －6 IKONOS 最终 DOM 影像 ( 1m)

图 4 －7 QuickBird 最终 DOM 影像 ( 0. 6m)

类似于采取卷及滤波的方法将较大类别中的虚假像元归到该类中，定理一个变换核尺寸，主要分析用变换核中占主要地位的像元类别代替中心像元类别；次要分析是将变换核中占次要地位的像元的类别代替中心像元的类别。

工具位置：Toolbox/Classification/Post Classification/

空间连续性不太好的可以使用

聚类处理解决了这个问题。首先将被选的分类用一个膨胀 *** 作合并到一块，然后用变换核对分类图像进行腐蚀 *** 作

工具位置：Toolbox/Classification/Post Classification/Clump Classes

解决图像中的孤岛问题，通过分析周围的4个或8个像元，判定一个像元是否与周围的像元同组。如果一类中被分析的像元小于输入的阈值，这些像元就会被从该类中删除，删除的像元归为未分类的像元。【用的很少】

工具位置：Toolbox/Classification/Post Classification/Sieve Classes

工具位置：Toolbox/Classification/Post Classification/class statistics

（1）打开分类结果和原始影像——"\13.分类后处理\数据\can_tmr_class.dat"和"can_tmr.img"；

（2）打开分类统计工具，路径为Toolbox/Classification/Post Classification/Class Statistics，在d出对话框中选择"can_tmr_class.dat"，点击OK；

（3）在Statistics Input File面板中，选择原始影像"can_tmr.img"，点击OK；

（4）在d出的Class Selection面板中，点击Select All Items，统计所有分类的信息，点击OK；

注：可根据需要只选择分类列表中的一个或多个类别进行统计。

（5）在Compute Statistics Parameters面板可以设置统计信息（如下图所示），按照图中参数进行设置，点击Report Precision…按钮可以设置输入精度，按默认即可。点击OK；

为显示统计结果的窗口，统计结果以图形和列表形式表示。从Select Plot下拉命令中选择图形绘制的对象，如基本统计信息、直方图等。从Stats for标签中选择分类结果中类别，在列表中显示类别对应输入图像文件DN值统计信息，如协方差、相关系数、特征向量等信息。在列表中的第一段显示的为分类结果中各个类别的像元数、占百分比等统计信息。

工具位置：Toolbox/Classification/Post Classification/Overlay Classes

将各种类别叠加在一幅RGB彩色合成图或灰度图像上，从而生产一幅RGB图像，如果要想得到较好的效果，在叠加前背景图像经过拉伸并保存为字节型图像

（1）打开拉伸工具（Toolbox/Raster Management/ Stretch Data ），在d出的对话框中选择"can_tmr.img"文件，然后点击下方的Spectral Subset（如下图所示），在d出面板中选择波段1、2、3，（真彩色）点击OK；

图：选择拉伸文件和波段选择。在Data Stretching面板中，按照下图进行参数设置，点击OK即可

工具位置：Toolbox/Classification/Post Classification/Classification to Vector

将一定范围内像元都并入其他一个类别中

（1）在Interactive Class Tool面板中，选择Edit >Mode: Polygon Add to Class；

（2）在Interactive Class Tool面板中，鼠标左键单击"Unclassified"前面的方型色块，让"Unclassified"类别处于激活状态；

（3）选择一个编辑窗口：Image，在Image窗口中绘制多边形，多边形以内的像元全部归于"Unclassified"一类。

将一定范围内某一类像元并入其他一类中

(1)在Interactive Class Tool面板中，选择Edit >Mode: Polygon Delete from Class；

（2）选择Edit >Set delete class Value，选择并入的目标类，

（3）在Interactive Class Tool面板中，鼠标左键单击"Unclassified"前面的方型色块，让"Unclassified"类别处于激活状态。

（4）选择一个编辑窗口：Image，在Image窗口中绘制多边形，多边形以内的类别"Unclassified"全部归于"林地"。

注：浏览修改效果，需要勾选林地，在视图中显示林地这一类。

（5）选择Edit >Undo Changes，可以取消修改，选择File >Save changes to File，可以将修改结果保存。

总体分类精度

等于被正确分类的像元总和除以总像元数。被正确分类的像元数目沿着混淆矩阵的对角线分布，总像元数等于所有真实参考源的像元总数，如本次精度分类精度表中的Overall Accuracy = (2258/2346) 96.2489%

Kappa 系数

它是通过把所有真实参考的像元总数（N）乘以混淆矩阵对角线（XKK）的和，再减去某一类中真实参考像元数与该类中被分类像元总数之积之后，再除以像元总数的平方减去某一类中真实参考像元总数与该类中被分类像元总数之积对所有类别求和的结果。Kappa计算公式如下：

错分误差

指被分为用户感兴趣的类，而实际属于另一类的像元，它显示在混淆矩阵里面。本例中，总共划分为林地有441个像元，其中正确分类418个，23个是其他类别错分为林地（混淆矩阵中林地一行其他类的总和），那么其错分误差为23/441=5.22%。

漏分误差

指本身属于地表真实分类，当没有被分类器分到相应类别中的像元数。如在本例中的林地类，有真实参考像元419个，其中418个正确分类，其余1个被错分为其余类（混淆矩阵中耕地类中一列里其他类的总和），漏分误差1/419=0.24%

制图精度

是指分类器将整个影像的像元正确分为A类的像元数（对角线值）与A类真实参考总数（混淆矩阵中A类列的总和）的比率。如本例中林地有419个真实参考像元，其中418个正确分类，因此林地的制图精度是418/419= 99.76 %。

用户精度

是指正确分到A类的像元总数（对角线值）与分类器将整个影像的像元分为A类的像元总数（混淆矩阵中A类行的总和）比率。如本例中林地有418个正确分类，总共划分为林地的有441，所以林地的用户精度是418/441=94.78%。

监督分类中的样本选择和分类器的选择比较关键。在样本选择时，为了更加清楚的查看地物类型，可以适当的对图像做一些增强处理，如主成分分析、最小噪声变换、波段组合等 *** 作，便于样本的选择；分类器的选择需要根据数据源和影像的质量来选择，比如支持向量机对高分辨率、四个波段的影像效果比较好。

ENVIMAG/VLF地球物理系统是由加拿大Scintrex公司生产的磁测系统，其外貌如图4-5-13所示。它具有多项测量功能，包括磁力仪、磁梯度仪测量和甚低频测量。磁力仪测量方式时，不但可定点测量还可以采用连续行走方式测量地磁总场强度。其探头可以安装在背架上以更快的速度测量。它还可以自身做日变校正。当采用磁梯度仪测量方式时，它可同时提供总场和梯度数据。将两个探头固定在背架上可以测量地磁场的垂直梯度，这样消除了日变的影响。甚低频测量方式观测来自三个不同甚低频发射点上感应出来的二次场和总场的甚低频信号相位和相位差。唯一的三线圈探头为一个全方位探头，因此不必选择探头方向，提高了测量速度。它也可以利用两个或三个探头来测量甚低频感应电场，最后给出视电阻率和相位角度。它还可将已有的测量结果以剖面图，数据表等形式回放检查。它是由探头、主机、探杆、充电器、电缆、背架等组成。ENVI型磁力仪是与MP-4质子磁力仪同类型新的一种带微处理机的高分辨率仪器，其测地磁总场的原理与CZM-3型相同。

图4-5-13 ENVIMAG型质子旋进磁力仪外貌

(一)仪器的结构原理及应用范围

ENVI系统使您能灵活地寻找越来越难以发现的异常体。一整套ENVI系统经济、轻便，是具有VLF功能的便携式质子旋进式磁力仪/梯度仪(图4-5-13)，可确保您进行大面积快速、准确地测量，并能以图标形式实时监控数据质量，避免无效工作能使您寻找越来越难以发现的异常体。它适用于:①环境场地特征②地下水勘探③矿产勘探④考古学⑤地质填图⑥地下水研究⑦土工技术研究⑧土木工程。

(二)主要技术指标

应用地球物理教学实习指导

在2s采样率时“行走”模式

0.5s循环率

本键盘共有17个键，其中有最常用的键分别安置在仪器的面板的左侧和右侧(图4-5-14)，即开始/停止“START/STOP”键和记录键“RECORD”。有些键还有三重功能，各种功能体现在不同工作进程之中。

图4-5-14 ENVIMAG型质子旋进磁力仪控制台面板

·ON/OFF开启和关闭仪器。

·STAPT/STOP/CLEAR开始或停止某种 *** 作。

*数据采集数据回放

*查看数据

*当仪器处于文字信息显示状态时，可删除输入的内容。

·SETUP/5/MNO(设置)进入各种状态设置显示，所显示的设置菜单随当前的显示屏内容而定。只要按动此键就会出现所处的设置状态。

·AUX/1/ARC(辅助)进入该辅助功能键完成下述功能。

*设定液晶显示屏的亮度

*数据的输出

*锁定所设置的参数

*重新编制主系统软件。

·←/2/DEF→/3/GHI上、下、左、右移动光标。

·ENTER/RATE(进入)键盘上有两个类似的键，分别位于面板左、右角。

*打开和关闭设置状态期间的各项参数的选择范围

*打开和关闭图形显示期间绘图比例因子的选择范围

*在连续行走测量方式下选择坐标记的速率。

(三)ENVI磁测系统 *** 作键盘

·ESCAPE/0(退回)可从当前菜单退回到上一级菜单。最上一级菜单总是主 *** 作菜单显示。

*按此键还可以中断数据输出。

·INFO/7/STU(信息)按此键进入信息显示，在磁场测量状态下它具有不同的功能:

*在基本状态下和搜索状态下可进入帮助菜单显示

*高级状态下，它具有下述功能

*设置年、月、日和小时、分、秒

*输入仪器号和工作号以及 *** 作员识别号

*查看内存量。

·NEXT(下一个)键盘上有两个类似的键，分别位于面板左右两侧，其功能如下:

*按页顺序滚动全部数字数据显示

*按页顺序滚动全部图形数据显示

*将光标从当前页进入到下一页

*进入文字显示期间，将光标依次顺着字符移动

*在步行连续测量期间通过点号增量选择点号。

NOTE/6/FQR(注释)进入文字显示，其功能如下:

*利用所选择的读数可进入五种普通的重复记录信息

*利用特殊的读数进入所记录的唯一信息。

·NUM/GRA/4/JKL数字图形。

*只限于数据采集期间对数字和图形数据显示时做标记用。

RECALL/9/YZ(回放)。

*回放浏览数据显示选择

*显示所要回放浏览数据的类型

*设置所要回放浏览数据的开始位置“点号”及时间。

·RECORD(记录)对所测量到的数据进行手动记录并提示内存容量。

·TIE-PT/8/VWX(连接点校正)在进行连接点测量时，该键的作用如同START键，即测量键。该键只有在进行测线观测过程中遇到连接点时测量用，并还在利用上述连接点的数据进行循环状态日变校正时用。

·+-具有向前或向后反复滚动菜单表的功能:

*给数字加符号

*在回访浏览所测数据期间可按某条测线显示数据

*在显示图形时可进行比例尺选择

*在进行逐点测量时可增加或减小点线号。

·1-9，A-Z对文本输入和功能设置时都可输入字母和数字。

·Beeper/Port仪器面板左上角的空白键是一个目前还没有实际功能的键，但该键表面柔韧性的膜进一步增强蜂鸣器的音量。

·ON/OFF+AUX/LCD/1/ABC同时按下这两个键可进行冷启动 *** 作，使仪器重新回到出厂时的参数设置。

·ON/OFF+SETUP/5/MNO同时按下这两个键，可以完成下述功能:

*在总磁场测量(MAG)方式、甚低频测量(VLF)方式和总磁场与甚低频同时测量(MAG＆VLF)这三种方式任选一项

*在总磁场测量(MAG)方式下，从三种测量状态中任选一种，例如基本状态(Basic)、搜索状态(Search)或高级状态(Advanced)。

(四) *** 作过程

(1)在每个不同的测量之前，按“SETUP”和“ON”两个键，初始化仪器，进入特征菜单(图4-5-15):按“1”，“2”或“3”键选择要进入的测量状态。

图4-5-15 ENVI测量主菜单

(2)按“1”键进入磁力仪工作状态菜单，该菜单占用两个显示屏，按“+”、“-”切换，如图4-5-16所示。

图4-5-16 ENVI总场状态下应用菜单

(3)基本状态。当按“1”，“2”或“3”键时进入基本测量状态，显示如图4-5-17所示。

图4-5-17 ENVI基本测量状态菜单

在任何时候，按“NUM/GRA”键显示当前测线的最后178个点的剖面图，按“NEXT”后显示最后一个数据的进动衰减信号，再按“NUM/GRA”回到数据显示。

(4)搜索状态。从(2)中选择“4”时，显示如图4-5-18所示。

图4-5-18 ENVI搜索状态菜单

设置完毕后，按“START”键开始测量。

在任何时候，按“NUM/GRA”键显示当前测线的最后178个点的剖面图，按“NEXT”后显示最后一个数据的进动衰减信号，再按“NUM/GRA”回到数据显示。

(5)高级状态。从(2)中选择“5”，“6”或“7”进入高级状态，显示如图4-5-19所示。

图4-5-19 ENVI高级状态菜单

按“SETUP”键进入仪器功能设置显示，如图4-5-20所示。

图4-5-20 ENVI高级状态下仪器设置菜单

设置完毕后，按“START”键开始测量。按“SETUP”键进入磁力仪设置显示，如图4-5-21所示。

图4-5-21 ENVI高级状态下磁力仪设置设置菜单

设置完毕后，按ESC回到(5)，按“START”键开始测量。

在任何时候，按“NUM/GRA”键显示当前测线的最后178个点的剖面图，按“NEXT”后显示最后一个数据的进动衰减信号，再按“NUM/GRA”回到数据显示。

测量完毕后，下载数据并连接基站和野外测量磁力仪，进行日变校正:

测量前或结束后，按“AUX/LCD”键，进入附加功能显示，在这里，你可以调节显示器亮度、输出数据、进行因子测试、锁定参数设定功能和重装系统。

测量前或结束后，按“NOTER”键，进入注释功能显示，在这里，你可以输入数据文件注释。

测量前或结束后，按“INFO”键，进入信息功能显示，在这里，你可以输入数据文件的信息:测量时间、设备序列号、测量日期、 *** 作员姓名和 *** 作号。

测量前或结束后，按“RECALL”键，进入数据回放显示。

(6)这是甚低频工作状态，从(1)中按“2”进入 *** 作显示，如图4-5-22所示。

图4-5-22 甚低频测量主菜单

按“SETUP”键进入仪器设置显示，如图4-5-23所示。

图4-5-23 ENVI甚低频模式下仪器设置菜单

按“ENTER”键进入(或从主菜单中按“ENTER”进入)甚低频设置菜单，如图4-5-24所示。

图4-5-24 ENVI甚低频模式下测量参数设置菜单

设置完毕后，按ESC回到(6)，按“START”键初始化显示，初始化完毕后自动进入测量状态。

测量过程中，如需要查看正在保存的单个频率的数据，按“1”，“2”或“3”键。

Q框显示 *** 作者质量数，显示在5-6为优。

S框显示信号噪声比数，显示在9为最优。

在任何时候，按“NUM/GRA”键显示当前测线的最后178个点的剖面图，按“NEXT”后显示最后一个数据的进动衰减信号，再按“NUM/GRA”回到数据显示。

数据的回放、下载与磁力仪测量类似。

(7)如果想同时测量VLF和磁场，从(1)中按“3”键，显示如图4-5-25所示。

图4-5-25 ENVI同时测量甚低频和磁场主菜单

按“SETUP”键进入仪器设置显示，如图4-5-26所示。

图4-5-26 ENVI同时测量甚低频和磁场模式下仪器设置菜单

设置完毕后，按ESC回到(7)，按“START”键开始测量

测量过程中或完毕后，数据的显示、回放、下载和校正与前面的单独测量类似。

---