基于地面实测的GF-4AOD产品高时间分辨率检验的制作方法

本发明涉及卫星遥感反演产品的真实性检验技术领域，特别涉及一种基于地面实测的gf-4aod产品高时间分辨率检验。

背景技术：

气溶胶在包括气候、大气环境与人类健康在内的各个领域具有重要作用，为了更好地研究气溶胶，需要对以气溶胶光学厚度(aerosolopticaldepth，aod)为主的物理参数进行长期观测。目前对气溶胶光学厚度观测主要有卫星遥感反演及地面站点观测两种方式，前者空间覆盖范围广，但精度较低；后者精度较高但空间覆盖较小，故可使用真实性检验站的观测数据对遥感观测结果进行验证。

对于地面上某一气溶胶观测站点来说，其观测具有一定的时间间隔，以太阳分光光度计ce-318为例，在日间晴朗条件下，其观测时间频率可高达3min/次；而卫星成像同样具有一定的重访周期，如太阳同步轨道卫星terra/aqua上搭载的modis传感器对中纬度地区的重访周期为0.5日。由于aod的时空变异程度较大，为了保证验证过程的可靠性，需要使得地面观测时间与卫星过境时间尽可能接近，同时选择空间上尽可能接近站点的像元，即需要确定合适的时空匹配窗口。当前相关研究中，较多研究者针对太阳同步轨道卫星传感器(如aqua/terra搭载的modis)气溶胶产品开展精度验证工作，为了获得较多的匹配数据对，这些研究多将卫星过境前后各30min内(时间窗口)的地面观测aod数据均值与以地面站点为中心的50km*50km空间邻域(空间窗口)内的aod反演均值进行比较。这种方法增加了匹配成功的数据量，但是客观上引入了由于时、空间窗口较大而导致的aod异质性误差。

技术实现要素：

本发明针对现有技术中的不足，提供一种基于地面实测的gf-4aod产品高时间分辨率检验，面向高分四号气溶胶光学厚度产品，利用其高时间分辨率的特点，采用更为精细化的时空匹配策略，即更短的时间窗口、更小的空间窗口，基于地面观测，实现对卫星遥感aod产品的真实性验证。

为实现上述目的，本发明采用以下技术方案：

基于地面实测的gf-4aod产品高时间分辨率检验，包括如下步骤：s1、判断gf-4aod产品的观测频率，如果卫星在同一个地点一段时间内进行了连续多帧成像，获得多个连续aod值，转至步骤s2；如果卫星在同一个地点一段时间内只进行了单景成像，转至步骤s3；

s2、获取以站点为中心，边长5～10km、地表类型均一的正方形空间窗口内gf-4aod像元平均值，同时获取与成像时间最相近的地面前后两次观测值的平均，将其与卫星像元值生成匹配对，进行真实性检验；转至步骤s4；

s3、采用站点50km*50km空间窗口内至少10个质量等级最高像元的aod值求平均，然后与地面实测的在卫星成像前后各0.5小时内所有观测值的平均进行匹配，并进行真实性检验；进入步骤s5；

s4、将匹配的散点绘制在二维坐标系之中，根据地面观测值、各站点处卫星反演aod值，进行线性拟合；转至步骤s6；

s5、将卫星高质量等级像元均值与地面相应时段内的观测平均值匹配生成数据对，对这些数据对进行线性拟合；进入步骤s6；

s6、计算精度评价指标，考察卫星反演结果与地面观测之间的一致性，评估aod反演精度，并分析误差来源。

为优化上述技术方案，采取的具体措施还包括：

进一步地，对于短时间段内获得多个连续aod值的情况采用精细时空匹配方法：对于每一帧aod影像，将其与各个地面真实性验证站点进行叠置，提取以每一个站点为中心，大小为5km*5km至10km*10km、且地表类型较为均一的空间窗口内对应的所有aod像元，求这些像元的均值，完成空间匹配；在时间上，在每一站点处获取与卫星成像时刻相距最近的前后各一条地面观测记录，对其求均值后完成时间匹配。

进一步地，设定验证时间段为0.5小时，将每0.5小时内的各帧aod影像生成的卫星与地面匹配数据对合并验证。

进一步地，对每一个站点建立缓冲区，缓冲区为以站点为中心的50km*50km范围，即50km*50km空间窗口。

进一步地，步骤s3包括将研究时段内每一幅高分四号卫星aod产品与缓冲区的图层进行叠置，统计完全落入各个缓冲区的像元数。

进一步地，对各个缓冲区内的aod像元值分别求均值，计算方法如下式所示：

其中，n表示高质量等级的像元数量，aodi表示第i个高质量aod像元，aodkj表示经过空间匹配操作(即求取空间窗口内高质量aod均值)之后的平均aod。

进一步地，对于每一景aod产品，选择其成像前0.5小时至成像后0.5小时时间段作为研究时段，选择各个站点在研究时段内的aod观测。

进一步地，步骤s4中，以地面数据为横坐标，卫星观测数据为纵坐标，绘制散点图。

进一步地，步骤s6中，计算精度评价指标包括均方根误差rmse、落入期望误差范围ee、反演数据的比例p、线性相关系数r，计算公式如下式所示：

δaod＝±(0.05+0.2·aodground)

其中，均方根误差rmse定量描述卫星反演aod(即aodgf4)较地面aod真值(即aodground)之间的差异的大小，于同一地点、同一时间段来说，其值越小，表明验证精度越高；n表示匹配数据对的数量。δaod定义了理论上卫星反演aod的绝对误差范围的最大值与最小值，而期望误差范围内匹配数据比例p则反映了卫星反演aod的总体精度。

本发明的有益效果是：首先，本专利采用更为精细化的时空匹配方法，面向具有高时间分辨率的高分四号卫星aod产品数据进行精度验证，可以在保证比对数据量的同时有效规避传统时空匹配方法容易引入的较高匹配误差，提升真实性检验结果的可靠性；其次，通过本发明可以分析得到高分四号卫星aod的潜在误差来源，有利于优化气溶胶反演算法；最后，该验证结果还可以为高分四号气溶胶产品的下游产品(如pm2.5空间分布)提供依据。

附图说明

图1为本发明检验方法的流程示意图。

图2为本发明检验方法的高分四号aod影像与加密验证站点分布示意图。

图3为本发明检验方法采用1*1km空间窗口方案时，各加密站处各幅影像aod的变化图。

图4为本发明检验方法采用10*10km空间窗口方案时，各加密站处各幅影像aod的变化图。

图5为本发明检验方法采用50*50km空间窗口方案时，各加密站处各幅影像aod的变化图。

图6为本发明检验方法中国各站点处使用1*1km方案与10*10km方案所得aod的差异示意图。

图7为本发明检验方法中国各站点处使用1*1km方案与50*50km方案所得aod的差异示意图。

具体实施方式

现在结合附图对本发明作进一步详细的说明。

需要注意的是，发明中所引用的如“上”、“下”、“左”、“右”、“前”、“后”等的用语，亦仅为便于叙述的明了，而非用以限定本发明可实施的范围，其相对关系的改变或调整，在无实质变更技术内容下，当亦视为本发明可实施的范畴。

对时间变化较快的遥感产品，如大气产品，在进行产品真实性检验时需同时进行时间与空间的匹配，制定合理的匹配策略。高分四号卫星对特征区域可实现高频次观测，使得产品检验技术与极轨卫星存在较大的差异。极轨卫星重返周期通常较长，且经过同一特征区域的过境时间相对固定，因此在进行产品检验时通常选择卫星过境前后一段时间内采集的样本数据；同时，为了增加匹配数据对从而有效提高统计质量，这些研究往往采用较长的时间窗口与较大的空间范围窗口，这将引入匹配过程误差。而由于高分四号可以实现全天候多频次的观测，因此在实验区的选择、实验数据的获取及检验方法上都较极轨卫星发生改变，必须开展时空匹配技术研究，结合长时间序列的样本数据，从而实现高精度的产品检验技术。

