u_8 22Z @&##c6\$ 使用FRED的图形用户界面和它的内置脚本语言，我们可以轻松地实现热辐射和
成像。尽管强力的光线追迹同样是可能的，FRED使用了应用标准
光学工程算法的高效运算器来实现热成像和辐射计算。使用源自辐射度量学的技术，用FRED追迹必要数量光线的可能需要的时间，我们可以高效并精确地完成热成像、冷反射、杂散光、热
照明均匀性和热自发辐射的计算。
A5ID I<a jt--w"|-r 1. 热辐射和热成像是什么？ o7XRa] O yZ$;O0f&& 热成像定义为产生一个场景的可视化二维图像的过程，该图像依赖于从场景到达成像仪器孔径的热辐射或红外辐射的差异。热成像系统通常会减去背景来增强在红外场景中变化的对比度。当背景不均匀时，由于冷反射的存在，可能产生杂散信号。对于国防和安全问题尤为重要，在其中我们可以发现具有不同热温度或辐射率的物体，此时可以从图像场景的剩余部分区分出它们。对于这个问题的主要应用是：探测、分类和追迹隐藏在个人身上、包裹中、车辆上或船运集装箱中的武器、人员、车辆、物品和材料。图1是一个非常好的案例，当在FRED中进行仿真时，一个日常用品：茶壶，通过一个具有热探测仪的摄像头成像。
sLp
LY1X YO0x68 热辐射是从一个光学仪器周围的环境或结构中发出的能量，它会引起杂散光问题。冷反射是一个热辐射问题，由于反射到探测器上的辐射，在一个红外系统中的热辐射表现为在一个显示图像中的黑色圆形区域。
|l(lrJ{ h_yR$H&tX 通常，这些系统通过探测叠加在大的背景上的小信号工作。在室温下，黑体辐射曲线的峰值大致在10μm处。因此世界在这个波长处“发光”，发光的微小变化表明了温度或辐射率的变化。特别的，当一个冷却的探测器图像反映了自身，那么就会产生一个局部背景的缺失。这通常表现为在图像中央的黑点。有人可能称之为“杂散黑”，而不是杂散光。
nF$)F?|| O Zm[i H 在测量绝对辐射而不是相对信号的红外辐射仪中，任何背景辐射是不可接受的。在这样一个仪器中，冷却整个仪器到低温度来消除由于自辐射导致的杂散光是必须的。
K<JP9t6Qd 
ss8v4@C 图1 此图演示了一个在茶壶表面具有不同辐射系数和温度分布的暖茶壶的简单问题。然后茶壶通过单透镜成像，探测器放置在单透镜后面。这种机械结构辐射到了探测器上类型的问题，可以在许多引起热问题的红外系统中找到。
_PXo'*j 7=A9E] : 最近，已经开始致力于红外遥感应用中，包括温度的测量和绘图、森林火势的感知和控制、监督和多
光谱地表成像等。
>8O=^7 N-YZ0/c 这些应用种有许多是经过长距离完成的，透过大气，在大气中IR能量的吸收是这些系统性能的一个影响因素。军事的和基于空间的应用一般来说可以通过探测器处理，探测器的工作波长落在8.0-15微米之间，在这个波段内大气的吸收是最小的。其他的应用的波带较宽，为0.9-300微米。
1>y=i +T/b G5J ZB7C 2. FRED如何进行光线追迹和显现热辐射和成像？ '|N4fbZd Bjj^!T/# 在FRED中追迹热辐射有几个方法。第一种方法是创建一个
光源，然后在
光学系统中对它进行强力光线追迹。第二种方法是通过光学系统从探测器后端进行光线追迹，这需要较少的光线。在两种方法中，能够显现二维和三维图的热成像是非常重要的。
D D;+& fe `Kw8rG\]: 事实上这里有两个问题：计算时间和精度。在一个复杂系统中，如果分析者尝试去了解设计中递增量的影响，并且想要实时这样做，光线追迹时间可能会特别长。反向光线追迹能够使计算变得几乎是交互式的。另外，由于幂指数收敛比均匀性快，即使每个微分区域只有几束光线到达热源，分析者几乎可以肯定保证精确的结果。
D@W3;T^ 
%E7.$Gj% 图2 两个辐照度计算的比较：一个是使用前向光线追迹，另一个是使用后向光线追迹。后者需要比前者少1/53倍的光线达到同样的准确度。
]u ~Fn2 4$GRCq5N; 热自发辐射可以简单描述：每个光学和机械结构像朗伯辐射体一样辐射能量，辐射出的能量是自身温度和辐射系数的函数。通过追迹光线来模拟发射出的能量；在传播过程中发生透射和反射时，它们遵循几何光学的规律。这些光线（它们表示热能）累积到了FPA上。
[<a%\:c m4 0a 6z "K} 根据这个计算，大多数软件会让用户设置物体的温度和辐射系数。从统计学的角度来看，这样做完全是错误的！在大多数“真实”的系统中，相对于光学和机械元件，FPA朝向一个非常小的固体角，因此当大量光线被追迹时，如果有的话，到达FPA的光线很少（图3）。结果是热自发辐射的错误的估计。
C,vc
aC? N: jiZ) 有一种更加高效但是非最佳的方法。直接或者间接，大多数软件允许用户指定一个优先的辐射方向；这在文献中被称为“重点采样”。使用该方法，用户为每一个光学和机械元件指定了一个重点采样。在光线追迹的过程中，光线散射到了这些优先的方向，这在引导光线到达FPA方面十分高效（图4）。这极大的改善了统计结果，并且产生了热自发辐射的一个准确评估。然而，在一个具有许多结构性元件的复杂、“真实生活”的系统，这样做单调乏味且耗时，特别是如果分析者需要设定温度和辐射系数。
.r %|RWs6W >1pH 91c' 图3 热发射主机筒的强力光线追迹。在这个例子中，入射光线没有到达FPA。
.x?zky^ 图4 发射朝向一个重要方向，计算更加高效，但是设置麻烦。
g4SYG)'R+ Y6? mY! 处理这个计算最好的方法是采用热辐射度量学的数学方法。计算热自发辐射（TSE）基本上是该形式的总和
[HiTR !o* QN&^LaB<T 
vP)~j1 因为热辐射是朗伯的，我们可以用由热出射度推导出的等值表达式替换L
*Q 120 R 4?M3#],'h 
:f Kz^@mY4 其中ε是辐射系数，f是黑体积分，σ是斯特藩-玻尔兹曼常数，T是温度（K）。问题是如何高效的计算Aobject和Ωdetector。
h]DE Cd{ #]a51Vss 反转光源和收集器的作用（使用方程3），我们可以写成
B%:9 P aC
Lg~g4 
jTUf4& b- "M0l; 注意到探测器的面积Adetector是一个定量，为了高效的计算每个发射物体的立体角，我们在辐射方程中引用了对称性概念。
#L=
eK8^e 53*, f 
15T[J%7f 从方程4中，我们注意到如果我们从探测器（L=1/πAdetector）发射光线，入射到物体上的功率数值上将等于它的投影立体角（回想投影立体角等于Ω/π）。发射的探测器功率因此等于
N z3%}6F: o``>sB ZOq 
)^>XZ*eK （事实上只有探测器辐射到了一个圆锥体中，这个公式才是正确的。如果我们希望辐射到一个矩形体中，正确的探测器功率是4/π因子）
!v4j` A;% AD%D ,l 知道如何高效的计算投影立体角是该方法的核心；这是第二个“高超的技巧”。由于立体角是常数，它们只需要计算一次。
{%~4RZA 
4'hcHdL9 图5 从FPA反向光线追迹来获得精确的投影立体角。在大多数例子中，杜尔窗口是相对于探测器最靠近的通光孔径，因此我们一般将光线发射到它的立体角内。热自发辐射的计算由方程6获得，这很容易的应用到一个表格中（图6）。出于完整性，热自发辐射的方程由下式给出：
)Z]y.W ) J[Yg]6 
`CEj 4 其中表达式（Ω/π）是由光线追迹计算出的投影立体角。
图6 计算热自发辐射使用的表格，如Exceltm。列“Incident Power”实际上是投射立体角，列“Contribution”实现了用方程6实现。
=<zSF\Zr_ h 4.=sbzZ 热自发辐射的计算使用反向光线追迹的好处有很多，包括
U;Ll.BFP 1. 立体角的精度由从探测器追迹的光线数目决定
D<3V#Opw 2. 为了“切合实际”，温度和辐射系数在表格中可以很容易的改变。
chMc(.cN0 3. 这里没有冗长的设置时间，设置时间不再是模型复杂性的函数。
^N2M/B|0 4. 通过把探测器划分为微分区域，光线追迹并将每个区域求和，我们可以获得在探测器上热自发辐射的图像。
Z*9]:dG:! 9C)3
b3 现代光学软件十分强大且适应性广，通过点击工具条上的按钮，为有经验和没有经验的用户提供了来完成许多复杂的任务的方法。然而，有许多地方用户需要以软件预定义的方式执行计算。应用“高超的技巧”常常是唯一的办法，在合理的期限内执行这些运算并到达所需要的精确度水平。
d J%Rk#?;A &< |-> *v 3. FRED对于热分析如何创建几何结构？ 5 -WRv; 1H,tP|s 我们可以在FRED中使用图形用户界面直接创建系统几何结构，或者从IGES或STEP CAD格式和
光学设计程序导入。程序有很多选项来创建表面，包括标准的平面、圆锥面、圆柱面、椭球面、双曲面、圆环面、多项式表面、Zernike、Nurb、Meshed、旋转曲面、挤压面、复合曲线、样条曲线和用户定义的表面。
V>j hGf 11y .z^ 由于FRED具有多文件用户界面，可以在文件之间剪切、复制和粘贴元件。实体可根据逻辑关系安排到对应于系统物理布局的装配体、组件和元件等的分类中。任何表面可能会被隐式表面或者下面定义的孔径集曲线裁剪。
6^IqSNn- X})Imk7&E 任何表面可以根据温度和辐射系数定义。另外，把每个表面定义为热辐射源来传播热能和并创建热成像图。
w Al}:|+n =i^<a7M~ 4. FRED如何记录热辐射路径？ e_~ fJ ^ ?7dOW FRED具有进行高级光线追迹的能力，并且记录所有在系统中追迹光线的所有路径，图7表示的是一个茶壶热成像的光线路径的列表，这是使用从探测器通过透镜然后到图1中所示茶壶的后向光线追迹形成的。
D@\;@(
| =1
BNCKT< 该光线历史是每个路径功率、沿着该路径的光线数、它们如何到达最后一个实体和它们通过了多少表面（事件计数）的一个完整报告。使用任何光线追迹路径并把它们复制到用户定义的路径列表中（选择该路径，点击鼠标右键，选择复制该路径到用户定义路径列表的选项）是可能的。该路径现在将会作为在高级光线追迹的可选路径和可使用的光线方法之一出现。然后在仅有的路径上，绘制点列图或辐照度分布函数是可行的。[/td][/tr][/table][/td][/tr][tr][td]
H&=n:'k^ z'oiyXEE3 }uI7 \\ S 通过使用该方法，可以查看相对于信号路径，每个鬼像、直射、单个或多个散射路径贡献了多少功率。
Z 5)_B,E:X 
:IVk_[s 图7 如图1所示的茶壶案例的光线路径。注意到有三条路径到达了茶壶主体（.teapot7.Body.B-Spline Surface 4），它们是该系统中的路径1,2和8，显示为在第2行到最后一行。路径数2是从探测器通过透镜然后传播到茶壶本体的直线路径。注意一下这两个路径每个的功率差异，第1条路径的0.25448和第2条路径的0.000227487。第2条路径是该路径的散射部分；第1条路径是直接热贡献。还可以查看该路径的定义；示于图9.
图8 通过透镜追迹的第二条路径
Gkv~e?Kc~^ 图9 示于图8中路径的光线路径信息 Gl8&Fr R
7{An@hNh 5. FRED如何显示辐照度或热辐射图？ =0PRAc 默认条件下，FRED使用3图颜色面板显示辐照度或热成像图。左上图是一个等距的伪彩色图，显示了在选择的分析面上单位面积的功率。右边有一个比例尺显示了图中的功率等级。右上角和右下角的面板是左上角面板的横截面。点击左上图的任何部位会在水平和垂直方向创建三维图的一个横截面；该位置处的位置和辐照度信息显示在左上角面板的左下角处。茶壶系统的辐照度图如图10所示。
?)'
2l6 图10 茶壶系统的辐照度图
{ 8f+h 只考虑一个面板，在那个面板上双击鼠标左键可以放大该图。上面板正处于透视图模式，为了创建一个二维图，点击鼠标右键弹出绘图选项，选择perspective view来创建一个如图11所示的二维图。
"7 yNKO;W 图11茶壶的热影像二维图模式
t `Y!"l 6. FRED在热应用中如何定义散射表面？ 
ran^te^Ks( 图13重要性采样的菜单选项应用到一个特定表面。如图所示，多个散射属性可以应用于多个重要性采样目标。注意到在该图中指定了一个米散射散射体属性和一个Harvey Shack抛光表面散射体，另外定义了两个重要性采样目标，一个是指定了一个方向，另一个朝向一个表面-焦平面。
GjeUUmr 7. FRED如何记录热散射路径 .5Knb c 在高级光线追迹完成后，当选择保持光线历史选项，FRED就具有了一个可得到的杂散光报告。然后我们就可以从工具菜单请求一个详述鬼像和散射路径如何到达任何表面的杂散光报告。图14显示了具有创建/使用光线历史文件选项和确定选择的光线路径选项的高级光线追迹对话框。图15显示了一个简单Cassegrain望远镜的杂散光报告，该报告详述了从一个光源以离轴5度角进入
望远镜的杂散光。从报告中我们可以看到，在离轴5度角处的最大的杂散光光源是偏离第二个支架的侧面，因为该路径的功率的73%到达了探测器。总功率如第4列所示，表明只有0.00197290w从该路径到达了探测器。
a4T~\\,dZ> 图14 具有创建/使用光线历史文件选项和确定选择的光线路径选项的高级光线追迹对话框
V@v1a@=W ,'C30 A*p 图15 显示了杂散光报告电子表格，可以用于追迹任何级次的散射和鬼像路径。所示的报告部分显示了详述每个报告路径的路径数、光线数、功率百分比和总功率。
ss`P QN 8. FRED如何在探测器上穿过多个像素点迭代来创建热图？ I%9bPQ FRED具有一个内置的可编译的Basic脚本语言。从Visual Basic脚本语言里，几乎所有用户图形界面（GUI）命令是可用这里的。FRED同样具有自动的客户端和服务器能力，它可以被调用和并调用其他可启动程序，如Excel。因此可以在探测器像素点上定义多个离轴光源，及在FRED Basic脚本语言里的For Next loops语句沿着探测器像素点向上和向下扫描来反向追迹光线。注意到在图16中写了功率信息到数据文件中，这样可以使用三维图表查看器调用和查看数据。
xEVLE,*?> 图16 FRED Basic脚本语言创建茶壶的热成像
rps(Jos _~ 'script for calculating thermal image map'edited rnp 4 november 2005
os1?6 z~ 'declarationsDim op As T_OPERATION
WDE e$k4. Dim trm As T_TRIMVOLUMEDim irrad(22,22) As Double 'make consistent with sampling
% Cu.u)/+ Dim temp As DoubleDim emiss As Double
JAlU%n?R 'setupnx = 21
!8Z2X!$m{< ny = 21numRays = 1000
6X7s 4 minWave = 7 'micronsmaxWave = 11 'microns
-]Q(~'a sigma = 5.67e-14 'watts/mm^2/deg k^4Print ""
n$XdSh/ Print "THERMAL IMAGE CALCULATION"detnode = FindFullName( "Geometry.detector.Surf 1" )
IX< r5!
Print "found detector array at node " & detnodesrcnode = FindFullName( "Optical Sources.DiffDetArea" )
_)KY Print "found differential detector area at node " & srcnodeGetTrimVolume detnode, trm
dU]i -NF detx = trm.xSemiApedety = trm.ySemiApe
aqs%m ( area = 4 * detx * detyPrint "detector array semiaperture dimensions are " & detx & " by " & dety
v<;: 0 Print "sampling is " & nx & " by " & ny'reset differential detector area dimensions to be consistent with sampling
JBR[;
zM pixelx = 2 * detx / nxpixely = 2 * dety / ny
* me,(C SetSourcePosGridRandom srcnode, pixelx / 2, pixely / 2, numRays, FalsePrint "resetting source dimensions to " & pixelx / 2 & " by " & pixely / 2
l#D-q/k? 'zero out irradiance arrayFor i = 0 To ny - 1
J F M"ii{8 For j = 0 To nx - 1irrad(i,j) = 0.0
9<|m4 Next jNext i
Ys-Keyg 'main loopEnableTextPrinting( False )
e(yQKwVD ypos = dety + pixely / 2For i = 0 To ny - 1
z,{e]MB)M xpos = -detx - pixelx / 2ypos = ypos - pixely
# (%t*"IY; EnableTextPrinting( True )Print i
~{L.f94N EnableTextPrinting( False )For j = 0 To nx - 1
'h-3V8m^e xpos = xpos + pixelx'shift source
m|pTn#*` LockOperationUpdates srcnode, trueGetOperation srcnode, 1, op
CQdBf 3q op.val1 = xposop.val2 = ypos
KNg8HYFW\ SetOperation srcnode, 1, opLockOperationUpdates srcnode, false
}B y)y;~ 'ray traceCreateSource srcnode
aje^Z=] TraceExisting 'draw'radiometry
6* ZU}xT For k = 0 To GetEntityCount()-1If IsSurface( k ) Then
Fr-[UZ~V temp = AuxDataGetData( k, "temperature" )emiss = AuxDataGetData( k, "emissivity" )
U~aWG\h#X If ( temp <> 0 And emiss <> 0 ) ThenProjSolidAngleByPi = GetSurfIncidentPower( k )
[tUv*jw % frac = BlackBodyFractionalEnergy ( minWave, maxWave, temp )irrad(i,j) = irrad(i,j) + frac * emiss * sigma * temp^4 *
- $U@By<SJ ProjSolidAngleByPiEnd If
mZwi7s&u End IfNext k
&Hv;< Next jNext i
/ZHO>LNN| EnableTextPrinting( True )'write out file
D" L|"qJ Open "C:\Documents and Settings\teapotimage.dat" For Output As #1Print #1, "GRID " & nx & " " & ny
xE%O:a?S Print #1, "1e+308"Print #1, pixelx & " " & pixely
% drJ p6n% Print #1, -detx+pixelx/2 & " " & -dety+pixely/2maxRow = nx - 1
}wG,BB %N maxCol = ny - 1For rowNum = 0 To maxRow ' begin loop over rows (constant X)
A j, ]n>{ row = ""For colNum = maxCol To 0 Step -1 ' begin loop over columns (constant Y)
eY
T8$ row = row & irrad(colNum,rowNum) & " " ' append column data to row stringNext colNum ' end loop over columns
mA& =q_gS Print #1, rowNext rowNum ' end loop over rows
+%P t_ Close #1Print "all done!!"
