爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 17162|回复: 4

[其他] 【已解决】提问:NCL绘图-兰波托(Lambert)投影数据

[复制链接]
发表于 2021-3-25 19:42:53 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册

x
本帖最后由 Julia12138 于 2021-3-26 15:19 编辑

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

                               
登录/注册后可看大图


ncl_filedump后数据内部:
微信图片_20210325194812.png


                               
登录/注册后可看大图


数据用panoply打开:
微信图片_20210325193358.png


--------------------------------------------------------------
使用NCL画图后:风场分布旋转了90度,如下:↓ 微信图片_20210325193614.png
--------------------------------------------------
NCL脚本如下:
微信图片_20210325193915.png


                               
登录/注册后可看大图

脚本文字版:

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)
密码修改失败请联系微信:mofangbao
发表于 2021-3-25 20:34:32 | 显示全部楼层
你用panoply画图的时候选择经纬度的,不要选择XY的,ncl的结果应该是没有问题的
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

 楼主| 发表于 2021-3-25 21:47:38 | 显示全部楼层
尽头的尽头 发表于 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,您是这个意思嘛?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

发表于 2021-3-26 08:43:34 | 显示全部楼层
using lon for X axis and lat for Y axis ,这个才是正确的
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

 楼主| 发表于 2021-3-26 15:17:32 来自手机 | 显示全部楼层
尽头的尽头 发表于 2021-3-26 08:43
using lon for X axis and lat for Y axis ,这个才是正确的

嗯好哒!谢谢您。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表