关于雷达资料画图的问题

静雅武

Mild_Alice 发表于 2014-3-3 17:31

                               
登录/注册后可看大图

雷达型号：CINRADC/CC雷达，数据说明就是数据区有三种数据（Z，V，W）。总共有14层，这些都已经解析出来， ...

是哪种格式啊？rec还是txt？rec是二进制的，你可以用MATLAB试试能不能打开。fopen。txt就直接用记事本打开，你就能看见数据是什么样的了。

Mild_Alice

莫名╰☆╮殇 发表于 2014-3-3 17:54

                               
登录/注册后可看大图

要画图要么借助已有的雷达画图软件，如果自己fortran编程的话就得看看雷达坐标系和直角坐标系的转换方法。将 ...

请问你做过这方面的没？已有的雷达画图软件有些什么？我现在虽然知道了14个层中每一层的Z  V  W数据（存在三维数组中），后面图形化显示就不知如何着手？请不吝赐教

Mild_Alice

静雅武 发表于 2014-3-3 21:24

                               
登录/注册后可看大图

是哪种格式啊？rec还是txt？rec是二进制的，你可以用MATLAB试试能不能打开。fopen。txt就直接用记事本打开 ...

原始文件是二进制文件（.bin文件），现在我已经解析出来了，就是已经得到14层中每一层的回波强度（Ｚ），速度（Ｖ），谱宽（Ｗ）数据并且存在数组中了。
　　现在我是要将每层中的三种数据分别绘制出图形，也就是绘制１４×３张图。
　　问题就是不知如何绘制雷达图？？？？

莫名╰☆╮殇

Mild_Alice 发表于 2014-3-3 22:19

                               
登录/注册后可看大图

请问你做过这方面的没？已有的雷达画图软件有些什么？我现在虽然知道了14个层中每一层的Z  V  W数据（存在 ...

bin文件就是雷达的原始文件吧，呃有一个画图软件叫PUP，但是得要购买的，大部分处理雷达数据的程序貌似都是C++转换成grd之后在用fortran处理的

莫名╰☆╮殇

Mild_Alice 发表于 2014-3-3 22:19

                               
登录/注册后可看大图

请问你做过这方面的没？已有的雷达画图软件有些什么？我现在虽然知道了14个层中每一层的Z  V  W数据（存在 ...

你说的知道了每一层的数据，是已经可以用fortran读取到了么，如果读取到了，就直接把雷达坐标系转换到直角坐标系上，转换的方法比较多，你可以参考VRADS里面的程序

Mild_Alice

莫名╰☆╮殇 发表于 2014-3-4 10:03

                               
登录/注册后可看大图

你说的知道了每一层的数据，是已经可以用fortran读取到了么，如果读取到了，就直接把雷达坐标系转换到直角 ...

我用c#读到的数据。能具体说说如何坐标转换不？然后比如说是第一层仰角为0.5度的512*500的回波强度数据，它投影之后是不是表示不同经度、纬度的回波强度值？然后将各个回波强度值按经度纬度值标注到地图上？？

静雅武

Mild_Alice 发表于 2014-3-3 22:26

                               
登录/注册后可看大图

原始文件是二进制文件（.bin文件），现在我已经解析出来了，就是已经得到14层中每一层的回波强度（Ｚ）， ...

一般是用MATlab画，也有很多人用c++，我把我的一个m文件的一段给你参考，你也只能参考，我是画rec的，所以你自己研究下，
% filepath
clear;close all
Rec_filepath='H:\data';
fns = dir([Rec_filepath '1*.rec']);
savefilepath='D:\My_HMBQ\Figures\';

% read
mm=length(fns);
height=(0:499)*30/1000;

for ff=1:mm;
    fid=fopen([Rec_filepath '\' fns(ff).name],'r');
    [reflectivity,velocity_v,spectwidth_v,depolarized_ratio,time_rec]=readREC(fid);
   
   
    % figure-reflectivity
    figure(1)
    imagesc(1:size(reflectivity,2),height',reflectivity);colorbar;
    set(gca,'YDir','normal');
    ylabel('Height (km)','FontName','Times New Roman','FontSize',16);
    ylim([0,15]);
    xlabel('Time','FontName','Times New Roman','FontSize',16);
    text(10,14,['Reflectivity ' num2str(time_rec(1)) '-' num2str(time_rec(2)) '-' num2str(time_rec(3)) ' ' num2str(time_rec(4)) ':' num2str(time_rec(5)) ':' num2str(time_rec(6))],'FontName','Times New Roman','FontSize',16);
    Colormap=[1 1 1;0.8571 0.8571 1;0.7143 0.7143 1;0.5714 0.5714 1;0.4286 0.4286 1;0.2857 0.2857 1;0.1429 0.1429 1;0 0 1;0 0.0625 1;0 0.125 1;0 0.1875 1;0 0.25 1;0 0.3125 1;0 0.375 1;0 0.4375 1;0 0.5 1;0 0.5625 1;0 0.625 1;0 0.6875 1;0 0.75 1;0 0.8125 1;0 0.875 1;0 0.9375 1;0 1 1;0.0625 1 0.9375;0.125 1 0.875;0.1875 1 0.8125;0.25 1 0.75;0.3125 1 0.6875;0.375 1 0.625;0.4375 1 0.5625;0.5 1 0.5;0.5625 1 0.4375;0.625 1 0.375;0.6875 1 0.3125;0.75 1 0.25;0.8125 1 0.1875;0.875 1 0.125;0.9375 1 0.0625;1 1 0;1 0.9375 0;1 0.875 0;1 0.8125 0;1 0.75 0;1 0.6875 0;1 0.625 0;1 0.5625 0;1 0.5 0;1 0.4375 0;1 0.375 0;1 0.3125 0;1 0.25 0;1 0.1875 0;1 0.125 0;1 0.0625 0;1 0 0;0.9375 0 0;0.875 0 0;0.8125 0 0;0.75 0 0;0.6875 0 0;0.625 0 0;0.5625 0 0;0.5 0 0];
    set(gca,'clim',[-40,20]);
    saveas(gcf,[savefilepath 'emf\' num2str(time_rec(1)) '-' num2str(time_rec(2)) '-' num2str(time_rec(3)) ' ' num2str(time_rec(4)) '-' num2str(time_rec(5)) '-' num2str(time_rec(6)) ' Reflectivity.emf'], 'emf');
    close;

莫名╰☆╮殇

Mild_Alice 发表于 2014-3-4 10:55

                               
登录/注册后可看大图

我用c#读到的数据。能具体说说如何坐标转换不？然后比如说是第一层仰角为0.5度的512*500的回波强度数据， ...

嗯哈，是的你说的是对的，画出来的图就是格点图。

wowiwf

Mild_Alice 发表于 2014-3-3 16:27
我也在研究这个，跟楼主的问题完全一样。我们要解析这个雷达文件Z_RADR_I_Z9931_20131003035900_O_DOR_CC_C ...

我还不会解析

liudage1983

一个像素一个像素的画图