CINRAD-SA雷达系统主要产品与算法简介之五

winter6789

超强，版主到底是谁？

meifhaisxiac

谢谢楼主分享哈O(∩_∩)O~

碧水寒

很好的帖子，期待楼主下次再介绍雷达数据处理方面的介绍啊，比如将SA型号的数据处理成MDV格式之类的

ayzqs

给力呀，这个收了！谢谢

pang

nice thank u

智慧de苹果

谢谢楼主，学习了!

liangyuli

讲的不错 谢谢啦 学习了

liyuting72

好东西，感觉巩固了很多知识

gumblear

版主，您好，我按照您的公式编写了matlab程序，但是算出来的VIL和产品差异很大，麻烦您帮我看看错在哪儿，万分感谢{:eb511:}
以下是程序：
clear;clc;
nwg=301;%格点数
data=load(['d:\data\ref',num2str(t),'.txt']);
E=data(1,:);%%每层的仰角值，已经转换为弧度
c=length(E);%一共有c层仰角

vil=zeros(nwg,nwg);
reff=zeros(nwg,nwg,c);
for k=1:c
     ref=load(['d:\data\ref1_',num2str(k),'.txt']);%雷达反射率，单位dBZ
     reff(:,:,k)=ref；
end
    for m=1:nwg
        for n=1:nwg
               if max(reff(n,m,:))>0
                   s1=((n-151)^2+(m-151)^2)^0.5;%反射率到圆点的投影距离               
                   %-----------------垂直液态水含量
                   for i=1:c
                       if reff(n,m,i)>55
                           reff(n,m,i)=55;%将大于55dBz的值取为55
                       end
                   end
                   z=zeros(c,1);
                   M=zeros(c,1);
                   for j=1:c
                       z(j)=10^(reff(n,m,j)/10);%转化为反射率因子，单位是mm^6/m^3
                   end                  
                   dh=zeros(c-1,1);
                   for j=1:c-1
                       dh(j)=s1.*(tan(E(j+1))-tan(E(j)))*1000;%相邻仰角的高度差，单位m         
                       M(j)=(3.44*10^(-6)).*((z(j)+z(j+1))/2).^(4/7).*dh(j);%相邻仰角的液态水含量，单位kg/m^2
                   end
                   vil(n,m)=sum(M);                  
               end
        end
    end
图片不知道怎么放上来，反正两个差异很大

wangjin0919

所以垂直累计液态水为什么分辨率是4km*4km，却是以极坐标形式表示的呢