Changeset 356 for trunk/src/tps_eval.m

Timestamp:

Jan 3, 2012, 12:58:52 AM (12 years ago)

Author:

sommeria

Message:

civ updated with new functions for opening files, consistently with uvmat
Bugs to be expected (use previous version then)

File:

• trunk/src/tps_eval.m (modified) (1 diff)

trunk/src/tps_eval.m

@@ -1 @@
-%   EM:     MxN matrix whose i,j position contains the Euclidean
-%              distance between the i-th data site and j-th center
-  function EM = tps_eval(dsites,ctrs)
-  [M,s] = size(dsites); [N,s] = size(ctrs);
-  EM = zeros(M,N);
-  
-  % calculate distance matrix: accumulate sum of squares of coordinate differences
-  % The ndgrid command produces two MxN matrices:
-  %   Dsite, consisting of N identical columns (each containing
-  %       the d-th coordinate of the M data sites)
-  %   Ctrs, consisting of M identical rows (each containing
-  %       the d-th coordinate of the N centers)
-  for d=1:s
-     [Dsites,Ctrs] = ndgrid(dsites(:,d),ctrs(:,d));
-     EM = EM + (Dsites-Ctrs).^2;%EM=square of distance matrices
-  end
+%'tps_eval': calculate the thin plate spline (tps) interpolation at a set of points
+% see tps_ceff for more information
+%------------------------------------------------------------------------
+% function EM = tps_eval(dsites,ctrs)
+%------------------------------------------------------------------------
+% OUPUT:
+% EM:  Mx(N+s) matrix representing the contributions at the M sites 
+%   from unit sources located at each of the N centers, + (s+1) columns
+%   representing the contribution of the linear gradient part.
+%
+%INPUT:
+%dsites:  Nxs matrix representing the postions of the N 'observation' sites, with s the space dimension
+%ctrs: Mxs matrix  representing the postions of the M centers, sources of the tps,
+%
+% related functions:
+% tps_coeff, tps_eval_dxy
 
-  % calculate tps
-  np=find(EM~=0);
-  EM(np) = EM(np).*log(EM(np))/2;%= tps formula r^2 log(r) (EM=r^2)
-  
-  % add linear gradient part:
-  EM = [EM ones(M,1) dsites];
+function EM = tps_eval(dsites,ctrs)
+[M,s] = size(dsites); [N,s] = size(ctrs);
+EM = zeros(M,N);
+
+% calculate distance matrix: accumulate sum of squares of coordinate differences
+% The ndgrid command produces two MxN matrices:
+%   Dsite, consisting of N identical columns (each containing
+%       the d-th coordinate of the M data sites)
+%   Ctrs, consisting of M identical rows (each containing
+%       the d-th coordinate of the N centers)
+for d=1:s
+ [Dsites,Ctrs] = ndgrid(dsites(:,d),ctrs(:,d));
+ EM = EM + (Dsites-Ctrs).^2;%EM=square of distance matrices
+end
+
+% calculate tps
+np=find(EM~=0);
+EM(np) = EM(np).*log(EM(np))/2;%= tps formula r^2 log(r) (EM=r^2)
+
+% add linear gradient part:
+EM = [EM ones(M,1) dsites];