Changeset 1139 for trunk/src/filter_tps.m

Timestamp:

Apr 25, 2024, 5:56:32 PM (6 months ago)

Author:

sommeria

Message:

series/test patch added + some bug corrections

File:

• trunk/src/filter_tps.m (modified) (6 diffs)

trunk/src/filter_tps.m

@@ -1 +1 @@
 %'filter_tps': find the thin plate spline coefficients for interpolation-smoothing
 %------------------------------------------------------------------------
-% [SubRange,NbCentre,Coord_tps,U_tps,V_tps,W_tps,U_smooth,V_smooth,W_smooth,FF] =filter_tps(Coord,U,V,W,SubDomainSize,Rho,Threshold)
+% [SubRange,NbCentre,Coord_tps,U_tps,V_tps,W_tps,U_smooth,V_smooth,W_smooth,FF] =filter_tps(Coord,U,V,W,SubDomainSize,FieldSmooth,Threshold)
 %
 % OUTPUT:
@@ -16 +16 @@
 % U,V, possibly W: set of velocity components of the initial data
 % SubdomainSize: estimated number of data points in each subdomain
-% Rho: smoothing parameter
+% FieldSmooth: smoothing parameter
 % Threshold: max diff accepted between smoothed and initial data 
 
@@ -38 +38 @@
 %=======================================================================
 
-function [SubRange,NbCentre,Coord_tps,U_tps,V_tps,W_tps,U_smooth,V_smooth,W_smooth,FF] =filter_tps(Coord,U,V,W,SubDomainSize,Rho,Threshold)
+function [SubRange,NbCentre,Coord_tps,U_tps,V_tps,W_tps,U_smooth,V_smooth,W_smooth,FF] =filter_tps(Coord,U,V,W,SubDomainSize,FieldSmooth,Threshold)
 
 %% adjust subdomain decomposition
@@ -61 +61 @@
 
 %% smoothing parameter
-rho=Siz(1)*Siz(2)*Rho/1000;%optimum rho increase as the area of the subdomain (division by 1000 to reach good values with the default GUI input)
+smoothing=Siz(1)*Siz(2)*FieldSmooth/1000;%optimum smoothing increase as the area of the subdomain (division by 1000 to reach good values with the default GUI input)
 
 %% default output
@@ -100 +100 @@
         % subdomain includes enough vectors, perform tps interpolation
         else
-            [U_smooth_sub,U_tps_sub]=tps_coeff(Coord(ind_sel,:),U(ind_sel),rho);
-            [V_smooth_sub,V_tps_sub]=tps_coeff(Coord(ind_sel,:),V(ind_sel),rho);
+            [U_smooth_sub,U_tps_sub]=tps_coeff(Coord(ind_sel,:),U(ind_sel),smoothing);
+            [V_smooth_sub,V_tps_sub]=tps_coeff(Coord(ind_sel,:),V(ind_sel),smoothing);
             UDiff=U_smooth_sub-U(ind_sel);% difference between interpolated U component and initial value
             VDiff=V_smooth_sub-V(ind_sel);% difference between interpolated V component and initial value
@@ -132 +132 @@
             % else interpolation-smoothing is done again with the selected vectors
             else
-                [U_smooth_sub,U_tps_sub]=tps_coeff(Coord(ind_sel(ind_ind_sel),:),U(ind_sel(ind_ind_sel)),rho);
-                [V_smooth_sub,V_tps_sub]=tps_coeff(Coord(ind_sel(ind_ind_sel),:),V(ind_sel(ind_ind_sel)),rho);
+                [U_smooth_sub,U_tps_sub]=tps_coeff(Coord(ind_sel(ind_ind_sel),:),U(ind_sel(ind_ind_sel)),smoothing);
+                [V_smooth_sub,V_tps_sub]=tps_coeff(Coord(ind_sel(ind_ind_sel),:),V(ind_sel(ind_ind_sel)),smoothing);
                 x_width=(SubRange(1,2,isub)-SubRange(1,1,isub))/pi;
                 y_width=(SubRange(2,2,isub)-SubRange(2,1,isub))/pi;