基于面元網格化的空間目標光學特性計算方法

鮑文卓，叢明煜，張 偉，程 軍，曹移明

(哈爾濱工業大學空間光學工程研究中心，哈爾濱 150001，latermouse@126.com)

基于面元網格化的空間目標光學特性計算方法

鮑文卓，叢明煜，張 偉，程 軍，曹移明

(哈爾濱工業大學空間光學工程研究中心，哈爾濱 150001，latermouse@126.com)

為實現復雜結構空間目標光學特性的工程化計算，基于計算幾何的面元網格化思想，提出了一種空間目標光學特性計算方法.建立了空間目標標準幾何結構、面元網格劃分、面元遮擋、面元照度計算的數學模型，設計開發了空間目標光學特性計算軟件.針對不同的衛星目標，進行了光學特性的計算，并將計算結果分別與解析法計算結果和試驗測量結果進行比較分析，結果表明面元網格化的計算精度優于解析方法，且與試驗結果相比其計算偏差小于一個視星等，滿足工程應用的要求.

空間目標;光學特性計算;幾何體建模;面元網格化

近距離導航測量、空間目標天基監視等領域中，可見光載荷具有質量輕、體積小、定位測量精度高等特點，已得到廣泛應用.在空間目標光學特性計算方面［1－7］，主要采用2種計算方法:標準幾何體解析法、蒙特卡羅光線追跡法.標準幾何體解析法計算模型簡單、精度較低，主要用于光學載荷的總體方案設計階段;蒙特卡羅光線追跡法考慮了追跡光線在幾何體面型之間的多次反射現象，計算模型復雜、計算速度慢，主要用于目標特性典型與邊界觀測條件下的計算中.

本文應用計算幾何面元網格化思想，提出并實現了一種基于面元網格化的空間目標光學特性計算方法，建立了空間目標標準幾何結構、面元網格劃分、面元遮擋、面元照度計算的數學模型，實現了空間目標光學特性的高效計算.

1 面元網格化計算方法

1.1 計算方法原理

空間目標光學特性可等效為目標離散化的表面面元光學特性的疊加.圖1給出了將空間目標表面分解為一系列相互獨立的面元集合的示意圖，通過求取各個面元的反射特性，進行整體疊加，可求得空間目標整體的光學特性.

圖1 空間目標面元網格化示意圖

式中:Pk為面元中心坐標(xk，yk，zk)，Rk為面元法線方向矢量，Ak為面元面積，εk為面元表面漫反射系數，Sk為面元光照加權因子，Vk為面元觀測加權因子.

空間目標反射光譜輻照度為

式中:E(λ)為空間目標反射光譜輻照度，Ek(λ)為面元Fk的反射光譜輻照度，m為面元數目.

計算方法流程如圖2所示，主要包括:1)目標標準幾何體建模，將目標描述為二次曲面幾何體的形式;2)目標網格面元生成，將各二次曲面幾何體離散化為面元的集合;3)面元遮蔽性判斷，對離散化后的面元進行可見性條件判斷;4)面元特性計算;5)目標整體特性計算.

圖2 目標光學特性計算流程

1.2 坐標系變換關系

空間目標光學特性計算需要定義下列坐標系:衛星本體坐標系Rsat(OXYZ)、子結構體參考坐標系 Rref(oxyz)、觀測視線坐標系 R′(Ox′y′z′)、光照方向坐標系 R″(Ox″y″z″)，各坐標系之間的坐標變換關系如圖 3，圖4所示(坐標系原點O、o分別位于目標與子結構體質心位置，?，σ分別為觀測方向Robs與光照方向Rsun在衛星本體坐標系中的方位角、高低角)［8］.

圖4中，矩陣BR、BL、BO分別為衛星本體坐標系到子結構體參考坐標系、觀測視線坐標系、光照方向坐標系的變換矩陣.各個坐標系之間的空間位置矢量變換滿足:

式(2)中的變換矩陣為

式中:l，m，n分別為子結構體參考系坐標方向在本體坐標系下的單位矢量，RY()，RZ()分別為繞Y、Z軸的旋轉變換矩陣，Rs為子結構體參考坐標系坐標原點在本體坐標系中的位置矢量.

圖3 坐標系的相互關系示意圖

圖4 坐標系的變換關系

1.3 目標標準幾何體與面元網格建模

為了提高計算效率，忽略對空間目標光學特性影響較小的結構體，將空間目標等效為長方體、球體、圓柱體與圓錐體的組合.圖5給出了目標結構等效示例，可看出，簡化后的目標整體結構保持了原有目標的主要形狀特征.

