MATLAB: Removing part of a dimension

image processingindexinglogical?maskMATLABmrmri

Hi all,

I have a 4D MRI dataset of 80 x 80 x 16 x 23 from which I calculate movement outliers based upon standard deviation. I receive as output an logical array of 1×23 showing outliers with a value of 0. Example: OutlierArray = [1 1 1 1 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]

I want to remove the outlier located at, in this case, the 5th volume using the 'OutlierArray' variable and concatenate the remaining volumes.

The resulting 4D dataset would then look like 80 x 80 x 16 x 22.

Can someone please help me? Thanks in advance!

Best Answer

Hi Mathijs,

If your data is called MRI, then the following should do the job:

MRI(:,:,:,logical(OutlierArray))=[];

Related Solutions

MATLAB: How to loop a DICOM and corresponding ROI folder and write results to a .csv

diFolderName = '.\DICOM';
roiFolderName = '.\ROI';
numfiles = 60;
filenames = cell(numfiles, 1);
meanHU = cell(numfiles, 1);
for K = 1 : numfiles
  basename = sprintf('F%06d.dcm', K);
  folderFileName = fullfile(diFolderName, basename );
    dinfo = dicominfo(folderFileName);
    % Read the dicom image data and rescale to the correct values
    myImage = double(dicomread(folderFileName))*dinfo.RescaleSlope + dinfo.RescaleIntercept;
    ROIname = fullfile(roiFolderName, sprintf('%d.nii', K));
    ROI = niftiread(ROIname);      %requires R2017b or later!
    values = myImage(ROI);
    filenames{K} = basename;
    meanHU{K} = mean(values);
  end
outdata = [filenames, meanHU];
xlswrite('meanHU_Summary.csv', outdata);

Note: niftiread() is new in R2017b. For earlier versions it is necessary to use something from the File Exchange.

MATLAB: Index variable based on time points

output = zeros(size(mri,1), size(mri,2), 35, 40, class(mri));
for t = 1 : 40
    output(:,:,:,t) = mri(:,:,35*(t-1) + (1:35));
end

However if this was your purpose then just use

output = reshape(mri, size(mri,1), size(mri,2), 35, 40);

Best Answer

Related Solutions

MATLAB: How to loop a DICOM and corresponding ROI folder and write results to a .csv

MATLAB: Index variable based on time points

Related Question