平码五不中公式规律
  • / 9
  • 下载费用:30 金币  

一种立体测绘卫星三线阵CCD相机的大气订正方法.pdf

关 键 ?#21097;?/dt>
一种 立体 测绘 卫星 三线 CCD 相机 大气 订正 方法
  专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
摘要
申请专利号:

CN201410228861.5

申请日:

2014.05.27

公开号:

CN104133203A

公开日:

2014.11.05

当前法律状态:

授权

有效性:

有权

法?#19978;?#24773;: 授权|||实质审查的生效IPC(主分类):G01S 7/497申请日:20140527|||公开
IPC分类号: G01S7/497 主分类号: G01S7/497
申请人: 中国资源卫星应用中心
发明人: 韩启金; 王爱春; 傅俏燕; 闵祥军; 潘志强
地址: 100094 北京市海淀区永丰产业基地丰贤东路5号
优?#28909;ǎ?/td>
专利代理机构: 中国航天科技专利中心 11009 代理人: 庞静
PDF完整版下载: PDF下载
法律状态
申请(专利)号:

CN201410228861.5

授权公告号:

||||||

法律状态公告日:

2016.08.31|||2014.12.10|||2014.11.05

法律状态类型:

授权|||实质审查的生效|||公开

摘要

一种立体测绘卫星三线阵CCD相机的大气订正方法,(1)计算每个像元的入瞳处辐射亮度;并确立每个像元的地表反射率计算模型;(2)建立不同气溶胶类型和气溶胶光学厚度下的立体测绘卫星三线阵CCD相机每个波段的反演参数查?#20918;恚?3)依据Rahman?BRDF模型,建立立体测绘卫星前视、后视CCD相机与正视CCD相机地表反射率关系;并利用地表反射率计算模型计算三个相机不同气溶胶类型和光学厚度下的地表反射率;(4)建立最小化判定函数,确定函数的最优解,根据最优解确定三线阵CCD相机的气溶胶类型和光学厚度;(5)对步骤(4)的结果进行修正,获得三线阵CCD相机修正后的气溶胶光学厚度;(6)利用地表反射率计算模型计算得到最终的地表反射率。

权利要求书

权利要求书
1.  一种立体测绘卫星三线阵CCD相机的大气订正方法,所述的三线阵CCD相机包括前视、正视、后视三台CCD相机;其特征在于步骤如下:
(1)根据立体测绘卫星三线阵CCD相机的遥?#22411;?#20687;数据,计算每个像元的入瞳处辐射亮度;并依据太阳-大气-地表-大气-相机相互作用的大气辐射传输方程确立每个像元的地表反射率计算模型;
(2)依据步骤(1)中建立的地表反射率计算模型,设置气溶胶类型和气溶胶的光学厚度,建立不同气溶胶类型和气溶胶光学厚度下的立体测绘卫星三线阵CCD相机每个波段的反演参数查?#20918;恚?BR>(3)依据Rahman BRDF模型,建立立体测绘卫星前视CCD相机与正视CCD相机地表反射率关系,以及后视CCD相机与正视CCD相机地表反射率关系;并利用上述反演参数查?#20918;懟?#27599;个像元入瞳处辐射亮度以及地表反射率计算模型计算三个相机不同气溶胶类型和光学厚度下的地表反射率;
(4)根据步骤(3)的结果,建立最小化判定函数,确定函数的最优解,根据最优解确定三线阵CCD相机的气溶胶类型和光学厚度;
(5)根据气溶胶光学厚度与波长的关系对步骤(4)的结果进行修正,获得三线阵CCD相机修正后的气溶胶光学厚度;
(6)根据步骤(4)中确定的气溶胶类型以及步骤(5)中获得的气溶胶光学厚度,代入上述反演参数查?#20918;恚?#21033;用地表反射率计算模型计算得到最终的地表反射率。

