[發明專利]基于GIS的X波段面雨量雷達數據提取方法

在審

• 申請號: 201710094219.6
• 申請日: 2017-02-21
• 公開（公告）號: CN106908797A
• 公開（公告）日: 2017-06-30
• 發明（設計）人: 牛睿平;李輝;唐躍平;張永兵;杜紅娟;王巖
• 申請（專利權）人: 水利部南京水利水文自動化研究所
• 主分類號: G01S13/95
• 分類號: G01S13/95;G01W1/14;G06F17/30
• 代理公司: 南京眾聯專利代理有限公司32206
• 代理人: 葉涓涓

1.基于GIS的X波段面雨量雷達數據提取方法，其特征在于，包括如下步驟：

步驟一，數據的存儲

以雷達所在位置為圓心，按極坐標方式存儲，共360個徑向，從0度至359度排列，每徑向數據存儲若干雨強值，每個雨強值占用兩個字節；

步驟二，極坐標數據轉換為柵格數據

步驟二-1，生成柵格數據

(1)生成空柵格

將雷達所在位置的經緯度通過投影轉換為坐標，調用函數生成空柵格數據；

(2)極坐標數據寫入到空柵格

①定義一個二維數組，數組中存儲有極坐標方位角angel，庫數n，利用二重循環將雷達極坐標數據寫入到二維數組中；

②從用來存儲極坐標數據的空柵格抽象獲得像素塊，用于臨時存儲雨量數據；

③通過下式計算柵格點到雷達中心的距離：

$dis=\sqrt{{(x-centx)}^2-{(y-centy)}^2}$

式中dis為柵格點到雷達中心距離，x，y為柵格點坐標；

④通過下式計算柵格點到雷達中心的方位角：

式中azim為方位角，θ＝arctan(y/x)；

⑤根據柵格點到雷達中心的距離dis、方位角azim在極坐標中找到相應的雨量值，從而獲得該柵格點的雨量值，從像素塊左上角開始，先X方向在Y方向，依次遍歷，賦予柵格點雨量值，將像素塊重新寫回空柵格，完成雷達極坐標數據到柵格數據的轉換。

2.根據權利要求1所述的基于GIS的X波段面雨量雷達數據提取方法，其特征在于，還包括以下步驟：

步驟三，雷達雨量圖像展示

步驟三-1，對柵格圖進行分級渲染；

步驟三-2，利用色階渲染形成雨量顏色分級；

步驟三-3，對雨量圖進行透明度設置。

3.根據權利要求1所述的基于GIS的X波段面雨量雷達數據提取方法，其特征在于，還包括步驟四，雷達雨量數據提取：用點、面批量化方式提取雷達雨量值，根據站點和/或行政邊界和/或小流域提取。