批量读取.nc转为.grd格式后画图出现问题：偏移+数值错误

菟子君

新手求助：我将2天的TRMM资料（.nc格式）批量读取后转为.grd格式画图，发现图像出现偏移以及有数值上的差错。
                            
上为直接用nc文件画的第一个时刻，下为用批量读取再转换后的grd画的第一个时刻的图。

批量读取的ctl:
dset E:/2016/3B42.%y4%m2%d2.%h2.7.nc
undef -9999.9
title TRMM rainfall - 3 hourly
options little_endian template
xdef 1440 linear -179.875 0.25
ydef 400 linear -49.875 0.25
zdef 1 linear 0 1
tdef 16 linear 00Z01JAN2016 3hr
vars 6
pcp=>pcp  0 t,y,x   precipitation [mm hr**-1]
err=>err  0 t,y,x   relative error [mm hr**-1]
source=>source   0 t,y,x   source
hqp=>hqp   0 t,y,x   high quality precipitation [mm hr**-1]
irp=>irp    0 t,y,x   IR precipitation [mm hr**-1]
satobstime=>satobstime    0 t,y,x   observation time [minutes]
ENDVARS

转换用的gs:
'reinit'
'open E:/trmm.ctl'
'set gxout fwrite'
'set fwrite E:/2016.grd'
'set x 1 1440'
'set y 1 400'
'set z 0'
it=1
while (it<=16)
'set t 'it''
'd pcp'
it=it+1
endwhile
'disable fwrite'
;

用grd画图用的ctl和gs:
dset E:/2016.grd
undef -9999.9
title TRMM rainfall - 3 hourly
xdef 1440 linear -179.875 0.25
ydef 400 linear -49.875 0.25
zdef 1 linear 0 1
tdef 16 linear 00Z01JAN2016 3hr
vars 1
pcp=>pcp  0 t,y,x   precipitation [mm hr**-1]
endvars


'reinit'
'open e:/16/trmm.ctl'
'set x 1 1440'
'set y 1 400'
'set z 0'
'set t 1'
'd pcp'
;

用nc画图时的gs:
'reinit'
'sdfopen e:/2016/3B42.20160101.00.7.nc'
'set x 1 1440'
'set y 1 400'
'd pcp'
;

附nc文件的描述：

pearal

转换后的缺省值应该为-9.99e+08吧~试试看？

菟子君

pearal 发表于 2016-11-15 19:55
转换后的缺省值应该为-9.99e+08吧~试试看？

改了缺省值图像还是没有变化，错误还在

river

你先用你的批量描述文件画画图，看看和原始的对的上不

376394482

谢谢楼主，楼主辛苦了

菟子君

river 发表于 2016-11-15 21:38
你先用你的批量描述文件画画图，看看和原始的对的上不

已解决。是批量描述的时候没有dtype的缘故，谢谢您的建议~