一种基于Himawar-8 AOD的PM2.5遥感监测反演方法专利检索-土球作物管理专利检索查询-专利查询网

一种基于Himawar-8 AOD的PM2.5遥感监测反演方法

阅读：204发布：2023-12-31

专利汇可以提供一种基于Himawar-8 AOD的PM2.5遥感监测反演方法专利检索，专利查询，专利分析的服务。并且本发明公开了一种基于Himawar-8AOD的PM2.5遥感监测反演方法，包括以下步骤：S1、获取Himawari-8AOD产品数据，PM2.5地面监测站点数据，气象数据以及土地利用数据；并进行时空配准，建立遥感估计近地面PM2.5 质量浓度的混合效应回归模型；S2、分别对获取到的数据进行预处理，提取出研究区域的AOD地理信息，各气象因子地理信息，以及地理信息；通过Matlab将预处理得到的信息以.mat格式存储；S3、编写Matlab代码读取步骤S2中预处理得到的AOD地理信息、各气象因子地理信息以及地理信息，并应用于步骤S1中建立的混合效应回归模型，得到遥感估计近地面PM2.5质量浓度结果。本发明能够对PM2.5浓度进行大面积、连续性时间的动态监测，并且该方法的监测稳定性好，准确性高，具有较好的应用前景。，下面是一种基于Himawar-8 AOD的PM2.5遥感监测反演方法专利的具体信息内容。

权利要求

