begin ncol = 6 ;臺風(fēng)參數(shù) nrow = 31 ;時次總數(shù) nbin = 6 ;已知該該氣旋共經(jīng)歷了6個等級的演變 ;讀入臺風(fēng)資料 data = asciiread("NEPARTAK.txt",(/nrow,ncol/),"integer") ;/31,6/ 31行6列，integer整數(shù)類型 ;;數(shù)據(jù)讀取函數(shù)總結(jié)，asciiread讀取cav、txt文件，addfile讀取nc、hdf、grib、shp文件，fbindrread讀取bin、ctl、dat文件 amp = data(:,1) ; 強(qiáng)度 lat = data(:,2)/10. ; 緯度 lon = data(:,3)/10. ; 經(jīng)度 wks = gsn_open_wks("pdf","plot-typhoon") gsn_define_colormap(wks, "rainbow") cmap = read_colormap_file("rainbow") ;讀取色板中的顏色 ncol = dimsizes(cmap(:,0)) ;色板的顏色總數(shù) res = True res@gsnFrame = False res@gsnDraw = False res@mpMinLatF = 5 res@mpMaxLatF = 35 res@mpMinLonF = 105 res@mpMaxLonF = 150 res@mpOutlineOn = True res@pmTickMarkDisplayMode = "Always" ;坐標(biāo)軸上添加度符號 map = gsn_csm_map(wks,res) ;繪制出地圖底圖 ;--添加任意折線，表示出臺風(fēng)路徑 lnres = True lnres@gsLineColor = "black" lnres@gsLineThicknessF = 1.5 lnres@gsLineDashPattern = 0 str = unique_string("string") map@$str$ = gsn_add_polyline(wks,map,lon,lat,lnres) ;繪制Nepartak txres = True txres@txFontHeightF = 0.01 str = unique_string("string") ;字符串形式 map@$str$ = gsn_add_text(wks,map,"1601Nepartak",135,22,txres) ;添加標(biāo)識，表示出每個時次的臺風(fēng)中心位置 mkres = True mkres@gsMarkerSizeF = 0.01 ;標(biāo)識大小，在0~1范圍內(nèi)設(shè)置 mkres@gsMarkerIndex = 3 ;標(biāo)識形狀 mkres@gsMarkerThicknessF = 1.5 do i = 0, nrow-1 mkres@gsMarkerColor = toint(floor(ncol/nbin*amp(i)-1)) ;在整個色板中準(zhǔn)等距離的挑選顏色 str = unique_string("string") ;字符串形式 map@$str$ = gsn_add_polymarker(wks,map,lon(i),lat(i),mkres);繪制臺風(fēng)中心位置標(biāo)識 end do ;在圖形正下方添加等級圖例 x_ndc = (/0.25,0.35,0.45,0.55,0.65,0.75/) ;相對坐標(biāo)x方向的位置 y_ndc = (/0.14,0.14,0.14,0.14,0.14,0.14/) ;相對坐標(biāo)y方向的位置 t_ndc = tostring(ispan(1,6,1)) ;假定已知等級為1-6 labels = (/"10.8-17.1m/s","17.2-24.4m/s","24.5-32.6m/s","32.7-41.4m/s","41.5-50.9m/s",">=51.0m/s"/) mkres@gsMarkerSizeF = 0.01 ;標(biāo)識的大小，在 0~1范圍內(nèi)設(shè)置。同上 mkres@gsMarkerThicknessF = 1.5 ;標(biāo)識的粗細(xì)。同上 do i = 0,5 ; 繪制每個等級的圖例,它由標(biāo)識和文本構(gòu)成 mkres@gsMarkerColor = toint(floor(ncol/nbin*(i+1)-1)) gsn_polymarker_ndc(wks,x_ndc(i),y_ndc(i),mkres) gsn_text_ndc(wks,t_ndc(i),x_ndc(i)+0.02,y_ndc(i),txres) ;將文本向右邊移動0.02 gsn_text_ndc(wks,labels(i),x_ndc(i),y_ndc(i)-0.02,txres);將文本向下邊移動0.02 end do ;根據(jù)經(jīng)緯度位置繪制多邊形 x_lonlat = 117 ;中心位置 y_lonlat = 25 box_lon = (/x_lonlat-6,x_lonlat+6,x_lonlat+6,x_lonlat-6,x_lonlat-6/) box_lat = (/y_lonlat+5,y_lonlat+5,y_lonlat-5,y_lonlat-5,y_lonlat+5/) gonres = True gonres@gsFillIndex = 3 ;用形狀填充，默認(rèn)值是0表示顏色填充，3是陰影 gonres@gsFillColor = "red" dum = gsn_add_polygon(wks,map,box_lon,box_lat,gonres) ;以下繪制第23時次(二次登陸前)的大致風(fēng)圈范圍：以其中心為圓心，以0.05長度為半徑的圓形區(qū)域 xout_ndc = lon yout_ndc = lat datatondc(map,lon,lat,xout_ndc,yout_ndc) ;將經(jīng)緯度數(shù)值轉(zhuǎn)換威單位坐標(biāo)系中的坐標(biāo)值 degrad = 0.017453292519943 ;Π/180=3.1415926/180 degrees = ispan(0,360,1) ;創(chuàng)建361點(diǎn) xcos = cos(degrad*degrees) ;各個角度的余弦值 xsin = sin(degrad*degrees) ;各個角度的正弦值 ;圓心點(diǎn)位置及半徑 xcenter = xout_ndc(22) ;第22時次中心所在x位置 ycenter = yout_ndc(22) ;第22時次中心所在y位置 radius = 0.05 ;計(jì)算出各個點(diǎn)在單位坐標(biāo)系中的坐標(biāo)位置 xc = xcenter + (radius * xcos) yc = ycenter + (radius * xsin) ;設(shè)置多邊形的繪圖參數(shù)，根據(jù)單位坐標(biāo)系中的位置繪制多邊形 lnres = True lnres@gsFillColor = 46 lnres@gsFillOpacityF = 0.5 ;50%透明 lnres@gsFillLineThicknessF = 2. ;多邊形邊框粗細(xì) gsn_polygon_ndc(wks,xc,yc,lnres) ?? ;;添加汕頭位置 lon_st = 116.7 lat_st = 23.5 mres = True mres@gsMarkerColor = "black" mres@gsMarkerSizeF = 0.02 mres@gsMarkerIndex = 3 dum1 = gsn_add_polymarker(wks,map,lon_st,lat_st,mres) ;添加SWATOW txtres = True txtres@txFontHeightF = 0.005 txtres@txFontColor = "black" dum2 = gsn_add_text(wks,map,"SWATOW",lon_st,lat_st+1.5,txtres) ngezlogo(wks) draw(map) frame(wks)end

