How to read only Axial view of lung CT DICOM image in matlab?

I have a database of lung Ct scans in Dicom format. In ITK snap 3 views of image are opened. I just want to open axial view of CT in matlab. How can I do that?

[harjeet singh]

dear zubda, do upload an example dicom format

[Zubda Butt]

This DCM image is opened in ITK SNAP . I want to process this image in matlab using only its axial view which is in top left corner.

[harjeet singh]

now do you want automatic clip top left corner image using matlab?

[Zubda Butt]

yes.

Answer 1

[harjeet singh]

dear zubda do use this code

clear all
close all
clc

a=imread('ct.png');
figure(1)
imshow(a)

b=~(a(:,:,1)==240 & a(:,:,2)==240 & a(:,:,3)==240);
figure(2)
imshow(b);
[m,n]=size(b);
[lab,num]=bwlabel(b);
for i=1:num
    [r,c]=find(lab==i);
    loco(i,:)=[i mean(r) mean(c)];
end
loco(:,4)=sqrt( (loco(:,2)-m*0.2885).^2 + (loco(:,3)-n*0.7794).^2);
loco=sortrows(loco,4);
[r,c]=find(lab==loco(1,1));

c=a(min(r):max(r),min(c):max(c),:);

figure(3)
imshow(c)
drawnow

[Zubda Butt]

Thank you Sir for your help. Code is working. :) Can you please guide me to some good matlab tutorial source where I can learn and understand code and get help for errors in my code?

[harjeet singh]

try to solve as many tasks as you can for basics of matlab, do evolve your own little tasks whatever you think and do try to solve that with matlab. rest matlab help is the perfect for knowing commands and fuctions.