Lecture d'image DICOM et définition de ROI

Fermé
kaloskagathos - 8 avril 2013 à 17:20
Bonjour,

Je suis étudiant en M2 de physique médicale et je dois écrire un programme qui doit me permettre d'ouvrir plusieurs images DICOM d'IRM (.ima) en même temps puis de définir x région d'intérêt (ROI) sur chaque image afin d'en ressortir des mesures sous formes graphiques. Pensez-vous pouvoir m'aider? J'ai déjà commencer à programmer mais je n'arrive pas à ouvrir plusieurs images en même temps, ça ne marche que pour une image:
clear all;clc;

setappdata(0,'UseNativeSystemDialogs',false)

[FileName,PathName] = uigetfile('.ima','Select the picture','MultiSelect','on');

addpath(PathName);%path(path,'PathName');
% garde le nom dans les preference
% N=length(FileName);
%
% for j=1:N
% img(j)=dicomread(FileName(j));
% end
I=dicomread(FileName);
% lis l'image

R=double(I)+double(edge(I));%pas obligatoire je peux rester uniquement avec I

% les crochets sont pour recaler les niveaux de gris
figure;
imshow(R,[])%comme dit plus pareil avec I
caxis([min(min(R)) max(max(R))]);%juste pour centrer mieux mais pas nécessaire

% info = dicominfo('FANTOME_DQPRM_.MR.NEURO_CRANE_TYPE.0003.0001.2013.03.27.12.23.15.750000.73042279.IMA');
% Y = dicomread(info);
% figure, imshow(Y, 'DisplayRange',[]);

% ROI = roipoly; % sélectionne ROI sur l'image activée, clique droit
% edge; % calcul du gradient de l'image