一种基于H8/AHI的白天雾快速提取的方法与流程

本发明涉及环境监测与气象预报领域，特别涉及一种基于h8/ahi的白天雾快速提取的方法。

背景技术：

雾是一种灾害性天气现象，它会降低能见度、恶化空气质量，对交通安全尤其是航海、航空和高速公路运输等方面有较严重影响。随着卫星遥感技术的快速发展，相对传统雾检测，遥感技术雾检测具有明显的优势，主要体现在数据更新快、时效性高、检测范围大和技术成本低等方面。基于遥感技术雾检测主要使用极轨卫星数据(如terra/modis、noaa/avhrr)，极轨卫星数据具有光谱信息丰富的特点，但过境时间(白天上午10:30、白天下午2:30、晚上22：30和凌晨2：30)与雾生成期有一定的延迟，时间分辨率(一天获取2-4幅数据)较低，难以满足雾检测近实时快速响应的要求。静止卫星时间分辨率高，可对研究区域连续观测，尤其是生命周期短、变化快的雾,但以往的静止卫星(gms-5、meteosat8和mtsat-1r等)空间分辨率和光谱分辨率较低，难以进行雾的精准检测。日本新一代静止气象卫星葵花8号(himawari-8,简称h8)主载荷(ahi)具有全盘和区域扫描能力，能在10分钟内完成全盘扫描，成像仪通道从0.46μm-13.3μm共16个通道，空间分辨率最高达500m，相比于早期静止卫星，其时间分辨率(10分钟)，光谱分辨率(16个通道)和空间分辨率(500m)均有较大提高。h8/ahi数据为雾变化检测提供了一个良好的数据源^【1】,同时也对雾近实时快速提取技术提出了更高的要求。

由于受到遥感数据时间分辨率、资金、设备等因素的影响，纵观国内外雾检测研究现状，雾近实时准确快速提取仍有一定困难。相比以往静止气象卫星，h8/ahi数据独特的特性为雾近实时快速提取提供了良好的数据源，因此发明一种基于h8/ahi的白天雾检测算法减少因雾造成的财产和生命损失是十分必要和迫切的。

技术实现要素：

本发明所要解决的技术问题是，针对现有技术不足，提供一种基于

h8/ahi的白天雾快速提取的方法，能快速、准确地检测白天雾，减少因雾造成的财产和生命损失。

为解决上述技术问题，本发明所采用的技术方案是：步骤1：利用多天可见光数据合成晴空底图ir，利用反射率差异，阈值法去除晴空地表，获取去除地表后的影像iscf；

步骤2：归一化雪指数法去除影像iscf中的雪和不透明冰云，获取去除后的影像icf；

步骤3：去除影像icf中纹理复杂且快速运动的云，获取雾检测影像iitwcf；

步骤4：去除影像iitwcf中冰云和薄卷云，获取去除后影像iwcf；

步骤5：去除影像iwcf中水云，获取去除后影像ic2f；

步骤6：利用多天红外数据合成晴空底图ibt，阈值去除影像ic2f中的低层云,获取卫星雾检测结果if；

所述的利用多天可见光数据合成晴空底图ir，阈值法去除晴空地表，获取去除地表后的影像iscf，包括以下步骤：

利用多天可见光数据合成晴空底图，使用当天及前9天当前时刻和前一时刻可见光ib3数据，结合可见光波段晴空地表反射率较云雾最低这一特征，分别合成晴空底图ir-现和ir-前；

对当前时刻ib3-现数据与晴空合成底图ir-现做差，前一时刻b3数据ib3-前与晴空合成底图ir-前做差，分别获得差值影像i差值-现和i差值-前；

对差值影像采用固定阈值(阈值设置为0.05)快速去除地表，分别与原始影像ib3-现和ib3-前相乘获取去除地表后的影像iscf-现和iscf-前；其中iscf-现表示使用当前时刻数据获取的去除晴空地表后的影像iscf；iscf-前表示使用前一时刻数据获取的去除地表后影像iscf；

