windows下ncl產生tiff圖（案例）

標籤：

一：安裝軟體和準備資料

1.需要安裝Vapor（注意安裝路徑不要存在空格）

註：版本2.4.2及以後

2.安裝NCL，方法見http://www.cnblogs.com/striver-zhu/p/4574819.html

註：

　　a.得包含這兩個可執行檔：convert && psplit，一般安裝時會自行安裝進去

　　b. NCL_5.1.0版本存在一個問題：wrf2geotiff.ncl中的492行的str_lower需要去掉，否則運行時報錯

3.下載資料：jamingWRF2

4.3個ncl檔案（NCAR\VAPOR\share\examples\NCL目錄下）：
　　wrf_Height.ncl （this plots humidity, temperature, pressure and wind at a fixed elevation）
　　wrf_Precip.ncl （this plots precipitation tendency with pressure isobars）
　　wrf_crossSection2.ncl （this is a vertical plot of relative humidity）

二：資料轉換
1.將WRF資料轉換成VDC資料，從Vapor2.0開始，便不需要轉換，但是大資料還是轉換以下比較好

　　通過cmd執行以下命令(進入命令所在目錄下)：

　　wrfvdfcreate：掃描WRF資料，產生用於描述WRF資料的Vapor中繼資料
　　　　（wrfvdfcreate wrfout_d02_2008-09-28.nc jangmi-09-28.vdf）
　　wrf2vdf：將所有的WRF資料中的2D和3D變數轉換成Vapor使用的資料格式，即VDC格式
　　　　（wrf2vdf jangmi-09-28.vdf wrfout_d02_2008-09-28.nc   =》 產生的變數會存在jangmi-09-28_data目錄下）

三：
1.將NCL繪製的資料圖形轉換成georeferenced格式（在工作目錄下，需要添加.hluresfile檔案，該檔案只能通過cygwin的vi建立（：wq退出儲存））

　　wrf_Height.ncl：產生在一個特定平面上的圖形（）
　　wrf_Precip.ncl：產生一個呈現降水趨勢和海平面壓力的圖形
　　wrf_CrossSection2.ncl：產生溫度和相對濕度的垂直圖形，沿著一個特定的平面上，平行於XZ軸
　　

　　我們將圖形轉換成.PS輸出檔案和將這些.PS檔案轉換成地理參考TIFF檔案

　1）.修改指令碼，使其僅產生一個映像，在每個時間步長。 （如果想在相同的時間步長，使用多個映像，則需要為每個不同的映像的單獨運行指令碼）

　2）.修改這個指令碼，遍曆WRF檔案中的所有步長。這可能涉及外迴圈（迴圈多個檔案）和內迴圈（迴圈每個檔案中的時間步）

　3）.修改指令碼，輸出為.PS檔案

　4）.將下面的新行到NCL指令碼：

　　a.在頂部，插入一行來載入wrf2geotiff.ncl：load “wrf2geotiff.ncl”（存在C:\NCAR\VAPOR\share\examples\NCL）
　　（這個指令碼和其他NCL的例子是預設安裝子目錄,unix下在share/vapor-x.x.x/examples/NCL，Windows下在$(VAPOR_HOME)/share/examples/NCL/）

　　b.在NCL工作站建立之後，調用wrf2geotiff=wrf2geotiff_create(wks)

　　c.如果這是一個垂直圖形，禁止地理參考：
　　　　wrf2geotiff_disableGeoTags(wrf2gtiff)

　　d.設定圖形屬性[email protected]=False

　　e.圖形每產生一次，插入兩行：
　　　　wrf2geotiff_write()
　　　　frame(wks)

　　f.在檔案末尾調用wrf2geotiff_close()，將會合并圖形產生tiff格式圖形

四：擷取地理檔案

1.為類比地區獲得地形圖片：需要經緯度範圍，同時需要連網，進入getWMSImage.sh所在檔案夾

cygwin下：getWMSImage.sh –o jangmiTerrain.tiff 115 15 130 30

 

windows下ncl產生tiff圖（案例）