MATLAB的滑动空间百分位法

kongfeng0824

1. <font size="3">function StaPrctile = StaFiltPrctile( location, stavalues, prc, winsize )
   
2.     [~,StaCount] = size(stavalues);
   
3.     StaPrctile = zeros(StaCount,1);
   
4.     %% the distance of each station
   
5.     s = length(location);
   
6.     DisMat = zeros(s,s);
   
7.     for m = 1:s
   
8.         for n = m:s
   
9.             tempdis = sqrt(( location(1,m) - location(1,n) )^2 + ( location(2,m) - location(2,n) )^2 );
   
10.             DisMat(m,n) = tempdis;
    
11.             DisMat(n,m) = tempdis;
    
12.             clear temp*;
    
13.         end
    
14.     end
    
15.    
    
16.     %%
    
17.     for m = 1:s           
    
18.         tempvalues = stavalues(:, DisMat(m,:) < winsize);
    
19.         StaPrctile(m) = prctile(tempvalues(tempvalues > 0),prc);
    
20.         clear tempvalues;
    
21.     end
    
22. 
    
23. end</font>
    

复制代码

来自群组: 北京师范大学地学-气象

zhchlue

滑动空间百分位法是什么意思？