NCL 读取 融合CMORPH后的出图问题

hhumuer

大家好，我在中国气象数据网上下载了中国自动站与CMORPH降水产品融合的逐时降水量网格数据，翻遍论坛后，用NCL写了个脚本读取这种二进制文件，并把结果绘出来。但是，同一个时间段的数据，我的出图结果（右图）与官网上公布的（左图）大不相同， 如图。把数据按Y轴反转后，图片也还是不对。还希望各位高手能帮忙找找问题，在下不甚感激！

读取数据部分的脚本如下：
tim = "1.grd"
setfileoption("bin","ReadByteOrder","LittleEndian")
rain = fbindirread(tim,0,(/nlat,nlon/),"double")

nlat = 440
nlon = 700
r=new((/nlat,nlon/),"double")
k=0
do i  = 0, nlat-1
  do j = 0, nlon-1
   r(i,j) = rain(nlat-1-i,j)
        k = k+1
  end do
end do

lon = fspan(70.05,139.95,nlon)
lon!0 = "lon"
lon@units = "degrees_east"
lon@long_name = "longitude"

lat = fspan(15.05,58.95,nlat)
lat!0 = "lat"
lat@long_name = "latitude"
lat@units = "degrees_north"

r!0 = "lat"
r!1 = "lon"
r&lat = lat
r&lon = lon
r@units = "mm/hour"
r@long_name = "merged CMORPH data "
r@_FillValue = -999.00

权威

字节顺序是指占内存多于一个字节类型的数据在内存中的存放顺序，通常有小端、大端两种字节顺序。小端字节序指低字节数据存放在内存低地址处，高字节数据存放在内存高地址处；大端字节序是高字节数据存放在低地址处，低字节数据存放在高地址处。基于X86平台的PC机是小端字节序的（原文参考）。在跨系统处理二进制数据流时，要注意这个问题。我就是在处理++服务端的BML（Binary Markup Language）二进制标记数据流时碰到了这个问题。
所以不需要setfileoption("bin","ReadByteOrder","LittleEndian")

freekiller

print 你的降水出来看看

诸事安康

我也出图了下，也这样，感觉是数据反了

好久不见!

请问您解决了吗  ？我的降水都出现在新疆和东北地区，也不知道为什么

好久不见!

zlove璐 发表于 2018-4-16 20:54
我也遇到这个问题了，求问如何解决的？

您好 ，你的cmorph降水出图问题解决了吗

cqlqlq

好人啊！谢谢大家！‘’