[發明專利]一種基于柵格化GIS 空間關系判斷方法及其系統

無效

• 申請號: 201110222077.X
• 申請日: 2011-08-04
• 公開（公告）號: CN102270236A
• 公開（公告）日: 2011-12-07
• 發明（設計）人: 程振林;董慧;趙艷偉;方金云
• 申請（專利權）人: 中國科學院計算技術研究所
• 主分類號: G06F17/30
• 分類號: G06F17/30
• 代理公司: 北京律誠同業知識產權代理有限公司 11006
• 代理人: 祁建國;梁揮

1.一種基于柵格化GIS空間關系判斷方法，其特征在于，包括：

步驟1，輸入矢量方式表示的幾何體表達的圖形，由柵格單元對所述圖形進行壓蓋，柵格單元以對其壓蓋所述圖形的要素ID為要素索引，柵格單元以坐標值為位置索引；

步驟2，對形成索引的圖形進行分塊，每個柵格單元具有對其壓蓋的所述圖形在所述柵格單元的實際占用面積的信息；

步驟3，兩個經過柵格單元壓蓋的所述圖形中，坐標相同的柵格單元相互對應，將相對應的兩個柵格單元的實際占用面積的信息進行比較，通過地理信息空間關系判斷方法實現誤差分析和置信度分析，得出所述兩個經過柵格單元壓蓋的以矢量方式表示的圖形的GIS空間關系判斷結果。

2.如權利要求1所述的基于柵格化GIS空間關系判斷方法，其特征在于，所述步驟1進一步為，

步驟21，輸入矢量方式表示的圖層中地理要素的圖形的矢量點，按顯示屏幕的分辨率對地理要素的矢量點坐標進行坐標轉換，按顯示屏幕的像素點進行柵格劃分，柵格單元以坐標值為位置索引，柵格單元以對其壓蓋的圖形的要素ID為要素索引；

步驟22，根據線段斜率及增量依次判斷該線段與所有橫縱柵格線的交點，輪廓由N個矢量點組成，每兩個相鄰矢量點組成一條線段，輪廓掃描過程就是逐個線段的掃描，這里的線段是指其中兩個點的坐標組成的線段，交點的計算需要轉換成亞像素精度，增量的計算公式描述如下：

delta=(256-fy)×ΔxΔy]]>

其中delta為增量，fy為該線段壓蓋柵格點的實數坐標分量Y部分經過坐標轉換通道所得的小數部分，設線段的起始點坐標為(x1，y1)，終止點坐標為(x2，y2)，Δx＝x₂-x₁，Δy＝y₂-y₁，可計算出受該直線所影響的整數柵格單元的cover和area，計算公式描述如下：

cov?er＝fy2-fy1

area＝(fx2+fx1)×cov?er

其中(fx1，fy1)為線段起始點的小數坐標部分，(fx2，fy2)為線段終止點的小數坐標部分；

步驟23，對于每個圖形，繪制控制器遍歷圖形的輪廓經過的柵格單元，依據柵格單元的覆蓋面積判斷所述柵格單元是否被圖形完全填充，對完全填充的柵格單元和未完全填充的柵格單元分別進行標記，將輪廓內的柵格單元進行跨度填充，將這些柵格單元標記為完全填充；

步驟24，將完全填充的柵格單元的圖形在柵格單元的實際占用面積的比值設置為1；對于未完全填充的柵格單元，根據所述柵格單元的權值和覆蓋面積計算所述柵格單元的圖形在柵格單元的實際占用面積的比值，保存柵格單元的位置索引、要素索引和實際占用面積的比值，進而生成所述圖形對應的柵格底圖。

3.如權利要求1所述的基于柵格化GIS空間關系判斷方法，其特征在于，所述步驟2進一步為，

步驟31，對輸入的圖層進行柵格索引生成，存放在基于圖片像素與矢量要素對應算法的柵格單元中；

步驟32，根據輸入的圖層，找到壓蓋的圖層片段。

4.如權利要求1所述的基于柵格化GIS空間關系判斷方法，其特征在于，在柵格化GIS空間關系判斷中兩個圖層中的圖形是否有相交時，所述步驟3進一步為，

步驟41，返回被搜索圖層中與搜索圖層交集不為空的所有對象；

步驟42，滿足下面任何一種情況都為相交：

其中，上圖層內部與下圖層內部交集不為空；

上圖層內部與下圖層邊界交集不為空；

上圖層邊界與下圖層內部交集不為空；

上圖層邊界與下圖層邊界交集不為空。