所述的归一化雪指数法去除影像iscf中的雪和不透明冰云，获取去除后的影像icf，包括以下步骤：

分别使用当前时刻ib3-现、ib5-现和前一时刻ib3-前、ib5-前数据，获取当前时刻和前一时刻归一化雪指数ndsi现和ndsi前，ndsi现＝(ib3-现-ib5-现)/(ib3-现+ib5-现)；ndsi前＝(ib3-前-ib5-前)/(ib3-前+ib5-前)；

对去除地表后的影像iscf-现和iscf-前使用归一化雪指数法，设置阈值0.4去除影像iscf-现和iscf-前中雪和不透明冰云，分别与影像iscf-现和iscf-前相乘得到影像icf-现和icf-前；其中icf-现表示使用当前时刻数据获取的去除雪和不透明冰云的影像icf；icf-前表示使用前一时刻数据获取的去除雪和不透明冰云的影像icf。

所述的去除影像icf中纹理复杂且快速运动的云，获取雾检测影像itwcf，包括以下步骤：

对当前时刻icf-现数据和前一时刻icf-前数据的b3波段做比值运算，获取比值影像i比值，i比值＝icf-现/icf-前；

对比值影像i比值采用一阶概率纹理滤波法获取datarange(数据范围)纹理数据idr；

设置阈值0.3去除影像icf-现中运动的及类型变化较大的云,与影像icf-现相乘得到影像iitwcf；

所述的去除影像iitwcf中冰云和薄卷云，获取去除后影像iwcf，包括以下步骤：

使用当前时刻ib13-现(中心波段：10.4um)数据，设置阈值230k去除影像iitwcf中冰云，与影像iitwcf相乘得到影像itwcf；

对当前时刻ib11-现(中心波段：8.7um)数据和ib14-现(中心波段：11.2um)数据做差值运算，获取差值影像i差值1，i差值1＝ib11-现-ib14-现；

设置阈值0k去除影像itwcf中薄卷云，与影像itwcf相乘得到影像iwcf；

所述的去除影像iwcf中水云，获取去除后影像ic2f，包括以下步骤：

对当前时刻ib15-现(中心波长：12.3um)数据和ib11-现数据做差值运算，获取差值影像i差值2，i差值2＝ib15-现-ib11-现；

其阈值依赖太阳高度角θ的变化，所以归一化处理获取动态阈值vt，

使用归一化阈值vt去除影像iwcf中水云，与影像iwcf相乘得到影像ic2f；

所述的利用多天红外数据合成晴空底图ibt，阈值法去除影像ic2f中的低层云,获取卫星雾检测结果if，包括以下步骤：

利用热红外波段中，晴空地表亮度温度大于云雾的亮度温度这一特征，使用当天及前9天当前时刻ib14-现数据合成晴空底图ibt；

对当前时刻ib14-现数据与晴空底图ibt做差，获取差值的绝对值影像i差值3,i差值3＝ib14-现–ibt；

采用阈值法(设置阈值11)快速去除影像ic2f中低层云，与影像ic2f相乘获取最终卫星雾检测结果if；

与现有技术相比，本发明所具有的有益效果为：本发明是当前具有高时间分辨率的静止卫星数据的条件下的一种可准确、高效快速实现白天雾提取的方法，检测结果可靠，精度高。

附图说明

图1示出了基于h8/ahi的白天雾检测算法流程；

图2示出了卫星雾检测结果与地面观测数据叠加图；其中，(a)8:00卫星雾检测结果与地面观测数据叠加图；(b)14:00卫星雾检测结果与地面观测数据叠加图；注：黑色区域为卫星雾检测结果；不同形状表示地面站点不同观测结果：▲：特浓雾浓雾、●：雾、×：非雾；

