一种基于鲁棒性卫星数据分析技术的地震红外辐射异常特征多角度遥感提取方法
【专利摘要】本发明提出一种基于鲁棒性卫星数据分析技术的地震红外辐射异常特征多角度遥感提取方法。步骤1)剔除遥感影像上那些受到云影响的观测区域；步骤2)对卫星传感器接收到的大气层顶的信号进行大气纠正；步骤3)利用水陆分离掩膜分别对陆地和海洋的红外多角度遥感图像进行处理，计算单幅红外影像像元差值；步骤4)构建红外背景参考场作为地震红外辐射异常特征多角度遥感提取的参考背景；步骤5)针对单一角度红外数据，提出地震红外辐射异常特征多角度遥感异常鲁棒性测算指标，用于单一角度地震红外辐射特征的提取；步骤6)基于红外多角度遥感数据，提出多角度亮温合成地震红外辐射分析方法，用于地震红外辐射特征多角度遥感的合成和分析。
【专利说明】一种基于鲁棒性卫星数据分析技术的地震红外辐射异常特 征多角度遥感提取方法
【技术领域】
[0001] 本发明属于地震红外辐射异常特征多角度遥感研究领域，涉及一种基于鲁棒性卫 星数据分析技术的地震红外辐射异常特征多角度遥感提取方法。
【背景技术】
[0002] 卫星红外观测以其大视野、高精度、周期短和不受地面条件限制等诸多优势成为 研究地震红外辐射异常特征很有前途的遥感观测手段。近20年来，随着卫星遥感技术的发 展，红外遥感数据越来越丰富，为地震红外遥感研究提供了保证，地震红外遥感研究工作取 得了很大的进展。与此同时，红外多角度遥感机理、建模以及多角度遥感技术方面的研究取 得了很大的进步，红外多角度遥感技术与应用研究相互促进、相互发展，已成为重要的前沿 研究方向，但是红外多角度遥感在地震科学的应用研究仍属空白，因此，发展地震红外辐射 异常特征多角度遥感识别与提取技术，探索红外多角度在地震红外辐射特征研究中的适用 性，将为地震红外遥感的研究提供了一种新的技术和方法。
[0003] 目前，尽管国内外专家开展了大量的震例研究工作，不断探索震前热红外遥感的 异常征兆，也取得了一定进展，但某些主要技术环节尚需深入研究。例如，现有的地震红外 异常信息提取方法主要有"涡度"计算法、"距平"计算法和小波时频分析方法。这些方法在 提取和突出地震异常信息上有一定的效果，但这些方法本身都是数学处理方法，存在着一 些缺陷，计算时用到的数据不同就会得到不同的结果，不确定性大。因此，有必要引进或研 究新的方法，解决如何有效地在热红外辐射数据中提取地震异常信息。
[0004] 在过去的几年里，鲁棒性卫星数据分析技术（Robust Satellite data analysis Technique)已经被成功的应用于监测重大的自然灾害和环境灾害，例如：火山、水灾、森林 火灾等等。鲁棒性卫星数据分析技术是基于多年多时相热红外卫星遥感数据的初步分析， 专门用于表征热红外遥感卫星图像中像素信息，因此对热红外遥感地震异常信息的提取更 加有效。
[0005] 本研究发展了一种基于鲁棒性卫星数据分析技术的地震红外辐射异常特征多角 度遥感提取方法，该方法考虑了云、陆地与水体之间不同辐射特性和气候、季节变化等影响 因素，能够消除由于季节交替或气候突变所引起的地面增温对异常判别和提取的影响；并 基于多角度多年背景参考场，研究利用红外异常鲁棒性测算指标检测地震红外辐射异常信 息，而且检测出来的红外辐射异常信息可以利用信噪比来刻画，与地震相关的噪声信息越 大，地震信息检测指标值越小，反之亦然。所以，红外异常鲁棒性测算指标能够很好避免非 震引起的热红外异常，并且在地震红外辐射异常特征多角度遥感识别和提取上具有很好的 鲁棒性。
【发明内容】
[0006] 针对现有技术中存在的问题，本发明的目的在于提出和发展一种基于鲁棒性卫星 数据分析技术的地震红外辐射异常特征多角度遥感提取方法，用于准确提取和识别地震红 外辐射多角度遥感异常特征。
[0007] 为了实现上述目的，本发明提出了一种地震红外辐射异常特征多角度遥感异常鲁 棒性测算指标，用来检测地震红外辐射异常信息，而且检测出来的红外辐射异常信息可以 利用信噪比来刻画，与地震相关的噪声信息越大，地震信息检测指标值越小，反之亦然。所 以，该指标能够很好避免非震引起的热红外异常，并且在地震红外辐射特征检测上具有很 好的鲁棒性。具体包括如下步骤：
[0008] 步骤1)剔除影像上受到云影响的观测区域，以确保是在完全晴空条件下对影像 进行处理；
[0009] 步骤2)对卫星传感器接收到的大气层顶的信号进行大气纠正；
[0010] 步骤3)利用水陆分离掩膜分别对陆地和海洋的红外多角度遥感图像进行处理， 计算单幅红外影像像元差值；
[0011] 步骤4)构建红外背景参考场作为地震红外辐射异常特征多角度遥感提取的参考 背景（明确如何构建红外背景参考场），以便于排除非震因素引起的红外辐射异常现象；
[0012] 步骤5)针对单一角度红外数据，提出地震红外辐射异常特征多角度遥感异常鲁 棒性测算指标（补充具体测算方法），用于单一角度地震红外辐射特征的提取；
[0013] 步骤6)基于红外多角度遥感数据，提出多角度亮温合成地震红外辐射分析方法 (补充具体测算方法），用于地震红外辐射特征多角度遥感的合成和分析。
[0014] 进一步，所述步骤4)具体方法为
[0015] 步骤a)将多年同一区域的大量红外多角度影像数据投影到统一的标准网格上， 生成标准网格化数据，保证数据的统一性；步骤b)利用无云的，完成重采样的红外多角度 遥感影像，来建立红外背景均值参考场和红外背景标准差参考场。
[0016] 以下结合附图与实施例对本发明作进一步的说明。
【专利附图】
【附图说明】
[0017] 图1 2010年3月份玉树地区红外背景参考场。
[0018] 图2 2010年4月份玉树地区红外背景参考场。
[0019] 图3 2010年5月份玉树地区红外背景参考场。
[0020] 图4基于AATSR星下点方向的玉树地震时空演化异常特征。
[0021] 图5基于AATSR前向视角方向的玉树地震时空演化异常特征。
[0022] 图6多角度亮温玉树地震红外辐射时空演化异常特征。
【具体实施方式】
[0023] 下面结合附图，以红外多角度AATSR卫星数据为例，并且以玉树地震为例，对本发 明基于鲁棒性卫星数据分析技术的地震红外辐射异常特征多角度遥感提取方法作进一步 说明。
[0024] -、单幅红外影像像元差值计算
[0025] 红外影像像元差值计算的目的是为了消除由于气候突变或者季节交替所引起的 地表温度升高对红外辐射异常判别提取分析的影响。说明以AATSR数据为例对红外图像 进行了云消除处理和大气纠正，并利用水陆分离掩膜分别对陆地和海洋的红外图像进行处 理，具体流程如下：
[0026] (1)构建云掩膜影像
[0027] 利用前向和星下点云测试算法（nadir-forward cloud tests):此云检测算法与 双角度AATSR亮度温度之间的差异相关，在星下点和前向视角方法11 μ m处测量的亮度温 度差值，与同一视角方向的11 μ m到12 μ m通道亮温差值不同，因为这主要是因为两个方向 视角的大气传输路径不一致，还有大气水气的影响。我们在这里给出一个云识别的阈值，如 果实际像素的差值T 11-T12超过阈值，则认为此图像像素有云。通过云探测算法，构建云掩膜 影像CM (x，y)，对于影像的每个像元点（x，y):无云的情况，设置像素标记CM= 1 ;在有云的 情况，设置像素标记CM = 0。
[0028] (2)大气纠正
[0029] AATSR数据在经过云检测处理后，利用分窗算法被用来得到图像中的大气水汽含 量。得到水汽含量后加入AATSR热红外波段数据可以计算出地表的两个角度的亮度温度。
[0030] (3)构建水陆分离影像
[0031] 利用AATSR LlB的辅助产品，结合像元的经纬度，对像元进行下垫面土地覆盖类型 识别，进行水陆识别，并结合土地覆被分类图和AATSR LlB的辅助产品，构建水陆分离影像 LS (X，y)，对于影像像素的每个像元点（x，y)，当LS = 1时，标记像素为海表，当LS = 0时， 标记像素为陆表。
[0032] (4)单景AATSR图像基准值计算
[0033] 对于无云的海水像素，即影像的每个像元点（X，y)，当CM = 1和LS = 1时，求和 y
【权利要求】
1. 一种基于鲁棒性卫星数据分析技术的地震红外福射异常特征多角度遥感提取方法， 该方法具体包括如下步骤： 步骤1)剔除影像上受到云影响的观测区域，W确保是在完全晴空条件下对影像进行 处理； 步骤2)对卫星传感器接收到的大气层顶的信号进行大气纠正； 步骤3)利用水陆分离掩膜分别对陆地和海洋的红外多角度遥感图像进行处理，计算 单幅红外影像像元差值； 步骤4)构建红外背景参考场作为地震红外福射异常特征多角度遥感提取的参考背 景，W便于排除非震因素引起的红外福射异常现象； 步骤5)针对单一角度红外数据，提出地震红外福射异常特征多角度遥感异常鲁棒性 测算指标，用于单一角度地震红外福射特征的提取； 步骤6)基于红外多角度遥感数据，提出多角度亮温合成地震红外福射分析方法（补充 具体测算方法），用于地震红外福射特征多角度遥感的合成和分析。
2. 如权利要求1所述的方法，其特征在于，所述步骤4)具体方法为： 步骤a)将多年同一区域的大量红外多角度影像数据投影到统一的标准网格上，生成 标准网格化数据，保证数据的统一性；步骤b)利用无云的，完成重采样的红外多角度遥感 影像，来建立红外背景均值参考场和红外背景标准差参考场。
【文档编号】G06T7/00GK104463835SQ201410444133
【公开日】2015年3月25日 申请日期:2014年9月3日 优先权日:2014年9月3日
【发明者】熊攀, 孟庆岩, 顾行发, 申旭辉, 卢显 申请人:中国科学院遥感与数字地球研究所