主权项 |
一种基于高光谱载荷的多光谱载荷无场地交叉定标方法,其特征在于包括以下步骤:(1)分别获取高光谱载荷和多光谱载荷在某试验区域的遥感图像,并在两幅遥感图像上选择同一个感兴趣区域,记为第一感兴趣区域;(2)根据该感兴趣区域的高光谱载荷遥感图像数据得到高光谱载荷各波段的等效入瞳处辐亮度值;根据该感兴趣区域的多光谱载荷遥感图像数据得到多光谱载荷每个波段在感兴趣区域的像素值;(3)根据高光谱载荷和多光谱载荷各波段的波长范围,确定多光谱载荷的每个波段与高光谱载荷各波段的对应关系;(4)根据高光谱载荷和多光谱载荷各波段的光谱响应函数,结合热红外光谱仪对该感兴趣区域实测的热红外光谱,计算得到在该感兴趣区域多光谱载荷每个波段相对于高光谱载荷对应波段群的光谱匹配因子;设在某感兴趣区域,多光谱载荷的第i个波段与高光谱载荷第i<sub>1</sub>,i<sub>2</sub>,…,i<sub>n</sub>个波段相对应,则在该感兴趣区域多光谱载荷第i个波段相对于高光谱载荷对应波段群的光谱匹配因子的计算方式为:(4.1)根据多光谱载荷第i个波段的光谱响应函数,结合热红外光谱仪对该感兴趣区域实测的热红外光谱,利用公式<maths num="0001"><math><![CDATA[<mrow><msub><mi>L</mi><mrow><mi>B</mi><mi>i</mi></mrow></msub><mo>=</mo><mfrac><mrow><mo>∫</mo><msub><mi>L</mi><mi>M</mi></msub><mrow><mo>(</mo><mi>λ</mi><mo>)</mo></mrow><msub><mi>f</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>λ</mi><mo>)</mo></mrow><mi>d</mi><mi>λ</mi></mrow><mrow><mo>∫</mo><msub><mi>f</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>λ</mi><mo>)</mo></mrow><mi>d</mi><mi>λ</mi></mrow></mfrac></mrow>]]></math><img file="FDA0001145515230000011.GIF" wi="491" he="174" /></maths>计算得到多光谱载荷的第i个波段在该感兴趣区域的模拟入瞳处辐亮度值L<sub>Bi</sub>;其中f<sub>B</sub>(λ)为多光谱载荷第i个波段的光谱响应函数,L<sub>M</sub>(λ)为热红外光谱仪实测的第i个波段在该感兴趣区域的热红外光谱,λ为第i个波段中的波长;(4.2)根据高光谱载荷第i<sub>1</sub>,i<sub>2</sub>,…,i<sub>n</sub>个波段的光谱响应函数,结合热红外光谱仪对该感兴趣区域实测的热红外光谱,利用公式<maths num="0002"><math><![CDATA[<mrow><msub><mi>L</mi><mi>A</mi></msub><mrow><mo>(</mo><msub><mi>λ</mi><mrow><mi>i</mi><mi>m</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mo>∫</mo><msub><mi>L</mi><mi>M</mi></msub><mrow><mo>(</mo><msub><mi>λ</mi><mrow><mi>i</mi><mi>m</mi></mrow></msub><mo>)</mo></mrow><msub><mi>f</mi><mi>A</mi></msub><mrow><mo>(</mo><msub><mi>λ</mi><mrow><mi>i</mi><mi>m</mi></mrow></msub><mo>)</mo></mrow><mi>d</mi><mrow><mo>(</mo><msub><mi>λ</mi><mrow><mi>i</mi><mi>m</mi></mrow></msub><mo>)</mo></mrow></mrow><mrow><mo>∫</mo><msub><mi>f</mi><mi>A</mi></msub><mrow><mo>(</mo><msub><mi>λ</mi><mrow><mi>i</mi><mi>m</mi></mrow></msub><mo>)</mo></mrow><mi>d</mi><mrow><mo>(</mo><msub><mi>λ</mi><mrow><mi>i</mi><mi>m</mi></mrow></msub><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0001145515230000021.GIF" wi="682" he="175" /></maths><maths num="0003"><math><![CDATA[<mrow><msub><mi>L</mi><mrow><mi>A</mi><mi>i</mi></mrow></msub><mo>=</mo><munderover><mo>Σ</mo><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>L</mi><mi>A</mi></msub><mrow><mo>(</mo><msub><mi>λ</mi><mrow><mi>i</mi><mi>m</mi></mrow></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001145515230000022.GIF" wi="323" he="126" /></maths>计算得到高光谱载荷第i<sub>1</sub>,i<sub>2</sub>,…,i<sub>n</sub>个波段组成的波段群在该感兴趣区域的模拟入瞳处辐亮度值L<sub>Ai</sub>;其中L<sub>A</sub>(λ<sub>im</sub>)为高光谱载荷第i<sub>m</sub>个波段在该感兴趣区域的模拟入瞳处辐亮度值,f<sub>A</sub>(λ<sub>im</sub>)为高光谱载荷第i<sub>m</sub>个波段的光谱响应函数,L<sub>M</sub>(λ<sub>im</sub>)为热红外光谱仪实测的第i<sub>m</sub>个波段在该感兴趣区域的热红外光谱,λ<sub>im</sub>为第i<sub>m</sub>个波段中的波长,i∈自然数,n∈自然数;(4.3)利用公式ρ<sub>i</sub>=L<sub>Bi</sub>/L<sub>Ai</sub>计算得到多光谱载荷第i个波段相对于高光谱载荷第i<sub>1</sub>,i<sub>2</sub>,…,i<sub>n</sub>个波段组成的波段群的光谱匹配因子ρ<sub>i</sub>;(5)根据步骤(2)得到的高光谱载荷各波段在该感兴趣区域的等效入瞳处辐亮度值,结合步骤(4)得到的多光谱载荷每个波段相对于高光谱载荷对应波段群的光谱匹配因子,计算得到多光谱载荷每个波段在该感兴趣区域的等效入瞳处辐亮度值;(6)在步骤(1)得到的两幅遥感图像上同时选择另一个感兴趣区域,记为第二感兴趣区域,重复步骤(2)—(5),得到多光谱载荷每个波段在另一个感兴趣区域的像素值和等效入瞳处辐亮度值;(7)根据多光谱载荷每个波段在两个不同感兴趣区域的像素值和等效入瞳处辐亮度值,计算得到多光谱载荷每个波段的辐射定标系数,从而完成多光谱载荷的无场地交叉定标。 |