由传递函数分析光电成像系统的分辨率_图文

1、第竺纛景1:期激光与光电子学进展2003年11月”。、V014咀Nnll Non 20032000.Mdch.GemaII.2000.1:119一12024MikkelserI B.Auopacal wavelengmcoHve岫scheI岫for m曲一speedRz涮Fc帅ats.皿打咖.船o,1997,38:2137213925W札fson n AH一0p廿cal 2RRegene功妇at40Gbi妇jIlanSoAb粕edMa曲一姗der InterfemmeteL OF1999.S趾Diego,cA.1999.PD0626k11恤0ld J.100Gbh/s Auopdcalwave

2、leng山cor帆r8ionw札hanhlte鲫地d s(ADelayedinterferencecoIl丘gu姻虹0n.郇bcalAmp】访er龃dApp吐cadon CorIference2000.Qu曲ec,锄ada_2000.OWB327Leu帅ld丑Novel3RRegener砒0r Basedonse咖c蚰ductor Op赴al AIllp】蚯er De蛔ed一血teIferencecorI丘gLlra舡on脚鼢细knBL厶ff,200L AI璃86086228nyshevPv.Allop廿calD北IRegene枷0rIB丛edse】f_phaseModIlla廿0IlE地cL

3、ECOCl99&Madddsp血L199&475476Realisa廿on舡Id Development of l刀缸丑一Speed optical 1阻me DivisionMllltiphe】iIIgzHU珏enC衄NarI舭lGU0F、lsheIIgnAI GuoHangWANGH0ng咖0.面069,D1,伽删471000Abs虹a肚To鼬tis如the de咖dformore位瑚sfni鼹ioncapaci睁,bom忙n咖ber 0f w删e砸:tlls(wDMarIdmernemIll一plexed b址ra土ecIDhDperwDMcha蚰elareirm勰edIn鸿p印ewep

4、Msent仕Ie scateof_art w】_lich enable me TDM bn ra把of160Gbi以by mea地of 0p廿cal廿】memvisiOnm山ple】drtg,w曲舭鹏on seIrIicondu咖rdevic船.Key wordsop廿cal胁edivislon嘶lexmgde删埘岫ledIIgdock recovery由传递函数分析光电成像系统的分辨率达争尚陈良益(中国科学院西安光机所水下光电探测技术室,西安710068一IN 2月提要光电成像系统的空间分辨是表征系统性能的一个最主要指标,通过分析光电成像系统中的传递函数,用特征频率法和传递函数对频率的积分方

5、法分析了系统的空间分辨,并针对一组典型的光电成像系统参数进行了计算,结果为:系统在Nyqujst频率处显示图像的wrF=0.017;传递函数对频率税分法计算表明:系统实际的分辨率为18.6cyck妇吼。分析过程表明:显示器是限制实时成像系统空间分辨的主要因素,其次是探测器的像元尺寸,苒次是电子系统的带宽。同时表明,分辨率只有在蕊统中讨论才有意义,仅甩rl斌频率表示是欠妥的。关键词光电成像 空问分辨wrF 系统带宽岬jst频率通常.在表示光电成像系统的分辨率时.总是用光电器件的像素表示。例如,对数码相机,资料显示可以达到400万的像素数,假设光电器件的宽高比为4:3,光学格雹旺巫塑雯囹式为1/2

6、英寸(6.4mm4.8Inm,则其像素的大小应为2.8叫lz,就成像器件目前的发展水平讲.这是不可能的。实际上.400万是个电子数据概念,是经过软件插值以后得到的,并不表示真实像素的多少:光电采样成像系统分辨率的另一种表示方法是用探测器的Nyquist频率表示m目,而实际系统能否达到这一分辨也是一个需要讨论的问题。作为一个复杂的成像系统.从输入的光学信息到人眼在cRT上观察第40卷.第ll期20年11月到最终的图像,其间包含有众多的信息传递环节,仅仅用这一数值表示其分辨是欠妥的。对比于传统的胶片摄影系统,用户关心的是系统的物理分辨,亦即系统单位长度所能分辨的空间线对数空间分辨率,这就是本文所要

7、阐述的主要内容。本文从传递函数的角度,分别用日1特征频率法(取为Nyq陋st 频率;2传递函数的平方对空间频率的积分方法f类似于能量集中度指标对光电系统的空间分辨予以分析,得出限制系统空间分辨的主要因素。为系统设计提供一定的参考,同时也与胶片式摄影系统作以比较。首先分析光电成像系统中的传递函数。图1为光电成像系统的组成环节。整个系统包括以下几个环节:1光学系统,2探测器,3信号预处理.4量化,5图像处理,6图像重建,7显示。在信息传输过程中有两次图像抽样:一是探测器对输人光学图像的二维空间抽样,=是在图像量化过程中AD转换器对图像的一维时间抽样。以上的每一个环节和过程都对信息的传递产激光与光电

8、子学进展生影响,有各自的传递函数,对照图l,下面对各个环节的传递函数(进行分析,为便于表述,只写出一维方向的MTF表达式。2.1光学系统MIF在可见光波段,对像差矫正较好的系统可以近似认为是衍射极限系统。圆形孔径的光学系统其衍射斑强度的分布为一阶贝塞尔函数,调制传递函数为MTF日00=导c酬(专卜鲁、/一等】,当时=O,(1当其他情况时,%。为截止频率;。=D棚,地是像方空间分辨。Do是光瞳直径,是焦距,A是波长,通常取o.55岬。2.2探测器肼宽度为dH的探测器姗F(调制传递函数为H蚴o=掣|=l曲M讲0l(2上式是一个典型的siIlc函数,第一个零点地=l肠:所对应的频率即为探测器本身的截

9、止频率。而探测器阵列的IF为:图l光电成像系统组成图V01.40llNov.2003lsm(mol(3表示相邻像元的中心距,阵列的Nyqllist频率为仳州=l/2吐唧。当dH_d。时,(2式和(3式相等。采样时,当探测器相对于信号的相位改变时,输出的幅值改变,从而使采样的调制度是一个变化的值,如图2所示。在任一位置采样阵列与图像相对相位对探测器MrF的影响可表示为nE,。.”MT=。傩(专8(4自是信号和探测器阵列的相对相位角,在地=铭w时,若e=o,lnFlI=l【图2(a】;口_90。,m日=o【图2】;平均的日=矾扫肛E一=0.707,这一数值是计算电视分辨率时常用的一个系数。而相位对

10、探测器mF的平均影响可表示为“M仉一4shl。(嚣J(5考虑到相位影响后,探测器阵列的MTF为:MTFa。0=fs妇of f咖c(等l(6实际上影响探测器rF的因素还有电荷转移效率以及电荷包形成时电子在相邻势阱边沿扩散的影墨基嗅蔷蓦信导一际,噫L卜m囊 忏婴二第40卷,第1l期2003年11月激光与光电子学进展茸目l目目.万日*探测器输出n(a口口口口口口广1(1图2相位对探测器MTF的影响响。但是只有当像素数很多时电荷转移效率对rrF的影响才明显起来,而电子扩散的影响只有在红外波长以后才明显.在可见光范围可不计。2.3信号预处理模拟滤波器MTF滤波器处理的电子信号.信号的频率正与空间图像的频

11、率地有内在的对应关系忙旦娶盟羔她(7BHmBPOV为视场,b,。为读出一行信号的时间,与视频带宽有关。在Nyq陋st频率时,上对应于视频带宽。理想的滤波器是不存在的,实际滤波器的模型可以用N一阶巴特奥斯ut塘rworm滤波器近似,传递函数为16lH一=/、/,+(点“(8 ,棚是当信号幅值衰减为输人的O.707倍时的频率.当_+。时.蜀一(国接近理想的滤波器,截止频率为J二。,对应于(7式中探测器的Nyqllist频率。2.4采样一保持器的A册图像量化过程中的采样一保持器产生的MIF为Hm昧s疵争(9当刚好满足采样定理时。名对应于探测器的空间采样频率。2.5数字图像处理环节的MTF 墨数字图像

12、处理的算法很多.但是只有少数算法的传递函数可以定量描述,这些传递函数可以包括在系统模型中,其他算法的效果只可通过系统对典型输入的实验结果得到。在图像处理算法中,最基本的是数字滤波.而数字滤波算法又可细分为好多种_。算法作用的对象是存储器中的图像数据,当将像素和这些数据一一对应时,滤波器的时间采样频率对应于探测器的空间采样频率,这样滤波器实际是对每一个像素进行处理。一种常用的加权数字平均滤波算法的MTF为“册。I刮%舢f学1fMTF。(D=lA心osf兰警吐1k-0、J(10为平均数值个数,A。为加权系数。当正=O时,MrFm:,联=0=1,所以有乞A产l。2.6图像重建滤波器的脚nA转换后的重

13、建信号轮廓显现为台阶状,为消除轮廓状(即重建形成的高频成分,需用滤波器滤除重建信号中的高频成分,恢复原始图像信号的频谱,同样用N一阶巴特奥滤波器近似,传递函数为”如一d/、/+(丢“(112.7显示器的MTT假设显示器上的亮点强度为vol蚰.11N唧.2003高斯分布,并且水平和垂直方向的分布一致。其MTF为HM哪d=e冲f_2晌粕-e印【一2斧(者叫2J(12其中。赢是高斯分布的标准差,S为高斯分布相对强度为一半时的亮点宽度,S=2.3勋。锄是CRT显示的空间频率。它和像方的频率乱。有如下关系【同(7式1”甓孚铒F罨警%I!氍堕l f13【I【n Juv、vFuv分别为像方水平和垂直方向的视

14、场,W。、E。为显示器的宽度和高度。工业标准要求髓。m一2.3勋0.7x每1.645M,为扫描线数,和探测器垂直方向的像素数对应.由此得册F岛为MTFc=e印(一2稍瑶础:一唧掰(缶%2】(t4基于姗F的分辨评价方法很多H,介绍本文拟采用的两种方法:1特征频率法,根据系统的应用要求确定特征频率。对应于特征频率的姗F不低于某一数值(例如O.1。对光电采样系统,人们通常用Nyqujst频率表征系统的分辨率,这里取特征频率为Nyquist频率进行分析。2传递函数平方对空间频率的积分法,亦即shade的等儿分辨玩。方法,sImdeI司教授将图像的锐度与下式联系起来f忙(【mF(%圳d毗(15O(15式

15、的含义实际是对各个空间频 率用其MTF的平方加权求和,M第帅卷,第ll期20年11月即表示了系统的空间分辨率。e越大,图像质量越好。uoyd将shade 的公式变通为,/i扁;=寿=1/2J【MTF(蚴】2(蚍一叩f i(16 jo嘶表示分辨力,岛驯、,则分辨率大,(15式和(16式是一个问题的两种表示。风。并不能直接测量得出,而是一个计算的数值,由于各个分系统是彼此独立互不相关,整个系统的jk。可以由各个分系统的等价兄。按下式求出局辩湃=V倚品一。+鼬_2+佃南吨(17 &mde用MTF的平方来强调册F相对大的空间频率成分.该方法对于评价肼逐级递减系统的性能来说是一个好方法,并且假设整个系统

16、是模拟的,忽略采样的影响。作为一个整体性能指标,品。较之其他的分辨方法更能表示系统的分辨性能。在计算R。4。时,若某一子系统的酣rF在要求频段上为1,则其R。fO。将前面讨论的各分系统的M耶表达式用(16式积分计算其置。,结果如表1所示。光电成像系统在宏观上可以认为是线性不变系统日,系统的MrFk激光与光电子学进展可以由各个分系统的mF相乘得到MTFk=MTF:州M晒l。MTFI llTrF矗函ITFE曲MTFR IMTT;目(18分析光电成像系统的空间分辨时.为使问题具体化.作以下假设:光学系统焦距.厂=100mm,/ Do=5.6,=O.55岬,探测器的幽= d。=d。=lO畔l,低通滤波

17、器是理想的(一*1.数字滤波器的MTFj。=1,并假设系统的带宽、采样满足采样定理的要求。1NyqIIist频率时的MkNyqllist频率为乱i=l/2d吲= 50呵des,mm代人以上各个分系统M7rF表达式有BFO神=O.805;雌阿。0曲=0.405;盯IE。0=O,707,Mkn啪=O.637;MI艮0=1; MTFR山曲=O.707;MTF(m(钆曲=0.16l(19则MIm曲=0.017;这一调制度还达不到人眼能够分辨的调制度阈值要求.即系统的分辨实际上达不到Nyqllist频率。2shade等效分辨率参照表l,忽略采样的影响,仅考虑模拟环节的子系统船%1.845护光学系统【os

18、、7虬v、/J2风探测器sincldlli d一一* N阶低通滤波器、四L当_l一CRT e如(一2枇32订噼2订(蛊V01.40.1lN叫2003R暾dk=1.845妒,DFO.006,置%“=0.0lO,局dr=O.01仉兄*玎=0.022,(20由(17式得:甬目!f、厢面霭骊而面飘丽而霜而万= 0.027(21由(16式、(16式得系统的分辨率为:他2赤2186cycl8如“,此即为显示图像的空间分辨率。若01式中不计五。,则“电子图像”的分辨率为。=、/丽丽颓瓜面面丽矿=0.015(22设想用高分辨率的打印机(MTF1将图像打印出来,则其分辨率为:赤=33cydes, mm,与胶片摄

19、影系统的分辨率相当。但不幸的是,正是由于打印设备性能欠佳.才使数码照片质量与摄影照片质量有一定的差距。需要说明的是在用传递函数建模对光电系统分辨率分析时,没有考虑电学系统噪声对图像的影响。实际上人眼在观察有噪声的图像时.所要求的调制度阚值将更大一些日.相应的系统的分辨率减小。而由于建模时的理想化处理,实际系统分辨率应较理论计算还要小些。通过以上对一组光电成像系统典型参数的分析计算,得出结论如下:1(19式、(21式表明,对实时成像系统而言,显示器是限制系统MTF和空间分辨的主要因素,其次是探测器像元尺寸.再次是电子系统的带宽,而光学系统对整体分辨影响相对来说要小一些,但这是一般情况。当探测器的

20、尺寸相对较大 墨第40卷,第ll期2003年11月时.它对系统分辨率的影响更多一些.例如红外成像系统;当系统的信息速率高时.电子带宽对分辨的影响可能更多,例如高速视频成像系统。2对Nyquist频率wrF的分激光与光电子学进展析说明.光电成像系统实际上很难达到这一分辨。作为一个复杂的信息传递系统。其间包含有光学信息向时间信息的空间一时问变换、图像的探测器空间采样、图像重建的时间采样以及显示等过程,图像传参考文献V014n.1lNdv.2008输的每过程都会对最终图像的分辨造成影响.因此必须从系统的角度讨论图像的分辨率,探测器只是系统图像信息传递的一个环节,仅用探测器的Nyqu龇频率表示光电采样

21、成像系统的分辨率是欠妥的。1杨桦.焦村文ccD相机在奈奎斯特频率处的传递函数.光学学报,2002,22(3:3133162郝彩线阵ccD相机细分采样成像的像质研究.光学学报.2000,20(10:140714113沙定国.高精度光学传递函数的计算.光学学报.18(4:4374414庄松林.钱振邦.光学传递函数.北京:机械工业出版社,1981_63695I蛐t G C.CCD Arrays,C枷eI鹤,锄d Di印la涔.Beuin曲am,w盘sPPr螂,2000.7389,1561856胡广书数字信号处理.北京:清华大学出版社,1997.3l,507董思白模拟电子技术北京:清华大学出版社.19

22、92.2372508陈自宽,母国光.光电成像系统中的光学传递函数.光学技术,1998,l(1:5763A啦rsis of Spa廿al Resolution of optoElectroIIic h岫ge System111ongh tlIe Mbdula_don叻瓢nsfer F珊c廿佃Da Zl他rIgs叫曜Chen UalIgyi(巩蝴伽椰。,咖翻蝴幻一脚咖耙脑咖脚池州佩倒姗胁玩据qr却妇n州A诚t帆埘醣概tcso,oiS皿船710068Absh喊皿mdugh血e MTF ana】ysjs of0p协一ele咖nic irnage町砒em,t【Ie印们al r郫0lu廿on js cal

23、anated based叩me f0(月血喀t啊o MrF metllods:mecharac把枷c freqI】即cy of佗印咖m arId仕Ie iIlte鲫n0f姗F叽t|Ie d咒rall sp删自刈岫n呵.The calcl】la抽n re卯lt of a typical s粥础nis:the册is0.017at仕le N硼【llist nquency,me systems ac机al resohlnon isl86cycles眦.The卸a】ysis shows CRTis血e丘rstfactor仕m de鲫dest|le sy8tem resolu虹0n;山e ne啦isme

24、pixelpnch0f 血e detect0L t|l朗廿Ie circu址b皿d而mh.Also,me author eonside瑙me resolll廿on should be卸alyzed h廿le sy8teI珥n is not p彤per to璐e Nyqllist盘equency0rdy t08Iow t王Ie resolu廿0n.Key words opt0一electmIlic ir【Iage syst哪spa廿al Iesol曲叽MTF systemb龃d-wm NW删龇行equen呵盈*:l:木:l:木:|:I:术:l:l:木水:l:木木术;l:l:水木;l=:l=术:l

25、:术术枣木;l:木:l:l:木木术术术术* *羹月刊,代号4179,每期定价20元羹*:l:|:l:木:I:术木:l:木木,Ic:l:l:木木木;l:木:l:l:术:l:木木水:术:I:l:l:I:I:水木半术木*量享年娘雌进O眨矿7加学阅子订电局艺隋光驯崩蒙光激由传递函数分析光电成像系统的分辨率 作者： 作者单位： 刊名： 英文刊名： 年，卷(期： 被引用次数： 达争尚， 陈良益 中国科学院西安光机所水下光电探测技术室,西安,710068 激光与光电子学进展 LASER & OPTOELECTRONICS PROGRESS 2003，40(11 0次 参考文献(8条 1.杨桦.焦村文 CCD

26、相机在系统奈奎斯特频率处的调制传递函数期刊论文-光学学报 2002(03 2.郝彩 线阵CCD相机细分采样成像的像质研究期刊论文-光学学报 2000(10 3.沙定国 高精度光学传递函数的计算 4.庄松林.钱振邦 光学传递函数 1981 5.HOLST G C CCD Arrays, Cameras, and Displays 2000 6.胡广书 数字信号处理 1997 7.董思白 模拟电子技术 1992 8.陈自宽.母国光 光电成像系统中的光学传递函数 1998(01 相似文献(10条 1.期刊论文 朱华征.范大鹏.马东玺.张文博.ZHU Hua-zheng.FAN Da-peng.MA

27、Dong-xi.ZHANG Wen-bo 载体运动对 光电成像系统性能的影响分析 -红外技术2008,30(10 首先根据光电成像系统与目标的几何关系,分析了积分时间内载体运动在焦平面上产生的像移;然后从调制传递函数的角度分析了运动对成像系统性 能的影响,最后分析了运动对光电成像系统分辨率的影响.结果表明:载体的角运动会严重影响光电成像系统性能,必须将积分时间内载体的角运动控制在 探测器瞬时视场内. 2.会议论文 金伟其.王霞.苏秉华.刘扬阳.张楠.徐超 亚像元与超分辨力成像技术及其应用进展 2005 目前,光电成像系统分辨力的提高主要受体积、重量以及探测器水平等因素的限制,结合现代图像处理与

28、光学成像技术已成为国际上突破以上限制的 主要技术途径.亚像元成像技术通过光学微扫描或探测器错位拼接等方法以及后续信号处理,获得小于探测器几何尺寸的空间分辨;超分辨力成像则利用图 像的先验知识,提高图像信息的高频信息.亚像元处理和超分辨力图像处理技术已被证明是在探测器和光学系统限制下提高成像系统分辨力的有效方法,在 国外的一些遥感和军用成像系统中获得成功的应用.本文介绍了国内外亚像元和超分辨力成像技术及其进展,并对发展对地观测系统提出了一些建议. 3.期刊论文 朱华征.范大鹏.马东玺.张文博.ZHU Hua-zheng.FAN Da-peng.MA Dong-xi.ZHANG Wen-bo 动载

29、体光电 成像系统视轴稳定精度研究 -应用光学2009,30(4 载体扰动会对光电成像系统的性能产生影响.分析载体扰动的特点及产生的原因,主要分析载体角运动对光电成像系统性能的影响.在此基础上从空间 分辨率、人眼观察性能以及跟踪精度3个方面分析光电成像系统对视轴稳定精度的要求,指出对于探测器限制的光电成像系统,应将积分时间内载体的角扰 动控制在探测器限制的空间分辨率内,或将像移控制在一个探测单元尺寸内,对于衍射限制的光电成像系统,应将积分时间内载体的角扰动控制在衍射限制 的空间分辨率内,得出了视轴稳定精度必须在满足成像系统空间分辨率的基础上,根据实际功能需求进行确定的结论.最后通过实例验证了上述

30、结果. 4.期刊论文 薛玲玲.陈波.石晓光.薛国俊 30.4 nm光电成像系统分辨率的初步实验研究 -量子电子学报2003,20(4 本文基于微通道板MCP(Microchannel Plate探测器件设计一套成像系统,用于对波长为30.4 nm的极紫外EUV(Extreme Ultraviolet光进行成像.结 果获得了一宽度为3 mm的狭缝的像,实验测得490m的成像系统的空间分辨率,并分析了影响系统分辨率的各种因素及为提高系统分辨率所应采取的措施. 5.期刊论文 薛玲玲.陈波 19.5 nm波段光电成像系统分辨率的实验研究 -光电子激光2003,14(1 基于微通道板(MCP探测器件设计

31、一套成像系统,用于对波长为19.5 nm的极紫外(EUV光进行成像.结果获得了宽度为3 mm的狭缝的像,实验测得175 m的成像系统的空间分辨率.并分析了影响系统分辨率的因素及为提高系统分辨率所应采取的措施. 6.期刊论文 张斌.李朝晖.ZHAN Bin.LI Zhaohui 光电成像系统的分辨率鉴定与测量技术 -现代电子技术 2010,33(1 论述了光电成像系统中广泛使用的分辨率指标及分类,对空间分辨率模拟度量法的原理和测量方法进行了论述和分析.通过研究指出用空间分辨率指 标来描述成像系统的质量,具有较好的直观性和归一性.由于单一的空间分辨率测量指标还不可能给出总的图像系统的性能,仅仅基于

32、分辨率指标的图像评 估不可能同时保证系统灵敏度设计的技术要求.因此,结合模拟度量法研究光电成像系统的分辨率测量法,给出成像分辨率测量准则. 7.期刊论文 盛亮.黑东炜.邱孟通.魏福利.王培伟.王奎禄.SHENG Liang.HEI Dong-wei.QIU Meng-tong.WANG Kuilu.WANG Pei-wei.WEI Fu-li 纳秒级时间选通成像系统动态空间分辨率 -核电子学与探测技术2008,28(2 成像系统动态空间分辨率的下降与静态相比主要发生在光阴极与微通道板输入面之间.本文采用时空传递函数理论,分析了动态下空间分辨率的下降 因子,计算得到的结果与实验结果相比基本一致,