1. dicom 파일 2개 여는 방법 -> 2개 파일의 size 는 같아야함 (ex: 512*512)
InputDicomA = dicomread('A.dcm');
InputDicomB = dicomread('B.dcm');
2. B영상에서 A영상을 빼는 방법 -> matrix 연산과 동일함
OutputDicomC = InputDicomB- InputDicomA ;
3. 처리가 완료된 영상을 dicom 파일로 저장하는 방법
dicomwrite(OutputDicomC, 'output.dcm');
4. 정확한 영상 출력을 위해 각 영상의 최소/최대값을 확인 (생략해도 무방함)
min_valueA = min(InputDicomA(:));
max_valueA = max(InputDicomA(:));
min_valueB = min(InputDicomB(:));
max_valueB = max(InputDicomB(:));
min_valueC = min(OutputDicomC(:));
max_valueC = max(OutputDicomC(:));
5. 여러 개의 영상을 하나의 window (창) 에 출력하는 방법 -> subplot 함수
6. 영상 출력 방법 -> imshow 함수
7. 각 영상의 제목 -> title 함수
subplot(1,3,1)
imshow(InputDicomA,[min_valueA, max_valueA]);
title(' Input A ');
subplot(1,3,2)
imshow(InputDicomB,[min_valueB, max_valueB]);
title(' Input B ');
subplot(1,3,3)
imshow(OutputDicomC, [min_valueC, max_valueC]);
title(' Output C ');
8. 의료 영상인 Dicom 에 들어가있는 모든 정보를 확인하는 방법 -> dicominfo 함수
InfoA = dicominfo('A.dcm');
InfoB = dicominfo('B.dcm');
InfoC = dicominfo('output.dcm');
==============================================================
ⓒ 2016. All rights reserved.
Created by: E.V.Joo.
Website: evjoo.tistory.com
==============================================================