hi. i need to automatically crop the interested region of the image. The images are the scanned MRI images. i want to segmentation prostate and i need to remove other organs .i want to use level set method . plz help me to solve this project.
MATLAB: How to crop automatically the interested region of an image
croppingImage Processing Toolbox
Related Solutions
OK, try this:
clc; % Clear the command window.
close all; % Close all figures (except those of imtool.)
clear; % Erase all existing variables. Or clearvars if you want.
workspace; % Make sure the workspace panel is showing.
format long g;format compact;fontSize = 25;%===============================================================================
% Get the name of the image the user wants to use.
% Have user browse for a file, from a specified "starting folder."
% For convenience in browsing, set a starting folder from which to browse.
startingFolder = pwd;if ~exist(startingFolder, 'dir') % If that folder doesn't exist, just start in the current folder.
startingFolder = pwd;end% Get the name of the file that the user wants to use.
defaultFileName = fullfile(startingFolder, '*.*');[baseFileName, folder] = uigetfile(defaultFileName, 'Select a file');if baseFileName == 0 % User clicked the Cancel button.
return;end% Get the full filename, with path prepended.
fullFileName = fullfile(folder, baseFileName);%===============================================================================% Read in a demo image.
grayImage = imread(fullFileName);% Get the dimensions of the image.
% numberOfColorChannels should be = 1 for a gray scale image, and 3 for an RGB color image.
[rows, columns, numberOfColorChannels] = size(grayImage)if numberOfColorChannels > 1 % It's not really gray scale like we expected - it's color.
% Use weighted sum of ALL channels to create a gray scale image.
grayImage = rgb2gray(grayImage); % ALTERNATE METHOD: Convert it to gray scale by taking only the green channel,
% which in a typical snapshot will be the least noisy channel.
% grayImage = grayImage(:, :, 2); % Take green channel.
end% Display the image.
subplot(1, 2, 1);imshow(grayImage, []);axis on;axis image;caption = sprintf('Original Gray Scale Image');title(caption, 'FontSize', fontSize, 'Interpreter', 'None');drawnow;hp = impixelinfo();% Set up figure properties:
% Enlarge figure to full screen.
set(gcf, 'Units', 'Normalized', 'OuterPosition', [0 0 1 1]);% Get rid of tool bar and pulldown menus that are along top of figure.
% set(gcf, 'Toolbar', 'none', 'Menu', 'none');
% Give a name to the title bar.
set(gcf, 'Name', 'Demo by Mark Hayworth', 'NumberTitle', 'Off')drawnow;uiwait(msgbox('Drag out a 128x128 box'));% Let user draw the box.
k = waitforbuttonpress;point1 = get(gca,'CurrentPoint'); % button down detected
finalRect = rbbox; % return figure units
point2 = get(gca,'CurrentPoint'); % button up detected
point1 = point1(1,1:2); % extract x and y
point2 = point2(1,1:2);x1 = point1(1);y1 = point1(2);x2 = point2(1);y2 = point2(2);% Find the coordinates of the box.
xUserBox = [x1, x2, x2, x1, x1];yUserBox = [y1, y1, y2, y2, y1];% The box from rbbox() disappears after drawing.
% Redraw what they drew.
hold onaxis manualplot(xUserBox, yUserBox, 'r-', 'LineWidth', 2); % Draw box.
% Now we need to get a 128x128 box centered at that location.
% Compute the center
xCenter = (x1 + x2)/2; % Get average.
yCenter = (y1 + y2)/2; % Get average.plot(xCenter, yCenter, 'r+', 'MarkerSize', 30, 'LineWidth', 2); % Plot cross at center.
% Get box edges for a 128x128 box centered at the center of the box that they drew.
x1 = round(xCenter - 64)y1 = round(yCenter - 64)x2 = x1 + 127y2 = y1 + 127% Find the coordinates of the box.xBox128 = [x1, x2, x2, x1, x1];yBox128 = [y1, y1, y2, y2, y1];% Draw the 128x128 box
plot(xBox128, yBox128, 'b-', 'LineWidth', 2);% Display the image.croppedImage = grayImage(y1:y2, x1:x2);[rows, columns, numColors] = size(croppedImage) % Show values just to make sure.
subplot(1, 2, 2);imshow(croppedImage, []);axis on;title('Cropped 128x128 Image', 'FontSize', fontSize, 'Interpreter', 'None');
Make it integer like it's asking:
rect = getrect(gcf);% Might be floating point, so round:
rect = int32(rect);% Make sure it's never zero because an index of 0 is not allowed.
rect(rect == 0) = 1;Related Question
- Imrect function how I can show the image
- Converting a rectangle position on an image to image matrix index
- Select particular region of image without imcrop
- How to draw a free hand box on an image to select that region of interest
- Imrect not visible after updating an image
- Hi, I am using Matlab version 7.9.0(R2009b). I need to select a rectangular region of interest for the tif image. I tried imrect but the rectangular appeared on an empty coordinat system, not in the tif image. Could you please help me on this
- Making a rectangle with one mouse click
- Editing images using MATLAB
Best Answer