1.一种基于Himawar-8 AOD的PM2.5遥感监测反演方法，其特征在于，该方法包括以下步骤：
S1、获取研究区域内Himawar-8卫星采集到的Himawari-8 AOD产品数据，地面监测站点采集到的PM2.5地面监测站点数据，气象数据以及土地利用数据；对获取到的数据进行时空配准，建立遥感估计近地面PM2.5质量浓度的混合效应回归模型；
S2、分别对Himawari-8 AOD产品数据进行预处理，提取出研究区域的AOD地理信息；对气象数据进行预处理，提取出研究区域的各气象因子地理信息；对土地利用信息进行预处理，提取出研究区域的地理信息；通过Matlab将预处理得到的信息以.mat格式存储；
S3、编写Matlab代码读取步骤S2中预处理得到的AOD地理信息、各气象因子地理信息以及地理信息，并应用于步骤S1中建立的混合效应回归模型，得到遥感估计近地面PM2.5质量浓度结果。
2.根据权利要求1所述的基于Himawar-8 AOD的PM2.5遥感监测反演方法，其特征在于，步骤S2中建立遥感估计近地面PM2.5质量浓度的混合效应回归模型的方法具体为：
S21、提取某一时间段内Himawari-8 AOD产品数据、PM2.5地面监测站点数据、气象数据与土地利用数据；
S22、对Himawari-8 AOD产品数据、PM2.5地面监测站点数据、气象数据进行空间坐标上的匹配获得样本数据；
S23、加入土地利用信息对样本数据进行回归分析，建立遥感估计近地面PM2.5质量浓度的混合效应回归模型，该模型为：
-3
其中，PM2.5代表近地面PM2.5质量浓度的小时均值，单位为ug·m ；AOD是指Himawari-
8提供的气溶胶光学厚度值；BLH是指欧洲中期天气预报中心ECMWF资料获取的边界层高度，单位m；RH为地基测量的近地表相对湿度，单位为％；TEM为地基测量的近地表温度，单位为K；PRES为近地表压强，单位为Pa；为考虑了地表风速的大小及方向的分量，单位为m/s；DEM为研究区域的数字高程模型信息，单位为m；NDVI为月均值归一化植被指数，无单位；β0为公式中固定的截距；β1～β8为公式中相应的回归系数，与预测变量相关。
3.根据权利要求1所述的基于Himawar-8 AOD的PM2.5遥感监测反演方法，其特征在于，步骤S1中获取数据的方法具体为：
利用Himawari-8卫星对大气进行遥感监测提供的气溶胶二级产品，获得Himawari-8 AOD产品数据；
利用欧洲中期天气预报中心ECMWF中的再分析数据集，获取气象数据；
利用美国航空航天局NASA提供的MODIS Level 3月均值归一化植被指数和空间分辨率为90米的DEM产品，获得土地利用信息数据；
利用PM2.5地面监测站点获取近地表的PM2.5地面监测站点数据。
4.根据权利要求1所述的基于Himawar-8 AOD的PM2.5遥感监测反演方法，其特征在于，气象数据包括边界层高度、相对湿度、近地面温度、近地面压强、近地面风速；气象数据从欧洲中期天气预报中心再分析数据集资料库中下载获得。
5.根据权利要求2所述的基于Himawar-8 AOD的PM2.5遥感监测反演方法，其特征在于，步骤S22中空间坐标上的匹配方法为：
气象数据采用与地面监测站点最近的4个格点气象数据进行加权最小二乘匹配；
Himawari-8 AOD数据采用以PM2.5地面监测站点为中心周围5公里半径范围内的非零象元均值进行匹配。
6.根据权利要求1所述的基于Himawar-8 AOD的PM2.5遥感监测反演方法，其特征在于，步骤S2中对Himawari-8 AOD产品数据进行预处理的方法为：利用LeapFTP 软件下载NetCDF格式的Himawari-8 AOD产品数据，并在Matlab中编写代码提取其地理信息数据并转成.mat格式存储。
7.根据权利要求1所述的基于Himawar-8 AOD的PM2.5遥感监测反演方法，其特征在于，步骤S2中对气象数据进行预处理的方法为：通过欧洲中期天气预报中心ECMWF官网下载气象数据，在Matlab中编写代码提取其地理信息数据并转化成.mat格式存储。
8.根据权利要求1所述的基于Himawar-8 AOD的PM2.5遥感监测反演方法，其特征在于，步骤S2中对土地利用信息进行预处理的方法为：通过美国航空航天局NASA下载MODIS Level 3月均值归一化植被指数NDVI和空间分辨率为90米的DEM产品，获得土地利用信息数据，在Matlab中编写代码提取其地理信息数据并转化成.mat格式存储。
9.根据权利要求1所述的基于Himawar-8 AOD的PM2.5遥感监测反演方法，其特征在于，该方法中Himawari-8 AOD产品数据是Himawari-8卫星采集到的Himawari-8 AOD Level 2产品数据，时间分辨率为10min，空间分辨率为10km；气象数据是欧洲中期天气预报中心全球再分析数据资料ERA-Interim中获取的，空间分辨率是0.125°；土地利用信息数据是美国航空航天局NASA下载得到的，月均值归一化植被指数NDVI时间分辨率是16days，空间分辨率是0.05°，DEM空间分辨率是90m。

说明书全文

一种基于Himawar-8 AOD的PM2.5遥感监测反演方法

技术领域

[0001] 本发明涉及大气遥感技术领域，尤其涉及一种基于Himawar-8 AOD的PM2.5遥感监测反演方法。

背景技术

[0002] 暴露于空气中的细颗粒物(PM2.5)，指空气动力学直径小于等于2.5μm的颗粒，是大气污染的重要组成成分之一,也是造成雾霾的主要原因。最新的研究表明，自1990年以来，全球环境的空气污染负担一直在持续增加。PM2.5的影响主要体现在对人体呼吸道健康的危害、地表温度和大气降水的变化等多个方面。流行病学研究表明，PM2.5与心血管呼吸系统疾病发病率和死亡率的增加有关，当亚微米大小的颗粒物渗透到肺部，可能诱发呼吸道疾病。因此，长期、大面积的PM2.5监测和准确的PM2.5浓度预测对于评估空气质量以及保障公众的健康至关重要。

