Changeset 1122 for trunk/src/transform_field/ima_filter.m

Timestamp:

Apr 21, 2023, 11:01:06 AM (13 months ago)

Author:

sommeria

Message:

various updates

File:

• trunk/src/transform_field/ima_filter.m (modified) (2 diffs)

trunk/src/transform_field/ima_filter.m

@@ -1 @@
-% 'ima_filter': example of image transform with input parameters: low-pass filter of an image
+% 'ima_filter': low-pass filter of an image or other 2D fields defined on a regular grid
+% the size of the filtering window in x and y is interactivement defined
 
 %------------------------------------------------------------------------
 %%%%  Use the general syntax for transform fields with a single input and parameters %%%%
-% OUTPUT: 
-% DataOut:   output field structure 
+% OUTPUT:
+% DataOut:   output field structure
 %
 %INPUT:
 % DataIn:  input field structure
+%
 % Param: matlab structure whose field Param.TransformInput contains the filter parameters
+% DataIn_1: variables possibly introduced as a second input field
 %-----------------------------------
 
@@ -78 +81 @@
 if exist('DataIn_1','var')
     [CellInfo,NbDim,errormsg]=find_field_cells(DataIn_1);
-for icell=1:numel(CellInfo)
-    if isfield(CellInfo{icell},'CoordType')&& strcmp(CellInfo{icell}.CoordType,'grid')
-        for ivar=1:numel(CellInfo{icell}.VarIndex)
-            VarName=DataIn_1.ListVarName{CellInfo{icell}.VarIndex(ivar)};
-            Atype=class(DataIn_1.(VarName));% detect integer 8 or 16 bits
-            if numel(size(DataIn_1.(VarName)))==3
-                DataOut.(VarName)=filter2(Mfiltre,sum(DataIn_1.(VarName),3));%filter the input image, after summation on the color component (for color images)
-                DataOut.(VarName)=uint16(DataOut.(VarName)); %transform to 16 bit images
-            else
-                DataOut.(VarName)=filter2(Mfiltre,DataIn_1.(VarName));
-                DataOut.(VarName)=feval(Atype,DataOut.(VarName));%transform to the initial image format
+    for icell=1:numel(CellInfo)
+        if isfield(CellInfo{icell},'CoordType')&& strcmp(CellInfo{icell}.CoordType,'grid')
+            for ivar=1:numel(CellInfo{icell}.VarIndex)
+                VarName=DataIn_1.ListVarName{CellInfo{icell}.VarIndex(ivar)};
+                Atype=class(DataIn_1.(VarName));% detect integer 8 or 16 bits
+                if numel(size(DataIn_1.(VarName)))==3
+                    DataOut.(VarName)=filter2(Mfiltre,sum(DataIn_1.(VarName),3));%filter the input image, after summation on the color component (for color images)
+                    DataOut.(VarName)=uint16(DataOut.(VarName)); %transform to 16 bit images
+                else
+                    DataOut.(VarName)=filter2(Mfiltre,DataIn_1.(VarName));
+                    DataOut.(VarName)=feval(Atype,DataOut.(VarName));%transform to the initial image format
+                end
             end
         end
     end
 end
-end