[發(fā)明專利]適用于瓦片地圖服務的海洋環(huán)境預報產(chǎn)品并行處理方法

在審

• 申請?zhí)枺?/td>: 201410319012.0
• 申請日: 2014-07-04
• 公開（公告）號: CN104063496A
• 公開（公告）日: 2014-09-24
• 發(fā)明（設計）人: 劉仁義;杜震洪;張豐;林賢輝;王葉晨梓
• 申請（專利權）人: 浙江大學
• 主分類號: G06F17/30
• 分類號: G06F17/30
• 代理公司: 杭州求是專利事務所有限公司 33200
• 代理人: 張法高

1.一種適用于瓦片地圖服務的海洋環(huán)境預報產(chǎn)品并行處理方法，其特征在于包括如下步驟：

1)讀取海洋環(huán)境預報產(chǎn)品的元數(shù)據(jù)，提取元數(shù)據(jù)中的經(jīng)緯度，并將WGS84坐標下的經(jīng)緯度轉換為Web Mercator投影坐標下的經(jīng)緯度；

2)海洋環(huán)境預報產(chǎn)品具有多個深度層，將每個深度層根據(jù)并行機器的核數(shù)切分為n份數(shù)據(jù)，設計MATLAB并行策略；

3)根據(jù)實際需求設置地圖分級數(shù)；

4)將步驟2)切分得到的數(shù)據(jù)根據(jù)當前地圖級別插值，利用MATLAB繪圖指令并行生成分級切片文件；

5)重復步驟4)，完成當前地圖級別下海洋環(huán)境預報產(chǎn)品所有參數(shù)的繪制；

6)將步驟4)生成的切片文件采用LOD方式組織，重復步驟3)～步驟5)，完成所有級別下海洋環(huán)境預報產(chǎn)品所有參數(shù)的繪制；

7)等待所有MATLAB并行策略處理完畢，生成元數(shù)據(jù)。

2.根據(jù)權利要求1所述的一種適用于瓦片地圖服務的海洋環(huán)境預報產(chǎn)品并行處理方法，其特征在于所述的步驟1)為：讀取海洋環(huán)境預報產(chǎn)品的元數(shù)據(jù)，提取元數(shù)據(jù)中的經(jīng)緯度，海洋環(huán)境預報產(chǎn)品中的經(jīng)緯度是WGS84經(jīng)緯坐標，瀏覽器客戶端需要加載的產(chǎn)品用的投影坐標系是Web Mercator投影，所以對多源地理數(shù)據(jù)處理需要涉及投影變換，利用公式(1)～公式(6)對WGS84坐標和Web Mercator投影坐標進行轉換；

N0=6378137.0Sqrt(1-Pow(0.081819190843,2)*Pow(sin(RefLat*PI180),2))---(1)]]>

q1=Log(Tan((180.04.0+Lat2.0)*PI180.0))---(2)]]>

q2=0.0818191908432*Log(1+0.081819190843*Sin(Lat*PI180.0)1-0.081819190843*Sin(Lat*PI180.0))---(3)]]>

q＝q1-q2(4)

X=N0*Cos(RefLat*PI180.0)*(Long-CenMeri57.29577951)---(5)]]>

Y=N0*Cos(RefLat*PI180.0)*q---(6)]]>

式中Lat代表WGS84緯度；Long代表WGS84經(jīng)度；RefLat代表基準緯度；CenMeri代表中央子午線；X代表Web Mercator橫坐標；Y代表Web Mercator縱坐標；N0、q1、q2、q為中間過程參數(shù)。