nc4格式数据处理求助

Linux · 发表于 2021-2-7 12:19:29

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

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

x

在论坛里拜读了诸位大神的帖子，然小生愚钝，还是无法弄清楚nc4格式的数据该如何处理。跪求大神指导！
我的想法如下：首先在grads官网查了xdfopen命令，然后写出相应描述文件。但是失败了！！！
附件里是数据文件和最新的ncdump.exe。

下面是利用ncdump命令获取的数据文件头文件信息：

netcdf \3B42_Daily.19980101.7 {
dimensions:
lon = 1440 ;
lat = 400 ;
variables:
float precipitation(lon, lat) ;
string precipitation:units = "mm" ;
string precipitation:long_name = "Daily accumulated precipitation (combined microwave-IR) estimate with gauge calibration over land" ;
string precipitation:coordinates = "lat lon" ;
precipitation:_FillValue = -9999.9f ;
byte precipitation_cnt(lon, lat) ;
string precipitation_cnt:units = "count" ;
string precipitation_cnt:long_name = "Count of valid 3-hr (combined microwave-IR) precipitation retrievals for the day" ;
string precipitation_cnt:coordinates = "lat lon" ;
float IRprecipitation(lon, lat) ;
string IRprecipitation:units = "mm" ;
string IRprecipitation:long_name = "Daily accumulated precipitation (IR) with HQ calibration" ;
string IRprecipitation:coordinates = "lat lon" ;
IRprecipitation:_FillValue = -9999.9f ;
byte IRprecipitation_cnt(lon, lat) ;
string IRprecipitation_cnt:units = "count" ;
string IRprecipitation_cnt:long_name = "Count of valid 3-hr IR precipitation retrievals for the day" ;
string IRprecipitation_cnt:coordinates = "lat lon" ;
float HQprecipitation(lon, lat) ;
string HQprecipitation:units = "mm" ;
string HQprecipitation:long_name = "Daily accumulated High Quality (microwave) precipitation" ;
string HQprecipitation:coordinates = "lat lon" ;
HQprecipitation:_FillValue = -9999.9f ;
byte HQprecipitation_cnt(lon, lat) ;
string HQprecipitation_cnt:units = "count" ;
string HQprecipitation_cnt:long_name = "Count of valid 3-hr HQ (microwave) precipitation retrievals for the day" ;
string HQprecipitation_cnt:coordinates = "lat lon" ;
float randomError(lon, lat) ;
string randomError:units = "mm" ;
string randomError:long_name = "Daily total error of combined microwave-IR precipitation estimate" ;
randomError:_FillValue = -9999.9f ;
byte randomError_cnt(lon, lat) ;
string randomError_cnt:units = "count" ;
string randomError_cnt:long_name = "Count of valid 3-hr error retrievals for the day" ;
string randomError_cnt:coordinates = "lat lon" ;
float lon(lon) ;
string lon:units = "degrees_east" ;
string lon:long_name = "Longitude" ;
float lat(lat) ;
string lat:units = "degrees_north" ;
string lat:long_name = "Latitude" ;
// global attributes:
:BeginDate = "1998-01-01" ;
:BeginTime = "01:30:00.000Z" ;
:EndDate = "1998-01-02" ;
:EndTime = "01:29:59.999Z" ;
:FileHeader = "StartGranuleDateTime=1998-01-01T01:30:00.000Z;\nStopGranuleDateTime=1998-01-02T01:29:59.999Z" ;
:InputPointer = "3B42.19980101.03.7.HDF;3B42.19980101.06.7.HDF;3B42.19980101.09.7.HDF;3B42.19980101.12.7.HDF;3B42.19980101.15.7.HDF;3B42.19980101.18.7.HDF;3B42.19980101.21.7.HDF;3B42.19980102.00.7.HDF" ;
:title = "TRMM (TMPA) L3 Daily 0.25 x 0.25 degree (TRMM_3B42_Daily)" ;
:ProductionTime = "2016-07-20T13:48:20.000Z" ;
}

复制代码

根据上述信息，我写出了描述文件如下：
DSET ^3B42_Daily.19981213.7.nc4
TITLE TRMM (TMPA) L3 Daily 0.25 x 0.25 degree (TRMM_3B42_Daily)
UNDEF -9999.9
XDEF lon 1440 LINEAR -179.875 0.25
YDEF lat 400 LINEAR -59.875 0.25
TDEF 1 LINEAR 0Z13DEC1998 6hr
VARS 1
precipition=>pcp 1 99 rain
ENDVARS

然后xdfopen 命令打开描述文件，出现了错误如下，请问该如何处理呢？

苏拉苏拉 · 发表于 2021-2-7 16:09:20

本帖最后由苏拉苏拉于 2021-2-7 16:55 编辑

tdef后面是不是少time

Linux · 发表于 2021-2-7 20:28:50

苏拉苏拉发表于 2021-2-7 16:09
tdef后面是不是少time

加上了time,虽然没有上面的问题了，但出现了一堆not supported。最后一行是sdf descriptor file test. ctl was not successfully opened&parsed

苏拉苏拉 · 发表于 2021-2-8 01:16:32

Linux 发表于 2021-2-7 20:28
加上了time,虽然没有上面的问题了，但出现了一堆not supported。最后一行是sdf descriptor file test. ct ...

endvars 后面或另起一行加上；（英文分号）试试

Linux · 发表于 2021-2-9 09:11:47

苏拉苏拉发表于 2021-2-8 01:16
endvars 后面或另起一行加上；（英文分号）试试

现在搞定了，能打开。
但是因为文件的dimensions只有两个： lon = 1440 ; lat = 400 。这样是不是就不能定义tdef了，也就不能进行批处理了？

苏拉苏拉 · 发表于 2021-2-9 13:29:40

Linux 发表于 2021-2-9 09:11
现在搞定了，能打开。
但是因为文件的dimensions只有两个： lon = 1440 ; lat = 400 。这样是不是就不能 ...

时间维度也不能算没有吧，只是只有1个，tdef在ctl里还是要写的。批处理的话也是可以的，文件名有规律就行。

Linux · 发表于 2021-2-10 06:47:03

苏拉苏拉发表于 2021-2-9 13:29
时间维度也不能算没有吧，只是只有1个，tdef在ctl里还是要写的。批处理的话也是可以的，文件名有规律就行 ...

我要是加上时间维度，ctl就打不开了。有没有time都是，各种打不开。请问您是怎么操作的？能针对附件数据给个描述文件吗

苏拉苏拉 · 发表于 2021-2-10 13:49:25

Linux 发表于 2021-2-10 06:47
我要是加上时间维度，ctl就打不开了。有没有time都是，各种打不开。请问您是怎么操作的？能针对附件数据 ...

我没弄过，只是根据你给的错误提示猜问题出在哪。TDEF time 1 LINEAR 13DEC1998 1dy 这样试试呢。你多试试啊，根据报错。加油。

		自动登录	找回密码
密码			立即注册

[脚本编辑] nc4格式数据处理求助

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