CN106706132A - 一种用于海面太阳亮带内目标识别的红外探测装置及方法 - Google Patents
一种用于海面太阳亮带内目标识别的红外探测装置及方法 Download PDFInfo
- Publication number
- CN106706132A CN106706132A CN201611100992.0A CN201611100992A CN106706132A CN 106706132 A CN106706132 A CN 106706132A CN 201611100992 A CN201611100992 A CN 201611100992A CN 106706132 A CN106706132 A CN 106706132A
- Authority
- CN
- China
- Prior art keywords
- matrix
- middle infrared
- infrared camera
- infrared
- temperature control
- 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.)
- Granted
Links
- 241001662443 Phemeranthus parviflorus Species 0.000 title claims abstract description 42
- 238000000034 method Methods 0.000 title claims abstract description 18
- 238000001514 detection method Methods 0.000 claims abstract description 54
- 230000010287 polarization Effects 0.000 claims abstract description 48
- 229910052732 germanium Inorganic materials 0.000 claims abstract description 30
- GNPVGFCGXDBREM-UHFFFAOYSA-N germanium atom Chemical compound [Ge] GNPVGFCGXDBREM-UHFFFAOYSA-N 0.000 claims abstract description 29
- 238000007789 sealing Methods 0.000 claims abstract description 5
- 239000011159 matrix material Substances 0.000 claims description 170
- 230000005855 radiation Effects 0.000 claims description 84
- 230000010354 integration Effects 0.000 claims description 24
- 230000004044 response Effects 0.000 claims description 14
- 238000003384 imaging method Methods 0.000 claims description 12
- 235000012149 noodles Nutrition 0.000 claims description 12
- 230000004313 glare Effects 0.000 claims description 11
- 210000001747 pupil Anatomy 0.000 claims description 8
- 238000005516 engineering process Methods 0.000 claims description 7
- 230000001360 synchronised effect Effects 0.000 claims description 7
- 230000003287 optical effect Effects 0.000 claims description 5
- 229910000838 Al alloy Inorganic materials 0.000 claims description 4
- 239000000956 alloy Substances 0.000 claims description 4
- 230000003647 oxidation Effects 0.000 claims description 4
- 238000007254 oxidation reaction Methods 0.000 claims description 4
- 230000005457 Black-body radiation Effects 0.000 claims description 3
- 239000005357 flat glass Substances 0.000 claims description 3
- 239000004065 semiconductor Substances 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims description 2
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 8
- 230000002452 interceptive effect Effects 0.000 description 6
- 238000005259 measurement Methods 0.000 description 4
- 229910052757 nitrogen Inorganic materials 0.000 description 4
- 238000007747 plating Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 239000007788 liquid Substances 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000007613 environmental effect Effects 0.000 description 2
- 230000006641 stabilisation Effects 0.000 description 2
- 238000011105 stabilization Methods 0.000 description 2
- 241001270131 Agaricus moelleri Species 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000001035 drying Methods 0.000 description 1
- 238000007667 floating Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000003331 infrared imaging Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 238000000545 stagnation point adsorption reflectometry Methods 0.000 description 1
- 235000019640 taste Nutrition 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V8/00—Prospecting or detecting by optical means
- G01V8/10—Detecting, e.g. by using light barriers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J2005/0077—Imaging
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
- Radiation Pyrometers (AREA)
Abstract
本发明公开一种用于海面太阳亮带内目标识别的红外探测装置及方法,红外探测装置是两台中红外相机的镜头成一定角度,其中一台中红外相机的镜头前方设有红外偏振片,并在两台中红外相机的镜头前方的光路交叉处设置线栅偏振器,线栅偏振器的前方光路设置有光阑,光阑前方设有锗封窗;锗封窗和控温罩壳构成密封控温腔体,中红外相机、红外偏振片、线栅偏振器、光阑均置于控温腔体中。探测时先对红外探测装置的系统定标,采集海面太阳亮带图像,经图像处理后最终实现目标的识别。本发明提高了目标和周围海域的图像灰度差异,减少了亮带海域出现的等灰度块状海域,提高了海面太阳亮带内的目标识别能力。
Description
技术领域
本发明属于海面目标识别技术领域,特别涉及一种用于海面太阳亮带内目标识别的红外探测装置及方法。
背景技术
晴日,海面光学成像探测设备迎着太阳观测海面目标时,太阳辐射会经海面反射至探测设备。探测图像上表现为一条耀眼闪烁的带状区域,称之为海面太阳亮带。海面太阳亮带严重影响光学成像探测设备对亮带海域内的目标探测,导致其目标识别能力大打折扣,如图1和图2所示,所述海面目标包括海平面附近的岛屿、岩礁和各种海面漂浮物等。
国内刘燕、沈国土、靳兰海等人研究了海面太阳亮带对海面红外目标识别的影响,发现海面太阳亮带对中红外波段的海面目标识别有很大干扰,目标会被湮没在太阳亮带海域中。罗姣采用目标轨迹预测方法预估目标在亮带内的定位。目标进入亮带海域时,依靠目标进入亮带前的运动轨迹和速度,预测更新目标位置。轨迹预测方法可应用于遵循既定航线的目标跟踪,无法跟踪故意沿着亮带海域行进的敌方目标。
现有海面太阳亮带内目标识别技术无法准确识别目标的原因主要包括以下3个方面:
(1)部分亮带海域的辐射亮度与目标接近。相机像元对应的视场海域内会包含许多微小的太阳辐射反射区域,导致部分亮带海域对应的相机像元的AD输出值与目标相差不大。
(2)部分亮带海域出现大量等灰度的块状海域,成为了干扰识别的假目标。
(3)海面太阳亮带造成红外成像设备的红外相机饱和。
美国蒙特利市海军研究院物理系的Alfred W.Cooper和美国SPAR宇航公司的Dennis J、Gregoris等人,采用红外相机镜头前方加装偏振片,研究旋转偏振片对海面太阳亮带内目标识别的改善效果。美国海军武器中心下属的红外与可见光部的J.L.Beard也尝试使用一片红外偏振片抑制海面亮带的干扰,研究了最佳的偏振透振方向。相机镜头前方加装偏振片可消除部分海面反射的太阳辐射,一定程度上提高图像信噪比,但剩余海面反射的太阳辐射亮度仍远大于海面目标辐亮度,目标依旧被隐藏在海面亮带。
综上所述,目前公开文献上,暂未见准确识别太阳亮带海域内目标的研究报道,因此研发一种探测装置,能消除海面的强烈反射太阳辐射干扰,是实现海面太阳亮带内的目标识别的关键。
发明内容
为克服现有海面目标识别设备不能识别海面太阳亮带内目标,本发明提供一种用于海面太阳亮带内目标识别的红外探测装置及方法。
本发明通过以下技术方案实现:一种用于海面太阳亮带内目标识别的红外探测装置,包括中红外相机和红外偏振片,其特征在于:还包括线栅偏振器、光阑、同步触发器、控温罩壳和锗封窗,所述中红外相机设置为两台,两台中红外相机的镜头成一定角度,其中一台中红外相机的镜头前方设有红外偏振片,以保证进入该中红外相机的海面辐射都经过红外偏振片;并在两台中红外相机的镜头前方的光路交叉处设置线栅偏振器,线栅偏振器的前方光路设置有光阑,光阑用以限制非被测海域的环境辐射杂散光,光阑前方设有锗封窗;锗封窗和控温罩壳构成密封控温腔体,中红外相机、红外偏振片、线栅偏振器、光阑均置于控温腔体中;两台中红外相机和同步触发器通过信号线连接。
所述红外偏振片的透振方向与线栅偏振器反射的偏振辐射的偏振方向相同。
所述成一定角度是要满足下列条件:入射至线栅偏振器的海面辐射通过透射和反射被分成两束偏振方向正交且传播方向不同的偏振辐射,两束偏振辐射的偏振方向一个平行于海平面,一个垂直海平面,且两束线偏振辐射的中心线和两台中红外相机的光轴分别重合;两束偏振辐射分别在两台中红外相机成像,两幅图像对应相同的成像区域。
所述控温罩壳上设有线孔,控温罩壳内的中红外相机成像控制电路通过信号线穿过线孔连接控温罩壳外的同步触发器。
所述同步触发器内置有触发信号发生芯片,同步触发器发出的触发信号一分为二,分别输入给两台中红外相机。
所述两台中红外相机的性能配置相同,中红外相机的成像控制电路中设有用于同步触发功能的芯片,可以接受同步触发器输入的触发信号,触发控制两台中红外相机同步拍摄。
所述控温罩壳由铝合金材质制成,控温罩壳内表面阳极氧化发黑且设有消光沟槽,用于减少仪器热辐射的干扰;控温罩壳外表面设有基于半导体控温技术的控温层。
所述锗封窗为锗平板玻璃,镀增透膜以提高中红外相机工作波段辐射的透过率。
所述线栅偏振器为金属线栅偏振器,线栅偏振器的工作波段覆盖中红外相机的工作波段,线栅偏振器反射光的偏振方向平行于线栅偏振器的反射面。
本发明的另一目的在于提供一种用于海面太阳亮带内目标识别的方法,采用如上所述的红外探测装置,具体包括以下步骤:
步骤1、红外探测装置的系统定标
分别对两台中红外相机的自身响应误差、自身热辐射误差和辐射分束误差进行定标:
步骤1.1、两台中红外相机响应误差的辐射定标:
中红外相机的辐射定标用于获取相机的绝对辐亮度响应度α和相机光机结构热辐射,偏置电压及暗电流等引起的偏移量G0;采用大面源黑体覆盖中红外相机入瞳的定标方法,对两台中红外相机进行辐射定标,中红外相机的像元输出AD值即模数转换器的输出值,用Gt表示为:
Gt=αL+G0 (1)
式中,α为中红外相机的绝对辐亮度响应度;G0是中红外相机的光机结构热辐射、偏置电压、暗电流噪声等引起的偏移量;L为中红外相机的某像元所对应的海域在该像元位置处的总辐亮度;
两台中红外相机的积分时间为t,处于固定温度Ta下,通过测量不同温度的面源黑体,获取不同温度时每个像元的L和Gt值,进而对两台中红外相机辐射定标,分别获得两台中红外相机的二维矩阵G01、α1和二维矩阵G02、α2;
步骤1.2、红外探测装置的光机结构热辐射杂散光定标:
红外探测装置的中红外相机、线栅偏振器、红外偏振片和光阑置于控温罩壳内,控温罩壳的控温值与中红外相机定标时的温度Ta相同;控温后,控温罩壳、线栅偏振器、红外偏振片、光阑、中红外相机的温度稳定,热辐射也稳定;锗封窗镀增透膜后,热发射率很低,锗封窗受外界环境温度起伏影响引起的锗封窗热辐射起伏很微弱,红外探测装置自身光路内的热辐射杂散光基本稳定不变;
定标方法:用一个冷面完全覆盖锗封窗端部,冷面是指特殊物体表面,该特殊物体需满足下列条件:冷面自身热辐射在中红外相机的响应小于中红外相机的噪声,引起中红外相机的像元输出AD值可忽略不计;将两台中红外相机的积分时间设为t,控温罩壳的控温值设为Ta,启动红外探测装置拍摄图像,两台中红外相机分别获得的图像利用MATLAB、IDL(Interactive Data Language)等计算机语言读取得到二维数字矩阵,即为两台中红外相机在积分时间t下对应的热辐射杂散光定标矩阵R1和R2,矩阵R1和R2具有相同的行数和列数,且等于中红外相机的像元行、列数;
步骤1.3、线栅偏振器的辐射分束定标:
设置两台中红外相机的积分时间为t,控温罩壳的控温值为Ta,面源黑体置于红外探测装置入瞳前且覆盖红外探测装置入瞳,黑体辐射经线栅偏振器和红外偏振片分光后,两台中红外相机分别成像,所得图像分别利用MATLAB、IDL(Interactive Data Language)等计算机语言读取得到二维数字矩阵,即矩阵B1和B2,矩阵B1和B2具有相同的行数和列数且等于中红外相机像元的行、列数;
矩阵B1的每个元素值,分别减去步骤1.1所得矩阵G01的对应值,再除以矩阵α1的对应值,然后减去步骤1.2所得矩阵R1相对应的元素值,即得到矩阵C1;同理,矩阵B2的每个元素值,分别减去步骤1.1所得矩阵G02的对应值,再除以矩阵α2的对应值,然后减去步骤1.2所得矩阵R2相对应的元素值,即得到矩阵C2;矩阵C1和C2的行、列数等于中红外相机像元的行、列数;
与C1具有相同的行数和列数的全1矩阵H为其中一台中红外相机在积分时间t下的分束定标矩阵,以矩阵Ca为另一台中红外相机在积分时间t下的分束定标矩阵,矩阵Ca的每个元素值等于矩阵C2相对应元素值除以C1相对应元素值;
步骤2:红外探测装置进行海面太阳亮带内的目标探测
步骤2.1同步采集图像:
设置两台中红外相机的积分时间为t,控温罩壳的控温值为Ta,选取待测海域,红外探测装置的锗封窗一端指向待测海域,启动同步触发器,触发信号一分为二同步触发两台中红外相机拍摄,获取两台中红外相机对待测海域的两幅同时刻图像,将图像分别利用MATLAB、IDL(Interactive Data Language)等计算机语言读取得到二维数字矩阵,即矩阵D1和D2;
矩阵D1的每个元素值,分别减去步骤1.1所得矩阵G01的对应值,再除以矩阵α1的对应值,然后减去步骤1.2所得矩阵R1相对应的元素值,再乘以步骤1.3生成的全1矩阵H,上述运算消除了红外探测装置本身对海面辐射测量精度的影响,获得系统定标后AD输出值矩阵A1;
同理,矩阵D2的每个元素值,分别减去步骤1.1所得矩阵G02的对应值,再除以矩阵α2的对应值,然后减去步骤1.2所得矩阵R2相对应的元素值,再乘以步骤1.3生成的矩阵Ca,获得系统定标后AD输出值矩阵A2;矩阵A1和A2具有相同的行数和列数;
步骤2.2图像处理:
通过以下公式获得矩阵A3和A4:
A3=0.5(A1-A2-|A1-A2|) (2)
A4=A1+A2 (3)
式中,|A1-A2|表示矩阵,该矩阵的每个元素值等于A1和A2相对应元素值之差的绝对值,矩阵A3的每个元素值等于矩阵A1相对应元素值减去矩阵A2相对应元素值,再减去矩阵|A1-A2|相对应元素值,所得值乘以0.5;矩阵A4的每个元素值等于A1和A2相对应元素值的加和值;
式中,|A3|表示矩阵,该矩阵的每个元素值等于矩阵A3相对应元素值的绝对值;max(|A3|)为矩阵|A3|中的最大元素值,N为中红外相机的模数转换器(ADC)的位数,矩阵A5是矩阵|A3|与的数乘所获得的矩阵,矩阵A6获取方法与A5相同;
A5和A6为二维数字矩阵,利用MATLAB、IDL等计算机编程语言转化为人眼可识别的图像;
步骤2.3目标识别:
海面太阳亮带反射的太阳辐射具有很强的偏振特性,海面目标辐射包括自身辐射和太阳辐射在其表面的漫反射,偏振特性较弱;若海面目标辐射为无偏振,理论情况下,A5对应图像的目标区域为黑色区域;实际的各种海面目标都有一定偏振性,但偏振性很弱,图像中的目标区域接近黑色区域;海面亮带反射的太阳辐射有很强的偏振特性,越靠近亮带中心区域,偏振特性越强,对应的数值越大,越接近白色区域;A5对应图像中白色区域内的暗区域(即黑色区域或接近黑色区域)即为目标区域。
与现有技术相比,本发明具有如下优点:
(1)红外探测装置设有整体控温罩壳,抑制了探测设备中各器件的性能参数随温度变化产生的测量误差。结合系统定标,消除了辐射分束误差和红外探测装置自身热辐射杂散光的干扰。
(2)海面反射的太阳辐射具有较强的偏振特性,本发明装置生成的A5提高了目标和周围海域的图像灰度差异,减少了亮带海域出现的等灰度块状海域,提高了海面太阳亮带内的目标识别能力。线栅偏振器具有减光功能,有利于避免相机饱和。
(3)实现有、无太阳亮带两种情况下的海面目标识别。红外探测装置同时获得对应A5和A6的两幅图像,A5用于目标处于海面太阳亮带内的目标识别,A6用于没有太阳亮带时的海面目标识别。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其它的附图。
图1为海面太阳亮带的可见光图像;
图2为海面太阳亮带的中红外图像;
图3为本发明提供红外探测装置的结构示意图;图中,1-锗封窗,2-光阑,3-中红外相机Ⅰ,4-控温罩壳,5-线栅偏振器,6-红外偏振片,7-中红外相机Ⅱ;
图4为本发明提供红外探测装置的外观示意图;
图5为本发明提供红外探测装置的工作原理示意图;
图6为中红外相机Ⅰ采集的待测海域图像;
图7为中红外相机Ⅱ采集的待测海域图像;
图8为A5转化为的人眼可识别的图像;
图9为A6转化为的人眼可识别的图像;
图10为传统中红外相机拍摄的海面亮带图像。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有付出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图3、4所示,用于海面太阳亮带内目标识别的红外探测装置,包括中红外相机和红外偏振片,还包括线栅偏振器5、光阑2、同步触发器、控温罩壳4和锗封窗1,所述中红外相机设置为两台(中红外相机Ⅰ3和中红外相机Ⅱ7),两台中红外相机的镜头成一定角度,即满足下列条件:入射至线栅偏振器5的海面辐射通过透射和反射被分成两束偏振方向正交且传播方向不同的偏振辐射,两束偏振辐射的偏振方向一个平行于海平面,一个垂直海平面,且两束线偏振辐射的中心线和两台中红外相机的光轴分别重合;两束偏振辐射分别在两台中红外相机成像,两幅图像对应相同的成像区域。
所述中红外相机Ⅱ7的镜头前方设有红外偏振片6,以保证进入中红外相机Ⅱ7的海面辐射都经过红外偏振片6;并在两台中红外相机的镜头前方的光路交叉处设置线栅偏振器5,线栅偏振器5的前方光路设置有光阑2,光阑2用以限制非被测海域的环境辐射杂散光,光阑2前方设有锗封窗1;锗封窗1和控温罩壳4构成密封控温腔体,中红外相机Ⅰ3、中红外相机Ⅱ7、红外偏振片6、线栅偏振器5、光阑2均置于控温腔体中;控温罩壳4上设有线孔,控温罩壳4内的中红外相机成像控制电路通过信号线穿过线孔连接控温罩壳4外的同步触发器;同步触发器内置有触发信号发生芯片,同步触发器发出的触发信号一分为二,分别输入给两台中红外相机。
所述红外偏振片6的透振方向与线栅偏振器5反射的偏振辐射的偏振方向相同。所述两台中红外相机的性能配置相同,中红外相机的成像控制电路中设有用于同步触发功能的芯片,可以接受同步触发器输入的触发信号,触发控制两台中红外相机同步拍摄。所述控温罩壳4由铝合金材质制成,控温罩壳4内表面阳极氧化发黑且设有消光沟槽,用于减少仪器热辐射的干扰;控温罩壳4外表面设有基于半导体控温技术的控温层。
所述锗封窗1为锗平板玻璃,镀增透膜以提高中红外相机工作波段辐射的透过率。所述线栅偏振器5为金属线栅偏振器,线栅偏振器5的工作波段覆盖中红外相机的工作波段,线栅偏振器5反射光的偏振方向平行于线栅偏振器的反射面。
用于海面太阳亮带内目标识别的方法,采用如上所述的红外探测装置,具体包括以下步骤:
步骤1、红外探测装置的系统定标
分别对两台中红外相机的自身响应误差、自身热辐射误差和辐射分束误差进行定标:
步骤1.1、两台中红外相机响应误差的辐射定标:
中红外相机的辐射定标用于获取相机的绝对辐亮度响应度α和相机光机结构热辐射,偏置电压及暗电流等引起的偏移量G0;采用大面源黑体覆盖中红外相机入瞳的定标方法,对两台中红外相机进行辐射定标,中红外相机的像元输出AD值即模数转换器的输出值,用Gt表示为:
Gt=αL+G0 (1)
式中,α为中红外相机的绝对辐亮度响应度;G0是中红外相机的光机结构热辐射、偏置电压、暗电流噪声等引起的偏移量;L为中红外相机的某像元所对应的海域在该像元位置处的总辐亮度;
两台中红外相机的积分时间为t,处于固定温度Ta=25℃下,通过测量不同温度的面源黑体,获取不同温度时每个像元的L和Gt值,中红外相机各像元对应的AD输出值矩阵Gti(k,m);中红外相机的红外探测器为红外焦平面阵列探测器,其中k,m表示像元在红外焦平面阵列探测器中的行/列值;中红外相机探测温度Ti的面源黑体时,像元(k,m)对应的AD输出值为Gti(k,m),在横/纵坐标分别为辐亮度和像元AD值的坐标系中对应一个点,绘出不同黑体温度对应的点图,按照最小二乘法把点图拟合成一条直线,直线的斜率与纵轴的交点即为像元(k,m)的α(k,m)和G0(k,m)值;计算出每个像元对应α和G0值,进而分别获得中红外相机Ⅰ3在积分时间t下的二维矩阵G01、G02以及中红外相机Ⅱ7在积分时间t下的二维矩阵α1、α2;
步骤1.2、红外探测装置的光机结构热辐射杂散光定标:
红外探测装置的中红外相机、线栅偏振器、红外偏振片和光阑置于控温罩壳内,控温罩壳的控温值与中红外相机定标时的温度Ta相同;控温后,控温罩壳、线栅偏振器、红外偏振片、光阑、中红外相机的温度稳定,热辐射也稳定;锗封窗镀增透膜后,热发射率很低,锗封窗受外界环境温度起伏影响引起的锗封窗热辐射起伏很微弱,红外探测装置自身光路内的热辐射杂散光基本稳定不变;
定标方法:用一个冷面完全覆盖锗封窗端部,冷面是指特殊物体表面,该特殊物体需满足下列条件:冷面自身热辐射在中红外相机的响应小于中红外相机的噪声,引起中红外相机的像元输出AD值可忽略不计;本例采用铝合金材质制成中空圆柱薄壳体,圆柱体上表面进行阳极氧化发黑,作为冷面,下表面有液氮注入口,注入液氮后,用干燥氮气冲刷冷面;由于液氮的温度极低,冷面热辐射引起的红外探测装置的中红外相机AD输出值可忽略不计;红外探测装置测量冷面时,中红外相机AD输出值来源于红外探测装置自身热辐射;
将两台中红外相机的积分时间设为t,控温罩壳的控温值设为Ta=25℃,启动红外探测装置拍摄图像,两台中红外相机分别获得的图像利用MATLAB、IDL(Interactive DataLanguage)等计算机语言读取得到二维数字矩阵,即为中红外相机Ⅰ3在积分时间t下对应的热辐射杂散光定标矩阵R1以及中红外相机Ⅱ7在积分时间t下对应的热辐射杂散光定标矩阵R2,矩阵R1和R2具有相同的行数和列数,且等于中红外相机的像元行、列数;
步骤1.3、线栅偏振器的辐射分束定标:
设置两台中红外相机的积分时间为t,控温罩壳的控温值为Ta=25℃,面源黑体置于红外探测装置入瞳前且覆盖红外探测装置入瞳,黑体辐射经线栅偏振器和红外偏振片分光后,两台中红外相机分别成像,所得图像分别利用MATLAB、IDL(Interactive DataLanguage)等计算机语言读取得到二维数字矩阵,即矩阵B1和B2,矩阵B1和B2具有相同的行数和列数且等于中红外相机像元的行、列数;
矩阵B1的每个元素值,分别减去步骤1.1所得矩阵G01的对应值,再除以矩阵α1的对应值,然后减去步骤1.2所得矩阵R1相对应的元素值,即得到矩阵C1;同理,矩阵B2的每个元素值,分别减去步骤1.1所得矩阵G02的对应值,再除以矩阵α2的对应值,然后减去步骤1.2所得矩阵R2相对应的元素值,即得到矩阵C2;矩阵C1和C2的行、列数等于中红外相机像元的行、列数;
与C1具有相同的行数和列数的全1矩阵H为中红外相机Ⅰ3在积分时间t下的分束定标矩阵,以矩阵Ca为中红外相机Ⅱ7在积分时间t下的分束定标矩阵,矩阵Ca的每个元素值等于矩阵C2相对应元素值除以C1相对应元素值;
步骤2:红外探测装置进行海面太阳亮带内的目标探测
步骤2.1同步采集图像:
设置两台中红外相机的积分时间为t,控温罩壳的控温值为Ta=25℃,选取待测海域,红外探测装置的锗封窗一端指向待测海域,启动同步触发器,触发信号一分为二同步触发两台中红外相机拍摄,获取两台中红外相机对待测海域的两幅同时刻图像,如图6和图7,将图像分别利用MATLAB、IDL(Interactive Data Language)等计算机语言读取得到二维数字矩阵,即矩阵D1和D2;
矩阵D1的每个元素值,分别减去步骤1.1所得矩阵G01的对应值,再除以矩阵α1的对应值,然后减去步骤1.2所得矩阵R1相对应的元素值,再乘以步骤1.3生成的全1矩阵H,上述运算消除了红外探测装置本身对海面辐射测量精度的影响,获得系统定标后AD输出值矩阵A1;
同理,矩阵D2的每个元素值,分别减去步骤1.1所得矩阵G02的对应值,再除以矩阵α2的对应值,然后减去步骤1.2所得矩阵R2相对应的元素值,再乘以步骤1.3生成的矩阵Ca,获得系统定标后AD输出值矩阵A2;矩阵A1和A2具有相同的行数和列数;
步骤2.2图像处理:
通过以下公式获得矩阵A3和A4:
A3=0.5(A1-A2-|A1-A2|) (2)
A4=A1+A2 (3)
式中,|A1-A2|表示矩阵,该矩阵的每个元素值等于A1和A2相对应元素值之差的绝对值,矩阵A3的每个元素值等于矩阵A1相对应元素值减去矩阵A2相对应元素值,再减去矩阵|A1-A2|相对应元素值,所得值乘以0.5;矩阵A4的每个元素值等于A1和A2相对应元素值的加和值;
式中,|A3|表示矩阵,该矩阵的每个元素值等于矩阵A3相对应元素值的绝对值;max(|A3|)为矩阵|A3|中的最大元素值,N为中红外相机的模数转换器(ADC)的位数,矩阵A5是矩阵|A3|与的数乘所获得的矩阵,矩阵A6获取方法与A5相同;
A5和A6为二维数字矩阵,利用MATLAB、IDL等计算机编程语言转化为人眼可识别的图像,如图8和图9;
步骤2.3目标识别:
海面太阳亮带反射的太阳辐射具有很强的偏振特性,海面目标辐射包括自身辐射和太阳辐射在其表面的漫反射,偏振特性较弱;若海面目标辐射为无偏振,理论情况下,A5对应图像的目标区域为黑色区域;实际的各种海面目标都有一定偏振性,但偏振性很弱,图像中的目标区域接近黑色区域;海面亮带反射的太阳辐射有很强的偏振特性,越靠近亮带中心区域,偏振特性越强,对应的数值越大,越接近白色区域;A5对应图像中白色区域内的暗区域(即黑色区域或接近黑色区域)即为目标区域,如图8中能明显可见,提高了海面太阳亮带内的目标识别能力。A6表示两个正交方向偏振辐射分量的和,A6和传统中红外相机拍摄的海面亮带图像(如图10)基本相同。
Claims (10)
1.一种用于海面太阳亮带内目标识别的红外探测装置,包括中红外相机和红外偏振片,其特征在于:还包括线栅偏振器、光阑、同步触发器、控温罩壳和锗封窗,所述中红外相机设置为两台,两台中红外相机的镜头成一定角度,其中一台中红外相机的镜头前方设有红外偏振片,并在两台中红外相机的镜头前方的光路交叉处设置线栅偏振器,线栅偏振器的前方光路设置有光阑,光阑前方设有锗封窗;锗封窗和控温罩壳构成密封控温腔体,中红外相机、红外偏振片、线栅偏振器、光阑均置于控温腔体中;两台中红外相机和同步触发器通过信号线连接。
2.根据权利要求1所述的用于海面太阳亮带内目标识别的红外探测装置,其特征在于:所述红外偏振片的透振方向与线栅偏振器反射的偏振辐射的偏振方向相同。
3.根据权利要求1所述的用于海面太阳亮带内目标识别的红外探测装置,其特征在于:所述成一定角度是要满足下列条件:入射至线栅偏振器的海面辐射通过透射和反射被分成两束偏振方向正交且传播方向不同的偏振辐射,两束偏振辐射的偏振方向一个平行于海平面,一个垂直海平面,且两束线偏振辐射的中心线和两台中红外相机的光轴分别重合;两束偏振辐射分别在两台中红外相机成像,两幅图像对应相同的成像区域。
4.根据权利要求1所述的用于海面太阳亮带内目标识别的红外探测装置,其特征在于:所述控温罩壳上设有线孔,控温罩壳内的中红外相机成像控制电路通过信号线穿过线孔连接控温罩壳外的同步触发器。
5.根据权利要求1所述的用于海面太阳亮带内目标识别的红外探测装置,其特征在于:所述同步触发器内置有触发信号发生芯片,同步触发器发出的触发信号一分为二,分别输入给两台中红外相机。
6.根据权利要求1所述的用于海面太阳亮带内目标识别的红外探测装置,其特征在于:所述两台中红外相机的性能配置相同,中红外相机的成像控制电路中设有用于同步触发功能的芯片。
7.根据权利要求1所述的用于海面太阳亮带内目标识别的红外探测装置,其特征在于:所述控温罩壳由铝合金材质制成,控温罩壳内表面阳极氧化发黑且设有消光沟槽;控温罩壳外表面设有基于半导体控温技术的控温层。
8.根据权利要求1所述的用于海面太阳亮带内目标识别的红外探测装置,其特征在于:所述锗封窗为锗平板玻璃,镀增透膜。
9.根据权利要求1所述的用于海面太阳亮带内目标识别的红外探测装置,其特征在于:所述线栅偏振器为金属线栅偏振器,线栅偏振器的工作波段覆盖中红外相机的工作波段,线栅偏振器反射光的偏振方向平行于线栅偏振器的反射面。
10.一种用于海面太阳亮带内目标识别的方法,采用权利要求1至9中任一所述的红外探测装置,其特征在于包括以下步骤:
步骤1、红外探测装置的系统定标
分别对两台中红外相机的自身响应误差、自身热辐射误差和辐射分束误差进行定标:
步骤1.1、两台中红外相机响应误差的辐射定标:
中红外相机的像元输出AD值即模数转换器的输出值,用Gt表示为:
Gt=αL+G0 (1)
式中,α为中红外相机的绝对辐亮度响应度;G0是中红外相机的光机结构热辐射、偏置电压、暗电流噪声等引起的偏移量;L为中红外相机的某像元所对应的海域在该像元位置处的总辐亮度;
两台中红外相机的积分时间为t,处于固定温度Ta下,通过测量不同温度的面源黑体,获取不同温度时每个像元的L和Gt值,进而对两台中红外相机辐射定标,分别获得两台中红外相机的二维矩阵G01、α1和二维矩阵G02、α2;
步骤1.2、红外探测装置的光机结构热辐射杂散光定标:
用一个冷面完全覆盖锗封窗端部,冷面是指特殊物体表面,该特殊物体需满足下列条件:冷面自身热辐射在中红外相机的响应小于中红外相机的噪声,引起中红外相机的像元输出AD值可忽略不计;将两台中红外相机的积分时间设为t,控温罩壳的控温值设为Ta,启动红外探测装置拍摄图像,两台中红外相机分别获得的图像利用计算机语言读取得到二维数字矩阵,即为两台中红外相机在积分时间t下对应的热辐射杂散光定标矩阵R1和R2;
步骤1.3、线栅偏振器的辐射分束定标:
设置两台中红外相机的积分时间为t,控温罩壳的控温值为Ta,面源黑体置于红外探测装置入瞳前且覆盖红外探测装置入瞳,黑体辐射经线栅偏振器和红外偏振片分光后,两台中红外相机分别成像,所得图像分别利用计算机语言读取得到二维数字矩阵,即矩阵B1和B2;
矩阵B1的每个元素值,分别减去步骤1.1所得矩阵G01的对应值,再除以矩阵α1的对应值,然后减去步骤1.2所得矩阵R1相对应的元素值,即得到矩阵C1;同理,矩阵B2的每个元素值,分别减去步骤1.1所得矩阵G02的对应值,再除以矩阵α2的对应值,然后减去步骤1.2所得矩阵R2相对应的元素值,即得到矩阵C2;
与C1具有相同的行数和列数的全1矩阵H为其中一台中红外相机在积分时间t下的分束定标矩阵,以矩阵Ca为另一台中红外相机在积分时间t下的分束定标矩阵,矩阵Ca的每个元素值等于矩阵C2相对应元素值除以C1相对应元素值;
步骤2:红外探测装置进行海面太阳亮带内的目标探测
步骤2.1同步采集图像:
设置两台中红外相机的积分时间为t,控温罩壳的控温值为Ta,选取待测海域,红外探测装置的锗封窗一端指向待测海域,启动同步触发器,触发信号一分为二同步触发两台中红外相机拍摄,获取两台中红外相机对待测海域的两幅同时刻图像,将图像分别利用计算机语言读取得到二维数字矩阵,即矩阵D1和D2;
矩阵D1的每个元素值,分别减去步骤1.1所得矩阵G01的对应值,再除以矩阵α1的对应值,然后减去步骤1.2所得矩阵R1相对应的元素值,再乘以步骤1.3生成的全1矩阵H,获得系统定标后AD输出值矩阵A1;
同理,矩阵D2的每个元素值,分别减去步骤1.1所得矩阵G02的对应值,再除以矩阵α2的对应值,然后减去步骤1.2所得矩阵R2相对应的元素值,再乘以步骤1.3生成的矩阵Ca,获得系统定标后AD输出值矩阵A2;
步骤2.2图像处理:
通过以下公式获得矩阵A3和A4:
A3=0.5(A1-A2-|A1-A2|) (2)
A4=A1+A2 (3)
式中,|A1-A2|表示矩阵,该矩阵的每个元素值等于A1和A2相对应元素值之差的绝对值,矩阵A3的每个元素值等于矩阵A1相对应元素值减去矩阵A2相对应元素值,再减去矩阵|A1-A2|相对应元素值,所得值乘以0.5;矩阵A4的每个元素值等于A1和A2相对应元素值的加和值;
式中,|A3|表示矩阵,该矩阵的每个元素值等于矩阵A3相对应元素值的绝对值;max(|A3|)为矩阵|A3|中的最大元素值,N为中红外相机的模数转换器的位数,矩阵A5是矩阵|A3|与的数乘所获得的矩阵,矩阵A6获取方法与A5相同;
A5和A6为二维数字矩阵,利用计算机编程语言转化为人眼可识别的图像;
步骤2.3目标识别:
A5对应图像中白色区域内的暗区域即为目标区域。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611100992.0A CN106706132B (zh) | 2016-12-05 | 2016-12-05 | 一种用于海面太阳亮带内目标识别的红外探测装置及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611100992.0A CN106706132B (zh) | 2016-12-05 | 2016-12-05 | 一种用于海面太阳亮带内目标识别的红外探测装置及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106706132A true CN106706132A (zh) | 2017-05-24 |
CN106706132B CN106706132B (zh) | 2023-04-07 |
Family
ID=58935676
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611100992.0A Active CN106706132B (zh) | 2016-12-05 | 2016-12-05 | 一种用于海面太阳亮带内目标识别的红外探测装置及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106706132B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107860478A (zh) * | 2017-11-17 | 2018-03-30 | 北京长峰科威光电技术有限公司 | 一种红外成像系统冷反射黑斑的修复方法 |
CN108151887A (zh) * | 2017-12-25 | 2018-06-12 | 湖南航天诚远精密机械有限公司 | 一种微波实验炉 |
CN108761556A (zh) * | 2018-06-15 | 2018-11-06 | 北京理工大学 | 一种用于探测近岸水面太阳亮带内目标的方法 |
WO2020002837A1 (fr) | 2018-06-28 | 2020-01-02 | Safran Electronics & Defense | Dispositif d'imagerie infrarouge |
CN111609937A (zh) * | 2020-06-08 | 2020-09-01 | 北京环境特性研究所 | 一种红外热像仪外场标定方法及装置 |
CN113697031A (zh) * | 2021-09-10 | 2021-11-26 | 肖昌美 | 一种多用途浅吃水直通甲板船型 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1107182A2 (en) * | 1999-11-29 | 2001-06-13 | Eastman Kodak Company | Determining orientation of images containing blue sky |
CN104268831A (zh) * | 2014-09-12 | 2015-01-07 | 西安电子科技大学 | 海上阳光亮带干扰下的红外图像补偿方法 |
CN105043558A (zh) * | 2015-06-06 | 2015-11-11 | 中国科学院云南天文台 | 一种用于高反面红外辐射测量的屏蔽方法及装置 |
US20160078309A1 (en) * | 2014-09-12 | 2016-03-17 | Research Development Foundation | Apparatus and methods for identifying and evaluating bright spot indications observed through optical coherence tomography |
DE102014220818A1 (de) * | 2014-10-14 | 2016-04-14 | Signum Bildtechnik GmbH | Vorrichtung und Verfahren zur Reduzierung der Blendwirkung |
US20160234489A1 (en) * | 2014-06-23 | 2016-08-11 | Unlimited Optics Corp. | Method for measuring performance parameters and detecting bad pixels of an infrared focal plane array module |
CN206556765U (zh) * | 2016-12-05 | 2017-10-13 | 中国科学院云南天文台 | 一种用于海面太阳亮带内目标识别的红外探测装置 |
-
2016
- 2016-12-05 CN CN201611100992.0A patent/CN106706132B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1107182A2 (en) * | 1999-11-29 | 2001-06-13 | Eastman Kodak Company | Determining orientation of images containing blue sky |
US20160234489A1 (en) * | 2014-06-23 | 2016-08-11 | Unlimited Optics Corp. | Method for measuring performance parameters and detecting bad pixels of an infrared focal plane array module |
CN104268831A (zh) * | 2014-09-12 | 2015-01-07 | 西安电子科技大学 | 海上阳光亮带干扰下的红外图像补偿方法 |
US20160078309A1 (en) * | 2014-09-12 | 2016-03-17 | Research Development Foundation | Apparatus and methods for identifying and evaluating bright spot indications observed through optical coherence tomography |
DE102014220818A1 (de) * | 2014-10-14 | 2016-04-14 | Signum Bildtechnik GmbH | Vorrichtung und Verfahren zur Reduzierung der Blendwirkung |
CN105043558A (zh) * | 2015-06-06 | 2015-11-11 | 中国科学院云南天文台 | 一种用于高反面红外辐射测量的屏蔽方法及装置 |
CN206556765U (zh) * | 2016-12-05 | 2017-10-13 | 中国科学院云南天文台 | 一种用于海面太阳亮带内目标识别的红外探测装置 |
Non-Patent Citations (3)
Title |
---|
IRINA KULYK: "Brightness and polarization opposition effects at low phase angles of the Saturnian satellites Tethys, Dione, and Rhea", 《PLANETARY AND SPACE SCIENCE》 * |
周冠华等: "基于偏振原理的水面反射光的剥离", 《水科学进展》 * |
赵乃卓等: "多角度偏振遥感中太阳耀光剥离的深入研究", 《光学技术》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107860478A (zh) * | 2017-11-17 | 2018-03-30 | 北京长峰科威光电技术有限公司 | 一种红外成像系统冷反射黑斑的修复方法 |
CN107860478B (zh) * | 2017-11-17 | 2019-11-22 | 北京长峰科威光电技术有限公司 | 一种红外成像系统冷反射黑斑的修复方法 |
CN108151887A (zh) * | 2017-12-25 | 2018-06-12 | 湖南航天诚远精密机械有限公司 | 一种微波实验炉 |
CN108761556A (zh) * | 2018-06-15 | 2018-11-06 | 北京理工大学 | 一种用于探测近岸水面太阳亮带内目标的方法 |
CN108761556B (zh) * | 2018-06-15 | 2019-10-29 | 北京理工大学 | 一种用于探测近岸水面太阳亮带内目标的方法 |
WO2020002837A1 (fr) | 2018-06-28 | 2020-01-02 | Safran Electronics & Defense | Dispositif d'imagerie infrarouge |
FR3083416A1 (fr) * | 2018-06-28 | 2020-01-03 | Safran Electronics & Defense | Dispositif d'imagerie infrarouge |
US11125624B2 (en) | 2018-06-28 | 2021-09-21 | Safran Electronics & Defense | Infrared imaging device |
CN111609937A (zh) * | 2020-06-08 | 2020-09-01 | 北京环境特性研究所 | 一种红外热像仪外场标定方法及装置 |
CN111609937B (zh) * | 2020-06-08 | 2021-07-20 | 北京环境特性研究所 | 一种红外热像仪外场标定方法及装置 |
CN113697031A (zh) * | 2021-09-10 | 2021-11-26 | 肖昌美 | 一种多用途浅吃水直通甲板船型 |
Also Published As
Publication number | Publication date |
---|---|
CN106706132B (zh) | 2023-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106706132A (zh) | 一种用于海面太阳亮带内目标识别的红外探测装置及方法 | |
US9906718B2 (en) | Biomimetic integrated optical sensor (BIOS) system | |
CN102853916B (zh) | 一种用于对煤垛表面进行远距离红外测温的方法及系统 | |
CN110006463A (zh) | 一种光学遥感卫星的在轨绝对辐射定标方法及系统 | |
CN102538983B (zh) | Ccd测温装置 | |
CN104101297B (zh) | 一种基于光电观测的空间物体尺度获取方法 | |
CN110285712A (zh) | 一种实现动态拟态伪装隐身的装置 | |
CN102567977A (zh) | 一种基于小波的红外偏振图像自适应融合方法 | |
CN108287350A (zh) | 基于多信息综合的天基空中目标探测关键参数确定方法 | |
CN103792009A (zh) | 地基大口径望远镜的红外辐射定标方法 | |
CN105043552A (zh) | 比色测温系统显示与标定方法 | |
CN108120510A (zh) | 一种光学遥感器基于反射镜阵列的在轨绝对辐射定标方法 | |
Swamidoss et al. | Systematic approach for thermal imaging camera calibration for machine vision applications | |
CN206556765U (zh) | 一种用于海面太阳亮带内目标识别的红外探测装置 | |
CN108163223B (zh) | 一种便携式飞行器红外隐身性能评估装置和方法 | |
Deng et al. | A compact mid-wave infrared imager system with real-time target detection and tracking | |
Frías et al. | The Atmospheric Monitoring System of the JEM-EUSO space mission | |
Zhao et al. | Polarization imaging in atmospheric environment based on polarized reflectance retrieval | |
Cao et al. | Attitude measurement method of spinning projectile based on infrared focal plane array | |
CN106769895A (zh) | 一种标定测量整层大气光谱透过率的方法 | |
CN103411896B (zh) | 低轨空间碎片超低色散光谱特征获取方法 | |
Huang et al. | Research on radiometric calibration for super wide-angle staring infrared imaging system | |
CN111089651B (zh) | 一种渐变多光谱复合成像导引装置 | |
CN106251319A (zh) | 一种基于Curvelet变换的彩色图像融合算法 | |
Wu et al. | Study on large area laser reflective precision target |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |