大气传输引起的红外图像退化分析

第39卷第4期

V01.39No.4

红外与激光工程

InfraredandLaserEngineering

2010年8月

Aug.2010

大气传输引起的红外图像退化分析

李咏豪1,姜会林2,周

娜2

(1.北京特种车辆研究所,北京100072;2.长春理工大学光电工程学院,吉林长春130022)摘

要:分析并仿真了大气传榆引起的红外图像退化,.从大气湍流和气溶胶散射、吸收的物理过

程入手.分析了大气调制传递函数(MTF)及其造成的图像退化。提出了大气作用效果的计算机仿真方法,在目标距离为ll㈣、离地高度2m、C:=1.Oxl0…m…、光学系统为fig.04=600mrfl的情况下,针对远距离水平传输方向的红外成像系统进行了大气作用效果仿真.结果表明:大气传输引起的红外图像退化主要体现在由湍流引起的程度不同的局部模糊及图像跳动,

关键词:大气调制传递函数;中图分类号:TN929.12

红外仿真;

湍流;

散射;

吸收

文献标志码:A

文章编号:1007—2276(2010)04—0589—04

Analysisofinfraredimagedegradationby

LIYong-ha01,JIANGHui—lin2,ZHOUNa2

aero

transmission

(1.BeijingSpecialVehicleInstitute。Beijing100072,China;2.SchoolofOptoelectronicEnigeering,Changchun

ofScience

and

University

Technology,Changchun13(X)22,China)

Abstract:ne

analysisandsimulationoftheinfrared

effects

imagedegradationby

aero

transmission

were

presented.TheatmosphereMTFandthedegradationphysical

oftheimagewerepointedbyanalyzingthe

processofatmosphericturbulence,aerosolscatteringandabsorption.Thesimulationmethodsof

theatmosphericeffectswerepresented,andthecomputersimulationofatmosphericeffectswasdonefortheinfraredimaging

was

systemoflong—distance

horizontaltransmission

一l≈

on

condition

thatthetargetdistance

一2/3

km;theheightabovegroundwas2m;C。=x1.OxlO

showthattheinfrared

,opticalsystem

by

aero

is弘.0

and

f--600

mm.

Thesimulationresults

on

imagedegradationcaused

transmissionaremainly

focused

thelocalfuzzyandimageditheringinvariousdegreecausedbyatmosphereturbulence.

Infraredsimulation;Turbulence;

Keywords:AtmosphereMTF;

Scattering;Absorption

O引言

热成像系统用于观测远距离目标时,红外辐射在到达热成像系统的光学系统之前必须通过大气,并被

射源这些现象使远距离目标的热图像质量变差。研

究表明:肖红外成像系统用于远距离目标获取时,系统的整机性能受到系统自身MTF和大气湍流MTF

的很大限制|lI。因此.研究红外辐射的大气传输对像

质的影响,并对其进行建模和仿真,对红外成像系统的性能评估具有重要意义。

衰减。衰减的同时,大气梯度和湍流使辐射在结构上

发生畸变.引起空气折射率不均。此外,大气是固有辐

收稿日期:2009—10—10:修订日期:2(X)O—12一05

基金项目:国家自然科学基金资助项目(60977001);围家高技术研究发展计划(2007AAl2ZIlO)

作者简介:李咏豪(1977-)。男.山东招远人。I:程师.硕{:.主要研究方向为光电仿真‘j检测技术,Email:liyh@163.tom

万方数据

红外与激光工程第30卷

1红外成像中大气传输作用效果分析

l-l大气衰减

大气是一个包含湍流、散射和吸收粒子的随机介

质,大气对系统T作的影响基本体现在辐射衰减。红外辐射衰减与3种现象有关:大气气体分子的吸收;大气中分子、气溶胶和微粒的散射;冈气象条件(云、

雾、雨、雪)引起的衰减等。红外辐射的大气衰减满足布盖尔一朗伯定律,即光谱透过率r。(A)和衰减系数Ot。(A)之间满足如下关系:

下。(A)=exp【一d。(A)D】

(1)

式中:D为被观测目标和热成像系统间的距离。,

1.2大气湍流对成像的影响

红外辐射在大气中传播会造成图像模糊。原因有两个:一一是由于经气溶胶和其他大气成分散射后进人

人瞳的散射光:二是在光路中由湍流引起的介质折射

率变化。大气气溶胶散射造成的模糊通常被称为“邻接效果”。对于远距离的水平光路,由于接近地面的风和空气与地面之间的热交换造成温度的显著变化,使湍流成为影响图像质量的重要因素。即使存在气溶胶模糊.湍流也常常是模糊的决定性因素。湍流模糊常

被作为一个线性移不变物理过程来处理。湍流引起的系统性能下降可以用湍流MTF来描述。

Kopeika等提出r湍流MTF模型I创。该模型由古德曼提出的大气相十长度模型发展而来。在湍流

MTF模型中,大气的折射率变化用折射率结构参数t

表示。。湍流MTF分为长曝光模型和短曝光模型。长曝光模型集中反映r焦平面上模糊点位置随时间的变

化.这种变化是由传感器入瞳处的等晕波前倾斜引起

的。大气湍流的长曝光MTF是:

513

一l/3

M7FLE=exp(一57.40管eA

R)(2)式中:a为波形常数,当波形是球面时,a为3/8,当波

面是平面时,a为l;f为空间频率,c/mrad;C:为折射

率结构参数,m划3;A为波长;R为光程。大气湍流的短

曝光MTF是:

S13

—1/3l/3

M阿sE=exp(一57.4嘭eA

R)[1-Iz(Af/D)】(3)

式中:D为通光孔径;肛为无量纲系数,在近场和远场

万方数据

时取不同值.近场时取l,远场时取0.5。

长曝光记录的模糊光斑通常比短曝光的大,这是

因为造成模糊光斑的入射光角度随着光路的弯曲一直在变。使得光斑位置在像面上随机跳动,长曝光记录的是光斑跳动在时间上的积分。对于凝视型红外相机和红外监视型相机,其积分时间通常为毫秒级,典型值是

16.7ms,帧频60

Hz。所以。通常采用短曝光模型。

折射率结构参数C用下式来确定:

c。2=等

I,3

(4)

式中:<d‘>为波前入射角的变化量:D为光学孔径;尺

为物像间的距离。

1.3气溶胶对成像的影响

除了大气湍流,光线在大气传输过程中还受到分子和气溶胶产生的散射和吸收作用的影响。

散射作用大大扩展了光源的辐射范围,造成光线入射角度的H大差异.从而也造成像平面上图像的模糊。与短曝光湍流不同的是,由气溶胶造成的模糊是非常广泛的。散射光的角度受到成像系统视场的限制,并非所有的散射光都会出现在像面。此外,传感器的动态范围也限制了记录在像面的散射光角度.其影

响往往比视场的限制还大。如果在入瞳处由目标辐射

的散射光的辐射水平低于背景辐射,其对大气模糊的

贡献将是轻微的,通常町以忽略不计。接收器有限的

空间频率带宽限制了实际记录在像面的非散射光的最高频率。所有这些实际仪器的限制应适用于处理气

溶胶MTF的经典理论m41。

与可见光波段内不同的是。在热红外波段内吸收

是一个非常明显的效应。到达成像系统的入射光包含

非散射和散射两个部分。在长波红外大气窗口、中等气候条件下。长波红外的散射系数比可见光小得多,

光子通常只能与粒子发生一次作用,即要么被粒子吸

收。要么被散射15l。

因此,吸收只适用于非散射光的部分.而不适用

于所有辐射。所有与散射光有关的模糊,都不会因为微粒吸收而衰减。另外,能在无限空间频率带宽上增强图像的非散射光会被大气吸收而衰减,这使得图像

模糊作用更显著。因此,吸收造成气溶胶MTF在高空

间频率处衰减。气溶胶MTF可表示为:

第4期李咏豪等:大气传输引起的红外图像退化分析

:exp!、-s。Ri考”

w班i唧H刮1.f"v2忙卟≤K∞’

{exp卜SR}exp(-A。尺)}:VⅫ。

式中:S。和A。分别为大气有效的散射和吸收系数;尺

为路径长度;y为气溶胶MTF的角空间截止频率,

将湍流MTF和气溶胶MTF相乘。得到大气MTF:

MTF=M.M。

(6)

根据大气的热MTF分析,闩标与背景之间的温差通过有微粒的大气传播后将受到大气的调制.从而

使物体对比度下降、网像模糊、图像灰度值起伏及图

像形变、闪烁等忙引.

2大气退化作用仿真

2.1仿真方法

大气退化作用仿真器本身町以看作足时间函数F(・)和尺(・)组成的一个串联的系统.该系统作用于静

态输入图像l(x,Y),并创建・个输ff{图像序列O(u,',,t):

0(“,v,f)=JR(F(1(x,y)))

(7)

X和Y表爪输入图像的学标.U和v表示输出图像的坐标,输I叶I罔像序列的时间依赖性体现在t上。

第一个操作符F(・)用来模拟大气凋制传递厢数对图

像的模糊作用.定义为时间和夺问移不变的窄间滤波器

用高斯函数的标准差dw表爪模糊程度,

D(“.p,,)=(F(,(工,y);f)

(8)

第二个操作符尺(・)用来模拟|冬j像跳动,,这也是时间函数.且是卒间移不变操作符.O(u.V,f)=D(“(J,Y,,),v(x,Y,f),t(x,Y,f))=

R(1(x,y);d石(,),dy(,))=月(“工,y);,)(9)

这个操作处理输入图像的坐标为输入图像的每个位置计算出位移矢量.分别用出、dv分量表示横向

和纵向位移,图l、图2分别为滤波函数■维视l冬1和

滤波函数参数dw、血、dy爪意图.

存开始仿真之前.必须指定相机的特性、物距、大

每个像素。dw、血、dv的一一阶统计特性由前述公式给m.

万方数据

网l滤波函数一:维视网

网2滤波函数参数dw.出,dv爪崽网

Fig.1

VisualizationoftheFig.2Sketchof

parameters

dw.dx

filterfunction

dyofthefilterfunction

并由利用多重【日l归模型设计的随机数产乍器产生。

然后对输入图像使用操作符F(・)。需要计算㈣K

n个(m和n分别指输入网像在水平和垂直方向上的像素数)不同宽度的内核低通滤波。由此产幸的图像显爪出不同的局部模糊,,

接着对低通滤波后的图像运用操作符R(・)。这需

要汁算另一个肌x_,z的结果.根据每个输出罔像像素

的出和dv位移量.采用向后映射的方法找到对应的

输入图像中采样点的位置。向后映射能够将输入平面上各点随机分配到输出平面L,因此,可以精确模拟高动荡时的图像扭曲。使用双线性插值的方法考虑最近

的4个像素来计算该点的值,这种方法比最近邻法具

有更好的效果,仿真过程呵以用图3来表,】;,、

陶3仿真过稃小意陶

Fig.3Diagram

of

simulation

选用下列参数作为仿真对象:目标距离为lkm:}j标离地高度为2

m;C:为1.o×1.0—3m瑚(延光路方

向);光学系统为f/4.0,f=600mm,w=O.92。xO.73。;探

I.Lm:响应波根据公式(3)、(5)、(6)得出大气MTF,该传递函数

2.2仿真结果

测器探测元为320x256。像元尺寸为30段为7.7~9.5¨m,,

气情况等再将物理描述向仿真环境参数转换,,对于一次具体的仿真计算.需要为输m图像的每个像素计算

模糊程度和所期望的在指定动荡下的位置分布。对于

仅考虑大气传输的退化效果i,曲线图如图4所示。

592

红外与激光工程

第39卷

曲■

图4大气调制传递函数

Fig.4Atmosphere

modulationtransferfunction

对不同频率的条形靶进行了大气退化作用仿真。目的在于考察大气的独立退化作用。仿真结果如图5所示。经红外成像系统得到的实际红外图像还应包括系统自身及运动模糊、噪声等对图像的退化作用。

图5经大气传输退化后不同空间频率的条形靶仿真结果

Fig.5Simulatedbartargetsofdifferentfrequenciesdegraded

byatmospheretransmission

为了比较仿真结果与原始图像的差别,将仿真结果与原始图像做互相关运算,求出其互相关函数。为了简化计算,仅考虑一维的情况。

rfg(m)=E㈣咖训}=熙击蝥咖(棚)

(一M≤m≤M)

(10)

由于r缸与rf之间满足如下关系:

rfg(m)=rf(m)’h(m)

(11)

对公式(1)两端取傅里叶变换得:

Pfg∞)=Pf啦)‘H∞)

(12)

由上述关系式可求出退化过程的传递函数。图6

为在50lp/mm处的归一化互相关函数。

分别计算各频率处的互相关函数.并根据公式(12)求出退化过程的传递函数。图7为退化过程的传递函数与大气调制传递函数比较的结果。

万方数据

Spatialfreqtmcy/Ip・mm‘’

图650lp/mm处归一化互相关函数

Fig.6Normalized

cross-correlation

function

at

50lp/mm

函一

Spatial

frequncy/Ip’蛐‘

图7退化过程的传递甬数

Fig.7Transferfunctionofdegradation

从图7中可以看出:由于大气调制传递函数具有

大气传输引起的红外图像退化是个复杂的物理(1:转第626页)

较高的数值,大气传输造成的图像退化主要体现在湍流引起的介质折射率变化,造成调制传递函数的局部波动,使图像显示出不同程度的局部模糊及图像跳动,频率越高退化效果越明显。随着频率的升高,大气调制传递函数引起的边缘模糊变得更加明显,增加了对远距离目标探测和识别的难度。

3结论

过程。在远距离目标获取中,红外侦查相机的性能受到大气传输的限制。文中主要针对大气调制传递函数

和图像跳动.利用随机数产生器产生具有特定统计特性的串联操作符的3个参数dw、出、dy。对不同空间频率的条形靶进行了大气传输效果仿真,并将仿真结

果与原始图像进行了互相关分析,证明退化过程的传递函数与理论结果相符。算法主要针对的远距离水平传输方向的红外成像系统大气传输退化也可嵌入到其他成像仿真的图像处理过程中,与系统自身及运动模糊、噪声等对图像的退化作用进行联合仿真。大气

626红外与激光工程第39卷

法,根据电场作用下液滴面形变化的原理,直接在半导体激光管封装窗口上制作非球面透镜来准直LD光束。实验制得具有两个聚焦的非球面透镜,对LD光束

相互垂直方向的发散角分别进行准直。效果良好。

integrationof

lens

edge-emittinglaserdiodewith

ion

micro-elliptical

using

focused

on

beamdirect

deposition【J】.IEEE

Transactions(】):2—8.

Semiconductor

Manufacturing,2002,15

【8】LITong-hal,WUGuo-jun,HU

Bao-wen,etand

a1.Fabricationof

in

laser

参考文献:

FoR

thehemi-ellipticalmicrolensesbeam

application

diode

shaping【J】.ScienceTechnologyandEngineering(李同

KNERJF.KUN’lZ

W.Characteristicsofemcient

海,吴围俊,胡保文.等.椭圆微透镜的制备及其在LD光束整形中的应用.科学技术与工程),2005,24:1890一1892.

【9】GOUZhi-yong,WANGJiang,WANGLei,eta1.Diffraction・

limitedaspheric

collimation

lens【J】.Journal

of

Applied

laserdiode

collimators[C]//SPm,1997。740:145-151.

【2】

CHENGICe-hu扎CHENGDuo-si.Researchofthecollimationapplications

for

selfoc

lenseslJl.ActaPhotonicaSinica(程科

华,程多思.自聚焦透镜准直应用的长度计算与实验.光子学报),1992,21(2):145—151.

【31

XU

Xiao-fang,LILi-ua,WUJin-hui,eta1.Fibercoupling

of

the

Optics(勾志勇,乇江,王磊,等.衍射极限非球面准直透镜.

应用光学),2006,27(6):467—469.

【101

HE

Jian—hui,ZHAN

on

Zhen—xian,WU

Yun—liang,et

a1.

module

hi曲outputpower【J】.Infrared

andLaserMeasurement

surfaceprofileofasphericalliquidmicro—lens

Engineering(许孝芳,李丽娜,吴金辉,等.高功率半导体激光器列阵光纤耦合模块.红外与激光工程),2006,35(1):86—88.

【4】

BOKORN,DAVIDSONN.Beamshaping

use

【11.InfraredandLaserEngineering(何建慧,詹珍贤,吴云良.等.非球面液滴微透镜的轮廓测量方法.红外与激光工

程),2008,37(4):330—333.

【ll】ZHANG

with

diffuselightby

of

singlereflection【J】AppliedOptics,2001,40:2132.

Qiang,SHIYong-chao,LV

Nai-guang,et

a1.

【5】

SMITH

