【已解决】提问：NCL绘图-兰波托(Lambert)投影数据

Julia12138

数据：浙东南沿海风力数据
投影方式：Lambert投影（经纬度分别为400*400的数组）

                               
登录/注册后可看大图

ncl_filedump后数据内部：

                               
登录/注册后可看大图

数据用panoply打开：



--------------------------------------------------------------
使用NCL画图后：风场分布旋转了90度，如下：↓
--------------------------------------------------
NCL脚本如下：

                               
登录/注册后可看大图

脚本文字版：

f     = addfile ("E:/2020wind/ZJCALMET_d02.2021012101.nc", "r")
  UHTGL     = f->UHTGL(time|:,lvl_h|:,lat|:,lon|:)
  lat2d     = f->lat2d(lat|:,lon|:)
  lon2d     = f->lon2d(lat|:,lon|:)
;-----This is necessary so that data is overlaid properly on map
  UHTGL@lat = lat2d
  UHTGL@lon = lon2d
;---Create plot
  wks = gsn_open_wks("pdf","E:/2020wind/lcnative.pdf")         ; open a workstation
  gsn_define_colormap(wks,"gui_default")       ; choose colormap
  cnres                   =  True
  cnres@gsnAddCyclic      = False
  cnres@gsnDraw           = False              ; Turn these off. We
  cnres@gsnFrame          = False              ; will overlay plots
  cnres@tiMainString           = "Native Lambert Conformal Grid"
  cnres@tiMainFontHeightF      = 0.020              ; smaller title
  cnres@tiMainOffsetYF         = -0.03             ; move title down
  cnres@cnFillOn               = True
  cnres@cnLinesOn              = False

  cnres@mpMinLatF              = min(UHTGL@lat)     ; Add a little bit of
  cnres@mpMinLonF              = min(UHTGL@lon)    ; a margin around
  cnres@mpMaxLatF              = max(UHTGL@lat)     ; the plot.
  cnres@mpMaxLonF              = max(UHTGL@lon)

  cnres@mpAreaMaskingOn            = True
  cnres@mpDataSetName              = "Earth..4"
  cnres@mpDataBaseVersion          = "MediumRes" ; or "Ncarg4_1"
  cnres@mpOutlineSpecifiers        = (/"Zhejiang"/)
  cnres@mpMaskAreaSpecifiers       = (/"Zhejiang"/)
  cnres@pmTickMarkDisplayMode = "Always"

  cnres@tfDoNDCOverlay = True

  contour_fill_plot = gsn_csm_contour_map_ce(wks,UHTGL(0,0,:,:),cnres)
  draw(contour_fill_plot)

                               
登录/注册后可看大图

请问，有哪位朋友知道如何解决这个数据和地图不对应的问题嘛？谢谢~
test.ncl (2.46 KB, 下载次数: 0)

2021-3-25 19:48 上传

点击文件名下载附件

尽头的尽头

你用panoply画图的时候选择经纬度的，不要选择XY的，ncl的结果应该是没有问题的

Julia12138

尽头的尽头 发表于 2021-3-25 20:34
你用panoply画图的时候选择经纬度的，不要选择XY的，ncl的结果应该是没有问题的

您好，我重新选了一下，只有用using lon for X axis and lat for Y axis 图才能和ncl出来的对应上!
之前我一直都是选默认的georeferenced longitude-latitude color contour plot，您是这个意思嘛？

尽头的尽头

using lon for X axis and lat for Y axis ，这个才是正确的

Julia12138

尽头的尽头 发表于 2021-3-26 08:43
using lon for X axis and lat for Y axis ，这个才是正确的

嗯好哒！谢谢您。