附數(shù)據(jù)文件：可以直接復(fù)制，保存為NEPARTAK.txt

2016070212 1 82 1453 1000 152016070218 1 85 1452 1000 152016070300 2 88 1450 998 182016070306 2 92 1448 998 182016070312 2 99 1442 995 202016070318 2 109 1435 995 202016070400 2 118 1423 995 202016070406 2 128 1410 992 232016070412 3 135 1396 990 252016070418 3 145 1383 985 282016070500 4 153 1366 975 332016070506 4 163 1351 960 402016070512 6 171 1336 940 522016070518 6 180 1319 915 622016070600 6 187 1302 905 682016070606 6 195 1285 895 722016070612 6 201 1270 895 722016070618 6 207 1256 895 722016070700 6 213 1244 895 722016070706 6 216 1234 895 722016070712 6 219 1227 895 722016070718 6 224 1218 915 622016070800 5 226 1207 950 452016070806 4 232 1201 965 382016070812 4 231 1198 975 332016070818 3 233 1196 985 252016070900 2 243 1191 990 232016070906 2 247 1187 992 202016070912 2 251 1181 995 182016070918 1 258 1169 996 152016071000 1 264 1165 998 13

總結(jié)

以上是生活随笔為你收集整理的cesium 经纬度绘制点_NCL绘制2016年1号台风(Nepartak)的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。