Extracting max intensity coordinates from image sequence

Hello, I'm a complete beginner and I'm currently trying to extract maximum intensity pixel coordinates from multiple images at one time (I have about 2500 images) and store them in a matrix with slice number that correspond. After doing some research here I came up with this code :

filedir = '...';
imds = imageDatastore(filedir);
imgs = readall(imds);

for k=1:numel(imgs)
    [x,y]=size(imgs(k));
    max_int = max(imgs(:));
    for i=1:x
        for j=1:y
            if (imgs(i,j)==max_int)
            k;
            i;
            j;
            end
        end
    end
end

I understand how to do it for one image but doing it in a loop for multiple images is quite difficult for me.

Can someone give me some hint or solution on how to modify this loop to have what I want ?

Thank you very much.

Answer 1

[yanqi liu]

yes，sir，may be use some index to transfer，such as

im = imread('rice.png');
[max_p, ~] = max(im(:));
disp(max_p)

   204

ind = find(im(:)==max_p);
% use ind2sub
[r,c] = ind2sub([size(im,1) size(im,2)], ind);
% display
for i = 1 : length(r)
    fprintf('(%d, %d) is max value %d\n', r(i), c(i), im(r(i), c(i)));
end

(51, 93) is max value 204
(51, 94) is max value 204
(52, 94) is max value 204
(52, 95) is max value 204
(70, 95) is max value 204
(53, 96) is max value 204
(69, 96) is max value 204
(70, 96) is max value 204
(71, 96) is max value 204
(69, 97) is max value 204
(70, 97) is max value 204
(68, 98) is max value 204
(69, 98) is max value 204
(68, 99) is max value 204
(68, 100) is max value 204
(83, 116) is max value 204