WJ.ModernOpticalEngineering[MI.NewYork:

Detectionofthe

positionoflaser—spotcenter【J1.Journalof

of

McGraw—Hill.2000.

【6】

ZHAN

Zhen-xian,HE

Beijing

Jian-hui,YAO

Hai—tao,et

a1.

Institute

MachineryIndustry(张强,施涌潮,吕

乃光.等.激光光斑中心位置的探测.北京机械工业学院学

报),1996,11(2):24—29.

Aspheficalliquidmicro—lensmanipulatedbyelectrostaticforces

【J】.ActaOpticaSinica(詹珍贤,何建慧,姚海涛,等.电场作用下的变焦非球面液滴微透镜.光学学报),2008,28(2):

361-364.

【7】

FU

Yong-qi,NGOI

K,BRYAN

A.A

【12】TAN

and

Zuo-jun,XUESong,KANGJing-ran,eta1.Collimation

ofsemiconductorlasers

in

test

laserfuze【J1.Journal

of

Applied

novel

one

Optics(谭佐军,薛松,康竞然,等.激光引信巾半导

step

体激光器的准直及其测试.应用光学),2007,28(4):454—457.

(上接第592页)

(3):354-358.

传输的非等晕效应仿真还有待进一步研究。

参考文献:

KRAPELS

K,DRIGGERSRG,VOLLMERHAUSENRH.ct

turbulence

【4】

SADOTDAN,DVIRA,BERGELLimagesdistortedtheoretical

by

eta1.Restorationofthermal

on

theatmosphere,based

modulation

measuredand

atmospheric

transfer

function【J】.

OpticalEngineering,1994,33:44—53.

【5】

SADOTD,KOPEIKANS.Imagingthroughtheatmosphere:practical

instrumentation・based

theory

and

verification

of

a1.Atmospheric

modulationtransferfunctionfor

infraredtargetacquisitionmodeling【J1.OptEng,2001,加(9):

1906—1913.

SADOTD,KOPEIKA

aerosolMTF【Cl//SPIE,1992,1688:59—64.

【61

S。ROTMANSR.Incorporationof

modeling[C]//

【7】

SADOTD,KITRONG,KITRONN,eta1.Thermalimagingthrough

the

atmosphere:Atmospheric

MTF

theory

and

atmosphericblurringeffectsintargetacquisitionSPIE,1994,2224:95—107.

verification【C]//SPIE,1993,1968:59—66.GAO

Si—feng,WUPing,HEMan—li,eta1.Estimation

on

the

ZHU

Wen-yue,MAXiao-shan,RAORui-zhong,eta1.Optical

effects

on

operating

rangeofinfraredsystemundercomplexatmosphere

turbulence

electro-optical

