Matlab读取天气雷达基数据实现反射率因子三维可视化

小其其格

rabialee 发表于 2020-11-28 16:31
j1cor就是z方向的坐标

按照你的思路画了立体包络图，但是效果出来怎么不对呀？数据应该没问题，因为我画slice图效果是没问题的
楼主能帮我看看哪里不对么？


立体包络图效果图：


垂直剖面图效果：



主要代码如下：

1.     Point=rdBz_xyz(:,1:3);%--x，y，z坐标
   
2.     V=rdBz_xyz(:,4);%--反射率
   
3.     F = scatteredInterpolant(Point,V);
   
4.     %   F = scatteredInterpolant(P,v) 以数组形式指定样本点坐标。P 的行包含 v 中值的 (x, y) 或 (x, y, z) 坐标。
   
5.     xx=-100:1:100;
   
6.     yy=-100:1:100;
   
7.     zz=0:0.5:15;
   
8.     [xq,yq,zq] = meshgrid(xx,yy,zz);
   
9.     vq = F(xq,yq,zq);
   
10.     vq=smooth3(vq);
    
11. 
    
12.     fw=5;
    
13.     % --立体包络数值fw
    
14.     set(h1,'name','dBz定值包络立体图','MenuBar','none','ToolBar','none');
    
15.     title(['最外层表面的值为： ' , num2str(fw)]);
    
16.     fv=isosurface(xq,yq,zq,vq,fw);
    
17.     hiso=patch(fv);
    
18.     set(hiso,'facecolor','b','edgecolor','none');
    
19.     patch(isocaps(xq,yq,zq,vq,fw), 'FaceColor', 'interp', 'EdgeColor', 'none');
    
20. 
    

复制代码

schliezer

你好，想问下那11个lev的txt文件时怎么生成的，是通过循环fid=fopen('C:\Users\Rabia Lee\Desktop\Z_RADR_I_Z9737_20200820093300_lev1.txt','wt'); 这一步实现的吗？但我没有找到相应的循环体，而且第58行的n=2似乎也没有起到作用，感觉可能是后面的程序哪里漏了它。请指教，谢谢。

rabialee

schliezer 发表于 2020-12-1 16:22
你好，想问下那11个lev的txt文件时怎么生成的，是通过循环fid=fopen('C:%users\Rabia Lee\Desktop\Z_RADR_I ...

这个11层是由另外一个软件直接输出的，我上面写的代码的的一部分里面可以将9个仰角的反射率因子分别输出来，没有写循环，只是把三维可视化的思路写出来了。

rabialee

小其其格 发表于 2020-11-28 20:24
按照你的思路画了立体包络图，但是效果出来怎么不对呀？数据应该没问题，因为我画slice图效果是没问题的
...

看了你的代码，可能是画完包络面图没有加hold on;必须加这个才能把包络面和切片图一起显示，
hiso=patch(isosurface(xq,yq,zq,vq,25)); hold on;           %绘制雷达回波三维包络面
isonormals(xq,yq,zq,vq,hiso);
h=slice(xq,yq,zq,vq,-50,50,2);

rabialee

小其其格 发表于 2020-11-26 09:15
r_1(:,1)=Xcor;
    r_1(:,2)=Ycor;
    r_1(:,3)=j1cor;

不好意思我上次回复错了，j1cor是反射率，z坐标确实没有，可以自己写，前面的代码已经把极坐标换成了直角坐标

小其其格

rabialee 发表于 2020-12-4 15:40
看了你的代码，可能是画完包络面图没有加hold on;必须加这个才能把包络面和切片图一起显示，
hiso=patch ...

好的，前几天出差了，我试试

atmosphereic

请问怎么把笛卡尔坐标下的X、Y转换为经纬度呀

babylon1986

酷炫！{:eb502:}

lzwbeifeng

厉害，楼主交流一下

richie727

点赞大神点赞大神点赞大神