NCL繪圖一些注意事項(xiàng)

了解繪圖需求
繪制線性圖gsn_csm_xy呢還是地圖gsn_csm_contour_map呢?
熟悉數(shù)據(jù)屬性
什么叫熟悉數(shù)據(jù)屬性?比如你手中有一組降水?dāng)?shù)據(jù),那么你首先要了解它是站點(diǎn)數(shù)據(jù)還是再分析數(shù)據(jù)?是nc文件還是grib文件又或者是txt文件?
然后根據(jù)屬性進(jìn)行讀取
nc文件:f=addfile(path,"r")
txt文件:?f= asciiread(filename,-1,"string")
了解數(shù)據(jù)變量信息
了解變量信息的方式很多,比如nc文件你可以直接利用三方軟件panoply

(軟件需要搭建Java環(huán)境)讀取打開(kāi),可以看到變量的維度、屬性等。其次,利用nc_filedump查看nc文件。


當(dāng)然也可以直接編寫(xiě)ncl代碼printVarSummary。

一般情況下txt文件,官方為配有說(shuō)明文檔,如果沒(méi)有說(shuō)明你的文件不是官方的,記得找給你文件的人。
為什么要了解變量信息?
就像畫(huà)畫(huà)一樣,無(wú)論是風(fēng)景畫(huà)還是人物肖像畫(huà),在動(dòng)筆之前,都需要了解整體風(fēng)景分布構(gòu)造,人物輪廓(當(dāng)然我不是專業(yè)畫(huà)畫(huà)的, 可能也有那種一上來(lái)就畫(huà)的神仙畫(huà)手)。
當(dāng)你了解的數(shù)據(jù)的基本信息,以三維nc數(shù)據(jù)舉例,X(T,lat,lon),首先你得知道X代表的是什么?降水,那是日降水還是月降水?也就是cd_calendar(T)是yyyymmdd還是yyyymm,?如果是月降水單位是mm/day還是mm/month?這些都是需要知道的,以便必要時(shí)進(jìn)行單位轉(zhuǎn)換。其次對(duì)于經(jīng)緯度的拿捏,如果只繪制一個(gè)變量圖,相對(duì)簡(jiǎn)單一點(diǎn)不會(huì)涉及分辨率不對(duì)應(yīng)的問(wèn)題;如果是多個(gè)變量進(jìn)行對(duì)比,比如相關(guān),合成分析等涉及多個(gè)變量同時(shí)出現(xiàn)在一張圖上時(shí),就需要考慮經(jīng)緯度對(duì)應(yīng),需要將高分辨率轉(zhuǎn)低分或低分辨率轉(zhuǎn)高分,可以互轉(zhuǎn)的前提下建議高轉(zhuǎn)低,這樣不會(huì)過(guò)分損失數(shù)據(jù)精確度。需要特別注意:正常情況下數(shù)據(jù)各個(gè)維度均是按順序排列,eg. T=19900101-20220101, lat = (-90,90), lon =(-180,180), 有些數(shù)據(jù)集例外,時(shí)間,經(jīng)緯度可能是倒敘排列,導(dǎo)致繪出來(lái)的圖奇奇怪怪,這個(gè)時(shí)候就有必要檢查一下經(jīng)緯度信息并加以調(diào)整。
這里提一下,各位再向別人請(qǐng)教之前,一定要確保自己已經(jīng)有效搜索過(guò)各個(gè)平臺(tái)了,實(shí)在無(wú)計(jì)可施了,再去請(qǐng)教。
貼一個(gè)NCL官方文檔,不會(huì)還有人不知道吧~https://www.ncl.ucar.edu/??如果有人不會(huì)到相關(guān)目錄下查找指令,進(jìn)入該網(wǎng)址直接Ctrl+F然后輸入你要了解的函數(shù)或程序。
氣象家園:http://bbs.06climate.com/? 你遇到的問(wèn)題這里都曾遇到過(guò)。