空间CCD相机辐射标定方法的应用研究
2024-04-13
来源:易榕旅网
维普资讯 http://www.cqvip.com 2004年第l8卷第4期 (总第5O期) 测 试 技 术 学 报 JOURNAL OF TEST AND MEASUREMENT TECHNOLOGY VOj.18 NO.4 2004 (Sum NO.50) 文章编号:1 671—7449(2004)04—0325—05 空间CCD相机辐射标定方法的应用研究 吴国栋 ,张晓辉 , 韩昌元 , 赵志忠 (1.中国科学院长春光学精密机械与物理研究所,吉林长春130022; 2.辽宁工学院材料与化学工程学院,辽宁锦州121000) 摘要: 为了使空间CCD相机对不同辐亮度的地面目标获得合适的推扫图像,根据空间相机实际工作环 境计算出不同使用条件下相机入口处辐亮度.在地面建立了CCD相机辐亮度标定装置,并对此装置的传递 可靠性指标进行了测试,总误差不大于5.3%,可以在地面对CCD相机进行辐射标定.实验结果表明:相机 入口处辐亮度计算准确,通过地面标定.可在实际工作环境下获得最佳的最终图像. 关键词:空间CCD相机;辐亮度;标定;TDI 中图分类号:TB871 文献标识码:A Radiometric Calibration of the Space CCD Camera WU Guo dong ,ZHANG Xiao—hui ,HAN Chang—yuan ,ZHAO Zhi—zhong (1.Changchun Institute of Optics,Fine Mechanics and Physics,Chinese Academy of Sciences・Changchun 130022,China;2.Faculty of Materials and Chemistry。Liaoning Institute of Technology,Jingzhou 121000,China) Abstract: In order to acquire precision pushbroom images at different radiance of the earth ob)ect,and according to the working environment of space camera.The radiance at camera entrance port in difference space condition was calculated.The radiometric calibration instrument is established and the transfer reliability of the instrument iS tested,and its uncertainty of radiance is less than 5.3 9/6 which meets the demand of calibrating radiance in the lab.The results of experiment show that the calculation of the radiance at camera entrance port is accurate and the pushbroom images with high quality can be obtained through the calibration on the ground. Key words:the space CCD camera;radiance;calibration;TDI 早在20世纪50年代初中国科学院仪器馆(长春光机所前身)就开展了光度基准和光辐射标准的研 究工作,在这一期间解决了红外辐射标准源问题,还为国家的辐射温标的建立打下了技术基础.进入2O 世纪70年代,以多光谱相机为开端,开展多种遥感器的研制,并进行有关的辐射测量与标定.从20世 纪80年代到90年代,长春光机所开始承担重大的光学遥感器课题,遥感器的辐射标定成为迫切需要解 决的问题.为此进行了相应的研究,建立了光学遥感器辐射标定系统L1].20世纪70年代后期,开展了短 波段光学研究,承担了重大的短波段空间光学仪器研制任务,并将光辐射测量与标定工作延伸到紫外与 真空紫外区域[ 引. 空间CCD相机在搭载飞船发射后,用于获得地面目标图像.为了对不同辐亮度的地面目标获得曝 光量合适的推扫图像,空间CCD相机具有96,48,24,12,6等5挡时间延迟积分级数可调TDI CCD, 收稿日期:2004—03—24 作者简介:吴国栋(1965一),男,高级工程师,主要从事精密检测与装校技术研究 维普资讯 http://www.cqvip.com 326 测 试技术 学 报 2004年第4期 并配有增益可调的视频放大器.合理的选择TDI CCD时间延迟积分级数以及放大器的增益是保证获得 最佳最终图像的必要手段.根据CCD相机的使用要求,通过对相机工作时不同条件下的像面曝光量进 行了估算,确定了相机标定时所需辐亮度源的范围.在地面,对CCD延时积分级数和放大器增益进行精 细的调整必须提供稳定的辐亮度光源,来模拟不同辐亮度的地面目标.为此在地面对辐射标定装置进行 了精细标定,建立起辐射标定装置出I:I辐亮度和探测器数字化输出之间定量关系,并对计量传递误差进 行了分析,标定结果说明该装置能提供满足要求的辐亮度,并在CCD相机地面辐射标定工作中起到了 很好的作用. 1 相机人口处辐亮度的计算 相机像面曝光量是由像面辐照度和像元的曝光(积分)时间的乘积.要获得合适的曝光量必须对像面 辐照度进行估算.要知道像面辐照度必须求得相机入I:I处辐亮度,相机入I:I处辐亮度为地物产生的辐亮 度和天空后向散射的辐亮度之和,下面对相机入El处辐亮度进行计算. 地面的景物是被太阳直射光和天空散射光照亮的,它们产生的地表辐照度E。(相机工作波段zl= 0.5 m~O.8 m积分值,相机对此波段光谱辐射有较好的响应)与太阳天顶角 有关.表l列出大气能 见度为23 km时,地表辐照度E。值,单位w/m。(用美国空军地球物理实验室AFGL的大气传输软件 LOWTRAN计算). 表1地表辐照度Eo Tab.1 Earth surface irradlance Eo 0/。 Eo/W・m一 O 420 10 410 20 f 30 I 40 { 50 I 60 l 70 { 75 l 80 390 f 360 l 31 0 l 250 l 180 l 110 { 72 l 39 相机对地面景物进行拍摄时,地面亮度在相机入I:I处产生的辐亮度 与地物的平均反射率JD和大 气平均透过率 有关,表2列出大气能见度为23 km,大气平均透过率 =O.75,在不同0和|D条件下 求得的B。=pE。G/zr的值( =0.5 tLm ̄0.8 tLm积分值),单位W/m。Sr. 表2地面亮度在相机入口处产生的辐亮度 Tab.2 Radiance Bo at camera entrance port 0/。 户 O 1O 20 30 40 5O 60 70 80 9O 1.0 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0.05 100 90.3 8O.2 70.2 6O.2 50.1 40.1 30.1 20.1 10.0 5.01 97.9 88.1 78.3 68.5 58.7 48.9 39.2 29.4 19.6 9.79 4.89 93.1 83.8 74.5 65.2 55.9 46.6 37.2 27.9 l8.6 9.31 4.66 86.0 77.4 68.8 6O.2 51.6 43.0 34.4 25.8 1 7.2 8.60 4.30 74.0 66.6 59.2 51.8 44.4 37.0 29.6 22.2 14.8 7.40 3.70 59.7 53.7 47.8 41.8 35.8 29.8 23.9 1 7.9 l1.4 5.97 2.98 43.0 38.7 34.4 30.1 25.8 21.5 17.2 12.9 8.60 4.30 2.15 26.3 23.7 21.0 18.4 1 5.8 l3.1 10.5 7.89 5.26 2.63 1.31 1 7.2 1 5.5 l3.8 1 2.0 10.3 8.60 6.88 5.16 3.44 1.72 0.86 9.13 8.38 7.45 6.52 5.59 4.66 3.72 2.79 】.86 0.93 0.46 相机入口处辐亮度还有天空后向散射产生的辐亮度B ,表3列出大气能见度为23 km,不同太阳天 顶角0下求得B ( 一0.5/ ̄m---O.8 m积分值),单位w/m。Sr(由LOWTRAN软件计算). 表3天空后向散射产生的辐亮度B Tab.3 Radiance B via atmosphere scattering 0/。 j 0 1O 10.1 20 9.9O 30 {40 9.62 l 9.24 50 j 60 j 70 j 75 j 80 8.74}8.06 I 7.01 I 6.00 I 4.85 B /W・m一 Sr l 10.1 维普资讯 http://www.cqvip.com (总第5O期) 空间CCD相机辐射标定方法的应用研究(吴函一海等) 327 从而求得在不同0和p条件下相机入口处总辐亮度B=B。+B ( 一0.5 m~0.8 m积分值)的值,单 位w/m。Sr.如表4所示. 表4相机入口处总辐亮度B Tab.4 Summary radiance B at camera entrance port 8 P O 1O 20 3O 40 50 60 70 80 90 1。0 0.9 0.8 0.7 0。6 0.5 0.4 0.3 0.2 0.1 110 100 90.3 80.3 70.3 60.2 50.2 40.2 30.2 20.1 108 98.2 88。3 78.6 68.8 59.0 49.3 39.5 29.7 19.9 103 93、7 84.4 75。1 65.8 56.5 47、1 37.8 28.5 19.2 95.6 87.0 78.2 69。8 61.2 52.6 43.9 35.4 26.8 18。2 83.2 75.8 68.4 61.0 53.6 46.2 38.8 31.4 24.0 16.6 68.4 62.4 56。5 5O.5 44.5 38.5 32.6 26.6 2O.1 14.7 51.1 46.8 42.5 38.2 33.9 29.6 25.3 21.0 16.7 12.4 33.3 30.7 28.0 25.4 22.8 20.1 17.5 14.9 12.3 9.64 23.2 21.5 19.8 18.0 16。3 14.6 12.9 11.2 9.44 7.72 14.2 13.2 12.3 11.4 10.4 9.51 8.57 7.64 6.71 5。78 0.05 15、1 15.0 14。6 13.9 12.9 11.7 10。2 8.32 6.86 5。31 通过以上分析计算,可以得知需要在地面模拟的辐亮度范围,为辐射标定装置的建立提供了依据. 2地面标定装置的建立和标定 2.1标定设备组成 标定设备主要由积分球、光谱辐亮度计、标准灯、标准反射率白板组成,积分球由球体照明灯和辐 亮度监测用硅探测器组成Ⅲ,当所有光源点亮时,积分球出口处辐亮度可达138.36 W/m。Sr( 一0.5 m ~0.8 m积分值),满足相机摄影时,入口处最大辐亮度110 w/m。Sr的要求.标定设备如图1所示. 硅探测器 积分球 白板 图1 CCD标定装置示意图 Fig.1 The cha rt of CCD calibrating instrument 2.2辐亮度的计量传递 辐亮度的传递是通过经过检定的标准灯标定辐亮度计,求出辐亮度与辐亮度计的输出电压值关系, 并与积分球内设置的硅探测器相比较,确定积分球内监控用硅探测器输出电压与积分球辐亮度之间的关 系,通过调整光源亮度实现积分球辐亮度变化. 2.2.1辐亮度计的标定L4 将经过检定的标准灯置于经过检定的标准白板前1 m处,在与标准白板成45。方向放置辐亮度计,按 照检定证书要求给定标准灯电压和电流,并测出此时辐亮度计的输出电压值为0.1316 V,再根据标准灯检 定数据E和标准白板反射率lD,求得辐亮度B—pE/ ̄r,其中E( 一0.5 m~0.8 m积分值)一∑ ×A2, lD一0.950,此时B一4.022(W/m Sr),同理分别测出0.5 m~2 m多个位置的辐亮度计的输出电压值,在 不同位置时B :B。(Lf/L。) ,其中B。一4.022(w/m Sr),L。一1 m.根据不同位置的辐亮度计的输出电压 值和辐亮度值,拟合出辐亮度B和辐亮度计的输出电压 f的关系:B一0.159+28.909 41 Vf. 维普资讯 http://www.cqvip.com 328 测 试技术学报 2004年第4期 2.2.2积分球出口处辐亮度标定 用标定过的辐亮度计测积分球出口处辐亮度,并与积分球内设置的硅探测器相比较,确定积分球内 监控用硅探测器输出电压 与积分球辐亮度 之间的关系:Bi 一0.445 33+27.325 95V . 因此,通过改变12个灯的点亮数目及改变灯的电压可得到由高到低的连续可调辐亮度,表5是积 分球输出辐亮度值与辐亮度计以及积分球内监控用硅探测器输出电压的关系. 表5输出辐亮度值与辐亮度计及探测器输出电压的关系 Tab.5 The relation of output of radiance and output voltage of luminance meter and output voltage of sensor 辐亮度B/W・m Sr一1 110 1O0 90 80 70 60 vl| 3.799 489 3.453 581 3.107 673 2.761 765 2.415 856 2.069 948 /V 4.041 774 3.675 822 3.309 869 2.943 917 2.577 964 2.212 O12 辐亮度B/W・m Sr 5O 40 3O 2O 1O 5 vl| 1.724 040 1.378 132 1.O32 224 0.686 316 O.340 408 0.167454 /V 1.846 059 1.480 107 1.114 154 O.748 2O1 O.382 249 0.199273 3 标定装置可靠性实验和误差分析[5] 为了对CCD相机进行辐射标定,建立了上述标定系统.以国家计量院检定过的光谱辐照度标准灯 作为标定依据,利用光谱辐照度标准灯和标准白板建立小面积光谱辐亮度过渡标准.为了将这样建立起 来的小面积光谱辐亮度标准传递到大开口积分球光源,应用小视场光谱辐亮度计来传递这一标准,建立 起光谱辐亮度计电压和积分球出口光谱辐亮度关系.对标定装置的可靠性和误差进行分析如下: 3.1光谱辐亮度计标定误差 用光谱辐照度标准灯和标准白板建立小面积光谱辐亮度过渡标准,灯丝电流调至8.500 A,用小视 场光谱辐亮度计采集5组数据,每组采集30个数据,计算每组输出电压数据的平均值,计算出每组输出 电压间的误差(光谱辐亮度计标定误差)为0.25 (峰一峰值). 3.2积分球光源稳定性 为了提高光源的稳定性,采用德国OSRAM公司的Halogen Photo Optic Lamp,寿命300 h,并对 灯泡进行风冷,由12路稳压电源供电,开机30 min预热,使积分球达到热平衡.稳定性实验中,每隔 20 min测一次,共计4 h,获得13组数据如表6所示. 表6积分球光源的稳定性 Tab.6 Stability of a spherical integrating source 时间间隔/min 0 vl| 3.556 620 时间间隔/min 140 vt| 3.585 776 20 40 3.567 117 3.581 413 160 18O 3.585 37O 3.584 571 6O 80 3.58O O33 3.576 623 200 220 3.581 606 3.581 522 1O0 12O 3.58O 176 3.585 512 240 3.581 900 求得积分球光源辐射不稳定性误差为0.8 (峰一峰值). 3.3积分球出口处亮度均匀性 为了检测积分球出口处的亮度均匀性,将辐亮度计置于积分球出口处的中心及左右5个位置分别测 试,获得5组数据如表7所示. 维普资讯 http://www.cqvip.com (总第5O期) 空间CCD相机辐射标定方法的应用研究(吴国栋等) 表7积分球出口亮度均匀性 Tab.7 Uniformity of outlet with a spherical integrating source 329 辐亮度计测试位置/mm 。 vIf ——200 3,577 O31 —1OO 3.577 024 lf 0 3.567 117 100 3.576 491 200 3.573 891 求得亮度均匀性误差为0.3 (峰一峰值). 此方法误差的主要来源是:光谱辐照度灯的不确定度4.8%;光谱辐亮度计标定误差0.25 (峰一峰 值);积分球的稳定性误差0.8 9/5(峰一峰值);亮度均匀性误差0.3 9/5(峰一峰值).除上述主要因素外,还有 辐亮度计线性等其它引起的不确定度各项因素的影响按2 计算,用这种方法标定CCD相机的总误差 (均方根值)△ 5.3 .此精度对于获取几何信息为主要任务的CCD相机是可以满足要求的. 4 实际应用 根据积分球输出辐亮度与硅探测器输出电压的关系,当在6O。太阳天顶角地面反射率0.7时,B一 38.2 w/m。sr,求得V =1.414 236 V.调整光源使硅探测器输出电压为1.414 236 V,设置相机CCD固 定的行周期为0.2 ms,增益为32 H,将级数由6,12,24,96级依次改变,找出对应图像平均灰度 130 DN~240 DN的级数.在这个级数下,其他条件不变,改变增益设置,使4块CCD平均灰度为 220 DN,从而确定此条件下的级数和增益,依据上述标定方法,确定对应4O。,2O。,1O。,对太阳高角地 面反射率为0.7,0.2等条件下级数和增益参数进行设置.这样确保了相机工作时设置合理的级数和增 益,也为获得最佳图像提供了保证. 除此之外,在地面处理阶段还可以根据相机像面照度的不均匀性,对CCD相机焦面组件 TDICCD2048像元的响应函数进行测试,以达到对图像进行精细校正. 为了克服CCD器件响应随时间的推移而发生的退化衰减,在相机中设置了一组标定灯,定期比较 CCD器件的响应与标定灯亮度之间的关系,通过标定可以获得CCD对标定灯辐亮度响应的初始值. 5 结束语 通过对相机入V1处辐亮度的计算以及空间CCD相机地面标定装置的建立,较好的模拟了相机工作 时入口处辐亮度,在地面通过对CCD相机进行精细的辐射标定,空间实际应用中获得了最佳图像.但对 辐射测量精度要求较高时,利用自校准型探测器 等绝对辐射测量探测器可以提高辐射标定精度. 6致谢 在空间CCD相机地面标定过程中,得到了长春光机与物理研究所任建伟、张立国、万志等老师的大 力帮助,在此表示感谢! 参考文献: [1]崔敦杰.光学遥感器的辐射定标[A].纪念建所五十周年论文及论文摘要汇编(一)Ec9.陈星旦主编,长春:长春光 学精密机械与物理研究所,2002.148一l54. Eel Tang Y G,1.i F T.Investigation on the characteristics of synchrotron radiation emitted by the Hefei 800 MeVsy— nchrotron radiation facility[J'].Nuclear Instrument and Methods in physics Research・l994・337(A):598—602. [3] 周胜利.积分球在实验室内用于空间遥感器的辐射定标[J].航天返回与遥感。l998,(1):29—34. [4]李振国,魏青.光辐射与测量原理与技术[M].北京:中国建材工业出版社.1 998.37—39. [5]费业泰.误差理论与数据处理[M].北京:机械工业出版社.2000.52—55. [6]Peter J.Radiometric Calibration Transfer Chain from Primary Standards to the End to—End Hyperion Sensor[J]. SPTE.2000.(4l35):254—263.