[發明專利]一種基于CPU和GPU混合的高速并行地形遮蔽計算方法

有效

• 申請號: 201810968397.1
• 申請日: 2018-08-23
• 公開（公告）號: CN109118583B
• 公開（公告）日: 2022-09-13
• 發明（設計）人: 王毅;俞信;廉海明;張培
• 申請（專利權）人: 中國科學院電子學研究所蘇州研究院
• 主分類號: G06T17/05
• 分類號: G06T17/05;G06T1/20
• 代理公司: 南京經緯專利商標代理有限公司 32200
• 代理人: 樓高潮

1.一種基于CPU和GPU混合的高速并行地形遮蔽計算方法，其特征在于，該方法包括以下步驟：

步驟1）建立網格化數字高程模型，以高程矩陣方式來表示地面高程，根據需要計算的區域和高程的讀取范圍，將該區域劃分成若干網格，對網格進行編碼，記錄每個網格的高程值，計算網格化的高程矩陣，使得區域內任一位置有對應的網格編碼，從而得到該位置的高程值；

步驟2）采用遮蔽分析算法對區域的地形數據進行分析，將兩點在高程矩陣中的編碼擬合成一條直線，通過計算直線與高程網格矩陣的交點，得到視點到觀察點的地形剖面，根據點與點之間的夾角判斷是否通視，計算出最遠通視點；

步驟3）高程數據校正，將影響高程數據的地球曲率和大氣衰減因素作為校正參數，構建校正公式對高程數據進行校正。

2.根據權利要求1所述的基于CPU和GPU混合的高速并行地形遮蔽計算方法，其特征在于，所述步驟1）中，將需要計算的區域劃分成網格m×n，網格編碼以起始點的地理坐標為網格的原點(0，0)，水平方向為0到n-1，垂直方向為0到m-1。

3.根據權利要求2所述的基于CPU和GPU混合的高速并行地形遮蔽計算方法，其特征在于，所述步驟1）中，計算網格化的高程矩陣，包括以下步驟：

步驟1.1）根據需要計算的區域的網格大?。╩×n），最小點的位置P_min(lon, lat),其中lon和lat分別表示經度和緯度，最大點的位置P_max(lon, lat),計算點P(lon, lat)對應的網格編碼（x, y）,其計算公式如下：

單位網格的經緯度分辨率為：

R (?lon, ?lat) = P_max- P_min/ (n, m)；

步驟1.2）根據單位網格的經緯度分辨率和起始點的經緯度計算出對應點的網格編碼，計算公式如下：

（x, y） = (P-P_min) / R。