Finding Min Max Pixels and Co-ordinates from DICOM images?

hey, i am new here, i have dicom images i need to get 4 max pixel values and their co-ordinates and aromatically crop the 128 by 128 4 patches from that image keeping the center pixel one of the max pixel that has been found? please how can i do it

[maida]

no ones replying :/ :'( :/ :'(

Answer 1

[Image Analyst]

To get the 4 max values, sort the image in descending order:

sortedValues = sort(grayImage, 'descend');
% Get the 4 max values and their coords
for k = 1 : 4
    thisValue = sortedValues(k);
    % Find where it occurs
    [rows, columns] = find(grayImage, thisValue);
    % Plot them over the image
    for k2 = 1 : length(rows)
        thisRow = rows(k2);
        thisColumn = columns(k2);
        plot(thisColumn, thisRow, 'r+');
        hold on;
        text(thisColumn, thisRow, num2str(k));
        % Crop into a new image
        row1 = thisRow - 64;
        row2 = row1 + 127;
        col1 = thisColumn - 64;
        col2 = col1 + 127;
        subImage = grayImage(row1:row2, col1:col2);
        % Now do something with subimage....
    end
end

[maida]

subImage get the empty matrix, why is that so?