sensors[J】.Infraredand

condition[J].InfraredandLaserEngineering(高思峰,吴平.何曼丽,等.复杂大气条件下红外系统作用距离的估算.红外与激光工程),2008,37(6):941—944.

La赞r

Engineering(朱文越.马晓珊.饶瑞中.等.大气光学

湍流对光电探测器性能的影响.红外与激光工程),2006,3S

万方数据

大气传输引起的红外图像退化分析

作者:作者单位:刊名:英文刊名:年,卷(期):

李咏豪, 姜会林, 周娜, LI Yong-hao, JIANG Hui-lin, ZHOU Na

李咏豪,LI Yong-hao(北京特种车辆研究所,北京,100072), 姜会林,周娜,JIANG Hui-lin,ZHOU Na(长春理工大学光电工程学院,吉林,长春,130022)红外与激光工程

INFRARED AND LASER ENGINEERING2010,39(4)

参考文献(7条)

1. 高思峰;吴平;何曼丽 复杂大气条件下红外系统作用距离的估算[期刊论文]-红外与激光工程 2008(06)2. SADOT D;KITRON G;KITRON N Thermal imaging through the atmosphere:Atmospheric MrrF theory andverification 1993

3. SADOT D;KOPEIKA N S Imaging through the atmosphere:practical instrumentation-based theory andverification of aerosol MTF 1992

4. SADOT DAN;DVIR A;BERGEL L Restoration of thermal images distorted by the atmosphere,based onmeasured and theoretical atmospheric modulation transfer function 1994

5. 朱文越;马晓珊;饶瑞中 大气光学湍流对光电探测器性能的影响[期刊论文]-红外与激光工程 2006(03)

6. SADOT D;KOPEIKA N S;ROTMAN S R Incorporation of atmospheric blurring effects in target acquisitionmodeling 1994

7. KRAPELS K;DRIGGERS R G;VOLLMERHAUSEN R H Atmospheric turbulence modulation transfer function forinfrared target acquisition modeling[外文期刊] 2001(09)

本文链接:http://d.g.wanfangdata.com.cn/Periodical_hwyjggc201004002.aspx

第39卷第4期

V01.39No.4

红外与激光工程

InfraredandLaserEngineering

2010年8月

Aug.2010

大气传输引起的红外图像退化分析

李咏豪1,姜会林2,周

娜2

(1.北京特种车辆研究所,北京100072;2.长春理工大学光电工程学院,吉林长春130022)摘

要:分析并仿真了大气传榆引起的红外图像退化,.从大气湍流和气溶胶散射、吸收的物理过

程入手.分析了大气调制传递函数(MTF)及其造成的图像退化。提出了大气作用效果的计算机仿真方法,在目标距离为ll㈣、离地高度2m、C:=1.Oxl0…m…、光学系统为fig.04=600mrfl的情况下,针对远距离水平传输方向的红外成像系统进行了大气作用效果仿真.结果表明:大气传输引起的红外图像退化主要体现在由湍流引起的程度不同的局部模糊及图像跳动,

关键词:大气调制传递函数;中图分类号:TN929.12

红外仿真;

湍流;

散射;

吸收

文献标志码:A

文章编号:1007—2276(2010)04—0589—04

Analysisofinfraredimagedegradationby

LIYong-ha01,JIANGHui—lin2,ZHOUNa2

aero

transmission

(1.BeijingSpecialVehicleInstitute。Beijing100072,China;2.SchoolofOptoelectronicEnigeering,Changchun

ofScience

and

University

Technology,Changchun13(X)22,China)

Abstract:ne

analysisandsimulationoftheinfrared

effects

imagedegradationby

aero

transmission

were

presented.TheatmosphereMTFandthedegradationphysical

oftheimagewerepointedbyanalyzingthe

processofatmosphericturbulence,aerosolscatteringandabsorption.Thesimulationmethodsof

theatmosphericeffectswerepresented,andthecomputersimulationofatmosphericeffectswasdonefortheinfraredimaging

was

systemoflong—distance

horizontaltransmission

一l≈

on

condition

thatthetargetdistance

一2/3

km;theheightabovegroundwas2m;C。=x1.OxlO

showthattheinfrared

,opticalsystem

by

aero

is弘.0

and

f--600

mm.

Thesimulationresults

on

imagedegradationcaused

transmissionaremainly

focused

thelocalfuzzyandimageditheringinvariousdegreecausedbyatmosphereturbulence.

Infraredsimulation;Turbulence;

Keywords:AtmosphereMTF;

Scattering;Absorption

O引言

热成像系统用于观测远距离目标时,红外辐射在到达热成像系统的光学系统之前必须通过大气,并被

射源这些现象使远距离目标的热图像质量变差。研

究表明:肖红外成像系统用于远距离目标获取时,系统的整机性能受到系统自身MTF和大气湍流MTF

的很大限制|lI。因此.研究红外辐射的大气传输对像

质的影响,并对其进行建模和仿真,对红外成像系统的性能评估具有重要意义。

衰减。衰减的同时,大气梯度和湍流使辐射在结构上

发生畸变.引起空气折射率不均。此外,大气是固有辐

收稿日期:2009—10—10:修订日期:2(X)O—12一05

基金项目:国家自然科学基金资助项目(60977001);围家高技术研究发展计划(2007AAl2ZIlO)

作者简介:李咏豪(1977-)。男.山东招远人。I:程师.硕{:.主要研究方向为光电仿真‘j检测技术,Email:liyh@163.tom

万方数据

红外与激光工程第30卷

1红外成像中大气传输作用效果分析

l-l大气衰减

大气是一个包含湍流、散射和吸收粒子的随机介

质,大气对系统T作的影响基本体现在辐射衰减。红外辐射衰减与3种现象有关:大气气体分子的吸收;大气中分子、气溶胶和微粒的散射;冈气象条件(云、

雾、雨、雪)引起的衰减等。红外辐射的大气衰减满足布盖尔一朗伯定律,即光谱透过率r。(A)和衰减系数Ot。(A)之间满足如下关系:

下。(A)=exp【一d。(A)D】

(1)

式中:D为被观测目标和热成像系统间的距离。,

1.2大气湍流对成像的影响

红外辐射在大气中传播会造成图像模糊。原因有两个:一一是由于经气溶胶和其他大气成分散射后进人

人瞳的散射光:二是在光路中由湍流引起的介质折射

率变化。大气气溶胶散射造成的模糊通常被称为“邻接效果”。对于远距离的水平光路,由于接近地面的风和空气与地面之间的热交换造成温度的显著变化,使湍流成为影响图像质量的重要因素。即使存在气溶胶模糊.湍流也常常是模糊的决定性因素。湍流模糊常

被作为一个线性移不变物理过程来处理。湍流引起的系统性能下降可以用湍流MTF来描述。

Kopeika等提出r湍流MTF模型I创。该模型由古德曼提出的大气相十长度模型发展而来。在湍流

MTF模型中,大气的折射率变化用折射率结构参数t

表示。。湍流MTF分为长曝光模型和短曝光模型。长曝光模型集中反映r焦平面上模糊点位置随时间的变

化.这种变化是由传感器入瞳处的等晕波前倾斜引起

的。大气湍流的长曝光MTF是:

513

一l/3

M7FLE=exp(一57.40管eA

R)(2)式中:a为波形常数,当波形是球面时,a为3/8,当波

面是平面时,a为l;f为空间频率,c/mrad;C:为折射

率结构参数,m划3;A为波长;R为光程。大气湍流的短

曝光MTF是:

S13

—1/3l/3

M阿sE=exp(一57.4嘭eA

R)[1-Iz(Af/D)】(3)

式中:D为通光孔径;肛为无量纲系数,在近场和远场

万方数据

时取不同值.近场时取l,远场时取0.5。

长曝光记录的模糊光斑通常比短曝光的大,这是

因为造成模糊光斑的入射光角度随着光路的弯曲一直在变。使得光斑位置在像面上随机跳动,长曝光记录的是光斑跳动在时间上的积分。对于凝视型红外相机和红外监视型相机,其积分时间通常为毫秒级,典型值是

16.7ms,帧频60

Hz。所以。通常采用短曝光模型。

折射率结构参数C用下式来确定:

c。2=等

I,3

(4)

式中:<d‘>为波前入射角的变化量:D为光学孔径;尺

为物像间的距离。

1.3气溶胶对成像的影响

除了大气湍流,光线在大气传输过程中还受到分子和气溶胶产生的散射和吸收作用的影响。

散射作用大大扩展了光源的辐射范围,造成光线入射角度的H大差异.从而也造成像平面上图像的模糊。与短曝光湍流不同的是,由气溶胶造成的模糊是非常广泛的。散射光的角度受到成像系统视场的限制,并非所有的散射光都会出现在像面。此外,传感器的动态范围也限制了记录在像面的散射光角度.其影

响往往比视场的限制还大。如果在入瞳处由目标辐射

的散射光的辐射水平低于背景辐射,其对大气模糊的

贡献将是轻微的,通常町以忽略不计。接收器有限的

空间频率带宽限制了实际记录在像面的非散射光的最高频率。所有这些实际仪器的限制应适用于处理气

溶胶MTF的经典理论m41。

与可见光波段内不同的是。在热红外波段内吸收

是一个非常明显的效应。到达成像系统的入射光包含

非散射和散射两个部分。在长波红外大气窗口、中等气候条件下。长波红外的散射系数比可见光小得多,

光子通常只能与粒子发生一次作用,即要么被粒子吸

收。要么被散射15l。

因此,吸收只适用于非散射光的部分.而不适用

于所有辐射。所有与散射光有关的模糊,都不会因为微粒吸收而衰减。另外,能在无限空间频率带宽上增强图像的非散射光会被大气吸收而衰减,这使得图像

模糊作用更显著。因此,吸收造成气溶胶MTF在高空

间频率处衰减。气溶胶MTF可表示为:

第4期李咏豪等:大气传输引起的红外图像退化分析

:exp!、-s。Ri考”

w班i唧H刮1.f"v2忙卟≤K∞’

{exp卜SR}exp(-A。尺)}:VⅫ。

式中:S。和A。分别为大气有效的散射和吸收系数;尺

为路径长度;y为气溶胶MTF的角空间截止频率,

将湍流MTF和气溶胶MTF相乘。得到大气MTF:

MTF=M.M。

(6)

根据大气的热MTF分析,闩标与背景之间的温差通过有微粒的大气传播后将受到大气的调制.从而

使物体对比度下降、网像模糊、图像灰度值起伏及图

像形变、闪烁等忙引.

2大气退化作用仿真

2.1仿真方法

大气退化作用仿真器本身町以看作足时间函数F(・)和尺(・)组成的一个串联的系统.该系统作用于静

态输入图像l(x,Y),并创建・个输ff{图像序列O(u,',,t):

0(“,v,f)=JR(F(1(x,y)))

(7)

X和Y表爪输入图像的学标.U和v表示输出图像的坐标,输I叶I罔像序列的时间依赖性体现在t上。

第一个操作符F(・)用来模拟大气凋制传递厢数对图

像的模糊作用.定义为时间和夺问移不变的窄间滤波器

用高斯函数的标准差dw表爪模糊程度,

D(“.p,,)=(F(,(工,y);f)

(8)

第二个操作符尺(・)用来模拟|冬j像跳动,,这也是时间函数.且是卒间移不变操作符.O(u.V,f)=D(“(J,Y,,),v(x,Y,f),t(x,Y,f))=

R(1(x,y);d石(,),dy(,))=月(“工,y);,)(9)

这个操作处理输入图像的坐标为输入图像的每个位置计算出位移矢量.分别用出、dv分量表示横向

和纵向位移,图l、图2分别为滤波函数■维视l冬1和

滤波函数参数dw、血、dy爪意图.

存开始仿真之前.必须指定相机的特性、物距、大

每个像素。dw、血、dv的一一阶统计特性由前述公式给m.

万方数据

网l滤波函数一:维视网

网2滤波函数参数dw.出,dv爪崽网

Fig.1

VisualizationoftheFig.2Sketchof

parameters

dw.dx

filterfunction

dyofthefilterfunction

并由利用多重【日l归模型设计的随机数产乍器产生。

然后对输入图像使用操作符F(・)。需要计算㈣K

n个(m和n分别指输入网像在水平和垂直方向上的像素数)不同宽度的内核低通滤波。由此产幸的图像显爪出不同的局部模糊,,

接着对低通滤波后的图像运用操作符R(・)。这需

要汁算另一个肌x_,z的结果.根据每个输出罔像像素

的出和dv位移量.采用向后映射的方法找到对应的

输入图像中采样点的位置。向后映射能够将输入平面上各点随机分配到输出平面L,因此,可以精确模拟高动荡时的图像扭曲。使用双线性插值的方法考虑最近

的4个像素来计算该点的值,这种方法比最近邻法具

有更好的效果,仿真过程呵以用图3来表,】;,、

陶3仿真过稃小意陶

Fig.3Diagram

of

simulation

选用下列参数作为仿真对象:目标距离为lkm:}j标离地高度为2

m;C:为1.o×1.0—3m瑚(延光路方

向);光学系统为f/4.0,f=600mm,w=O.92。xO.73。;探

I.Lm:响应波根据公式(3)、(5)、(6)得出大气MTF,该传递函数

2.2仿真结果

测器探测元为320x256。像元尺寸为30段为7.7~9.5¨m,,

气情况等再将物理描述向仿真环境参数转换,,对于一次具体的仿真计算.需要为输m图像的每个像素计算

模糊程度和所期望的在指定动荡下的位置分布。对于

仅考虑大气传输的退化效果i,曲线图如图4所示。

592

红外与激光工程

第39卷

曲■

图4大气调制传递函数

Fig.4Atmosphere

modulationtransferfunction

对不同频率的条形靶进行了大气退化作用仿真。目的在于考察大气的独立退化作用。仿真结果如图5所示。经红外成像系统得到的实际红外图像还应包括系统自身及运动模糊、噪声等对图像的退化作用。

图5经大气传输退化后不同空间频率的条形靶仿真结果

Fig.5Simulatedbartargetsofdifferentfrequenciesdegraded

byatmospheretransmission

为了比较仿真结果与原始图像的差别,将仿真结果与原始图像做互相关运算,求出其互相关函数。为了简化计算,仅考虑一维的情况。

rfg(m)=E㈣咖训}=熙击蝥咖(棚)

(一M≤m≤M)

(10)

由于r缸与rf之间满足如下关系:

rfg(m)=rf(m)’h(m)

(11)

对公式(1)两端取傅里叶变换得:

Pfg∞)=Pf啦)‘H∞)