圖5 典型衛星目標結構建模圖

復雜的空間目標經過結構簡化，可等效為4種二次曲面幾何體的組合.這4種幾何體由矩形面、圓形面、球面、圓柱面以及圓錐面5類二次曲面組成.

采用一種快速矩形、梯形面元劃分方法，對各標準二次曲面參考系坐標軸向尺度范圍作多等份劃分，可達到面元網格劃分均勻性的目的.

以矩形平面網格劃分為例(如圖6)，將其在x、y軸方向(子結構體參考坐標系)劃分為L×W個面元，面元 Fi，j的特征參數為

圖6 矩形平面網格劃分方法示意圖

依據上述原理，可將圓形面、圓柱面、圓錐面在極坐標系下進行劃分，球面在球坐標系下進行劃分，面元劃分數學模型如表1所示.

表1 不同面型網格面元特征參數建模

1.4 面元可見與光照遮蔽計算模型

面元觀測可見性與光照遮蔽性邏輯特性參數是目標光學特性的加權因子(取1或0)，采用光線追蹤法求取面元加權因子.

空間目標經標準幾何結構建模與網格劃分后，其特征量在子結構體參考坐標系下可表達為面元特征量的集合形式:

將式(9)由子結構體參考坐標系變換至衛星本體坐標系，面元特征可表達為

式中:M為空間目標特征量集合，m為空間目標面元網格劃分個數.

面元觀測可見性計算的判定步驟為:

1)將式(10)的空間目標面元集合由本體坐標系變換到觀測視線坐標系，面元特征可表達為

2)將M′集合分成若干子集，每個子集內的面元具有相同的 y′、z′坐標值，即有

式中:Ni為 M′分類子集，n為分類子集的數目(n

3)將各分類子集內的面元按x′方向進行深度排序，具有最大x′值的面元視為可視面元，集合Ni中的面元特征量F′k中的可視加權因子為

1.5 面元光學特性計算模型

空間目標在軌時，接收到的輻射主要來自太陽、其他天體以及地球反照.其中，太陽輻射為目標光學特性的主要輻射源，其他輻射源在目標光學特性計算時可忽略.

面元光學特性計算原理如圖7所示，任意面元反射的太陽光能量在傳感器入瞳處產生的光譜照度為

式中:R為面元到光學載荷入瞳處的距離，Θ1，Θ2分別為面元法線方向n與太陽方向矢量Rsun、觀測方向矢量R的夾角，Esun(λ)為太陽在大氣層外的可見光波段光譜輻照度，

式中:rsun為太陽參考半徑，Rsun為太陽到空間目標的距離，M(λ)為太陽的可見光波段光譜輻出度.

圖7 面元光學特性計算原理示意圖

2 計算軟件結構與功能設計

基于上述數學模型與計算方法，采用Visual C++6.0編譯環境，開發出用于空間目標光學特性分析計算的集成軟件.軟件功能結構與操作界面如圖 8，圖9所示.

集成軟件具有:1)空間目標標準幾何體結構建模;2)空間目標標準幾何體解析法特性計算;3)空間目標面元網格化方法特性計算;4)計算結果的存儲與管理.

圖8 集成軟件功能結構

圖9 集成軟件界面及計算結果示意圖

3 計算結果與分析

為考核與驗證基于面元網格化光學特性計算的精度，采用下述方法進行分析與驗證:

1)面元網格化方法與標準幾何體解析法計算結果進行比較.通過二次曲面標準幾何體計算結果，確定滿足網格化方法計算收斂條件的網格劃分尺度參數;通過整星計算結果，驗證網格化方法的計算精度;

2)面元網格化方法與縮比模型(衛星本體、整星模型)試驗結果進行比較，進一步驗證網格化方法的計算精度和實際應用的可行性.

3.1 網格尺度參數確定

計算條件采用觀測方向與光照方向均選定為幾何體本體的 －x方向，觀測距離設定為15 km.計算結果如圖10所示，從圖10中可看出由于面元網格化會帶來結構體表面面形的畸變(曲面面元等效為平面面元)，使得面元網格化方法與解析法(真值)的計算結果有一定的誤差，但是隨著面元網格劃分密度的提高，誤差趨近于0.表2給出了滿足計算收斂條件(計算精度)下的標準幾何體網格劃分尺度參數.

圖10 標準幾何體目標光學特性對比結果

表2 滿足收斂條件下的結構體網格劃分尺度

3.2 計算精度驗證

