基于环境一号卫星的陆地气溶胶光学厚度反演方法技术

本发明专利技术公开了一种基于环境一号卫星的陆地气溶胶光学厚度反演方法，该方法包括步骤：S1.模拟观测几何、大气状况，进行辐射传输计算，获得大气参数，形成查找表；S2.根据辅助数据，提取各像元的观测天顶角以及方位角，并计算太阳天顶角以及方位角，获得观测几何；S3.对HJ-1A/B星的CCD数据进行预处理；S4.利用暗目标算法获得各波段地表反射率信息，将查找表插值到各像元的观测几何，根据经过预处理得到的表观反射率计算不同陆地气溶胶光学厚度下的地表反射率，然后拟合得到陆地气溶胶光学厚度。本发明专利技术的方法能够快速、方便的实现大范围的陆地气溶胶监测，为颗粒物污染监测提供了数据源。

【技术实现步骤摘要】

本专利技术涉及大气遥感
，尤其涉及一种基于环境一号卫星的陆地气溶胶光学厚度反演方法。
技术介绍
卫星遥感气溶胶的研究始于上世纪七十年代中期，在上个世纪最后的30多年里，用于研究和形成业务化使用的卫星遥感气溶胶包括在火山喷发气溶胶监测、NOAA，美国国家海洋和大气管理局)系列卫星对深海上空来自沙尘暴和森林火灾造成的气溶胶层光学厚度遥感、红外遥感沙尘性气溶胶遥感、TOMS紫外波段对吸收性气溶胶的遥感等几个方面。美国地球观测系统计划(EOS)Terra和Aqua携带的中分辨率成像光谱仪(MODIS)对全球发布了10千米气溶胶光学厚度的全球分布产品；法国自1996年起推动了POLDER(偏振多角度)相机探测海洋和陆地气溶胶的研究。到目前为止，利用卫星遥感气溶胶已经形成了一定的全球气溶胶探测体系。目前的卫星遥感陆地气溶胶反演算法多为暗目标算法，该方法多利用短波红外波段识别暗目标并获得地表反射率进而反演得到气溶胶，但是短波红外波段数据的获得复杂，需对观测卫星的遥感传感器的要求很高，需进行特殊的波段设置，另一方便也导致其成本的增加。
技术实现思路
(一)要解决的技术问题本专利技术要解决的技术问题是：提供一种基于环境一号卫星的可快速、方便的实现大范围的陆地气溶胶监测的陆地气溶胶光学厚度反演方法。(二)技术方案为解决上述问题，本专利技术提供了一种基于环境一号卫星的陆地气溶胶光学厚度反演方法，该方法包括步骤：S1.根据环境一号HJ-1A/B星的CCD相机的特征，模拟各像元的观测几何、大气状况，进行辐射传输计算，获得相应的大气参数，存储于数据文件中形成查找表；S2.根据环境一号HJ-1A/B星的CCD数据的辅助数据，提取各像元的观测天顶角以及观测方位角，并计算太阳天顶角以及太阳方位角，获得各像元的观测几何；S3.对环境一号HJ-1A/B星的CCD数据进行如下预处理：辐射定标获得相机各波段的表观反射率、去除云覆盖像元、以及海陆分离提取陆地像元；S4.利用暗目标算法获得各波段地表反射率信息，将所述查找表插值到各像元的观测几何，，根据经过预处理得到的表观反射率计算不同陆地气溶胶光学厚度下的地表反射率，拟合得到陆地气溶胶光学厚度。其中，步骤S1进一步包括：S1.1设定不同的观测情况：9个太阳天顶角，12个观测天顶角，16个相对方位角，-->气溶胶模式为大陆型气溶胶，相对于0.55μm波长处的气溶胶光学厚度设为6个等级，查找表计算的波段为环境一号HJ-1A/B星的CCD的蓝波段以及红波段，海拔设置为0m，所述相对方位角为太阳方位角与观测方位角之差；S1.2调用辐射传输模式进行模拟运算，得到相应的大气程辐射、大气整层双程透过率、大气向下的半球反射率的三个参数，并存储在数据文件中形成查找表。其中，步骤S2进一步包括：S2.1根据环境一号HJ-1A/B星的CCD数据中的几何定位信息，计算每个像元的经纬度，得到太阳天顶角和太阳方位角，太阳天顶角θs为：cosθs＝sin(lat)sinδ+cos(lat)cosδcost其中，lat为地理纬度，δ为太阳光与地球赤道平面的夹角，t为太阳的时角，定义纬地方时12点的时间为0，6点时为-π/2，18点时为π/2；太阳方位角φ为：cos(φ)=tan(π/2-θ)tan(lat)-sinδcos(π/2-θ)*cos(lat);]]>S2.2根据环境一号HJ-1A/B星的CCD数据的辅助数据中的观测方位角以及所述太阳方位角，计算所述相对方位角。其中，在步骤S3中，所述辐射定标获得表观反射率的方法包括步骤：S3.11从所述辅助数据中读取辐射定标系数g和L0，将遥感影像像元亮度DN值转换为表观辐亮度L：L＝DN/g+L0S3.12根据所述太阳天顶角θs，将所述表观辐亮度归一化为各波段的表观反射率：ρ=πLEλcosθs]]>其中，Eλ为大气层顶太阳辐照度；其中，在步骤S3中，所述去除云覆盖像元时，选取0.2作为云去除阈值。其中，在步骤S3中，所述海陆分离提取陆地像元是通过借助归一化植被指数NDVI阈值实现的：NVID=ρnir-ρredρnir+ρred]]>其中，ρnir和ρred分别表示环境一号HJ-1A/B星的CCD近红外波段和红波段的地表反射率，判断标准为：水体的NDVI小于0，土壤的NDVI为0.1，植被的NDVI大于0.3。其中，步骤S4进一步包括：S4.1使用NDVI识别出暗目标浓密植被；S4.2将所述查找表插值到各像元的观测几何，得到不同陆地气溶胶光学厚度下插值后的大气参数；S4.3根据步骤S3.12中得到的蓝波段和红波段的表观反射率以及下式：-->μS＝cosθs，μV＝cosθv，与μv分别为太阳天顶角以及观测天顶角，r为朗伯体地表反射率，S为大气向下的半球反射率，T为大气整层双程透过率，ρ0为大气的路径辐射项等效反射率，为相对方位角；拟合得到的蓝波段和红波段地表反射率；S4.4使步骤S4.3得到的蓝波段和红波段地表反射率符合蓝波段和红波段地表反射率的线性关系，从而得到的气溶胶光学厚度；所述暗目标在蓝波段和红波段地表反射率的线性关系表示为：rredS=krblueS]]>其中，分别表示暗目标在红波段和蓝波段的地表反射率，k为红波段和蓝波段地表反射率比率。其中，所述环境一号HJ-1A/B星的CCD数据为经过几何校正的蓝、绿、红、近红外四个波段，所述环境一号HJ-1A/B星的CCD数据的辅助数据包括：包含标定系数、过境时间辅助信息的xml文件，以及包含观测几何的文本文件。(三)有益效果本专利技术的方法从环境一号HJ-1A/B星的CCD数据中反演得到陆地气溶胶光学厚度，能够快速、方便的实现大范围的陆地气溶胶监测，为颗粒物污染监测提供了数据源。附图说明图1为依照本专利技术一种实施方式的基于环境一号卫星的陆地气溶胶光学厚度反演方法流程图。具体实施方式本专利技术提出的基于环境一号卫星的陆地气溶胶光学厚度反演方法，结合附图及实施例详细说明如下。HJ-1A/B是我国的环境与灾害监测预报小卫星星座系统，具有较高空间、时间分辨率和宽覆盖特点，其数据不仅能为环境与减灾业务运行系统提供重要保障，还将成为很多部门日常业务的重要数据源。基于环境卫星数据建立的环境与减灾应用系统，对推动遥感卫星业务服务具有重要的示范作用，为我国环境空气质量遥感监测提供新的遥感平台。本专利技术的方法主要基于中国环境一号A/B星(本申请中统称HJ-1A/B星)CCD数据，利用HJ-1A/B星CCD数据从HJ-1A/B星CCD数据的辅助数据计算得到的太阳天顶角、相对方位角等信息，经过预处理得到红光和蓝光波段表观反射率，最后反演得到的气溶胶光学厚度分布图。从而，本专利技术能够快速、方便的实现大范围的陆地气溶胶光学厚度的监测，为颗粒物污染监测提供了数据源。如图1所示，依照本专利技术一种实施方式的基于HJ-1A/B的陆地气溶胶光学厚度反演方法包括步骤：S1.利用辐射传输模式(本实施方式中辐射传输模型为SecondSimulation of the Satellite Signal in the本文档来自技高网...

【技术保护点】
１．一种基于环境一号卫星的陆地气溶胶光学厚度反演方法，其特征在于，该方法包括步骤：Ｓ１．根据环境一号ＨＪ－１Ａ／Ｂ星的ＣＣＤ相机的特征，模拟各像元的观测几何、大气状况，进行辐射传输计算，获得相应的大气参数，存储于数据文件中形成查找表；Ｓ２．根据环境一号ＨＪ－１Ａ／Ｂ星的ＣＣＤ数据的辅助数据，提取各像元的观测天顶角以及观测方位角，并计算太阳天顶角以及太阳方位角，获得各像元的观测几何；Ｓ３．对环境一号ＨＪ－１Ａ／Ｂ星的ＣＣＤ数据进行如下预处理：辐射定标获得相机各波段的表观反射率、去除云覆盖像元、以及海陆分离提取陆地像元；Ｓ４．利用暗目标算法获得各波段地表反射率信息，将所述查找表插值到各像元的观测几何，根据经过预处理得到的表观反射率计算不同陆地气溶胶光学厚度下的地表反射率，拟合得到陆地气溶胶光学厚度。

【技术特征摘要】
1.一种基于环境一号卫星的陆地气溶胶光学厚度反演方法，其特征在于，该方法包括步骤：S1.根据环境一号HJ-1A/B星的CCD相机的特征，模拟各像元的观测几何、大气状况，进行辐射传输计算，获得相应的大气参数，存储于数据文件中形成查找表；S2.根据环境一号HJ-1A/B星的CCD数据的辅助数据，提取各像元的观测天顶角以及观测方位角，并计算太阳天顶角以及太阳方位角，获得各像元的观测几何；S3.对环境一号HJ-1A/B星的CCD数据进行如下预处理：辐射定标获得相机各波段的表观反射率、去除云覆盖像元、以及海陆分离提取陆地像元；S4.利用暗目标算法获得各波段地表反射率信息，将所述查找表插值到各像元的观测几何，根据经过预处理得到的表观反射率计算不同陆地气溶胶光学厚度下的地表反射率，拟合得到陆地气溶胶光学厚度。2.如权利要求1所述的基于环境一号卫星的陆地气溶胶光学厚度反演方法，其特征在于，步骤S1进一步包括：S1.1设定不同的观测情况：9个太阳天顶角，12个观测天顶角，16个相对方位角，气溶胶模式为大陆型气溶胶，相对于0.55μm波长处的气溶胶光学厚度设为6个等级，查找表计算的波段为环境一号HJ-1A/B星的CCD的蓝波段以及红波段，海拔设置为0m，所述相对方位角为太阳方位角与观测方位角之差；S1.2调用辐射传输模式进行模拟运算，得到相应的大气程辐射、大气整层双程透过率、大气向下的半球反射率的三个参数，并存储在数据文件中形成查找表。3.如权利要求2所述的基于环境一号卫星的陆地气溶胶光学厚度反演方法，其特征在于，步骤S2进一步包括：S2.1根据环境一号HJ-1A/B星的CCD数据中的几何定位信息，计算每个像元的经纬度，得到太阳天顶角和太阳方位角，太阳天顶角θs为：cosθs＝sin(lat)sinδ+cos(lat)cosδcost其中，lat为地理纬度，δ为太阳光与地球赤道平面的夹角，t为太阳的时角，定义纬地方时12点的时间为0，6点时为-π/2，18点时为π/2；太阳方位角φ为：cos(φ)=tan(π/2-θ)tan(lat)-sinδcos(π/2-θ)*cos(lat);]]>S2.2根据环境一号HJ-1A/B星的CCD数据的辅助数据中的观测方位角以及所述太阳方位角，计算所述相对方位角。4.如权利要求3所述的基于环境一号卫星的陆地气溶胶光学厚度反演方法，...

【专利技术属性】
技术研发人员：王桥，王中挺，厉青，周春艳，张丽娟，王子峰，毛慧琴，杨幸，陈辉，黄陆雄，段文举，
申请(专利权)人：环境保护部卫星环境应用中心，
类型：发明
国别省市：11