(12)

由上述关系式可求出退化过程的传递函数。图6

为在50lp/mm处的归一化互相关函数。

分别计算各频率处的互相关函数.并根据公式(12)求出退化过程的传递函数。图7为退化过程的传递函数与大气调制传递函数比较的结果。

万方数据

Spatialfreqtmcy/Ip・mm‘’

图650lp/mm处归一化互相关函数

Fig.6Normalized

cross-correlation

function

at

50lp/mm

函一

Spatial

frequncy/Ip’蛐‘

图7退化过程的传递甬数

Fig.7Transferfunctionofdegradation

从图7中可以看出:由于大气调制传递函数具有

大气传输引起的红外图像退化是个复杂的物理(1:转第626页)

较高的数值,大气传输造成的图像退化主要体现在湍流引起的介质折射率变化,造成调制传递函数的局部波动,使图像显示出不同程度的局部模糊及图像跳动,频率越高退化效果越明显。随着频率的升高,大气调制传递函数引起的边缘模糊变得更加明显,增加了对远距离目标探测和识别的难度。

3结论

过程。在远距离目标获取中,红外侦查相机的性能受到大气传输的限制。文中主要针对大气调制传递函数

和图像跳动.利用随机数产生器产生具有特定统计特性的串联操作符的3个参数dw、出、dy。对不同空间频率的条形靶进行了大气传输效果仿真,并将仿真结