图3示出了雾检测结果与h8/ahi假彩色影像叠加图；其中，(a)8:00雾检测结果；(b)9:00雾检测结果；(c)10:00雾检测结果；(d)11:00雾检测结果；(e)12:00雾检测结果；(f)13:00雾检测结果；(g)14:00雾检测结果；(h)15:00雾检测结果；(i)16:00雾检测结果；

图4式出了可见光波段雾和地表反射率差异图。

具体实施方式

以下结合附图对本发明进行详细说明。

1、利用多天可见光数据合成晴空底图，阈值法去除晴空地表，获取去除地表后的影像iscf；

实验使用2015年11月16日-2015年12月1日h8/ahi卫星2km空间分辨率10分钟时间分辨率的时序影像ib3数据，分别合成晴空底图ir-现和ir-前；

对当前时刻ib3-现数据与晴空合成底图ir-现做差，前一时刻ib3-前数据与晴空合成底图ir-前做差，分别获得差值影像i差值-现和i差值-前；

对差值影像采用固定阈值(阈值设置为0.05)快速去除地表，分别与原始影像ib3-现和ib3-前相乘获取去除地表后的影像iscf-现和iscf-前；其中iscf-现表示使用当前时刻数据获取的去除晴空地表后的影像iscf；iscf-前表示使用前一时刻数据获取的去除地表后影像iscf；

2、归一化雪指数法去除影像iscf中的雪和不透明冰云，获取去除后的影像icf；

分别使用2015年11月26日-12月1日当前时刻ib3-现、ib5-现和前一时刻ib3-前、ib5-前数据，获取当前时刻和前一时刻归一化雪指数ndsi现和ndsi前，ndsi现＝(ib3-现-ib5-现)/(ib3-现+ib5-现)；ndsi前＝(ib3-前-ib5-前)/(ib3-前+ib5-前)；

对去除地表后的影像iscf-现和iscf-前使用归一化雪指数法，设置阈值0.4去除影像iscf-现和iscf-前中的雪和不透明冰云，分别与影像iscf-现和iscf-前相乘得到影像icf-现和icf-前；其中icf-现表示使用当前时刻数据获取的去除雪和不透明冰云的影像icf；icf-前表示使用前一时刻数据获取的去除雪和不透明冰云的影像icf。

3、去除影像icf中纹理复杂且快速运动的云，获取雾检测影像iitwcf；

对当前时刻icf-现数据和前一时刻icf-前数据的b3波段做比值运算，获取比值影像i比值，i比值＝icf-现/icf-前；

对比值影像i比值采用一阶概率纹理滤波法获取datarange(数据范围)纹理数据idr；

设置阈值0.3去除影像icf-现中运动的及类型变化较大的云,与影像icf-现相乘得到影像iitwcf；

4、去除影像iitwcf中冰云和薄卷云，获取去除后影像iwcf；

使用2015年11月26日-12月1日当前时刻ib13-现(中心波段：10.4um)数据，设置阈值230k去除影像iitwcf中冰云，与影像iitwcf相乘得到影像itwcf；

对2015年11月26日-12月1日当前时刻ib11-现(中心波段：8.7um)数据和ib14-现(中心波段：11.2um)数据做差值运算，获取差值影像i差值1，i差值1＝ib11-现-ib14-现；

设置阈值0k去除影像itwcf中薄卷云，与影像itwcf相乘得到影像iwcf；

5、去除影像iwcf中水云，获取去除后影像ic2f；

对2015年11月26日-12月1日当前时刻ib15-现(中心波长：12.3um)数据和ib11-现数据做差值运算，获取差值影像i差值2，i差值2＝ib15-现-ib11-现；

其阈值依赖太阳高度角θ的变化，所以归一化处理获取动态阈值vt，

使用归一化阈值vt去除影像iwcf中水云，与影像iwcf相乘得到影像ic2f；

6、利用多天红外数据合成晴空底图ibt，阈值去除影像ic2f中的低层云，获取卫星雾检测结果if；

