vessel segmentation in brain DSA and MRA image


Vessel segmentation in brain 3D DSA(Dyna CT) and MRA image is the first step of the project.

It is quite fortunate that the image quality is ideal, so the auto grayscale in Matlab imadjuct function is good enough to segment vessels.

The effect of segmentation and code for segmenting listed below.

600.jpg

3D Dyna CT image segmentation result of AVM patient

geshi={'*.dcm','Dicom image (*.dcm)';...
       '*.bmp','Bitmap image (*.bmp)';...
       '*.jpg','JPEG image (*.jpg)';...
       '*.*','All Files (*.*)'};
[FileName FilePath]=uigetfile(geshi,'Load MR files','*.dcm','MultiSelect','on');
if ~isequal([FileName,FilePath],[0,0]);
    FileFullName=strcat(FilePath,FileName);
    if  ~ischar(FileFullName)
        FileFullName=FileFullName([2:end 1])';
    end
else
    return;
end
MRo=[];
MR=[];
n=length(FileFullName);
for i=1:n
    I=dicomread(FileFullName{i});
    metadata=dicominfo(FileFullName{i});
    serialno=metadata.InstanceNumber;
    %I=rgb2gray(I);
    %
    MRo(:,:,serialno)=I;
    MR(:,:,serialno)=im2bw(imadjust(mat2gray(I),[],[]));
end
%calculate the mask of outside bone
geshi={'*.dcm','Dicom image (*.dcm)';...
       '*.bmp','Bitmap image (*.bmp)';...
       '*.jpg','JPEG image (*.jpg)';...
       '*.*','All Files (*.*)'};
[FileName FilePath]=uigetfile(geshi,'Load DSA files','*.dcm','MultiSelect','on');
if ~isequal([FileName,FilePath],[0,0]);
    FileFullName=strcat(FilePath,FileName);
    if  ~ischar(FileFullName)
        FileFullName=FileFullName([2:end 1])';
    end
else
    return;
end
DSAo=[];
DSA=[];
n=length(FileFullName);
for i=1:n
    I=dicomread(FileFullName{i});
    metadata=dicominfo(FileFullName{i});
    serialno=metadata.InstanceNumber;
    %I=rgb2gray(I);
    %
    DSAo(:,:,serialno)=I;
    DSA(:,:,serialno)=im2bw(imadjust(mat2gray(I),[],[]));
end

Code also available in Github: https://github.com/clatfd/avm/blob/remote/find_vessel.m


This is a part of my Capstone Research

Last Article Next article

Comment 评论



Share 分享

New Users 最新加入

  • refrain

  • 鑫鑫

  • :)

New comments 最新评论

test123: aasdas Details Apr 13 16:39
admin: Thanks! Details Apr 09 11:46
admin: Google map api Details Apr 09 11:46
lqj12: cooooooooool Details Apr 08 21:34
Yunhan Huang: 这个功能是如何实现的? Details Apr 08 13:23