高分四号(gf-4)卫星是我国自主研发的一颗地球同步轨道卫星，于2015年发射成功。与其它地球静止卫星相比，高分四号卫星可以达到较高的空间分辨率(可见光与近红外波段可达50米)与时间分辨率(20秒)，可以为aod的高时空分辨率反演提供数据支撑。目前，已有的高分四号aod四级产品的空间分辨率达到了1km。

如附图1-7所示，本发明提供了一种基于地面实测的gf-4aod产品高时间分辨率检验，针对高分四号气溶胶光学厚度产品真实性检验技术开展了针对性的研究，我们采用几个观测站连续半年观测的数据进行两种验证方法分析：

(1)选取以站点为中心的5km*5km～10km*10km的地表类型均一的空间窗口内的所有卫星aod像元并求平均，接着获取与卫星像元成像时间最相近前后两个地面观测，对其求平均后与卫星像元均值匹配，进行真实性检验；

(2)采用卫星50km*50km空间窗口(以站点为中心的50km*50km范围)内至少10个质量等级最高像元的aod值求平均，然后与地面实测的在卫星成像前后各0.5小时内所有观测值的平均进行匹配，并进行真实性检验。

研究分析表明，如果高分四号卫星在同一个地点较短时间内进行了连续多帧成像，即获得多个连续aod值，此时采用第一种检验方法，检验精度更高。反之，如果在同一个地点一段时间内只进行了单景成像，此时采用第二种检验方法，检验精度更高。

最后，将匹配散点绘制在二维坐标系之中：对于方案(1)，根据上述综合地面观测值、各站点处卫星反演aod值，进行线性拟合，计算精度评价指标，考察卫星反演结果与地面观测之间的一致性，评估aod反演精度，并分析误差来源。对于方案(2)，对卫星高质量等级像元均值与地面相应时段内的观测平均值进行线性拟合，分析误差来源。