果与原始图像进行了互相关分析,证明退化过程的传递函数与理论结果相符。算法主要针对的远距离水平传输方向的红外成像系统大气传输退化也可嵌入到其他成像仿真的图像处理过程中,与系统自身及运动模糊、噪声等对图像的退化作用进行联合仿真。大气

626红外与激光工程第39卷

法,根据电场作用下液滴面形变化的原理,直接在半导体激光管封装窗口上制作非球面透镜来准直LD光束。实验制得具有两个聚焦的非球面透镜,对LD光束

相互垂直方向的发散角分别进行准直。效果良好。

integrationof

lens

edge-emittinglaserdiodewith

ion

micro-elliptical

using

focused

on

beamdirect

deposition【J】.IEEE

Transactions(】):2—8.

Semiconductor

Manufacturing,2002,15

【8】LITong-hal,WUGuo-jun,HU

Bao-wen,etand

a1.Fabricationof

in

laser

参考文献:

FoR

thehemi-ellipticalmicrolensesbeam

application

diode

shaping【J】.ScienceTechnologyandEngineering(李同

KNERJF.KUN’lZ

W.Characteristicsofemcient

海,吴围俊,胡保文.等.椭圆微透镜的制备及其在LD光束整形中的应用.科学技术与工程),2005,24:1890一1892.

【9】GOUZhi-yong,WANGJiang,WANGLei,eta1.Diffraction・

limitedaspheric

collimation

lens【J】.Journal

of

Applied

laserdiode

collimators[C]//SPm,1997。740:145-151.

【2】

CHENGICe-hu扎CHENGDuo-si.Researchofthecollimationapplications

for

selfoc

lenseslJl.ActaPhotonicaSinica(程科

华,程多思.自聚焦透镜准直应用的长度计算与实验.光子学报),1992,21(2):145—151.

【31

XU

Xiao-fang,LILi-ua,WUJin-hui,eta1.Fibercoupling

of

the

Optics(勾志勇,乇江,王磊,等.衍射极限非球面准直透镜.

应用光学),2006,27(6):467—469.

【101

HE

Jian—hui,ZHAN

on

Zhen—xian,WU

Yun—liang,et

a1.

module

hi曲outputpower【J】.Infrared

andLaserMeasurement

surfaceprofileofasphericalliquidmicro—lens

Engineering(许孝芳,李丽娜,吴金辉,等.高功率半导体激光器列阵光纤耦合模块.红外与激光工程),2006,35(1):86—88.

【4】

BOKORN,DAVIDSONN.Beamshaping

use

【11.InfraredandLaserEngineering(何建慧,詹珍贤,吴云良.等.非球面液滴微透镜的轮廓测量方法.红外与激光工

程),2008,37(4):330—333.

【ll】ZHANG

with

diffuselightby

of

singlereflection【J】AppliedOptics,2001,40:2132.

Qiang,SHIYong-chao,LV

Nai-guang,et

a1.

【5】

SMITH

WJ.ModernOpticalEngineering[MI.NewYork:

Detectionofthe

positionoflaser—spotcenter【J1.Journalof

of

McGraw—Hill.2000.

【6】

ZHAN

Zhen-xian,HE

Beijing

Jian-hui,YAO

Hai—tao,et

a1.

Institute

MachineryIndustry(张强,施涌潮,吕

乃光.等.激光光斑中心位置的探测.北京机械工业学院学

报),1996,11(2):24—29.

Aspheficalliquidmicro—lensmanipulatedbyelectrostaticforces

【J】.ActaOpticaSinica(詹珍贤,何建慧,姚海涛,等.电场作用下的变焦非球面液滴微透镜.光学学报),2008,28(2):

361-364.

【7】

FU

Yong-qi,NGOI

K,BRYAN

A.A

【12】TAN

and

Zuo-jun,XUESong,KANGJing-ran,eta1.Collimation

ofsemiconductorlasers

in

test

laserfuze【J1.Journal

of

Applied

novel

one

Optics(谭佐军,薛松,康竞然,等.激光引信巾半导

step

体激光器的准直及其测试.应用光学),2007,28(4):454—457.

(上接第592页)

(3):354-358.

传输的非等晕效应仿真还有待进一步研究。

参考文献:

KRAPELS

K,DRIGGERSRG,VOLLMERHAUSENRH.ct

turbulence

【4】

SADOTDAN,DVIRA,BERGELLimagesdistortedtheoretical

by

eta1.Restorationofthermal

on

theatmosphere,based

modulation

measuredand

atmospheric

transfer

function【J】.

OpticalEngineering,1994,33:44—53.

【5】

SADOTD,KOPEIKANS.Imagingthroughtheatmosphere:practical

instrumentation・based

theory

and

verification

of

a1.Atmospheric

modulationtransferfunctionfor

infraredtargetacquisitionmodeling【J1.OptEng,2001,加(9):

1906—1913.

SADOTD,KOPEIKA

aerosolMTF【Cl//SPIE,1992,1688:59—64.

【61

S。ROTMANSR.Incorporationof

modeling[C]//

【7】

SADOTD,KITRONG,KITRONN,eta1.Thermalimagingthrough

the

atmosphere:Atmospheric

MTF

theory

and

atmosphericblurringeffectsintargetacquisitionSPIE,1994,2224:95—107.

verification【C]//SPIE,1993,1968:59—66.GAO

Si—feng,WUPing,HEMan—li,eta1.Estimation

on

the

ZHU

Wen-yue,MAXiao-shan,RAORui-zhong,eta1.Optical

effects

on

operating

rangeofinfraredsystemundercomplexatmosphere

turbulence

electro-optical