使用2015年11月16日-12月1日当前时刻b14数据合成晴空底图ibt；

对2015年11月26日-12月1日当前时刻ib14-现数据与晴空底图ibt做差，获取差值的绝对值影像i差值3，i差值3＝ib14-现-ibt；

采用阈值法(设置阈值11)快速去除影像ic2f中的低层云，与ic2f相乘获取最终卫星雾检测结果if；

7、检测结果精度定量验证

选取中国气象局2015年11月26日-2015年12月1日白天卫星影像对应时刻8:00和14:00的地面观测数据，验证白天雾检测算法的精度，图2示出了8:00卫星雾检测结果大约占整个雾区20％左右，其余5天8：00-9：00不同时刻白天雾检测结果不超过整个雾区40％；说明使用实验提出的白天雾检测算法此时间段检测效果不佳，进而说明在太阳高度角较低情况下获取的数据此算法不适用；14：00遥感雾检测结果与地面站点数据中的特浓雾、浓雾和雾区基本一致，使用精度指标进一步验证算法的有效性。

为进一步验证算法的有效性，对检测结果进行定量精度指标评价。常用指标评价体系评价分类的精度，验证算法的有效性^【12】，其中检验指标包括正确率(probabilityofdetection，pod)，误警率(falsealarmratio，far)和可靠性因子(criticalsuccessindex，csi)。这些指标定义为：

式中：nx为检测个数，其下标x为检测指标种类，包括h,m和f(h表示卫星检测结果与地面观测结果一致，即正确检测；m表示卫星检测结果中没雾而地面数据显示有雾，即失误检测；f表示卫星检测结果中有雾而地面数据显示没雾，即错误检测)。上述检测指标范围均为0-1，pod越大表示检测精度越高；csi越大说明方法越有效。

表1式出6天中14:00卫星雾检测正确率平均值为96.6％，错误率平均值为9.4％，可靠性因子平均值为87.9％。分析其检测错误原因如下：2015年11月27日14：00的遥感影像中存在少量小面积雾处于消散阶段，已逐步抬升成为低层云，此时地面能见度升高，地面观测站判定此区域无雾；而遥感影像上易被当做雾检测出来造成虚假报警。因此若对大面积雾进行检测，算法的正确性和可靠性因子均较高，但对于抬升为低层云的情况，算法易虚假报警。

表1白天雾检测精度

8、检测结果精度定性验证

雾具有时间连续性和空间位置变化较小的特性，因此对地面观测站点8:00和14:00确定为雾的区域，可认为此时间段内其他时刻该区域也是雾，图3示出了2015年11月30号8:00-16:00雾检测结果，定性验证算法的精度。

由图3可知：8:00卫星雾检测面积大约占整个雾区20％左右，8:00-9:00其他时刻雾检测面积不超过整个雾区60％，同时15:00-16:00雾检测面积随时间逐渐减小，16:00雾检测面积约占整个雾区15％左右；说明使用本文提出的算法对太阳高度角较低的时段(8:00-9:00,15:00-16:00)数据检测效果不佳；而白天时间段(9:00-15:00)雾检测结果占整个雾区90％以上，雾检测结果较好。

为进一步分析算法的适用时间段，选取2015年11月30日8:00-16:00的遥感数据进行了雾反射率，地表反射率及其差值的统计分析，见图4。日出前后(8：00-9：00)，太阳高度角较低，地表和雾顶接收到的太阳辐射较少，此时地表和雾在可见光波段反射率差异较小(低于0.05)，日出后日落前(9:00-15:00)太阳高度角升高，雾和地表反射率持续升高，此时地表和雾在可见光波段反射率差异较大(高于0.05)，而日落前(15:00-16:00)太阳高度角降低，雾与地表的反射率的差值小于0.05。本发明从理论上分析了该算法得适用时间段(9:00-15:00)，与图3中雾检测结果较好的时段一致。