You should be able to do
minValue = min(thisImage(:))
[row, column] = find(thisImage == minValue);
for each time point (slice of your 3-D image volume). Here is a full demo:
imgfilt = randi(65535, 64, 64, 1913);
imgfilt(44:66) = nan;
[rows, columns, numberOfTimePoints] = size(imgfilt)
for slice = 1 : numberOfTimePoints
thisImage = imgfilt(32:end, 32:end, slice);
minValue = min(thisImage(:));
[row, column] = find(thisImage == minValue);
row = row + 31;
column = column + 31;
ca{slice, 1} = row;
ca{slice, 2} = column;
for k = 1 : length(row)
r = row(k);
c = column(k);
fprintf('For time point %d, the min value of %f occurs at row %d, column %d.\n', ...
slice, imgfilt(r, c, slice), r, c);
end
end
Best Answer