2.  根据权利要求1所述的一种立体测绘卫星三线阵CCD相机的大气订正方法,其特征在于:所述步骤(1)中每个像元的地表反射率计算模型如下:
ρx,λObj(i,j)=Lx,λTOA(i,j)-Lx,λPath(i,j)Tx,λ(i,j)+Sx,λ(i,j)[Lx,λTOA(i,j)-Lx,λPath(i,j)]]]>
其中,为x相机第λ波段像元(i,j)处的地表反射率;x=f,n,b代表立体测绘卫星的前视、正视、后视相机;i,j为每台相机的像元位置;
为x相机第λ波段像元(i,j)处的入瞳处辐射亮度;
为x相机第λ波段像元(i,j)处的路径散射;
T(θx,s)为太阳到地面方向的大气漫射散射因子;T(θx,v)为地面到相机方向的大气漫射散射因子;θx,s为x相机观测时的太阳天顶角;ES,x,?#23435;獂相机第λ波段的大气外太阳光谱辐照度ES;
Sx,λ(i,j)为x相机第λ波段像元(i,j)处的大气球面反照率。

3.  根据权利要求1所述的一种立体测绘卫星三线阵CCD相机的大气订正方法,其特征在于:所述步骤(4)中的最小化判定函数如下:
min|Yλ(i,j)|=Σx=f,b[ρx,λR-Obj(i,j)-ρx,λObj(i,j)]2]]>
其中,为利用步骤(3)中的地表反射率关系以及步骤(3)中计算的正视CCD相机不同气溶胶类型和光学厚度下的地表反射率,得到的x相机的地表反射率,x=f,b分别代表立体测绘卫星的前视、后视相机;
为步骤(3)中计算的x相机不同气溶胶类型和光学厚度下的地表反射率。

说明书

