simulate_points_dem_15

% generate a set of points for dem-interpolation

0001 %% generate a set of points for dem-interpolation
0002 
0003 function [points,BB,dx,sigma_k,sigma_h,dem]=simulate_points_dem_15
0004 
0005 
0006 %% Example PCV (do not change)
0007 % bounding box
0008 BB=[0,0,9,6.000]*2/3;
0009 % grid size
0010 factor=1;
0011 dx = 1/factor;
0012 
0013 sigma_h = 1.0;
0014 sigma_k = 4/factor^(3/2);
0015 
0016 points = [...
0017     1.2, 1.5, 2,sigma_h;...
0018     1.8, 1.2, 3,sigma_h;...
0019     5.4, 2.4, 4,sigma_h;...
0020     2.1, 4.8, 5,sigma_h;...
0021     8.7, 3.3, 6,sigma_h;...
0022     6.0, 5.4, 1,sigma_h ...
0023     ];
0024 
0025 points(:,1:2)=points(:,1:2)*2/3;
0026 
0027 %%
0028 
0029 xmin = BB(1);
0030 ymin = BB(2);
0031 xmax = BB(3);
0032 ymax = BB(4);
0033 Nr = ceil((xmax-xmin)/dx)+1;
0034 Mc = ceil((ymax-ymin)/dx)+1;
0035 %points(:,3)=(points(:,3)-mean(points(:,3)))*sqrt(Nr)/4;
0036 % interpolation kernel
0037 sigma = min(Nr,Mc)/5;
0038 
0039 dem = zeros(Nr,Mc);
0040 
0041 
0042 return