本发明公开了一种能够区分不同粗糙面的被动毫米波辐射模拟方法。现有被动毫米波辐射模拟方法只适用于深粗糙面，而无法适用于中等粗糙面或是微粗糙面。本发明的步骤如下：一、建立被测场景的三维模型。二、对步骤一中建立的三维模型进行剖分。三、定义a个射线组，一个射线组包括b条射线。四、计算沿a×b条射线到达辐射计的亮温，计算中考虑第二层反射。五、将所得的a×b个到达辐射计天线的亮温，转化为a×b分辨率的亮温成像模拟图。本发明在进一步精确地表征了粗糙面的被动辐射模型，使得本发明获取的亮温成像模拟图和实际探测情况愈加接近。本发明能够用于与辐射计拍摄出的亮温图像进行对比，从而判断辐射计的品质好坏。
1.一种能够区分不同粗糙面的被动毫米波辐射模拟方法，其特征在于：步骤一、建立被测场景的三维模型；三维模型中包含地面、辐射计安装点和多个目标；
步骤三、定义a个射线组，一个射线组包括b条射线；所有射线均以三维模型中的辐射计安装点为端点，且经过三维模型的地面；同一射线组内射线的俯仰角相等，且方位角依次等距增大；a个射线组内射线的俯仰角依次等距增大；a×b的数量即为所得图像的分辨率；
步骤五、定义第x个射线组内第y条射线的俯仰角为θ0，方位角为 第x个射线组内第y条射线经过的目标为第一目标；第x个射线组内第y条射线射入第一目标的点作为特征反射点；第一目标上含有特征反射点的侧面即为特征侧面；
将沿着第x个射线组内第y条射线到达辐射计天线的毫米波作为入射波；入射波由第一反射波、第一目标辐射波和空气辐射波组成；第一目标辐射波为第一目标的特征反射点自身辐射向辐射计安装点的毫米波；空气辐射波为空气背景在特征反射点到辐射计安装点方向上辐射出的毫米波；在特征反射点反射出第一反射波的毫米波作为第二反射波；第二反射波与地面的交点为地面反射点；
式(1)中，T11(θ2,φ2)的表达式如式(2)所示；r1(θ2)为第一目标的特征反射点在入射波方向上的反射率；T1为第一目标的温度；e1(θ2)为第一目标的特征反射点在在入射波方向上的发射率； 为空气辐射波传达的亮温；
式(2)中， 为空气背景向下到达地面的亮温；p为一种极化方式；rp(θ1)为地面反射点在第二反射波方向上的p极化的菲涅尔反射率；e为自然对数的底；k为波数的均方根高度；σ为地面的均方根高度；θ1为第二反射波在地面反射点的反射角； 的表达式如式(3)所示；Td为地面的温度；edp(θ1)为地面反射点在第二反射波方向上的p极化发射率，edp(θ1)的表达式如式(6)所示； 为地面与第一目标的特征侧面之间的极化旋转矩阵；
式(3)中，定义N个特征反射波组，一个特征反射波组包括M条特征反射波；所有特征反射波均射入地面反射点；同一特征反射波组内特征反射波的俯仰角相等，且方位角依次等距增大；N个特征反射波组内特征反射波的俯仰角依次等距增大；θi′为第i个特征反射波组内特征反射波的俯仰角； 为同一特征反射波组内第j个特征反射波的方位角；
温； 为第i个特征反射波组内第j个特征反射波射入地面反射点的反射率，其表达式如式(5)所示；
若第i个特征反射波组内第j个特征反射波来自一个目标，则以该特征反射波作为二级特征反射波；以该目标作为第二目标，将θi′赋值给θs，将 赋值给 代入式(4)计算，得到否则， 的值等于式(2)中的的表达式如下：
式(4)中， 为空气背景向下到达第二目标的亮温；r2(θ3)为第二目标反射出二级特征反射波的点在二级特征反射波方向上的反射率；T2为第二目标的物理温度；e2(θ3)为第二目标反射出二级特征反射波的点在二级特征反射波方向上的发射率； 第二目标上反射出二级特征反射波的侧面与地面之间的极化旋转矩阵；
式(5)中，q为与p不同的极化方式； 为第二层在方向 上的射线到第一层方向为 的射线的同极化散射系数； 为为第二层在方向 上的射线到第一层方向为 的射线的交叉极化散射系数；
式(6)中， 分别为相干反射系数、非相干反射系数； 的取值等于式(2)中的rp(θ1)； 的表达式如式(7)所示；
为第一层在方向(θ1)上的射线到第二层方向为 的射线的同极化散射系数； 为第一层在方向(θ1)上的射线到第二层方向为 的射线的交叉极化散射系数；
将式将式(2),(3),(4),(5),(6),(7)代入式(1)，得到入射波传达的亮温 之后进入步骤六；
步骤六，若y＜b，则将y增大1，并重复执行步骤五；若y＝b且x＜a，则将1赋值给y，并将x增大1，并重复执行步骤五；若y＝b且x＝a，则进入步骤七；
步骤七、将步骤四、五和六所得的a×b个到达辐射计天线的亮温，转化为a×b分辨率的亮温成像模拟图。
2.根据权利要求1所述的一种能够区分不同粗糙面的被动毫米波辐射模拟方法，其特征在于：所述的剖分单元选择四面体。
3.根据权利要求1所述的一种能够区分不同粗糙面的被动毫米波辐射模拟方法，其特征在于：步骤一中所述的多个目标中包括高低不同的草地。
一种能够区分不同粗糙面的被动毫米波辐射模拟方法
技术领域
[0001] 本发明属于计算机辅助分析与设计以及软件设计技术领域，具体涉及一种被动毫米波辐射模拟的草地高度探测分辨方法。
背景技术
[0002] 被动毫米波成像模拟(Passive Millimeter Wave，PMMW)是被动毫米波成像领域里一项重要的环节，它可以帮助认识目标的辐射特性，解释辐射现象，寻找辐射规律以及判断实际测量结果的优劣。以上所述优点使得毫米波成像模拟近年来愈加受到重视。目前为止，大部分的工作都集中在待测目标表面为光滑平面的场景，而对于待测目标表面为粗糙面的场景则鲜有提及。
[0003] 传统的射线追踪方法中，当射线遇到粗糙面时即停止追踪，然后将粗糙面的发射率作为一个定值来处理。此方法的缺点是只适用于深粗糙面，而对于中等粗糙面或是微粗糙面，此方法则无法适用，以至于同种介质不同粗糙度的粗糙面无法加以区分。对于深粗糙面平面而言，在反射方向上几乎没有能量，所有入射能量会较为均匀的向周围各方向散射，因此可以适用朗伯近似定理。而对于中等粗糙度的粗糙面，除了在反射方向上有能量反射外，在四周也有能量的散射，而散射能量大小则是越靠近反射方向的能量越强，越远离反射方向的能量越弱。因此，需要一种新的计算模型及方法对同一介质，不同粗糙度的粗糙面进行更为精确的分辨，从而扩大被动毫米波辐射模拟的适用范围，以及提高了其分辨能力发明内容
[0004] 本发明的目的在于提供一种被动毫米波辐射模拟的草地高度探测分辨方法。
[0005] 本发明的具体步骤如下：
[0006] 步骤一、建立被测场景的三维模型。三维模型中包含地面、辐射计安装点和多个目标。
[0007] 步骤二、对步骤一中建立的三维模型进行剖分。
[0008] 步骤三、定义a个射线组，一个射线组包括b条射线。所有射线均以三维模型中的辐射计安装点为端点，且经过三维模型的地面。同一射线组内射线的俯仰角相等，且方位角依次等距增大；a个射线组内射线的俯仰角依次等距增大；a×b的数量即为所得图像的分辨率。
[0009] 步骤四、将1赋值给x和y。
[0010] 步骤五、定义第x个射线组内第y条射线的俯仰角为θ0，方位角为 第x个射线组内第y条射线经过的目标为第一目标。第x个射线组内第y条射线射入第一目标的点作为特征反射点。第一目标上含有特征反射点的侧面即为特征侧面。
[0011] 将沿着第x个射线组内第y条射线到达辐射计天线的毫米波作为入射波。入射波由第一反射波、第一目标辐射波和空气辐射波组成。第一目标辐射波为第一目标的特征反射点自身辐射向辐射计安装点的毫米波；空气辐射波为空气背景在特征反射点到辐射计安装点方向上辐射出的毫米波。在特征反射点反射出第一反射波的毫米波作为第二反射波；第二反射波与地面的交点为地面反射点。
[0012] 建立入射波传达的亮温 的表达式如式(1)所示：
[0013]
[0014] 式(1)中，T11(θ2,φ2)的表达式如式(2)所示；r1(θ2)为第一目标的特征反射点在入射波方向上的反射率；T1为第一目标的温度；e1(θ2)为第一目标的特征反射点在在入射波方向上的发射率； 为空气辐射波传达的亮温。
[0015] T11(θ2,φ2)的表达式如下：
[0016]
[0017] 式(2)中， 为空气背景向下到达地面的亮温；p为一种极化方式；rp(θ1)为地面反射点在第二反射波方向上的p极化的菲涅尔反射率；e为自然对数的底；k为波数的均方根高度；σ为地面的均方根高度；θ1为第二反射波在地面反射点的反射角。 的表达式如式(3)所示；Td为地面的温度；edp(θ1)为地面反射点在第二反射波方向上的p极化发射率，edp(θ1)的表达式如式(6)所示； 为地面与第一目标的特征侧面之间的极化旋转矩阵。
[0018] 的表达式如下：
[0019]
[0020] 式(3)中，定义N个特征反射波组，一个特征反射波组包括M条特征反射波。所有特征反射波均射入地面反射点。同一特征反射波组内特征反射波的俯仰角相等，且方位角依次等距增大；N个特征反射波组内特征反射波的俯仰角依次等距增大；θ′i为第i个特征反射波组内特征反射波的俯仰角； 为同一特征反射波组内第j个特征反射波的方位角。为第i个特征反射波组内第j个特征反射波表达的亮
温。 为第i个特征反射波组内第j个特征反射波射入地面反射点的反射率，其表达式如式(5)所示。
[0021] 若第i个特征反射波组内第j个特征反射波来自一个目标，则以该特征反射波作为二级特征反射波；以该目标作为第二目标，将θ′i赋值给θs，将 赋值给 代入式(4)计算，得到 否则， 的值等于式(2)中的
[0022] 的表达式如下：
[0023]
[0024] 式(4)中， 为空气背景向下到达第二目标的亮温；r2(θ3)为第二目标反射出二级特征反射波的点在二级特征反射波方向上的反射率；T2为第二目标的物理温度；e2(θ3)为第二目标反射出二级特征反射波的点在二级特征反射波方向上的发射率； 第二目标上反射出二级特征反射波的侧面与地面之间的极化旋转矩阵。
[0025] 的表达式如下：
[0026]
[0027] 式(5)中，q为与p不同的极化方式； 为第二层在方向 上的射线到第一层方向为 的射线的同极化散射系数； 为为第二层在方向
上的射线到第一层方向为 的射线的交叉极化散射系数。
[0028] 式(2)中的edp(θ1)的表达式如下：
[0029]
[0030] 式(6)中， 分别为相干反射系数、非相干反射系数； 的取值等于式(2)中的rp(θ1)； 的表达式如式(7)所示。
[0031] 的表达式如下：
[0032]
[0033] 为第一层在方向(θ1)上的射线到第二层方向为 的射线的同极化散射系数； 为第一层在方向(θ1)上的射线到第二层方向为 的射线的
交叉极化散射系数。
[0034] 将式将式(2),(3),(4),(5),(6),(7)代入式(1)，得到入射波传达的亮温之后进入步骤六。
[0035] 步骤六，若y＜b，则将y增大1，并重复执行步骤五；若y＝b且x＜a，则将1赋值给y，并将x增大1，并重复执行步骤五；若y＝b且x＝a，则进入步骤七。
[0036] 步骤七、将步骤四、五和六所得的a×b个到达辐射计天线的亮温，转化为a×b分辨率的亮温成像模拟图。
[0037] 进一步地，所述的剖分单元选择四面体。
[0038] 进一步地，步骤一中所述的多个目标中包括高低不同的草地。
[0039] 本发明具有的有益效果是：
[0040] 1、本发明在进一步精确地表征了粗糙面的被动辐射模型，使得本发明获取的亮温成像模拟图和实际探测情况愈加接近。
[0041] 2、本发明能够用于与辐射计拍摄出的亮温图像进行对比，从而判断辐射计的品质好坏。
[0042] 3、本发明在模拟草地场景时，对于处理同一种材质，不同高度的草地时，可以更加精确地探测出高低之间的区分，从而大幅增加了被动毫米波成像模拟的分辨能力。
附图说明
[0043] 图1为本发明中作为一种示例的机场跑道的光学图像；
[0044] 图2为本发明中作为一种示例的机场跑道的三维模型示意图；
[0045] 图3为本发明中作为一种示例的机场跑道的三维模型剖分后的示意图；
[0046] 图4为本发明的单根射线追踪示意图；
[0047] 图5为本发明的单根射线反演示意图；
[0048] 图6为通过传统方法得到的机场跑道亮温成像模拟图；
[0049] 图7为通过本发明得到的机场跑道亮温成像模拟图。
具体实施方式
[0050] 以下结合附图对本发明作进一步说明。
[0051] 如图1所示，步骤一、建立被测场景的三维模型。三维模型中包含地面、辐射计安装点和多个目标。以机场跑道图像为例，利用本发明对草地高度进行探测并加以分辨。如图1所示，图中为一机场跑道的实际光学照片，从照片中可以看出，跑道旁边的草地大致可以分为两种，较高的浅色区域和较矮的深色区域。对于机场跑道而言，建立的三维模型如图2所示。三维模型中，跑道以及远处的仓库、两边的草地都做了区分。地面及所有目标的物理温度、各侧面粗糙度、各角度的反射率、发射率均为已知量。
[0052] 步骤二、利用现有的剖分软件对步骤一中建立的三维模型进行剖分，剖分单元选择四面体。机场跑道三维模型剖分后的模型图如图3所示。
[0053] 步骤三、定义a个射线组，一个射线组包括b条射线。所有射线均以三维模型中的辐射计安装点为端点，且经过三维模型的地面。同一射线组内射线的俯仰角相等，且方位角依次等距增大；a个射线组内射线的俯仰角依次等距增大；a×b的数量即为所得图像的分辨率。
[0054] 步骤四、将1赋值给x和y。
[0055] 步骤五、如图4和5所示，定义第x个射线组内第y条射线的俯仰角为θ0，方位角为第x个射线组内第y条射线经过的目标为第一目标，即视为第x个射线组内第y条射线的方向与第一目标反射向辐射计的毫米波的方向相反。第x个射线组内第y条射线射入第一目标的点作为特征反射点。第一目标上含有特征反射点的侧面即为特征侧面。
[0056] 将沿着第x个射线组内第y条射线到达辐射计天线的毫米波(即经特征反射点反射向辐射计安装点的毫米波)作为入射波。入射波由第一反射波、第一目标辐射波和空气辐射波组成。第一目标辐射波为第一目标的特征反射点自身辐射向辐射计安装点的毫米波；空气辐射波为空气背景在特征反射点到辐射计安装点方向上辐射出的毫米波。
[0057] 在特征反射点反射出第一反射波的毫米波作为第二反射波；第二反射波与地面的交点为地面反射点；所有从空中射向地面反射点的毫米波以及经目标反射向地面反射点的毫米波均作为第三反射波；第二反射波是由所有第三反射波在地面反射点反射向特征反射点的毫米波叠加得到。
[0058] 建立入射波传达的亮温 的表达式如式(1)所示：
[0059]
[0060] 式(1)中，T11(θ2,φ2)·r1(θ2)为第一反射波传达的亮温(即地面经过第一目标并沿第x个射线组内第y条射线反射在辐射计安装点上的亮温)；T11(θ2,φ2)的表达式如式(2)所示；r1(θ2)为第一目标的特征反射点在入射波方向上的反射率(即反射角为θ2时的反射率)；T1·e1(θ2)为第一目标辐射波传达的亮温(即第一目标的特征反射点在反射角θ2方向上的辐射亮温；)T1为第一目标的温度；e1(θ2)为第一目标的特征反射点在在入射波方向上的发射率(即反射角为θ2时的发射率)；θ2为沿第x个射线组内第y条射线反射在辐射计安装点上的毫米波在特征反射点上的反射角； 为空气辐射波传达的亮温(即空气背景在沿着特征反射点到辐射计安装点方向上的亮温)。
[0061] T11(θ2,φ2)的表达式如下：
[0062]
[0063] 式(2)中， 是相干部分，其反射率是在菲涅尔反射率的基础上乘以该粗糙面的粗糙因子。 为空气背景向下到达地面反射点的亮温
p为极化方式(垂直极化或水平极化)；rp(θ1)为光滑地面(粗糙度为0)的地面反射点在第二反射波方向上的p极化的菲涅尔反射率(即反射角为θ1时的p极化发射率)；e为自然对数的底；k为波数(k＝2π/λ，意为2π长度上出现的全波数目，λ是真空中电磁波的波长)的均方根高度；σ为地面的均方根高度；θ1为第二反射波在地面反射点的反射角(即第二反射波与地面法向量的夹角)。 的表达式如式(3)所示；Td为地面的温度；edp(θ1)表示地面反射点在第二反射波方向上的p极化发射率(即射角为θ1时的p极化发射率)，edp(θ1)的表达式如式(6)所示； 为地面与第一目标的特征侧面之间的极化旋转矩阵。
[0064] 的表达式如下：
[0065]
[0066] 式(3)中，定义N个特征反射波组，一个特征反射波组包括M条特征反射波。所有特征反射波均射入地面反射点。同一特征反射波组内特征反射波的俯仰角相等，且方位角依次等距增大；N个特征反射波组内特征反射波的俯仰角依次等距增大；M·N个特征反射波均为第三反射波；θ′i为第i个特征反射波组内特征反射波的俯仰角； 为同一特征反射波组内第j个特征反射波的方位角。 为第i个特征反射波组内第j个特征反射波表达的亮温。 为第i个特征反射波组内第j个特征反射波射入地面反射点的反射率，其表达式如式(5)所示。
[0067] 若第i个特征反射波组内第j个特征反射波来自一个目标，则以该特征反射波作为二级特征反射波；以该目标作为第二目标，将θ′i赋值给θs，将 赋值给 代入式(4)计算，得到 否则， 的值等于式(2)中的
[0068] 的表达式如下：
[0069]
[0070] 式(4)中， 为空气背景向下到达第二目标的亮温；r2(θ3)为第二目标反射出二级特征反射波的点在二级特征反射波方向上的反射率；T2为第二目标的物理温度；e2(θ3)为第二目标反射出二级特征反射波的点在二级特征反射波方向上的发射率； 第二目标上反射出二级特征反射波的侧面与地面之间的极化旋转矩阵。
[0071] 的表达式如下：
[0072]
[0073] 式(5)中，q为与p不同的极化方式(即若p为垂直极化，则q为水平极化；若p为水平极化，则q为垂直极化)； 为第二层在方向 上的射线到第一层方向为的射线的同极化散射系数(即第i个特征反射波组内第j个特征反射波散射向俯仰角为θ，方位角为 方向的同极化散射系数)； 为为第二层在方向 上的射
线到第一层方向为 的射线的交叉极化散射系数(即第i个特征反射波组内第j个特征反射波散射向俯仰角为θ，方位角为 方向的交叉极化散射系数)。
[0074] 式(2)中的edp(θ1)的表达式如下：
[0075]
[0076] 式(6)中， 分别为相干反射系数、非相干反射系数； 的取值等于式(2)中的rp(θ1)，即光滑地面的地面反射点在第二反射波方向上的p极化的菲涅尔反射率(即反射角为θ1时的p极化发射率)； 的表达式如式(7)所示。
[0077] 的表达式如下：
[0078]
[0079] 为第一层在方向(θ1)上的射线到第二层方向为 的射线的同极化散射系数(即第二反射波散射向第i个特征反射波组内第j个特征反射波方向的同极化散射系数)； 为第一层在方向(θ1)上的射线到第二层方向为 的射线的交叉
极化散射系数(即第二反射波散射向第i个特征反射波组内第j个特征反射波方向的交叉极化散射系数)。
[0080] 将式(2),(3),(4),(5),(6),(7)代入式(1)，得到入射波传达的亮温 之后进入步骤六。
[0081] 步骤六，若y＜b，则将y增大1，并重复执行步骤五；若y＝b且x＜a，则将1赋值给y，并将x增大1，并重复执行步骤五；若y＝b且x＝a，则进入步骤七。
[0082] 步骤七、将步骤四、五和六所得的a×b个到达辐射计天线的亮温，转化为a×b分辨率的亮温成像模拟图(即每个亮温决定亮温成像模拟图上一个像素点的灰度值)。
[0083] 以步骤一中提及的机场跑道作为模拟对象，通过传统方法模拟出的图像如图6所示，通过本发明模拟出的图像如图7所示，对比两幅图可以发现，本发明得到的亮温成像模拟图中可以明显区分跑道旁边草地的高度，高的地方颜色偏浅，而矮草的颜色则偏深。可见，本发明相对传统方法，能够模拟出更加真实的成像结果。以本发明所得的亮温成像模拟图，与真实的辐射计获取的亮温图进行对比，能够判断辐射计获取的亮温图是否准确，从而判断辐射计的质量好坏。
