发明名称 一种凝视卫星面阵CCD相机的相对辐射校正方法
摘要 一种凝视卫星面阵CCD相机的相对辐射校正方法。首先建立面阵CCD相机探元的辐射响应模型;其次利用凝视卫星面阵CCD在没有月光的深夜对远离陆地的大面积深海进行实时、定点观测,通过多次获得的遥感图像确立噪声大小,去噪后得到无噪声的遥感图像数据;然后在晨昏时,利用凝视卫星面阵CCD相机在微弱均匀散射的天空光下对远离陆地的大面积的深海进行连续成像,通过多次获得的去除噪声的遥感图像计算相对辐射定标系数;最终得到凝视卫星面阵CCD相机的相对辐射校正图像。本发明解决了凝视卫星面阵CCD相机在轨运行无相对辐射校正方法的技术问题。
申请公布号 CN104065892A 申请公布日期 2014.09.24
申请号 CN201410287058.9 申请日期 2014.06.24
申请人 中国资源卫星应用中心 发明人 王爱春;傅俏燕;闵祥军;陆书宁;潘志强;李晓进;韩启金;张学文;刘李;李照洲
分类号 H04N5/357(2011.01)I 主分类号 H04N5/357(2011.01)I
代理机构 中国航天科技专利中心 11009 代理人 安丽
主权项 一种凝视卫星面阵CCD相机的相对辐射校正方法,其特征在于步骤如下:1)建立面阵CCD相机探元的辐射响应模型:DN(i,j)={f[L(i,j)]+N<sub>DC</sub>(i,j)+N<sub>S</sub>(i,j)+N<sub>R</sub>(i,j)}×G+N<sub>Q</sub>(i,j);式中,i、j为面阵CCD相机的探元位置(i∈[1,K],j∈[1,L]);K为面阵CCD相机的行探元总数,L为面阵CCD相机的列探元总数;DN(i,j)为探元(i,j)的数字DN值;f[L(i,j)]为探元(i,j)接收的辐射亮度L(i,j)函数;N<sub>DC</sub>(i,j)为探元(i,j)的暗电流噪声;N<sub>S</sub>(i,j)为探元(i,j)的散粒噪声;N<sub>R</sub>(i,j)为探元(i,j)的读出噪声;G为电路输出的综合增益;N<sub>Q</sub>(i,j)为探元(i,j)的量化噪声;将面阵CCD相机探元的辐射响应模型整理简化得到:DN(i,j)=f<sub>G</sub>[L(i,j)]+N(i,j);式中,f<sub>G</sub>[L(i,j)]=f[L(i,j)]×G;N(i,j)=[N<sub>DC</sub>(i,j)+N<sub>S</sub>(i,j)+N<sub>R</sub>(i,j)]×G+N<sub>Q</sub>(i,j);2)利用凝视卫星面阵CCD在没有月光的深夜对远离陆地的大面积深海进行实时、定点观测,通过N次获得的遥感图像确立噪声大小;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>N</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>DN</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow><mi>N</mi></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000526022010000011.GIF" wi="547" he="199" /></maths>式中,DN<sub>n</sub>(i,j)为第n景深夜成像数据中探元(i,j)的数字DN值,n∈[1,N];3)将步骤2)中得到的噪声大小,在凝视卫星面阵CCD相机获得的遥感图像数字DN值中去除,即得到无噪声的遥感图像数据<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mover><mi>DN</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mi>DN</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><mi>N</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>f</mi><mi>G</mi></msub><mo>[</mo><mi>L</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>]</mo><mo>;</mo></mrow>]]></math><img file="FDA0000526022010000012.GIF" wi="870" he="87" /></maths>4)在晨昏时,利用凝视卫星面阵CCD相机在微弱均匀散射的天空光下对远离陆地的大面积的深海进行连续成像,通过M次获得的去除噪声的遥感图像计算相对辐射定标系数;41)求取M景去除噪声的所有探元数字DN值样本均值μ<sub>DN</sub>:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>&mu;</mi><mi>DN</mi></msub><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><msub><mi>&mu;</mi><msub><mi>DN</mi><mi>m</mi></msub></msub></mrow><mi>M</mi></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000526022010000021.GIF" wi="340" he="204" /></maths>式中,<img file="FDA0000526022010000022.GIF" wi="539" he="218" /><img file="FDA0000526022010000023.GIF" wi="241" he="88" />为第m景中探元(i,j)去除噪声的数字DN值,m∈[1,M];42)将M次获得的去除噪声的所有探元数字DN值样本均值μ<sub>DN</sub>作为凝视卫星面阵CCD相机探元的输入能量期望值f<sub>G</sub>[L(i,j)]=μ<sub>DN</sub>,则第m景中探元(i,j)的相对辐射定标系数为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>R</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>DN</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow><msub><mi>&mu;</mi><mi>DN</mi></msub></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000526022010000024.GIF" wi="446" he="143" /></maths>其中DN<sub>m</sub>(i,j)为第m景晨昏成像数据中探元(i,j)的数字DN值,m∈[1,M];43)对探元(i,j)获得的M个相对辐射定标系数进行平均,最终得到凝视卫星面阵CCD相机各探元的相对辐射定标系数为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>R</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>R</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow><mi>M</mi></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000526022010000025.GIF" wi="437" he="195" /></maths>5)对凝视卫星面阵CCD相机获得的遥感图像数据,利用步骤2)中得到的噪声大小和步骤4)中得到的相对辐射校正系数,得到凝视卫星面阵CCD相机的相对辐射校正图像<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>DN</mi><mi>R</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mi>R</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow></mfrac><mo>&times;</mo><mo>[</mo><mi>DN</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><mi>N</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>]</mo><mo>=</mo><msub><mi>f</mi><mi>G</mi></msub><mo>[</mo><mi>L</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>]</mo><mo>.</mo></mrow>]]></math><img file="FDA0000526022010000026.GIF" wi="1093" he="138" /></maths>
地址 100094 北京市海淀区永丰产业基地丰贤东路5号