华北雷达组合反射率拼图数据RADA_NCN_DOR_L2_CR-2020xxxxxxxx.zip头文件格式

几杯清澈

有每6分钟一次的华北雷达组合反射率拼图数据，文件名为RADA_NCN_DOR_L2_CR-2020xxxxxxxx.zip，压缩包里的文件为20xxxxxxxx.dat，想绘制华北雷达组合反射率拼图，请教大家该如何读取数据？或者文件头格式是什么？我先随便试了下，读出了字符串“diamond 131 ”，通过搜索只能搜到“SWAN系统中说明书中的三维拼图数据格式”，试了一下发现格式不符合，因此想请教大家。

loongth

diamond文件头应该是MICAPS里面的 不过diamond 13似乎才是雷达拼图用的文件头

几杯清澈

loongth 发表于 2025-1-20 13:25
diamond文件头应该是MICAPS里面的 不过diamond 13似乎才是雷达拼图用的文件头

是的，我搜索到了MICAPS的用户手册，里面的第13类数据格式包含了雷达拼图，格式如下：

文件头：（文件头为TEXT格式）
diamond （8个字符）13 （3个字符）数据说明（40个字符）
年 （5个字符）月 （3个字符）日（3个字符）时次 （3个字符）X方向图象大小 （5个字符）Y方向图象大小 （5个字符）图象左下角经度坐标（8个字符）图象左下角纬度坐标 （8个字符）投影方式（2个字符，1-lambert  2-mecator  3-北半球  4-南半球） 放缩系数（5个字符）图象类型（2个字符，1—红外云图 2—雷达拼图 3—地形图 4—可见光云图 5—水汽图）象素值与相应物理量对照表文件名（12个字符）中心经度（8个字符）、中心纬度（8个字符）
（文件头部分共128个字符）
注：目前9210工程通讯系统传输的GMS5云图左下角经纬度（86.4，－1.3）
数据：（数据格式为二进制数据）
一个象素点占一个字节，先沿X方向后Y方向。
例子：
diamond 13 97年2月10日08时红外云图 1997 02 10 08 1280 1024 86.4 -1.3  1 1.0 1 ir.dat cloud
下面是256级灰度表示值......

但我按上述格式读取也不行，仍然报错，所以我觉得我的数据应该不是这一类。

几杯清澈

我使用如下python代码读取了前42个字节：
with open(read_filename, 'rb') as file:
    header = file.read(42).decode('ascii')
    print(header)
输出结果如下：

edwardli

一． D131格点数据格式 131格式经过数次修订，目前大多数数据仍未1.0格式，但是显示引擎已经能够解析2.0格式的数据，具体应用可以根据情况使用相应版本的格式，主要是从后面的Reserved中修改为对应版本的大小。保证文件头格式总长度1024个字节。 char ZonName[12];// diamond 131 12个字节 char DataName[38];//数据说明(例如 2008年5月19日雷达三维拼图)38个字节 charFlag[8];// 文件标志，"swan" charVersion[8];// 数据版本号，"2.0" 2.0表示为版本为2.0 unsigned short int year;//2008 两个字节 unsigned short int month;//05  两个字节 unsigned short int day;//19    两个字节 unsigned short int hour;//14   两个字节 unsigned short int minute;//31 两个字节 unsigned short int interval ;  //两个字节         unsigned short int XNumGrids;//1300 两个字节 unsigned short int YNumGrids;//800 两个字节 unsigned short int ZNumGrids;//20  两个字节 int RadarCount; //拼图雷达数 四个字节 float StartLon; //网格开始经度（左上角） 四个字节 float StartLat; //网格开始纬度（左上角） 四个字节 float CenterLon;//网格中心经度 四个字节 float CenterLat;//网格中心纬度 四个字节 float XReso;//经度方向分辨率 四个字节 float YReso;//纬度方向分辨率 四个字节 float ZhighGrids[40];//垂直方向的高度（单位km）数目根据ZnumGrids而得（最大40层） 160个字节。 char RadarStationName[20][16];    //相关站点名称,20*16字节 float  RadarLongitude[20];      //相关站点所在经度，单位：度， 4*20字节 float  RadarLatitude[20];       //相关站点所在纬度，单位：度， 4*20字节 float  RadarAltitude[20];      //相关站点所在海拔高度，单位：米， 4*20字节 unsigned char    MosaicFlag[20];    //该相关站点数据是否包含在本次拼图中，未包含:0，包含:1, 20字节 short  m_iDataType;   //0-unsigned char 1-char 2-unsigned short 3-short 4 //每一层的向量数，版本=2.0 short m_iLevelDimension; float  offset; float  scale char reserved[160]  保持头有1024个字节 接下来是数据块，从底层到高层进行排列共ZnumGrids层。每层的数据从起始点（左上角）开始，按维向（纬度y）减小写每行的经向（经度x增大）数据。数据存储按照m_iDataTyp表示，分为不同的数据类型，存储值和真实值之间的换算关系为 data（真实值） =scale* value(存储值)+ offset     一般情况下，反射率因子采用 unsigned char存储 offset = -33 scale = 0.5 A

edwardli

f = cinrad.io.SWAN(nFiles)

MeteoInfo

忘了从哪里要到的diamond 131格式的说明文档了，共享在这里，正如edwardli所说，这个格式可能并没有定版，是否和你的数据对应不得而知。有可能的话你是否可以共享一个示例数据文件，可以共同探讨一下。 SWAN三维拼图格式说明.pdf (64.19 KB, 下载次数: 0)

11 分钟前 上传

点击文件名下载附件

根据上面的格式文档，在MeteoInfo中实现了对micaps 131格式的解码（MeteoInfo/meteoinfo-data/src/main/java/org/meteoinfo/data/meteodata/micaps/MICAPS131DataInfo.java at master · meteoinfo/MeteoInfo）。