- 积分
- 5156
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-4-18
- 最后登录
- 1970-1-1
|
NCL
系统平台: |
|
问题截图: |
- |
问题概况: |
NCL画图时不设置投影方式经纬度范围没有问题,但设置为兰伯特投影后就有偏移 |
我看过提问的智慧: |
看过 |
自己思考时长(天): |
3 |
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
begin
data=fbindirread("t-mon.dat", -1, (/239,163,283/), "float")
printVarSummary(data)
lon = fspan(69.75,140.25,283)
lat = fspan(14.75,55.25,163)
time = fspan(1,239,239)
;time = tostring(time1)
lat!0 = "lat"
lon!0 = "lon"
time!0 = "time"
lat@units = "degrees_north"
lon@units = "degrees_east"
time@units= "days since 1998-04-01 00:00:00"
time@long_name = "Time"
lat&lat = lat
lon&lon = lon
time&time = time
data!1 = "lat"
data!2 = "lon"
data!0 = "time"
data&time = time
data&lat = lat
data&lon = lon
data@_FillValue = -9999
data@description=""
data@units = ""
wks_type="eps"
wks=gsn_open_wks(wks_type,"obs-t")
sum_t=data(66,:,:)
;printVarSummary(sum_t)
sum_t = sum_t / 30
res= True
res@gsnAddCyclic = False
res@mpDataSetName = "Earth..4" ; This new database contains
; divisions for other countries.
res@mpDataBaseVersion = "MediumRes" ; Medium resolution database
res@mpOutlineOn = True ; Turn on map outlines
res@mpOutlineSpecifiers = (/"China:states","Taiwan"/) ;China:states
res@mpProjection = "LambertConformal" ; Lambert
res@mpLambertMeridianF = 104.06
res@mpLimitMode = "LatLon"
res@mpLambertParallel1F = 30 ;Default: .001
res@mpLambertParallel2F = 60 ;Default: 89.999
res@mpMinLatF = 30 ; Asia limits
res@mpMaxLatF = 31.5
res@mpMinLonF = 102.45
res@mpMaxLonF = 104.45
res@cnFillOn=True
res@cnLinesOn=True
res@cnLineLabelsOn = False
res@cnFillDrawOrder= "PostDraw"
res@cnInfoLabelOn=False
res@cnLineColor = "Black"
res@cnLineLabelBackgroundColor = -1
res@gsnContourLineThicknessesScale = 2.0
res@pmTickMarkDisplayMode = "Always"
res@tmXTOn = False
res@tmYROn = False
res@cnLevelSelectionMode="ExplicitLevels"
res@cnLevels = (/8,10,12,14,16,19,20,22,24,26,28/)
;res@cnLevels = (/100,150,200,250/)
res@lbLabelBarOn = True ;LabelBar显示
map = gsn_csm_contour_map(wks,sum_t,res)
end
上面是我的脚本,区域设置是30-31.5,102.45到104.45,不设置投影方式出来没问题,但是设为兰伯特投影后就有偏差了,请问各位,这该怎么解决呢?
|
|