说明书一种立体测绘卫星三线阵CCD相机的大气订正方法
技术领域
本发明涉及立体测绘卫星三线阵CCD相机的大气订正方法,特别涉及ZY-3立体测绘卫星三线阵CCD相机的大气订正方法。
背景技术
立体测绘卫星三线阵CCD采用3台线阵CCD(前视相机、正视相机、后视相机)在不同时刻对同一地物进行多角度?#19978;瘢?#19981;仅可以获得目标立体信息,同时还可以得到多角度的定量化信息,但是要获得高精度的定量化信息,就必须详细考虑大气对于?#19978;?#27874;?#26410;?#36755;过程的定量影响,并对其影响做出准确的大气订正,实现地物目标辐射特性真实再现,?#28304;?#21040;提高遥?#34892;?#24687;质量及定量化应用水平目的。
目前国内外大气订正方面的成果较多,但由于大气状况瞬息万变,当时、当地的大气对遥感数据影响程度和重要性不同,再加上各卫星的相机波段设置不同,特性不同,研究目的和要求也不同,即使相同地区,大气状况也是复杂多变的,每景图像?#19978;?#26102;的大气状况也存在明显差异,因此很难将某种大气订正方法直接应用于立体测绘卫星三线阵CCD相机。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提供一种易于业务化运行的大气订正方法。
本发明的技术解决方案是:一种立体测绘卫星三线阵CCD相机的大气订正方法,所述的三线阵CCD相机包括前视、正视、后视三台CCD相机;步骤如下:
(1)根据立体测绘卫星三线阵CCD相机的遥?#22411;?#20687;数据,计算每个像元的入瞳处辐射亮度;并依据太阳-大气-地表-大气-相机相互作用的大气辐射传输 方程确立每个像元的地表反射率计算模型;
(2)依据步骤(1)中建立的地表反射率计算模型,设置气溶胶类型和气溶胶的光学厚度,建立不同气溶胶类型和气溶胶光学厚度下的立体测绘卫星三线阵CCD相机每个波段的反演参数查?#20918;恚?
(3)依据Rahman BRDF模型,建立立体测绘卫星前视CCD相机与正视CCD相机地表反射率关系,以及后视CCD相机与正视CCD相机地表反射率关系;并利用上述反演参数查?#20918;懟?#27599;个像元入瞳处辐射亮度以及地表反射率计算模型计算三个相机不同气溶胶类型和光学厚度下的地表反射率;
(4)根据步骤(3)的结果,建立最小化判定函数,确定函数的最优解,根据最优解确定三线阵CCD相机的气溶胶类型和光学厚度;
(5)根据气溶胶光学厚度与波长的关系对步骤(4)的结果进行修正,获得三线阵CCD相机修正后的气溶胶光学厚度;
(6)根据步骤(4)中确定的气溶胶类型以及步骤(5)中获得的气溶胶光学厚度,代入上述反演参数查?#20918;恚?#21033;用地表反射率计算模型计算得到最终的地表反射率。
所述步骤(1)中每个像元的地表反射率计算模型如下:
ρx,λObj(i,j)=Lx,λTOA(i,j)-Lx,λPath(i,j)Tx,λ(i,j)+Sx,λ(i,j)[Lx,λTOA(i,j)-Lx,λPath(i,j)]]]>
其中,为x相机第λ波段像元(i,j)处的地表反射率;x=f,n,b代表立体测绘卫星的前视、正视、后视相机;i,j为每台相机的像元位置;
为x相机第λ波段像元(i,j)处的入瞳处辐射亮度;
为x相机第λ波段像元(i,j)处的路径散射;
T(θx,s)为太阳到地面方向的大气漫射散射因子;T(θx,v)为地面到相机方向的大气漫射散射因子;θx,s为x相机观测时的太阳天顶角;ES,x,?#23435;獂相机第λ波段的大气外太阳光谱辐照度ES;
Sx,λ(i,j)为x相机第λ波段像元(i,j)处的大气球面反照率。
所述步骤(4)中的最小化判定函数如下:
min|Yλ(i,j)|=Σx=f,b[ρx,λR-Obj(i,j)-ρx,λObj(i,j)]2]]>
其中,为利用步骤(3)中的地表反射率关系以及步骤(3)中计算的正视CCD相机不同气溶胶类型和光学厚度下的地表反射率,得到的x相机的地表反射率,x=f,b分别代表立体测绘卫星的前视、后视相机;
为步骤(3)中计算的x相机不同气溶胶类型和光学厚度下的地表反射率。
本发明与现有技术相比有益效果为:
(1)本发明利用立体观测卫星三线阵CCD相机的多角度观测特点,结合大气辐射传输模型和Rahman BRDF模型,实现了立体观测卫星三线阵CCD相机的大气订正。
(2)本发明提出的立体观测卫星三线阵CCD相机大气订正方法,实现了逐像元的大气订正,克服以往大气订正算法中多区域采用单一气溶胶光学厚度的缺点,有效地提高了大气订正的精度。
(3)本发明充分利用了立体观测卫星三线阵CCD相机的多角度观测特点,直接从多角度遥?#22411;?#20687;本身出发进行大气订正,不要任何实际测量参数或其他辅助参数,实施简易、快速、经济使得工程化的业务化大气订正得以实现。
附图说明
图1为本发明方法流程图。
具体实施方式
下面结合附图对本发明做详细说明,如图1所示,本发明步骤具体如下:
(1)相机入瞳处辐射亮度计算:读取立体测绘卫星三线阵CCD相机的遥?#22411;?#20687;数据,依据数据附带文件中的绝对辐射定标系数,按照公式①计算每个 像元的入瞳处辐射亮度:
Lx,λTOA(i,j)=Gainx,λ×DNx,λ(i,j)+Biasx,λ]]>   ①
式中:x=f,n,b代表立体测绘卫星的前视、正视、后视相机;
λ代表每台相机的不同波段;
i,j为每台相机的像元位置(i∈[1,K],j∈[1,L]);
为x相机第λ波段像元(i,j)处的入瞳处辐射亮度;
DNx,λ(i,j)为x相机第λ波段像元(i,j)处的数字DN值;
Gainx,λ和Biasx,?#23435;獂相机第λ波段的绝对辐射定标系数。
(2)地表反射率计算模型确立:依据太阳-大气-地表-大气-相机相互作用的大气辐射传输方程,按照公式②确立每个像元的地表反射率计算模型:
Lx,λTOA(i,j)=Lx,λpath(i,j)+ES,x,λcos(θx,s)T(θx,s)T(θx,v)π×ρx,λObj(i,j)1-ρx,λObj(i,j)Sx,λ(i,j)]]>   ②=Lx,λPath(i,j)+Tx,λ(i,j)ρx,λObj(i,j)1-ρx,λObj(i,j)Sx,λ(i,j)]]>
式中:为x相机第λ波段像元(i,j)处的路径散射;
ES,x,?#23435;獂相机第λ波段的大气外太阳光谱辐照度ES;
θx,s为x相机观测时的太阳天顶角
θx,v为x相机观测时的相机观测天顶角;
φx,s为x相机观测时的太阳方位角;
φx,v为x相机观测时的相机观测方位角;
T(θx,s)为向下方向(太阳到地面)的大气漫射散射因子;
T(θx,v)为向上方向(地面到相机)的大气漫射散射因子;
为x相机第λ波段像元(i,j)处的地表反射率;
Sx,λ(i,j)为x相机第λ波段像元(i,j)处的大气球面反照率;
Tx,λ(i,j)=ES,x,λcos(θx,s)T(θx,s)T(θx,v)π.]]>
对②进行变化得到地表反射率计算模型:
ρx,λObj(i,j)=Lx,λTOA(i,j)-Lx,λPath(i,j)Tx,λ(i,j)+Sx,λ(i,j)[Lx,λTOA(i,j)-Lx,λPath(i,j)]]]>   ③
由上知已确定,故地表反射率由Tx,λ(i,j)和Sx,λ(i,j)参数确定;这3个反演参数是相机波段、太阳天顶角、卫星观测天顶角、太阳方位角、卫星观测方位角、气溶胶类型及气溶胶光学厚度τ的函数即有{Lx,λPath(i,j),Tx,λ(i,j),Sx,λ(i,j)}=f(λ,θx,s,θx,v,φx,s,φx,v,τ),]]>卫星?#19978;?#26102;除气溶胶类型及光学厚度的参数外其他参数都唯一确定,那么在其他参数状态确定的前提下这3个反演参数由气溶胶类型及光学厚度确定。
(3)参数的查?#20918;?#24314;立:依据步骤(2)中确立的参数模型,设置气溶胶类型和气溶胶光学厚度,通过大气辐射传输模型为立体测绘卫星每台线阵CCD相机的每个波段建立以(λ,θx,s,θx,v,φx,s,φx,v,τ)为索引的Tx,λ(i,j)和Sx,λ(i,j)在不同参数状态下的反演参数查?#20918;?#20540;;
(4)气溶胶类型及气溶胶光学厚度确定和修正:依据Rahman BRDF模型,建立立体测绘卫星前后、正视、后视相机的地表反射率关系以确定判定函数,通过判定函数确定气溶胶类型及气溶胶光学厚度,最后利用公式修正气溶胶光学厚度值;
(4.1)基于Rahman BRDF模型计算地表反射率
依据Rahman BRDF模型,立体测绘卫星前视、正视、后视相机的地表反射率关系为:
ρf,λR-Obj(i,j)=ρn,λObj(i,j)(cosθf,s)k-1(cosθf,v)k-1(cosθf,s+cosθf,v)1-kFf(Θ)[1+Hf,λ(G)]]]>
ρb,λR-Obj(i,j)=ρn,λObj(i,j)(cosθb,s)k-1(cosθb,v)k-1(cosθb,s+cosθb,v)1-kFb(Θ)[1+Hb,λ(G)]]]>   ④
式中:k为常数k∈[0,1],表征目标表面各向异性的特性;
F(Θ)为气溶胶散射特性的相函数;
F(Θf)=1-gf2[1+gf2-2gfcos(Θf)]1.5]]>
F(Θb)=1-gb2[1+gb2-2gbcos(Θb)]1.5]]>
Θx为相函数的相位角
cos(Θf)=cos(θf,s)cos(θf,v)-sin(θf,s)sin(θf,v)cos(φf,s-φf,v)]]>
cos(Θb)=-cos(θb,s)cos(θb,v)-sin(θb,s)sin(θb,v)cos(φb,s-φb,v)]]>
gx为常数依据气溶胶类型设置,表征不对称因子;
gf∈[0,1],gb∈[-1,0]
Hx,λ(G)为BRDF的热点效应;
Hx,λ(G)=1-ρn,λObj1+Gx]]>
Gx为几?#25105;?#23376;;
Gx=[tan(θx,s)]2+[tan](θx,v)2-2tan(θx,s)tan(θf,v)cos(φx,s-φx,v)]]>
(4.2)判定函数确立
利用Rahman BRDF模型确立的立体测绘卫星前后、正视、后视相机地表反射率关?#26723;?#21040;的与大气辐射传输模型确立的地表反射率的等价关系,建立最小化判定函数为:
min|Yλ(i,j)|=Σx=f,b|[ρx,λR-Obj(i,j)-ρx,λObj(i,j)]2|]]>   ⑤
其中,为利用步骤(3)中的地表反射率关系以及步骤(3)中计算的正视CCD相机不同气溶胶类型和光学厚度下的地表反射率,得到的x相机的地表反射率,x=f,b分别代表立体测绘卫星的前视、后视相机;为步骤(3)中计算的x相机不同气溶胶类型和光学厚度下的地表反射率。
(4.3)气溶胶类型及气溶胶光学厚度确立
利用立体测绘卫星?#19978;?#26102;获得的θx,s,θx,v,φx,s,φx,v,通过步骤(3)反演参数查?#20918;?#20540;获得以(λ,θx,s,θx,v,φx,s,φx,v)索引的不同气溶胶类型及光学厚度下的 Tx,λ(i,j)和Sx,λ(i,j),并根据大气辐射传输模型确立公式③计算出立体测绘卫星三线阵CCD相机对应的地表反射?#35797;?#26681;据步骤(4)中确定的地表反射率关系④计算出不同气溶胶类型及光学厚度下正视相机地表反射率对应的前视、后视相机的地表反射率和将上述两种方式计算多种气溶胶类型和光学厚度条件下的立体测绘卫星三线阵CCD相机的地表反射率和代入方程⑤式获得各种气溶胶类型和光学厚度条件下的判定函数值,将满足判定函数最小值对应条件下立体测绘卫星三线阵CCD相机各地表反射率对应的气溶胶类型和光学厚度确定为函数最优解;
(4.4)气溶胶类型及气溶胶光学厚度确修正
根据步骤(4.3)确定的立体测绘卫星三线阵CCD相机对应的气溶胶光学厚度,通过气溶胶光学厚度与波长满足的公式τ(λ)=?#21147;?β进行最小二乘法拟合以确定公式系数α、β,最后通过已知参数α、β的公式获得修正后立体测绘卫星三线阵CCD相机不同波段的气溶胶光学厚度。
(5)大气订正处理:根据步骤(4)获得的气溶胶类型及修正的气溶胶光学厚度,结合立体观测卫星三线阵CCD相机的几何参数θx,s,θx,v,φx,s,φx,v,查?#39029;?#20197;(λ,θx,s,θx,v,φx,s,φx,v)为索引的Tx,λ(i,j)和Sx,λ(i,j)反演参数值,即可通过方程③式获得经过大气订正后的立体测绘卫星三线阵CCD相机对应的地表反射率。
本发明未详细说明部分属于本领域技术人员公知常识。

关于本文
本文标题:一种立体测绘卫星三线阵CCD相机的大气订正方法.pdf
链接地址:http://www.pqiex.tw/p-6124576.html
关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

[email protected] 2017-2018 zhuanlichaxun.net网站版权所有
经营许可证编号:粤ICP备17046363号-1 
 


收起
展开
平码五不中公式规律 体彩内蒙古十一选五 2008排列五走势图 高频彩彩民 广东十一选五走势 最近双色球选号 大乐透杀跨度公式 广东11选5人5遗漏 2011短线股票推荐 新疆十一选五爱彩乐 聪明组六组合