本发明的一个实施例：由于地面站点仅在晴朗天气产生aod观测记录，故选择某日北京时间9：00-11：00期间的所有高分四号卫星aod产品。同时，从aeronet官网(https：//aeronet.gsfc.nasa.gov/)下载同时段的地面观测aeronetaod数据，考虑到地面站点较高的观测精度，将其视为aod真值。同时考虑到高分四号的成像空间范围，由北向南、由东向西依次可以选取如下站点的数据：beijing、beijing-cams、xianghe、gangneung_wnu、chiba_university、xuzhou-cumt、chen-kung_univ、dibrugarh_univ、hong_kong_polyu、manila_observatory、singapore、pontianak。

1.首先，判断gf-4aod产品的观测频率。对于即时监测等高时间频率需求研究，转至步骤2；对于气候研究等其它长时间序列aod动态分析需求研究，转至步骤3。

2.对于高时间频率(20s/帧)aod动态研究采用精细空间窗口时空匹配方法：提取以每一个站点为中心，大小为5*5km至10*10km且地表类型较为均一的空间窗口内对应的所有aod像元，求这些像元的均值，完成空间匹配；在时间上，在每一站点处获取与卫星像元成像时间最近的前后两次观测记录值的平均，将其与卫星像元值生成匹配对，完成时间匹配。为了保证样本数量的充足，设定验证时间段为0.5小时，将每0.5小时内的各帧aod影像生成的数据对合并验证。

3.对于低时间频率aod动态研究：根据各个站点的地理坐标(可从aeronet数据下载页面获取)，分别建立每一个站点的缓冲区，缓冲区为以站点为中心的50km*50km范围。将研究时段内每一幅高分四号卫星aod产品与缓冲区的图层进行叠置，统计完全落入各个缓冲区的像元数。由于存在云、雪、内陆水体的影响，aod反演产品存在像元缺失现象，故每一个缓冲区内部的像元数不一；仅当缓冲区内像元数大于10个时，对各个缓冲区内的aod像元值分别求均值，计算方法如下式所示，完成空间匹配，得到aodkj(下标kj表示空间匹配)；

接下来，针对研究时段内的每一景aod产品，选择其成像前30分钟至成像后30分钟时间段作为研究时段，选择各个站点在研究时段内的aod观测。由于在高分四号卫星aod反演的过程中同时会生产反映反演质量的标记图层，每一个aod反演像元均对应1个质量标记值(qac，其值为1，2或3)，值越大表明质量等级越高。为了提高数据验证结果的可靠性，挑选各站点在上述研究时段之内的具有最高质量等级(qac等于3)的aod反演像元，统计这些像元的数量，若其数量不少于3，则分别对各站点该时段内的aod数据进行求均值操作，完成时间匹配，得到aodsj。计算方式如下式所示：

4.最后，采用计算机编程方式，基于日期与时刻对时间匹配操作与空间匹配操作后的aod数据结果进行关联，而后以地面数据为横坐标，卫星观测数据为纵坐标，绘制散点图，计算评价指标，完成精度验证。使用到的主要评价指标为均方根误差(rmse)、落入期望误差范围(expectederror，ee)反演数据的比例(p)、线性相关系数(r)，计算公式如下式所示。

δaod＝±(0.05+0.2·aodground)(4)

其中，rmse定量描述卫星反演aod较地面aod真值之间的差异的大小，对于同一地点、同一时间段来说，其值越小，表明验证精度越高；而期望误差范围内匹配数据比例p则反映了卫星反演aod的总体精度，假设误差统计分布呈正态分布，当p大于68％时，即可认为aod卫星反演精度高。

以下以高频率观测方案为例，比较在此情景下采用大小分别为1km*1km、10km*10km与50km*50km(下文分别简称1km，10km，50km)的空间窗口的验证效果之差异。

本实例中的高分四号aod影像空间分布如图2所示，该数据由2016.07.23至2016.08.01期间的共8景影像构成，由于各景影像空间覆盖并不完全一致，故提取各幅影像的交集区域，在此区域中人工设置7个加密站点。接下来，基于这些站点，分别使用三种不同大小的空间窗口对各景影像aod值进行提取，将提取结果绘制成折线图，如图3-图5所示，分别采用1*1km、10*10km与50*50km空间窗口方案时，各加密站处各幅影像aod的变化图。图6、图7分别展示了各站点处使用1*1km方案与10*10km、50*50km方案所得aod的差异。

将各站点使用10km与50km空间窗口方案得到的aod值与使用1km方案得到的aod值作差，分别得到图6、图7。由图可知，随着影像(即观测日期)的变化，各空间窗口方案下各站点处的aod值变化趋势接近；但是，根据图6、图7，随着空间窗口的增大，各站点处所得aod值较之于1km方案差别愈大，例如在50km方案下，在站点1处，其与1km方案aod差值绝对值最高可达0.2(于160801影像达到)，而站点6处两种空间窗口方案结果也呈现出显著差异。根据aod真实性验证相关文献，这种差异是由验证的尺度效应造成的，当空间窗口较大时，对窗口内aod求得的均值将会受到邻域内非中心aod像元的影响而出现偏高或偏低的现象，选择1km方案得到的卫星aod结果将更好地代表站点观测值，而排除相邻像元的影响。因此，在高频次验证模式下，相较于传统的50km方案，使用更为精细化的1km方案将会提高验证精度。

本发明针对不同应用需求之下的不同aod反演时间频率，设计了适应于两种aod产品生产周期的验证方案，有效地保证了样本数，较好提升了真实性验证的可靠性。

以上仅是本发明的优选实施方式，本发明的保护范围并不仅局限于上述实施例，凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理前提下的若干改进和润饰，应视为本发明的保护范围。