[0003] 传统的PM2.5浓度遥感监测方法大部分是基于MODIS AOD数据来进行估算的，估算结果受AOD数据精度影响较大，并且这种极轨卫星提供的数据受观测频率的限制，因此这仍然不足以让我们理解PM2.5质量浓度在连续性时间上的演变。Himawari-8(向日葵8号)气象卫星是日本宇宙航空研究开发机构设计制造的向日葵系列卫星之一，是日本于2014年10月7日发射的新一代气象卫星。Himawari-8是世界上第一颗可以拍摄彩色图像的静止气象卫星，以往的极轨卫星受观测频率的限制，而Himawari-8的观测频率提高到了至少每10分钟一次，覆盖地球的三分之一(西太平洋、东亚、东南亚和大洋洲)，对云层等动向的持续观测性能也得到提升。

[0004] 故而，在这种背景下，研究一种新的能够准确对PM2.5浓度进行大面积、连续性时间的动态监测方法，为科学研究中待解决的问题。

发明内容

[0005] 本发明要解决的技术问题在于针对现有技术中遥感反演时间分辨率低的缺陷，提供一种基于Himawar-8 AOD的PM2.5遥感监测反演方法。

[0006] 本发明解决其技术问题所采用的技术方案是：

[0007] 本发明提供一种基于Himawar-8 AOD的PM2.5遥感监测反演方法，该方法包括以下步骤：

[0008] S1、获取研究区域内Himawar-8卫星采集到的Himawari-8 AOD产品数据，地面监测站点采集到的PM2.5地面监测站点数据，气象数据以及土地利用数据；对获取到的数据进行时空配准，建立遥感估计近地面PM2.5质量浓度的混合效应回归模型；

[0009] S2、分别对Himawari-8 AOD产品数据进行预处理，提取出研究区域的AOD地理信息；对气象数据进行预处理，提取出研究区域的各气象因子地理信息；对土地利用信息进行预处理，提取出研究区域的地理信息；通过Matlab将预处理得到的信息以.mat格式存储；

[0010] S3、编写Matlab代码读取步骤S2中预处理得到的AOD地理信息、各气象因子地理信息以及地理信息，并应用于步骤S1中建立的混合效应回归模型，得到遥感估计近地面PM2.5质量浓度结果。

[0011] 进一步地，本发明的步骤S2中建立遥感估计近地面PM2.5质量浓度的混合效应回归模型的方法具体为：

[0012] S21、提取某一时间段内Himawari-8 AOD产品数据、PM2.5地面监测站点数据、气象数据与土地利用数据；

[0013] S22、对Himawari-8 AOD产品数据、PM2.5地面监测站点数据、气象数据进行空间坐标上的匹配获得样本数据；

[0014] S23、加入土地利用信息对样本数据进行回归分析，建立遥感估计近地面PM2.5质量浓度的混合效应回归模型，该模型为：

[0015]

[0016] 其中，PM2.5代表近地面PM2.5质量浓度的小时均值，单位为ug·m-3；AOD是指Himawari-8提供的气溶胶光学厚度值；BLH是指欧洲中期天气预报中心ECMWF资料获取的边界层高度，单位m；RH为地基测量的近地表相对湿度，单位为％；TEM为地基测量的近地表温度，单位为K；PRES为近地表压强，单位为Pa；为考虑了地表风速的大小及方向的分量，单位为m/s；DEM为研究区域的数字高程模型信息，单位为m；NDVI为月均值归一化植被指数，无单位；β0为公式中固定的截距；β1～β8为公式中相应的回归系数，与预测变量相关。

[0017] 进一步地，本发明的步骤S1中获取数据的方法具体为：

[0018] 利用Himawari-8卫星对大气进行遥感监测提供的气溶胶二级产品，获得Himawari-8 AOD产品数据；

[0019] 利用欧洲中期天气预报中心ECMWF中的再分析数据集，获取气象数据；

[0020] 利用美国航空航天局NASA提供的MODIS Level 3月均值归一化植被指数和空间分辨率为90米的DEM产品，获得土地利用信息数据；

[0021] 利用PM2.5地面监测站点获取近地表的PM2.5地面监测站点数据。

