add orthogonalProjection.m

Author: dlegland

src/@Image/Image.m

@@ -293,7 +293,7 @@ function copyFields(obj, that)
         nc = obj.DataSize(4);
         switch lower(typeName)
             case 'binary', if nc ~= 1, tf = false; end
-            case 'grayscale', if nc ~= 1, tf = false; end
+            case 'grayscale', if nc ~= 1 || isfloat(obj.Data), tf = false; end
             case 'intensity', if nc ~= 1, tf = false; end
             case 'color', if nc ~= 3, tf = false; end
             case 'label', if nc ~= 1, tf = false; end

src/@Image/max.m

@@ -1,7 +1,7 @@
 function res = max(obj, varargin)
 % Compute the maximal value within image.
 %
-%   M = mean(IMG);
+%   M = max(IMG);
 %
 %   Example
 %     img = Image.read('rice.png');
@@ -10,7 +10,7 @@
 %        204
 %
 %   See also
-%     mean, median, min
+%     mean, median, min, orthogonalProjection
 %
 
 if isempty(varargin)

src/@Image/mean.m

@@ -10,7 +10,7 @@
 %       111.2468
 %
 %   See also
-%     median
+%     median, min, max, orthogonalProjection
 %
 
 res = mean(double(obj.Data(:)));

src/@Image/min.m

@@ -1,16 +1,16 @@
 function res = min(obj, varargin)
 % Compute the minimal value within image.
 %
-%   M = mean(IMG);
+%   M = min(IMG);
 %
 %   Example
 %     img = Image.read('rice.png');
-%     max(img)
+%     min(img)
 %     ans =
 %        40
 %
 %   See also
-%     mean, median, max
+%     mean, median, max, orthogonalProjection
 %
 
 

src/@Image/orthogonalProjection.m

@@ -0,0 +1,38 @@
+function res = orthogonalProjection(obj, dir, opName)
+% Project image intensities along one of the main directions..
+%
+%   PROJ = orthogonalProjection(IMG, DIR, OP)
+%   Project the intensities of the image IMG in the direction given by DIR.
+%   The operation OP applied on the values can be one of 'Max' (default),
+%   'Min', or 'Mean'.
+%
+%   Example
+%     img = Image.read('mri.tif');
+%     projXY = orthogonalProjection(img, 3, 'Mean');
+%     figure; show(projXY);
+%
+%   See also
+%     slice, min, max, mean
+%
+ 
+% ------
+% Author: David Legland
+% e-mail: david.legland@inrae.fr
+% INRAE - BIA Research Unit - BIBS Platform (Nantes)
+% Created: 2021-01-27,    using Matlab 9.8.0.1323502 (R2020a)
+% Copyright 2021 INRAE.
+
+if strcmpi(opName, 'max')
+    data = max(obj.Data, [], dir);
+elseif strcmpi(opName, 'min')
+    data = min(obj.Data, [], dir);
+elseif strcmpi(opName, 'mean')
+    data = mean(obj.Data, dir);
+else
+    error('Unknown operation name: %s', opName)
+end
+
+% create projected image
+newName = createNewName(obj, '%s-proj');
+res = Image('Data', data, 'Parent', obj, 'Name', newName);
+        

tests/image/test_orthogonalProjection.m

@@ -0,0 +1,67 @@
+function tests = test_orthogonalProjection
+% Test suite for the file orthogonalProjection.
+%
+%   Test suite for the file orthogonalProjection
+%
+%   Example
+%   test_orthogonalProjection
+%
+%   See also
+%     orthogonalProjection
+
+% ------
+% Author: David Legland
+% e-mail: david.legland@inrae.fr
+% Created: 2021-01-27,    using Matlab 9.8.0.1323502 (R2020a)
+% Copyright 2021 INRAE - BIA-BIBS.
+
+tests = functiontests(localfunctions);
+
+
+function test_Proj_Z(testCase) %#ok<*DEFNU>
+% Test call of function without argument.
+
+img = createTestImage3d();
+
+projXY = orthogonalProjection(img, 3, 'max'); 
+
+assertEqual(testCase, size(projXY), [30 20]);
+
+
+
+function test_Proj_Z_mean(testCase) %#ok<*DEFNU>
+% Test call of function without argument.
+
+img = createTestImage3d();
+
+projXY = orthogonalProjection(img, 3, 'mean'); 
+
+assertEqual(testCase, size(projXY), [30 20]);
+
+
+function test_Proj_Y(testCase) %#ok<*DEFNU>
+% Test call of function without argument.
+
+img = createTestImage3d();
+
+projXZ = orthogonalProjection(img, 2, 'max'); 
+
+assertEqual(testCase, size(projXZ), [30 1 10]);
+
+
+function test_Proj_X(testCase) %#ok<*DEFNU>
+% Test call of function without argument.
+
+img = createTestImage3d();
+
+projYZ = orthogonalProjection(img, 1, 'max'); 
+
+assertEqual(testCase, size(projYZ), [1 20 10]);
+
+
+function img = createTestImage3d
+% Create an image with size 30x20x10 containing a 20x10x5 cuboid.
+
+data = zeros([30 20 10], 'uint8');
+data(6:25, 6:15, 3:7) = 100;
+img = Image('Data', data);