對圖11(a)所示衛星構型，分別采用解析法與面元網格化法對衛星在一個軌道周期內光學特性變化情況進行對比計算.計算條件:衛星為柱體本體+長方體雙翼帆板結構(衛星柱體本體:底面半徑0.5 m，高1 m，反射率設為0.8;雙翼帆板尺寸參數如下:3 m×0.1 m×1 m，反射率均設為0.8);衛星運行軌道設定為600 km太陽同步軌道(降交點地方時06:30 am);計算起始時刻設為2009年1月1日0時;目標星姿態繞軌道系y軸的角速度為1(°)/s;觀測星與目標星處于同一軌道，且位于目標星后方15 km.

圖11 衛星結構示意圖

計算結果如圖12所示，由圖12可看出，由于面元網格化方法考慮了衛星各子結構體之間的相互遮擋關系，其特性計算結果優于解析法.

圖12 衛星光學特性計算結果對比

3.3 試驗結果比較

對圖11(b)所示衛星構型中的衛星本體及整星縮比模型進行了地面光學特性試驗測量，將試驗結果與面元網格化計算結果進行比較(如圖13)，結果表明面元網格化方法計算結果與試驗測量結果的偏差小于一個視星等，可滿足光學載荷系統地面數字仿真、半實物仿真對目標特性計算精度的要求.

圖13 面元網格化方法計算結果與實測數據比較

4 結論

1)建立了用于空間目標光學特性計算的面元網格化方法，給出了相關的數學模型與計算方法，并實現了用于空間目標光學特性分析計算的集成軟件，實現了空間目標光學特性的高效計算.

2)面元網格化方法具有較好的收斂性，方法特性計算精度優于解析法;

3)面元網格化方法與試驗測量結果的偏差小于一個視星等，可滿足光學載荷系統地面數字仿真、半實物仿真對目標特性計算精度的要求.

［1］RORK E W，LIN S S，YAKUTIS A J.Ground-based E-lectro Optical Detection of Artificial Satellites in Daylight from Reflected Sunlight［D］.Massachusetts Avenue，Cambridge，USA:Massachusetts Institute of Technology，1982:3－35.

［2］RASK J D.Modeling of Diffuse Photometric Signatures of Satellites for Space Objects Identification［D］.Ohio State，USA:Air Force Institute of Technology，1983:20－119.

［3］ZISSIS G J.The Infrared＆ Electro-optical Systems Handbook vol.1-Sources of Radiation［M］.Bellingham，Wash:SPIE Optical Engineering Press，1993:151－157.

［4］HROVAT D.Hyperspectral Analysis of Space Objects:Signal to Noise Evaluation［D］.Ohio State，USA:Air Force Institute of Technologies，1993:88－91.

［5］MCCUE G A，WILLIAMS J G，MORFORD J M.Optical characteristics of artificial satellites［J］.Planetary and Space Science， 1971，19(8):851－868.

［6］吳振森，竇玉紅.空間目標的可見光散射與紅外輻射［J］.光學學報， 2003，23(10):1250－1254.

［7］KOENDERINK J J，van DOORN A J.Shading in the case of translucent objects［C］//Proc of SPIE.San Jose，CA:Human Vision and Electronic Imaging VI， 2001，4299:312－320.

［8］許世文，龍夫年，付苓.實時星場模擬器中的坐標變換［J］.哈爾濱工業大學學報， 1998，30(5):118－120.

An optical characteristics calculating method based on surface mesh-creation for space targets

BAO Wen-zhuo，CONG Ming-yu，ZHANG Wei，CHENG Jun，CAO Yi-ming

(Research Center for Space Optical Engineering，Harbin Institute of Technology，Harbin 150001，China，latermouse@126.com)

In order to realize the optical characteristics calculation of space targets with complex shape，an optical characteristics calculating method for space targets was presented based on the general idea of surface mesh-creation.The mathematical models of standard geometries for targets，mesh-creation methods，visibility of the target，and the luminance calculation were given.The software calculating optical characteristics of space targets was designed and compiled.For several satellite targets，optical characteristics calculations were conducted with this method，and the results were compared with those of analytic method and experiments.It is showed that the mesh-creation method is better than the analytic method in precision，and the errors are less than one apparent magnitude compared with experimental results，so this method can be applied in practice.

space target;optical characteristics calculation;target geometry modeling;surface mesh-creation

O432.2

A

0367－6234(2010)05－0710－06

2009－04－01.

國家高技術發展研究計劃資助項目(2006AA704215).

鮑文卓(1983—)，男，博士研究生;

叢明煜(1964—)，男，教授，博士生導師;

張 偉(1962—)，男，教授，博士生導師.

(編輯 張 紅)