3830doppler数据格式.doc

CINRAD/CC/CCJ雷达原始数据格式声明：本数据格式适用CINRAD/CC和CINRAD/CCJ安徽四创电子股份有限公司保留最终解释权本文档仅供内部交流，请不要发表文件名：平扫（PPI）：NNNNYYDDHH.MMP高扫（RHI）：NNNNYYDDHH.MMR体扫（VPPI）：NNNNYYDDHH.MMV单库：NNNNYYDDHH.MMF等高PPI（CAPPI）：NNNNYYDDHH.MMCN：年 Y：月 D：日 H：时 M：分数据组成:整个数据由文件头(1个)和基于极坐标系的原始数据(512个径向)组成文件头:/tagWEATHERRADAR雷达信息的结构typedef struct tagWEATHERRADARchar cFileType16;/3830数据标识(CINRADC)char cCountry30;/国家名char cProvince20;/省名char cStation40;/站名char cStationNumber10;/区站号char cRadarType20;/雷达型号char cLongitude16;/天线所在经度char cLatitude16;/天线所在纬度long lLongitudeValue;/具体经度long lLatitudeValue;/具体纬度long lHeight;/天线海拔高度short sMaxAngle;/地物阻挡最大仰角short sOptAngle;/最佳观测仰角unsigned char ucSYear1;/观测开始时间的年千百位(19-20)unsigned char ucSYear2;/观测开始时间的年十个位(00-99)unsigned char ucSMonth;/观测开始时间的月(1-12)unsigned char ucSDay;/观测开始时间的日(1-31)unsigned char ucSHour;/观测开始时间的时(0-23)unsigned char ucSMinute;/观测开始时间的分(0-59)unsigned char ucSSecond;/观测开始时间的秒(0-59)unsigned char ucTimeFrom;/时间来源 0-计算机时钟(1天内未对时)/ 1-计算机时钟(1天内已对时)/ 2-GPS/ 3-其它unsigned char ucEYear1;/观测结束时间的年千百位(19-20)unsigned char ucEYear2;/观测结束时间的年十个位(00-99)unsigned char ucEMonth;/观测结束时间的月(1-12)unsigned char ucEDay;/观测结束时间的日(1-31)unsigned char ucEHour;/观测结束时间的时(0-23)unsigned char ucEMinute;/观测结束时间的分(0-59)unsigned char ucESecond;/观测结束时间的秒(0-59)unsigned char ucScanMode;/扫描方式 1-RHI/ 10-PPI和ZPPI / 1XX=VPPI(XX为扫描圈数)unsigned long ulSmilliSecond;/以微秒为单位表示的秒的小数位unsigned short usRHIA;/RHI所在的方位角(0.01度为单位)/ PPI和VPPI时为FFFFshort sRHIL;/RHI所在的最低仰角(0.01度为单位) /PPI和VPPI时为FFFFshort sRHIH;/RHI所在的最高仰角(0.01度为单位) /PPI和VPPI时为FFFFunsigned short usEchoType;/回波类型 0x405a-Z 0x406a-V 0x407a-W / 0x408a-ZVW三要素unsigned short usProdCode;/数据类型 0x8001-PPI数据 0x8002-RHI数据 / 0x8003-VPPI数据 0x8004-单强度PPI数据 /0x8005-CAPPI数据unsigned char ucCalibration;/标校状态 0-无 1-自动 2-1星期内人工/ 3-1月内人工unsigned char remain13;/保留字unsigned char remain2660;/保留字,放VPPISCANPARAMETER数据/该结构的说明见后long lAntennaG;/天线增益(0.001dB)long lPower;/峰值功率(瓦)long lWavelength;/波长(微米)unsigned short usBeamH;/垂直波束宽度(秒)unsigned short usBeamL;/水平波束宽度(秒)unsigned short usPolarization;/极化状态 0-水平 1-垂直 2-双偏振 / 3-圆偏振 4-其它unsigned short usLogA;/对数动态范围(0.01dB)unsigned short usLineA;/线性动态范围(0.01dB)unsigned short usAGCP;/AGC延迟量(微秒)unsigned short usFreqMode;/频率方式1-单重复频率 2-双重复频率3:2 /3-双重复频率4:3unsigned short usFreqRepeat;/重复频率unsigned short usPPPPulse;/PPP脉冲数unsigned short usFFTPoint;/FFT间隔点数unsigned short usProcessType;/信号处理方式1-PPP2-全程FFT/3-单库FFTunsigned char ucClutterT;/杂波消除阀值(即STC)char cSidelobe;/第一旁瓣(dB)unsigned char ucVelocityT;/速度门限unsigned char ucFilderP;/地物消除方式0-无1-IIR滤波器1/2-IIR滤波器23-IIR滤波器3/4-IIR滤波器4unsigned char ucNoiseT;/噪声消除阀值(即强度门限)unsigned char ucSQIT;/SQI门限unsigned char ucIntensityC;/DVIP强度值估算采用的通道 / 1-对数通道 2-线性通道unsigned char ucIntensityR;/强度值估算是否距离订正/ 0-无(dB) 1-已订正(dBZ)unsigned char ucCalNoise;/噪声系数标定值unsigned char ucCalPower;/发射功率标定值unsigned char ucCalPulseWidth;/脉冲宽度标定值unsigned char ucCalWorkFreq;/工作频率标定值unsigned char ucCalLog;/对数斜率标定值char remain392;/保留字unsigned long int liDataOffset;/数据偏移地址WEATHERRADAR;/tagVPPISCANPARAMETER仰角层的结构typedef struct tagVPPISCANPARAMETERunsigned short usMaxV;/最大可测速度(厘米/秒)unsigned short usMaxL;/最大可测距离(10米)unsigned short usBindWidth;/库长(米)unsigned short usBinNumber;/每径向库数unsigned short usRecordNumber;/本圈径向数unsigned short usArotate;/本圈转速(0.01度/秒)unsigned short usPrf1;/本圈第一次重复频率(0.1Hz)对应单重频或双重频的高者unsigned short usPrf2;/本圈第二次重复频率(0.1Hz)对应双重频的低者unsigned short usSpulseW;/本圈脉宽(微秒)short usAngle;/仰角(0.01度)unsigned char cSweepStatus;/1=单要素2=三要素(单重频)3=三要素(双重频)unsigned char cAmbiguousp;/0=无软件退模糊1=软件退模糊VPPISCANPARAMETER;/注：在WEATHERRADAR结构的remain2处填写了的VPPISCANPARAMETER;结构表明本仰角层的参数(对PPI和VPPI而言)或其他含义(对RHI和单库FFT文件而言)VPPI有几层填几个，PPI RHI FFT数据只填写一层即一个原始数据：PPIRHI VPPI的数据排列如下：第0个径向的数据： Z(字)500个 V(字)500个 W(字)500个.第512个径向的数据：Z(字)500个 V(字)500个 W(字)500个注:Z V W均为16位的字型数据(在C语言中用 short int 有符整型)所读出的值除以10得到一个小数就是最终的值Z以dBz为单位 V W均以m/s为单位读出的数据如果是0x8000，说明此处无回波这样做的目的是为了必须区别回波强度为和无回波单库FFT数据排列如下:当为256点时：256个通道的功率值(字型 单位:dB)当为128点时：128个通道的功率值(字型 单位:dB)附函数段：/检查要打开的文件/filename:文件名CString CRadar:CheckData(CString filename)CFile *fp;CStringstr=INVALID,error,tempfilename=filename;tryfp=new CFile(filename,CFile:modeRead);catch(CFileException* e)ProcessFileError(e,1);return str;WEATHERRADAR tempinfo;fp-SeekToBegin();fp-Read(&tempinfo,sizeof(WEATHERRADAR);if(strcmp(tempinfo.cFileType,CINRADC)error.Format(%s 不是3830雷达数据文件!,tempfilename);AfxMessageBox(error,MB_OK|MB_ICONSTOP);delete fp;return str;if(tempinfo.usProdCode=0x8001)str=PPI;else if(tempinfo.usProdCode=0x8002) str=RHI;else if(tempinfo.usProdCode=0x8003) str=VPPI;else if(tempinfo.usProdCode=0x8004) str=ZPPI;else if(tempinfo.usProdCode=0x8005) str=CAPPI;else if(tempinfo.usProdCode=0x8006) str=FFT;elseerror.Format(文件 %s 不是3830雷达数据文件!,tempfilename);AfxMessageBox(error,MB_OK|MB_ICONSTOP);delete fp;return str;/读取数据/filename:要打开的文件名BOOL CPPI:LoadData(CString filename)CFile *fp;CString str;tryfp=new CFile(filename,CFile:modeRead);catch(CFileException* e)radar-ProcessFileError(e,1);return FALSE;fp-SeekToBegin();fp-Read(&info,sizeof(WEATHERRADAR);fp-Seek(218L,CFile:begin);fp-Read(&vppipara,sizeof(VPPISCANPARAMETER);rmax=int(vppipara.usMaxL/100.0+0.5);scanyj=vppipara.usAngle/100.0;/装地图文件long l;/申请数据存储区hData=:GlobalAlloc(GHND,512L*1500L*2L);if(!hData)AfxMessageBox(内存不足,MB_OK|MB_ICONSTOP);:GlobalFree(hData);delete fp;return FALSE;pData=(short*):GlobalLock(hData);for(l=0;lSeek(info.liDataOffset,CFile:begin);int i,BEGINPOS,angle;long pos;short tempData500;/我的程序要把每一个径向的连续的ZVW数据块读出后再按每个点排列成/第1个点的Z V W第2个点的Z V W ./第500个点的Z V W/可以不用考虑if(radar-scanmode=PPI)for(angle=0;angle512;angle+)for(BEGINPOS=0;BEGINPOSRead(tempData,500L*2);pos=long(angle*1500+BEGINPOS);for(i=0;iscanmode=ZPPI)for(angle=0;angleRead(tempData,500L*2);pos=long(angle*1500+0);for(i=0;iProcessFileError(e,1);return FALSE;fp-SeekToBegin();fp-Read(&info,sizeof(WEATHERRADAR);fp-Seek(218L,CFile:begin);fp-Read(&vppipara,sizeof(VPPISCANPARAMETER);rmax=int(vppipara.usMaxL/100);scanfw=info.usRHIA/100.0;beginyj=info.sRHIL/100.0;endyj=info.sRHIH/100.0;hData=:GlobalAlloc(GHND,512L*1500L*2L);if(!hData)AfxMessageBox(内存不足,MB_OK|MB_ICONSTOP);:GlobalFree(hData);delete fp;return FALSE;pData=(short*):GlobalLock(hData);long l;for(l=0;lSeek(info.liDataOffset,CFile:begin);int i,BEGINPOS,angle;long pos;short tempData500;for(angle=0;angle512;angle+)for(BEGINPOS=0;BEGINPOSRead(tempData,500L*2);pos=long(angle*1500+BEGINPOS);for(i=0;iIsZip) filename=radar-unzipfile;CFile *fp;CString str;tryfp=new CFile(filename,CFile:modeRead);catch(CFileException* e)radar-ProcessFileError(e,1);return FALSE;fp-SeekToBegin();fp-Read(&info,sizeof(WEATHERRADAR);fp-Seek(218L,CFile:begin);ScanPara-num=info.ucScanMode-100;int yjc,temprf1,temprf2;for(yjc=0;yjcnum;yjc+)fp-Read(&vppiparayjc,sizeof(VPPISCANPARAMETER);ScanPara-yjyjc=float(vppiparayjc.usAngle/100.0);temprf1=int(vppiparayjc.usPrf1/10.0+0.5);temprf2=int(vppiparayjc.usPrf2/10.0+0.5);ScanPara-dsp.frequency=temprf1;if(temprf1*1.0/temprf2)=(3.0/2.0) ScanPara-dsp.freqctrl=2;else if(temprf1*1.0/temprf2)=(4.0/3.0) ScanPara-dsp.freqctrl=3;else ScanPara-dsp.freqctrl=1;ScanPara-dsp.pulsewidth=vppiparayjc.usSpulseW;rmax=int(vppipara0.usMaxL/100);long l;/申请数据存储区hData=:GlobalAlloc(GHND,512L*1500L*2L*ScanPara-num);if(!hData)AfxMessageBox(内存不足,MB_OK|MB_ICONSTOP);:GlobalFree(hData);delete fp;return FALSE;pData=(short*):GlobalLock(hData);for(l=0;lnum*512*1500);l+) pDatal=short(0x8000);/因为pData中存放的是处理后的数据(有0),为了区分无回波和0，规