[0022] 进一步地，本发明的气象数据包括边界层高度、相对湿度、近地面温度、近地面压强、近地面风速；气象数据从欧洲中期天气预报中心再分析数据集资料库中下载获得。

[0023] 进一步地，本发明的步骤S22中空间坐标上的匹配方法为：

[0024] 气象数据采用与地面监测站点最近的4个格点气象数据进行加权最小二乘匹配；Himawari-8 AOD数据采用以PM2.5地面监测站点为中心周围5公里半径范围内的非零象元均值进行匹配。

[0025] 进一步地，本发明的步骤S2中对Himawari-8 AOD产品数据进行预处理的方法为：利用LeapFTP 软件下载NetCDF格式的Himawari-8 AOD产品数据，并在Matlab中编写代码提取其地理信息数据并转成.mat格式存储。

[0026] 进一步地，本发明的步骤S2中对气象数据进行预处理的方法为：通过欧洲中期天气预报中心ECMWF官网下载气象数据，在Matlab中编写代码提取其地理信息数据并转化成.mat格式存储。

[0027] 进一步地，本发明的步骤S2中对土地利用信息进行预处理的方法为：通过美国航空航天局NASA下载MODIS Level 3月均值归一化植被指数NDVI和空间分辨率为90米的DEM产品，获得土地利用信息数据，在Matlab中编写代码提取其地理信息数据并转化成.mat格式存储。

[0028] 进一步地，本发明的该方法中Himawari-8 AOD产品数据是Himawari-8卫星采集到的Himawari-8 AOD Level 2产品数据，时间分辨率为10min，空间分辨率为10km；气象数据是欧洲中期天气预报中心全球再分析数据资料ERA-Interim中获取的，空间分辨率是0.125°；土地利用信息数据是美国航空航天局NASA下载得到的，月均值归一化植被指数NDVI时间分辨率是16days，空间分辨率是0.05°，DEM空间分辨率是90m。

[0029] 本发明产生的有益效果是：本发明的基于Himawar-8 AOD的PM2.5遥感监测反演方法，能够对PM2.5浓度进行大面积、连续性时间的动态监测，并且该方法的监测稳定性好，准确性高，具有较好的应用前景。附图说明

[0030] 下面将结合附图及实施例对本发明作进一步说明，附图中：

[0031] 图1是本发明实施例的方法流程图。

具体实施方式

[0032] 为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。

[0033] 如图1所示，本发明实施例基于Himawari-8 AOD产品数据、ECMWF气象数据及土地利用信息数据进行回归分析，建立了遥感估计近地面PM2.5质量浓度的混合效应回归模型，最终得到通过遥感技术反演的近地面PM2.5质量浓度。

[0034] 本发明实施例的遥感反演近地面细颗粒物质量浓度的方法，首先是需要下载获取Himawari-8卫星提供的气溶胶二级产品以及欧洲中期天气预报中心提供的气象数据；然后构建PM2.5地面监测站点的空间数据集，为步骤S1的数据时空配准提供基础。

[0035] 图1示出本发明基于Himawari-8 AOD产品数据、ECMWF气象数据及土地利用信息数据的遥感反演近地面细颗粒物质量浓度的方法流程。参考图1，本发明的遥感反演估算方法进一步包括以下步骤：

[0036] 步骤S1.通过对Himawari-8 AOD数据、地面站点PM2.5监测数据、气象数据和土地利用数据进行时空配准，建立遥感估计近地面PM2.5质量浓度的混合效应回归模型；

[0037] 步骤S2.对Himawari-8 AOD产品数据进行预处理：提取出研究区域的AOD地理信息，并在Matlab中转化成.mat格式存储；对ECMWF气象数据进行预处理：提取出研究区域的各气象因子地理信息进行加权最小二乘匹配到AOD产品格网中，并以.mat格式存储；对土地利用信息进行预处理：提取出研究区域的地理信息，并以.mat格式存储；

[0038] 步骤S3.编写Matlab代码读取S2中的数据并应用于S1中建立的混合效应回归模型，得到遥感估计近地面PM2.5质量浓度结果。

