fseek(fidRadar, 40, 'bof');
radarNum = fread(fidRadar, 1, 'ushort');%雷达站数
startyear=fread(fidRadar,1,'ushort');
startmonth=fread(fidRadar,1,'ushort');
startday=fread(fidRadar,1,'ushort');
starthour=fread(fidRadar,1,'ushort');
startminute=fread(fidRadar,1,'ushort');
slongitude = fread(fidRadar, 1, 'int')/1000;%雷达起始经度
slatitude = fread(fidRadar, 1, 'int')/1000;%雷达起始纬度
elongi= fread(fidRadar, 1, 'int')/1000;%雷达结束经度
elatitude = fread(fidRadar, 1, 'int')/1000;%雷达结束纬度
sheight=fread(fidRadar,1,'ushort');
eheight=fread(fidRadar,1,'ushort');
xgrids=fread(fidRadar, 1, 'ushort');%经度方向格点数
ygrids=fread(fidRadar, 1, 'ushort');%纬度方向格点数
fseek(fidRadar, 76, 'bof');
levels = fread(fidRadar, 1, 'ushort');%层数
dx=fread(fidRadar,1,'uint32')/10000;
dy=fread(fidRadar,1,'uint32')/10000;
dz=fread(fidRadar,1,'ushort');
fseek(fidRadar,108+4,'bof');%跳到数据段处
for i=1:radarNum
    staindex=fread(fidRadar,1,'int');
    staname=fread(fidRadar,20,'char');
    statype=fread(fidRadar,12,'char');
    radarlon=fread(fidRadar,1,'int')/1000;
    radarlat=fread(fidRadar,1,'int')/1000;
    radarheight=fread(fidRadar,1,'int');
    radarcover=fread(fidRadar,1,'ushort');
end
for i=1:levels
    height(i)=fread(fidRadar,1,'ushort');
end
data=fread(fidRadar,levels*ygrids*xgrids,'uchar');

for k=1:levels
    for j=1:ygrids
        for i=1:xgrids
            echo(k,j,i)=data((k-1)*ygrids*xgrids+(j-1)*xgrids+i);           
        end
    end
end
dx=0.01;
dy=0.01;
dz=0.5;