simulate_points_dem_square

% 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,oi,dem]=simulate_points_dem_square
0004 
0005 
0006 %% square box, all observed
0007 N=5;
0008 % bounding box
0009 BB=[0,0,N-1,N-1];
0010 % grid size
0011 factor=1;
0012 dx = 1/factor;
0013 
0014 sigma_h = 1.0;
0015 sigma_k = 4/factor^(3/2);
0016 k=0;
0017 for i=1:N
0018     for j=1:N
0019         k=k+1;
0020         points(k,1:4)=[i,j,1,sigma_h];
0021     end
0022 end
0023 
0024 %%
0025 
0026 xmin = BB(1);
0027 ymin = BB(2);
0028 xmax = BB(3);
0029 ymax = BB(4);
0030 Nr = ceil((xmax-xmin)/dx)+1;
0031 Mc = ceil((ymax-ymin)/dx)+1;
0032 %points(:,3)=(points(:,3)-mean(points(:,3)))*sqrt(Nr)/4;
0033 % interpolation kernel
0034 sigma = min(Nr,Mc)/5;
0035 
0036 dem = zeros(Nr,Mc);
0037 oi=ones(N^2,1);
0038 
0039 
0040 return