[0039] 其中，步骤S1进一步包括以下步骤：

[0040] A)所述ECMWF气象数据包括边界层高度、相对湿度、近地面温度、近地面压强、近地面风速(包括大小和二维平面上的方向)；所述ECMWF气象数据从欧洲中期天气预报中心再分析数据集资料库中下载获得。

[0041] B)所述步骤S1中的空间配准上，气象数据采用与地面监测站点最近的4个格点气象数据进行加权最小二乘匹配；Himawari-8 AOD数据采用以PM2.5地面监测站点为中心周围5公里半径范围内的非零象元均值进行匹配；时间上以卫星观测时间为准。

[0042] C)对上述处理后的数据通过统计分析软件进行回归分析，建立如下述公式(1)所示的遥感反演估算PM2.5质量浓度的多元回归模型。

[0043]

[0044] 公式(1)中，PM2.5代表近地面PM2.5质量浓度的小时均值，单位为ug·m-3；AOD是指Himawari-8提供的气溶胶光学厚度值，无单位；BLH是指ECMWF资料获取的边界层高度，单位m；RH为地基测量的近地表相对湿度，单位为％；TEM为地基测量的近地表温度，单位为K；PRES为近地表压强，单位为Pa；为考虑了地表风速的大小及方向的分量，单位为m/s；DEM为研究区域的数字高程模型信息，单位为m；NDVI为月均值归一化植被指数，无单位；β0为公式中固定的截距；β1～β8为公式中相应的回归系数，与预测变量相关。

[0045] 其中，步骤S2进一步包括利用LeapFTP软件下载NetCDF格式的AOD产品数据，并在Matlab中编写代码提取其地理信息数据并转成.mat格式存储；通过ECMWF官网下载气象数据，在Matlab中编写代码提取其地理信息数据并转化成.mat格式存储；通过美国航空航天局(NASA)下载MODIS Level 3月均值归一化植被指数(NDVI)和空间分辨率为90米的DEM产品，获得土地利用信息数据，并在Matlab中编写代码提取其地理信息数据并转化成.mat格式存储。

[0046] 其中，步骤S3进一步包括编写Matlab代码读取步骤S1、步骤S2中预处理后的数据并应用于步骤S1中建立的混合效应回归模型，计算出遥感估计近地面PM2.5质量浓度结果。

[0047] 上述卫星为Himawari-8，所述AOD产品是指Himawari-8 AOD Level 2产品数据，时间分辨率为10min，空间分辨率为10km；所述气象数据是欧洲中期天气预报中心全球再分析数据资料ERA-Interim中获取的，空间分辨率是0.125°；所述土地利用信息数据是美国航空航天局(NASA)下载得到的，月均值归一化植被指数(NDVI)时间分辨率是16days，空间分辨率是0.05°，DEM空间分辨率是90m。

[0048] 应当理解的是，对本领域普通技术人员来说，可以根据上述说明加以改进或变换，而所有这些改进和变换都应属于本发明所附权利要求的保护范围。

标题	发布/更新时间	阅读量
一种造山型岩浆熔离钴矿的找矿方法	2020-05-08	326
一种高铝耐火砖及其制备方法	2020-05-08	191
一种建筑废渣处金属制品分离装置	2020-05-08	103
一种净化空气花盆的通风结构	2020-05-08	171
一种用于地质勘查的岩土类试样钻取装置	2020-05-11	285
一种耐高温耐腐蚀负温度系数热敏电阻的制造方法	2020-05-08	612
精品机制砂质量检测方法	2020-05-11	141
一种静音井盖	2020-05-08	101
自閉症および関連障害の治療および予防のための腸管選択的隔離剤	2020-05-08	51
熱伝導性複合材料および放熱シート	2020-05-08	232

一种基于Himawar-8 AOD的PM2.5遥感监测反演方法

一种基于Himawar-8 AOD的PM2.5遥感监测反演方法

技术领域

背景技术

发明内容

具体实施方式

该功能需要专业版企业版VIP权限，您可以：