sensors[J】.Infraredand

condition[J].InfraredandLaserEngineering(高思峰,吴平.何曼丽,等.复杂大气条件下红外系统作用距离的估算.红外与激光工程),2008,37(6):941—944.

La赞r

Engineering(朱文越.马晓珊.饶瑞中.等.大气光学

湍流对光电探测器性能的影响.红外与激光工程),2006,3S

万方数据

大气传输引起的红外图像退化分析

作者:作者单位:刊名:英文刊名:年,卷(期):

李咏豪, 姜会林, 周娜, LI Yong-hao, JIANG Hui-lin, ZHOU Na

李咏豪,LI Yong-hao(北京特种车辆研究所,北京,100072), 姜会林,周娜,JIANG Hui-lin,ZHOU Na(长春理工大学光电工程学院,吉林,长春,130022)红外与激光工程

INFRARED AND LASER ENGINEERING2010,39(4)

参考文献(7条)

1. 高思峰;吴平;何曼丽 复杂大气条件下红外系统作用距离的估算[期刊论文]-红外与激光工程 2008(06)2. SADOT D;KITRON G;KITRON N Thermal imaging through the atmosphere:Atmospheric MrrF theory andverification 1993

3. SADOT D;KOPEIKA N S Imaging through the atmosphere:practical instrumentation-based theory andverification of aerosol MTF 1992

4. SADOT DAN;DVIR A;BERGEL L Restoration of thermal images distorted by the atmosphere,based onmeasured and theoretical atmospheric modulation transfer function 1994

5. 朱文越;马晓珊;饶瑞中 大气光学湍流对光电探测器性能的影响[期刊论文]-红外与激光工程 2006(03)

6. SADOT D;KOPEIKA N S;ROTMAN S R Incorporation of atmospheric blurring effects in target acquisitionmodeling 1994

7. KRAPELS K;DRIGGERS R G;VOLLMERHAUSEN R H Atmospheric turbulence modulation transfer function forinfrared target acquisition modeling[外文期刊] 2001(09)

本文链接:http://d.g.wanfangdata.com.cn/Periodical_hwyjggc201004002.aspx


相关内容

  • 遥感填空题
  • 1. 可见光遥感的的探测波段在之间,红外遥感的探测波段在0.76--1000微米之间. 2. 当大气中粒子的直径比波长小得多时发生的散射是 瑞丽散射,当大气中粒子的直径比波长大得多时发生的散射是 无选择性散射. 3. 已知某河流的宽度为20m ,在像片上量得的宽度为0.5cm ,则像片的比例尺,如果 ...

  • 遥感课后答案
  • 遥感的基本概念是什么?应用探测仪器,不与探测目标相接触,从远处把目标的电磁波特性记录下来,通过分析,揭示出物体的特征性质及其变化的综合性探测技术. 遥感探测系统包括哪几个部分?被侧目标的信息特征.信息的获取.信息的传输与记录.信息的处理和信息的应用. 作为对地观测系统,遥感与常规手段相比有什么特点? ...

  • 光电子技术课后答案第三版
  • 习 题1 1. 设在半径为Rc的圆盘中心法线上,距盘圆中心为l0处有一个辐射强度为Ie的点源S,如图所示.试计算该点源发射到盘圆的辐射功率. Ie ded 第1题图 , d πRcl 2 2 deIedIe πRcl 20 2 2. 如图所示,设小面源的面积为As,辐射亮度为Le ...

  • 遥感整理最新版
  • 第一章 绪论 1.遥感定义?遥感分类?(遥感平台.电磁波段.辐射源(工作方式)?) 遥感即遥远感知,是在不直接接触的情况下,对目标或自然现象远距离探测和感知的一种技 术. 遥感分类有不同的标准.1. 按遥感平台 (1)航天遥感 (2)航空遥感 (3)地面遥感 2 .按辐射源分类(1)被动遥感 (2) ...

  • 光电子技术课后答案
  • 习 题1 1. 设在半径为Rc的圆盘中心法线上,距盘圆中心为l0处有一个辐射强度为Ie的点源S,如图所示.试计算该点源发射到盘圆的辐射功率. deπRc2Ie, d2 dl0第1题图 πRc2deIedIe2 l0 2. 如图所示,设小面源的面积为As,辐射亮度为Le,面源法线与 ...

  • 遥感概论复习题
  • 遥感概论复习题 第一章 一.填空: 1. 遥感的分类方法很多,按遥感平台分:地面遥感. 航空遥感 .航天遥感 .航宇遥感. 2. 遥感的分类方法很多,按工作方式分:主动遥感和 被动遥感.成像遥感与非成像遥感. 二.简答及综合题 1.何谓遥感?遥感技术系统主要包括哪几部分? 遥感:广义上指一切无接触的 ...

  • 遥感总复习题库(含答案)
  • 第一章 电磁波及遥感物理基础 名词解释: 1. 电磁波 (变化的电场能够在其周围引起变化的磁场,这一变化的磁场又在较远的区域内引起新的变化电场,并在更远的区域内引起新的变化磁场.) 变化电场和磁场的交替产生,以有限的速度由近及远在空间内传播的过程称为电磁波. 2. 电磁波谱 电磁波在真空中传播的波长 ...

  • 遥感技术重点
  • 第一章 绪论 一遥感 概念:即不直接接触物体本身,从远处通过仪器(传感器)探测和接收来自目标物体的信息(如电场.磁场:电磁波.地震波等),经过信息的传输及其处理分析,来识别物体的属性及其分布等特征. 二 遥感技术 含义:从不同的高度平台上,使用各种传感器,接收来自地球表层各类地物的各种电磁波信息,并 ...

  • 一种计算红外辐射大气透过率的数学模型
  • 〈系统与设计〉 一种计算红外辐射大气透过率的数学模型 周国辉,刘湘伟,徐记伟 (解放军电子工程学院,安徽 合肥 230037) 摘要:从分析红外辐射在大气中传播的影响因素出发,分别给出大气吸收.大气散射和气象衰减的衰减系数,并计算出倾斜路程的海平面水平等效路程,最后给出一种计算红外大气透过率的数学模 ...