CN104977084B - 一种提高aotf成像空间分辨率和光谱分辨率的方法 - Google Patents
一种提高aotf成像空间分辨率和光谱分辨率的方法 Download PDFInfo
- Publication number
- CN104977084B CN104977084B CN201510391690.2A CN201510391690A CN104977084B CN 104977084 B CN104977084 B CN 104977084B CN 201510391690 A CN201510391690 A CN 201510391690A CN 104977084 B CN104977084 B CN 104977084B
- Authority
- CN
- China
- Prior art keywords
- aotf
- diffraction
- lambda
- light
- resolution
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Spectrometry And Color Measurement (AREA)
Abstract
本发明涉及AOTF光谱成像技术领域,涉及一种提高AOTF成像空间分辨率和光谱分辨率的方法;提出一种采用多次迭代去逼近真值的优化算法,在某波长下,将CCD成像的每个象元得到的光强减去由于光谱和衍射角展宽的干扰光,从而获得光强真值,该真值是经过多次迭代获得,因此方法将AOTF获得的二维图像和一维光谱的数据立方经过多次迭代来提高最终的成像空间分辨率和光谱分辨率;首先对AOTF衍射角展宽进行测量,然后对AOTF衍射光谱展宽进行测量,经过上述两者的测量得出衍射效率η(xi+m,yj,λk,mdλ);通过AOTF光谱成像获得被测目标在CCD的象元(xi,yj)探测的光强I0(xi,yj,λk);通过公式进行多次迭代;本发明主要应用在AOTF光谱成像方面。
Description
技术领域
本发明涉及AOTF光谱成像技术领域,更具体而言,涉及一种提高AOTF成像空间分辨率和光谱分辨率的方法,是一种根据AOTF衍射特点采用多次迭代优化算法提高AOTF空间分辨率和光谱分辨率的方法。
背景技术
声光可调谐滤光器(Acousto-optictunable filter,AOTF)与传统分光元件相比,具有体积小、调谐速度快、光谱范围宽、衍射效率高等优点。近年来,基于AOTF的成像光谱技术越来越受到研究者的关注。基于AOTF的成像光谱仪已经在遥感、环境监测、生物医学和食品检测中得到应用。而AOTF成像光谱***的空间分辨率和光谱分辨率的高低直接影响着AOTF光谱成像性能的优劣。
现有的AOTF成像光谱***中为消除图像的漂移现象,在AOTF出射面多采用光楔和棱镜补偿的方式消除漂移。但是此方法只考虑了AOTF滤光后主峰中心波长的色散带来的图像偏移现象,没有考虑AOTF滤光后光谱存在旁瓣和一定的展宽,也就是AOTF衍射中还存在衍射角和波长的展宽,这样将导致AOTF在衍射方向的成像空间分辨率较低,并且也导致光谱分辨率下降。
发明内容
针对现有AOTF光谱成像测量中由于衍射后存在光谱展宽和衍射角展宽导致衍射方向成像空间分辨率下降及光谱分辨率下降,根据衍射后光谱展宽和衍射角展宽一一对应关系,提出一种采用多次迭代去逼近真值的优化算法,在某波长下,将CCD成像的每个象元得到的光强减去由于光谱和衍射角展宽的干扰光,从而获得光强真值,该真值是经过多次迭代获得,因此方法将AOTF获得的二维图像和一维光谱的数据立方经过多次迭代来提高最终的成像空间分辨率和光谱分辨率。
为了解决上述技术问题,本发明采用的技术方案为:
一种提高AOTF成像空间分辨率和光谱分辨率的方法,包括下述步骤:
a、首先对AOTF衍射角展宽进行测量,然后对AOTF衍射光谱展宽进行测量,经过上述两者的测量得出衍射效率η(xi+m,yj,λk,mdλ);
b、通过AOTF光谱成像获得被测目标在CCD的象元(xi,yj)探测的光强I0(xi,yj,λk);
c、通过公式进行多次迭代,所述公式为:
上式中,In(xi,yj,λk)为第n次迭代时,(xi,yj)处波长为λk时的光强;I0(xi,yj,λk)为上述b步骤中探测的光强;η(xi+m,yj,λk,mdλ)为峰值波长为λk时,目标(xi+m,yj)经过AOTF后波长为λk+mdλ的衍射效率,由上述a步骤中测量获得;In-1(xi+m,yj,λk+mdλ)为第n-1次迭代时,(xi+m,yj)处波长为λk+mdλ时的光强;η(xi,yj,λk,0)为峰值波长为λk时,目标(xi,yj)经过AOTF后波长为λk的衍射效率;
迭代次数停止条件为:
上式中,σ为迭代中前后两次各象元探测强度的平均变化率;
(xN2~xN1)(yN2~yN1)(λN2~λN1)为数据立方块;
其中,数据立方块(xN2~xN1)(yN2~yN1)(λN2~λN1)选择有明显明暗变化的区域,并且(xN2~xN1)(yN2~yN1)(λN2~λN1)和σ应根据实际中测量精度和运算时间的要求进行选择。
所述a步骤中对AOTF衍射角展宽进行测量具体为:将白光光源、分光仪、CCD相机依次排列,所述分光仪中依次排列有带有狭缝的毛玻璃、平行光管、AOTF和望远镜,所述分光仪中还设置有刻度盘,白光光源照射在带有狭缝的毛玻璃上,由平行光管将狭缝的光变为平行光,并进入被测AOTF中,AOTF衍射后的光再由望远镜将其成像于CCD相机上。
所述刻度盘的最小角度测量精度为1'=(1/60)°,所述CCD相机的分辨率为640×480Pix,角度测量***的角度测量精度可以达到0.0015°。
所述a步骤中对AOTF衍射光谱展宽进行测量是采用高光谱分辨率的光谱仪测量。
所述b步骤中将前置光学***、AOTF、后置光学***和AOTF光谱成像CCD相机依次排列进行测量。
与现有技术相比本发明所具有的有益效果为:
该方法通过多次迭代优化处理AOTF成像光谱数据立方,提高AOTF成像的空间分辨率和光谱分辨率;
该方法采用分光仪和CCD结合的角度测量***来测量AOTF衍射角展宽,角度测量精度可达0.0015°,并且采用高光谱分辨率的光谱仪测得AOTF衍射光谱测展宽,获得AOTF衍射角展宽和衍射光谱展宽的对应关系;
该方法测量AOTF衍射角展宽和衍射光谱展宽时,对应一个AOTF只需一次完整的测量后,在后续使用中只需调用即可,无需重复测量;
该方法也无需改***件设备,只需对后续的数据进行优化处理即可。
附图说明
下面通过附图对本发明的具体实施方式作进一步详细的说明。
图1为本发明的流程图;
图2为AOTF衍射角展宽测量装置结构图;
图3为本发明中AOTF光谱成像结构图;
图4为CCD象元探测实际光强示意图;
图5为多次迭代AOTF光谱成像数据立方示意图;
图6a为实施例中水平狭缝CCD成像图;
图6b为实施例中竖直狭缝CCD成像图;
图6c为实施例中水平狭缝与竖直狭缝成像带宽对比图;
图7a为当超声驱动频率fa=110MHz时,衍射角展宽测量结果图;
图7b为当超声驱动频率fa=110MHz时,衍射光谱展宽测量结果图;
图8a为n=0时,竖直狭缝的CCD成像图;
图8b为n=5时,竖直狭缝的CCD成像图;
图8c为水平狭缝的CCD成像图;
图8d为图8a、图8b和图8c的狭缝成像带宽对比图。
图中:1为白光光源、2为分光仪、3为CCD相机、4为带有狭缝的毛玻璃、5为平行光管、6为AOTF、7为望远镜、8为刻度盘、9为前置光学***、10为后置成像光学***、11为AOTF光谱成像CCD相机。
具体实施方式
下面实施例结合附图对本发明作进一步的描述。
如图1所示,对于在AOTF每一个中心峰值波长滤光下,首先对AOTF衍射角展宽测量S101,并对AOTF衍射光谱展宽测量S102,获得AOTF衍射角展宽与衍射光谱展宽的对应关系,为后续的数据处理提供必要的参数;将AOTF光谱成像S103CCD获得的数据通过计算机多次迭代优化处理S104,最终实现高AOTF成像空间分辨率和光谱分辨率的目的。
由于AOTF中的超声换能器为矩形结构,这样将导致超声进入声光晶体时存在声波衍射,因此导致超声波失方向存在一定展宽,并且不同方向上的超声声压不同。根据AOTF动量匹配条件,超声波失方向的展宽将导致衍射光谱的展宽,并且超声波声压的不同将导致衍射光谱展宽后不同光的衍射效率不同。
本实施方式中的AOTF出射面已采用光楔方式消除了图像漂移。
如图2所示,AOTF衍射角展宽测量S101包括白光光源1、分光仪2和CCD相机3,白光光源1照射在分光仪2前端带有狭缝的毛玻璃4上,由平行光管5将狭缝的光变为平行光,并进入被测AOTF6中,AOTF6衍射后的光再由望远镜7将其成像于CCD相机3上。分光仪2的刻度盘8最小的角度测量精度为1'=(1/60)°,CCD相机3的分辨率是640×480Pix,CCD相机3固定在望远镜7的后端。当转动望远镜7使狭缝像的中心极大位置从CCD相机3的最左侧移动到最右侧,对应分光仪2转过的刻度为58',因此该角度测量***的角度测量精度可以达到58'/480≈0.0015°,该测量测得的是AOTF6外的衍射角展宽,并不是AOTF6晶体内的衍射角展宽。当然,如果CCD相机3中采用象元尺寸更小的CCD测量,角度测量精度会更高。
AOTF衍射光谱展宽测量S102采用高光谱分辨率的光谱仪测得。本实施方式中采用光谱分辨率为0.25nm@λ=600nm的光纤光谱仪。
AOTF衍射角展宽测量S101和AOTF衍射光谱展宽测量S102中是以波长为0.1nm间隔设置AOTF扫频间隔,测量光谱范围为450nm-700nm。
CCD相机3获得同一宽度的水平狭缝和竖直狭缝的像对比如图6a和图6b,提取图6a的第300列和图6b的第250行对比狭缝成像展宽,水平狭缝与竖直狭缝成像宽带比较如图6c所示,由图6c可以看出,与水平狭缝成像相比,竖直狭缝成像不仅存在旁瓣,而且主极大的宽带也较水平狭缝宽,也就是在AOTF衍射方向存在展宽,在AOTF超声驱动频率fa=120MHz时,竖直狭缝成像主极大的半峰宽是水平狭缝成像的3.5倍。
AOTF衍射角展宽测量S101与AOTF衍射光谱展宽测量S102对比结果如图7a和图7b所示。当超声驱动频率fa=110MHz时,衍射角展宽测量结果如图7a所示,衍射光谱展宽测量结果如图7b所示,由图7a和图7b可知,衍射角展宽与衍射光谱展宽基本吻合,包括特征光谱也能很好吻合,如图7a和图7b椭圆区域。
AOTF衍射角展宽测量S101和AOTF衍射光谱展宽测量S102,主要是为获得每个AOTF滤光下AOTF衍射角展宽与衍射光谱展宽的对应关系,为后续的多次迭代提供必要的参数。虽然AOTF衍射角展宽测量S101和AOTF衍射光谱展宽测量S102比较复杂,数据量大,但是对于一个AOTF6来说,只需要一次全面的测量,在后续的光谱成像中只需直接应用这些参数即可。
如图3所示,AOTF光谱成像S103包括前置光学***9、AOTF6、后置成像光学***10和成像的AOTF光谱成像CCD相机11。为了便于比较,AOTF光谱成像S103所有的光学***采用AOTF衍射角展宽测量S101的光学***相同,并且波长间隔为0.1nm,测量光谱范围为550nm-700nm,但是采用与AOTF衍射角展宽测量S101不同的白光光源。
如图4所示,当AOTF6设置的滤光波长为λk时,如果AOTF6是理想光谱成像,应该是AOTF光谱成像CCD相机11的象元与空间目标是一一对应关系,即CCD的象元(xi,yj)探测光强为空间目标(xi,yj)处波长为λk的光强;但是由于AOTF6衍射角和衍射光存在展宽,因此AOTF6设置的滤光波长为λk时,CCD的象元(xi,yj)探测到的光强I0(xi,yj,λk)为:
其中,I(xi,yj,λk)是目标(xi,yj)处波长为λk时的真实光强;η(xi+m,yj,λk,mdλ)是峰值波长为λk时,目标(xi+m,yj)经过AOTF6后波长为λk+mdλ的衍射效率。
因此,目标(xi,yj)处的真实的强度I(xi,yj,λk)为:
其中,I0(xi,yj,λk)可由AOTF光谱成像CCD相机11获得,η(xi+m,yj,λk,mdλ)可通过AOTF衍射角展宽测量S101和AOTF衍射光谱展宽测量S102获得,由(2)式可知,要得到真实的强度I(xi,yj,λk)还需要知道I(xi+m,yj,λk+mdλ),而I(xi+m,yj,λk+mdλ)需要多次迭代去逼近。
如图5所示,具体推导如下:
其中,In(xi,yj,λk)为第n次迭代时,(xi,yj)处波长为λk时的光强;I(xi+m,yj,λk+mdλ)采用前一次的In-1(xi+m,yj,λk+mdλ)代替,通过多次迭代去逼近。其中m和dλ的取值与AOTF光谱成像中的光谱间隔(λk+1-λk)、衍射角展宽和后置成像光学***10的焦距f有关。迭代次数停止条件根据实际要求合理选择,如下式:
其中(4)式等号右边表示迭代中前后两次各象元探测强度的平均变化率,其值会随着迭代次数的增加越来越小,因此σ选择越小,计算结果越精确,但运算也就越大;其中(xN2~xN1)*(yN2~yN1)*(λN2~λN1)为选择迭代次数停止条件的数据立方块,其越大越精确,但运算量越大;在选择(xN2~xN1)*(yN2~yN1)时应选择有明显变化的区域。因此,(xN2~xN1)*(yN2~yN1)*(λN2~λN1)和σ应根据实际要求合理选择。本实施方式中迭代停止条件选为:xN1=250,xN2=390,yN1=220,xN2=260,λN1=550nm,λN2=620nm,σ=0.00001。最终的实验结果对比如图8a、图8b、图8c和图8d所示,由图8a、图8b、图8c和图8d可看出,多次迭代以后竖直狭缝像与水平狭缝像的带宽基本相同,这将提高AOTF成像空间分辨率。根据(3)式可知,由于该多次迭代优化算法减去了衍射光谱展宽的影响,因此AOTF的光谱分辨率也得到了提高。
Claims (5)
1.一种提高AOTF成像空间分辨率和光谱分辨率的方法,其特征在于,包括下述步骤:
a、首先对AOTF衍射角展宽进行测量,然后对AOTF衍射光谱展宽进行测量,经过上述两者的测量得出衍射效率η(xi+m,yj,λk,mdλ);
b、通过AOTF光谱成像获得被测目标在CCD的象元(xi,yj)探测的光强I0(xi,yj,λk);
c、通过公式进行多次迭代,所述公式为:
上式中,k表示AOTF光谱成像数据立方中对应的第k张图的标识,λk为AOTF滤光第k次所对应的峰值波长,光谱成像数据立方中衍射In(xi,yj,λk)为第n次迭代时,(xi,yj)处波长为λk时的光强;I0(xi,yj,λk)为上述b步骤中探测的光强;η(xi+m,yj,λk,mdλ)为峰值波长为λk时,目标(xi+m,yj)经过AOTF后波长为λk+mdλ的衍射效率,由上述a步骤中测量获得;dλ为相邻光谱成像中CCD象元x方向相邻两象元多对应的波长间隔,m为(xi,yj)象元与(xi+m,yj)象元在x方向的象元间隔数,In-1(xi+m,yj,λk+mdλ)为第n-1次迭代时,(xi+m,yj)处波长为λk+mdλ时的光强;η(xi,yj,λk,0)为峰值波长为λk时,目标(xi,yj)经过AOTF后波长为λk的衍射效率;
迭代次数停止条件为:
上式中,σ为迭代中前后两次各象元探测强度的平均变化率;
(xN2~xN1)(yN2~yN1)(λN2~λN1)为数据立方块;
其中,数据立方块(xN2~xN1)(yN2~yN1)(λN2~λN1)选择有明显明暗变化的区域,N1为提取数据立方块中x、y和λ数据的下界,N2为提取数据立方块中x、y和λ数据的上界,并且(xN2~xN1)(yN2~yN1)(λN2~λN1)和σ应根据实际中测量精度和运算时间的要求进行选择。
2.根据权利要求1所述的一种提高AOTF成像空间分辨率和光谱分辨率的方法,其特征在于,所述a步骤中对AOTF衍射角展宽进行测量具体为:将白光光源(1)、分光仪(2)、CCD相机(3)依次排列,所述分光仪(2)中依次排列有带有狭缝的毛玻璃(4)、平行光管(5)、AOTF(6)和望远镜(7),所述分光仪(2)中还设置有刻度盘(8),白光光源(1)照射在带有狭缝的毛玻璃(4)上,由平行光管(5)将狭缝的光变为平行光,并进入被测AOTF(6)中,AOTF(6)衍射后的光再由望远镜(7)将其成像于CCD相机(3)上。
3.根据权利要求2所述的一种提高AOTF成像空间分辨率和光谱分辨率的方法,其特征在于:所述刻度盘(8)的最小角度测量精度为1'=(1/60)°,所述CCD相机(3)的分辨率为640×480Pix。
4.根据权利要求1所述的一种提高AOTF成像空间分辨率和光谱分辨率的方法,其特征在于,所述a步骤中对AOTF衍射光谱展宽进行测量是采用高光谱分辨率的光谱仪测量。
5.根据权利要求1所述的一种提高AOTF成像空间分辨率和光谱分辨率的方法,其特征在于,所述b步骤中将前置光学***(9)、AOTF(6)、后置光学***(10)和AOTF光谱成像CCD相机(11)依次排列进行测量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510391690.2A CN104977084B (zh) | 2015-07-06 | 2015-07-06 | 一种提高aotf成像空间分辨率和光谱分辨率的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510391690.2A CN104977084B (zh) | 2015-07-06 | 2015-07-06 | 一种提高aotf成像空间分辨率和光谱分辨率的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104977084A CN104977084A (zh) | 2015-10-14 |
CN104977084B true CN104977084B (zh) | 2017-03-29 |
Family
ID=54273813
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510391690.2A Expired - Fee Related CN104977084B (zh) | 2015-07-06 | 2015-07-06 | 一种提高aotf成像空间分辨率和光谱分辨率的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104977084B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105737983B (zh) * | 2016-03-30 | 2017-10-31 | 首都师范大学 | 一种超分重建的斜模高光谱成像地面试验装置及方法 |
CN106596064B (zh) * | 2016-12-13 | 2023-04-11 | 中国科学院西安光学精密机械研究所 | 同步扫描条纹相机动态空间分辨力的测量装置及方法 |
CN106840397B (zh) * | 2017-01-10 | 2018-11-16 | 长春理工大学 | 一种用于评价可见光烟幕遮蔽效应的测试方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201041514Y (zh) * | 2007-06-14 | 2008-03-26 | 武汉大学 | 便携式成像光谱仪 |
CN101907513A (zh) * | 2010-07-23 | 2010-12-08 | 中国科学院上海技术物理研究所 | 声光可调谐滤光器衍射性能弱光测试***和方法 |
CN203881515U (zh) * | 2014-03-28 | 2014-10-15 | 中国科学院上海技术物理研究所 | 自参考声光可调谐滤光器衍射性能测试装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008124446A1 (en) * | 2007-04-06 | 2008-10-16 | University Of Miami | Dynamic spectral imaging device with spectral zooming capability |
-
2015
- 2015-07-06 CN CN201510391690.2A patent/CN104977084B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201041514Y (zh) * | 2007-06-14 | 2008-03-26 | 武汉大学 | 便携式成像光谱仪 |
CN101907513A (zh) * | 2010-07-23 | 2010-12-08 | 中国科学院上海技术物理研究所 | 声光可调谐滤光器衍射性能弱光测试***和方法 |
CN203881515U (zh) * | 2014-03-28 | 2014-10-15 | 中国科学院上海技术物理研究所 | 自参考声光可调谐滤光器衍射性能测试装置 |
Non-Patent Citations (1)
Title |
---|
声光可调谐滤波器(AOTF)消色散设计;常凌颖等;《光子学报》;20091130;第38卷(第11期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN104977084A (zh) | 2015-10-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9519737B2 (en) | Filter modules for aperture-coded, multiplexed imaging systems | |
US6763149B2 (en) | Method and apparatus for correcting crosstalk and spatial resolution for multichannel imaging | |
US8023724B2 (en) | Apparatus and method of information extraction from electromagnetic energy based upon multi-characteristic spatial geometry processing | |
CN103592108A (zh) | Ccd芯片调制传递函数测试装置及方法 | |
CN104501955A (zh) | 一种可实现光谱超分辨率还原的光栅光谱仪 | |
CN104977084B (zh) | 一种提高aotf成像空间分辨率和光谱分辨率的方法 | |
CN104266755B (zh) | 一种提高光通量的光谱测量方法及*** | |
US8860942B1 (en) | Apparatus for multi-spectral imaging of point event detection | |
US20230204418A1 (en) | High-speed and high-precision spectral video system and method for flame shooting | |
CN105890529A (zh) | 测量细丝直径的方法及装置 | |
CN107655571B (zh) | 一种基于色散模糊的光谱成像***及其光谱重建方法 | |
CN104931138A (zh) | 一种采用棱镜提高aotf光谱成像质量的方法及装置 | |
CN1766531A (zh) | 激光光束质量m2因子实时检测仪 | |
US9194798B2 (en) | Imaging based refractometer for hyperspectral refractive index detection | |
CN104280120B (zh) | 一种光谱带宽测量方法和装置 | |
CN106022354A (zh) | 基于svm的图像mtf测量方法 | |
Kim et al. | Aperture-encoded snapshot hyperspectral imaging with a lensless camera | |
CN106595859B (zh) | 鬼成像方法和应用其的鬼成像装置 | |
CN110031100B (zh) | 一种多维度短波红外光谱成像检测装置 | |
CN116972969A (zh) | 一种基于二维色散光谱重构的光谱探测方法及装置 | |
CN104111243A (zh) | 一种荧光比率测量***及方法 | |
CN103792005A (zh) | 宽谱段二维色散高速成像光谱仪 | |
US20160021305A1 (en) | Method and apparatus for measuring optical systems and surfaces with optical ray metrology | |
CN107664533B (zh) | 一种三狭缝高光谱运动目标探测方法及装置 | |
CN105136295B (zh) | 一种aotf同一幅图中光谱不均匀的解决方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170329 Termination date: 20170706 |