发明名称 |
一种干涉型高光谱成像仪的光谱漂移检测方法 |
摘要 |
一种干涉型高光谱成像仪的光谱漂移检测方法,(1)将该高光谱成像仪光谱数据转化为各波段入瞳处辐亮度,作为待检验光谱辐亮度;(2)获得卫星高度的入瞳处辐亮度;(3)通过高光谱成像仪响应函数得到各波段入瞳处辐亮度,作为参考光谱辐亮度;(4)比对待检验光谱辐亮度与参考光谱辐亮度在大气吸收峰的位置,通过判断函数确定待检验光谱辐亮度与参考光谱辐亮度的偏移量是否达到约束条件,如果达到转步骤(5);否则调整高光谱成像仪响应函数中波长和带宽的调整量,从步骤(3)开始循环判断;(5)将满足约束条件时确定的高光谱成像仪光谱漂移量,对高光谱成像仪干涉数据重新进行傅立叶逆变换得到高光谱成像仪去除光谱漂移后的光谱数据图。 |
申请公布号 |
CN102322957B |
申请公布日期 |
2013.03.20 |
申请号 |
CN201110194801.2 |
申请日期 |
2011.07.12 |
申请人 |
中国资源卫星应用中心 |
发明人 |
王爱春;傅俏燕;闵祥军;潘志强;康倩;林军 |
分类号 |
G01J3/45(2006.01)I |
主分类号 |
G01J3/45(2006.01)I |
代理机构 |
中国航天科技专利中心 11009 |
代理人 |
庞静 |
主权项 |
一种干涉型高光谱成像仪的光谱漂移检测方法,其特征在于步骤如下:(1)读取高光谱成像仪干涉数据图,对该高光谱成像仪干涉数据进行离散傅立叶逆变换得到高光谱成像仪光谱数据图;通过绝对辐射定标系数,将该高光谱成像仪光谱数据转化为高光谱成像仪的各波段入瞳处辐亮度,将其作为待检验光谱辐亮度;(2)依据卫星成像时地面实测地表参数和大气参数,利用辐射传输模型获得卫星高度的入瞳处辐亮度;(3)通过高光谱成像仪响应函数得到高光谱成像仪的各波段入瞳处辐亮度,将其作为参考光谱辐亮度;(4)比对卫星成像时待检验光谱辐亮度与参考光谱辐亮度在大气吸收峰的位置,通过判断函数确定待检验光谱辐亮度与参考光谱辐亮度的偏移量是否达到约束条件,如果达到则转步骤(5);否则,调整高光谱成像仪响应函数中波长和带宽的调整量,从步骤(3)开始循环判断;所述的判断函数如下:min g=min{(1‑κ)σ+κα}其中:κ为权重系数比,κ∈[0,1];n是高光谱成像仪波段的总数; <mrow> <mi>σ</mi> <mo>=</mo> <msqrt> <mfrac> <mn>1</mn> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </mfrac> <munderover> <mi>Σ</mi> <mrow> <mi>λ</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msup> <mrow> <mo>(</mo> <msub> <mi>L</mi> <mi>unchecked</mi> </msub> <mrow> <mo>(</mo> <mi>λ</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>L</mi> <mi>referrence</mi> </msub> <mrow> <mo>(</mo> <mi>λ</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> </msqrt> </mrow> <mrow> <mi>α</mi> <mo>=</mo> <msup> <mi>cos</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mrow> <mo>(</mo> <mfrac> <mrow> <munderover> <mi>Σ</mi> <mrow> <mi>λ</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>L</mi> <mi>unchecked</mi> </msub> <mrow> <mo>(</mo> <mi>λ</mi> <mo>)</mo> </mrow> <msub> <mi>L</mi> <mi>referrence</mi> </msub> <mrow> <mo>(</mo> <mi>λ</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msqrt> <munderover> <mi>Σ</mi> <mrow> <mi>λ</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msup> <mrow> <mo>(</mo> <msub> <mi>L</mi> <mi>unchecked</mi> </msub> <mrow> <mo>(</mo> <mi>λ</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> <munderover> <mi>Σ</mi> <mrow> <mi>λ</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msup> <mrow> <mo>(</mo> <msub> <mi>L</mi> <mi>referrence</mi> </msub> <mrow> <mo>(</mo> <mi>λ</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> </msqrt> <mo>/</mo> <mfrac> <mi>π</mi> <mn>2</mn> </mfrac> </mrow> </mfrac> <mo>)</mo> </mrow> </mrow>Lunchecked(λ)对应待检验光谱波段的入瞳处辐亮度Lreferrence(λ)对应参考光谱波段的入瞳处辐亮度;λ为高光谱成像仪波段的波长;(5)将满足约束条件时确定的高光谱成像仪光谱漂移量,对高光谱成像仪干涉数据重新进行傅立叶逆变换得到高光谱成像仪去除光谱漂移后的光谱数据图。 |
地址 |
100094 北京市海淀区永丰产业基地丰贤东路5号 |