求各位大佬解，已定义坐标，但作图还是经纬度坐标信息警告

395033193

警告信息如下：
See http://www.ncl.ucar.edu/ for more details.
(0)     check_for_y_lat_coord: Warning: Data either does not contain
(0)     a valid latitude coordinate array or doesn't contain one at all.
(0)     A valid latitude coordinate array should have a 'units'
(0)     attribute equal to one of the following values:
(0)         'degrees_north' 'degrees-north' 'degree_north' 'degrees north' 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg north'
(0)     check_for_lon_coord: Warning: Data either does not contain
(0)     a valid longitude coordinate array or doesn't contain one at all.
(0)     A valid longitude coordinate array should have a 'units'
(0)     attribute equal to one of the following values:
(0)         'degrees_east' 'degrees-east' 'degree_east' 'degrees east' 'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'



代码如下：

begin

raw_2021_summer_anomaly=asciiread("/mnt/e/research/ALL_PROJECTS/2023-1-lanina/ncl/data/2021_summer_anomaly_dropna.csv",(/2249,5/),"float")

plon=raw_2021_summer_anomaly(:,1)
plat=raw_2021_summer_anomaly(:,2)

LAT=fspan(15, 58, 132)
;nLAT = dimsizes(LAT)
LON=fspan(70, 138, 204)

grid_pre=natgrid(plon, plat, raw_2021_summer_anomaly(:,4), LON, LAT)

grid_pre@units="%"
grid_pre@long_name = "precipitation anomaly"
grid_pre!0        = "longitude"
grid_pre!1        = "latitude"
grid_pre&longitude=LON
grid_pre&latitude=LAT
grid_pre&longitude@units="degrees_east"
grid_pre&latitude@units="degrees_north"

  wks_type = "png"
  wks_type@wkWidth  = 6000
  wks_type@wkHeight = 6000
  wks  = gsn_open_wks(wks_type,"/mnt/e/research/ALL_PROJECTS/2023-1-lanina/ncl/pictures/ceshi")
  res = True
  res@gsnDraw        = False
  res@gsnFrame       = False
  res@gsnAddCyclic  = False      
  res@mpDataSetName         = "/mnt/e/research/files/NCL-Chinamap-master/database/Earth..4"   
  res@mpDataBaseVersion     = "MediumRes"  
  res@mpOutlineOn           = True         
  res@mpOutlineSpecifiers   = (/"China","China:Provinces"/)        
  res@mpMinLatF             =  4         
  res@mpMaxLatF             =  56
  res@mpMinLonF             = 72
  res@mpMaxLonF             = 136
  res@gsnMajorLonSpacing      = 10
  res@gsnMinorLonSpacing      = 10
  res@gsnMajorLatSpacing      = 10
  res@gsnMinorLatSpacing      = 10
  res@mpGeophysicalLineThicknessF= 2.     
  res@mpNationalLineThicknessF= 1.5         
  res@mpAreaMaskingOn = True   
  res@mpMaskAreaSpecifiers = (/"China","China:Provinces"/)   
  res@mpLandFillColor = "White"
  res@mpOceanFillColor = "White"   
  res@mpInlandWaterFillColor = "White"   
  res@cnFillOn      = True
  res@cnLinesOn     = False         
  res@cnLineLabelsOn = False
  res@cnFillDrawOrder = "PreDraw"         
  res@gsnLeftString = " "                     
  res@gsnRightString = " "                     
  res@cnLevelSelectionMode = "ExplicitLevels"      
  res@cnLevels    =ispan(-10,10,1)*0.1
  res@cnFillPalette   = "WhiteYellowOrangeRed"
  res@lbLabelBarOn = True      
  plot = gsn_csm_contour_map(wks,grid_pre({72:136},{16:56}),res)

  draw(plot)                                 
  frame(wks)                                 
end
exit

1277105864

建议你给数据赋了属性之后，printVarSummary出来看一下，数据的属性有没有赋成功

皮皮皮皮皮皮

你好，现在解决了吗？是怎么解决的呢？我也遇到了这个问题

395033193

1277105864 发表于 2023-6-22 22:48
建议你给数据赋了属性之后，printVarSummary出来看一下，数据的属性有没有赋成功

赋值属性都是对的，就是不行，换了插值函数就行了

jinj

是不是没给经纬度定义单位
lon@units="degrees_east"
lat@units="degrees_north"