圖片摘要

本專欄目的是將netcdf格式的雪覆蓋數(shù)據(jù)進行重新投影，需要使用的工具包括matlab和Arcgis，下面進入正題。

1.數(shù)據(jù)的下載與讀取---matlab

最近我需要讀取北半球的冰雪覆蓋數(shù)據(jù)，下載的是MODIS/Terra Snow Cover Monthly L3 Global 0.05Deg CMG, Version 61，文件格式為netcdf，在matlab中讀取發(fā)現(xiàn)其為polar_stereographic極坐標投影。

matlab讀取程序

繪圖顯示與官網(wǎng)提供的數(shù)據(jù)與投影方式一致。但是在matlab中讀取netcdf文件的經(jīng)緯度坐標沒有辦法進行后續(xù)的處理與計算，因此需要進行重新投影轉(zhuǎn)換。

2.數(shù)據(jù)的投影柵格---arcgis

下面我們將NetCDF文件導(dǎo)入arcgis中，使用arcgis工具箱中的“多維工具”下的“創(chuàng)建NetCDF柵格圖層”導(dǎo)入，得到以下的結(jié)果：

接下來選擇【數(shù)據(jù)管理工具】->【投影和變換】下的【柵格】->【投影柵格】，選擇輸出坐標系為GCS_WGS_1984，點擊確定，注意需要將結(jié)果導(dǎo)出后新建空白地圖，再把文件重新導(dǎo)入arcgis中，即可得到正確的投影，此時再把這個文件導(dǎo)出成NetCDF文件。選擇【多維工具】下的【柵格轉(zhuǎn)NetCDF】。

再在matlab中讀取即可得到正確的結(jié)果。

程序代碼：

得到最終的結(jié)果：

后續(xù)有新的方法將更新！歡迎交流學(xué)習(xí)！