add static method Image.importRaw

Author: dlegland

src/@Image/Image.m

@@ -94,6 +94,7 @@
 
     img = read(fileName, varargin)
     img = readSeries(fileName, varargin)
+    img = importRaw(fileName, varargin)
 
     info = fileInfo(fileName)
 

src/@Image/importRaw.m

@@ -0,0 +1,91 @@
+function img = importRaw(fileName, dims, format, varargin)
+% Read image data from raw data file.
+%
+%   output = importRaw(FNAME, DIMS, FORMAT, OPTIONS)
+%
+%   Example
+%     img = Image.importRaw(fileName, [512 512 27], 'uint8');
+%
+%   See also
+%     read
+%
+ 
+% ------
+% Author: David Legland
+% e-mail: [email protected]
+% Created: 2020-12-23,    using Matlab 9.7.0.1190202 (R2020b)
+% Copyright 2020 INRA - Cepia Software Platform.
+
+% default parameter values
+pixelType = format;
+byteOrder = 'l'; % assume little endian is default
+offset = 0;
+verbose = false;
+
+
+% parse input options
+while length(varargin) > 1
+    name = varargin{1};
+    value = varargin{2};
+    
+    if strcmp(name, 'offset')
+        offset = value;
+    elseif strcmpi(name, 'byteOrder')
+        byteOrder = value;
+    elseif strcmp(name, 'verbose')
+        verbose = value;
+    else
+        error(['Unknown parameter name: ' name]);
+    end
+    
+    varargin(1:2) = [];
+end
+
+%% Allocate
+
+if verbose
+    disp('allocate memory');
+end
+% img = zeros(dims, pixelType);
+
+
+%% Open
+
+if verbose
+    disp('open data file');
+end
+
+% open for for binary reading
+f = fopen(fileName, 'rb');
+if f == -1
+    error(['Unable to open data file: ' fileName]);
+end
+
+% look for beginning of data
+fseek(f, offset, 'bof');
+
+
+%% Fill array
+
+if verbose
+    fprintf('read data...');
+end
+tic;
+data = fread(f, prod(dims), ['*' pixelType], byteOrder);
+t = toc;
+if verbose
+    fprintf(' elapsed time: %7.3f s\n', t);
+end
+
+%% Cleanup
+
+% close file
+fclose(f);
+
+% permute dims 1 and 2
+if verbose
+    disp('permute');
+end
+data = reshape(data, dims);
+
+img = Image('Data', data);