医学成像原理实验报告 (2)
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验一DICOM图像的读取和显示
DICOM（Digital Imaging and Communications in Medicine）即医学数字成像和通信，是医学图像和相关信息的国际标准。它定义了质量能满足临床需要的可用于数据交换的医学图像格式，利用不同的灰度值实现成像。
DICOM被广泛应用于放射医疗，心血管成像以及放射诊疗诊断设备（X射线，CT，核磁共振，超声等），并且在眼科和牙科等其它医学领域得到越来越深入广泛的应用。当前大约有百亿级符合DICOM标准的医学图像用于临床使用。
I = dicomread('CT-MONO2-16-ankle.dcm');
info = dicominfo('CT-MONO2-16-ankle.dcm');
I = dicomread(info);
imshow(I,'DisplayRange',[]);
dicomwrite(I,'ankle.dcm');
info = dicominfo('CT-MONO2-16-ankle.dcm');
I = dicomread(info);
dicomwrite(I,'ankle.dcm',info);
（2）图像读取
程序如下：
I = dicomread('CT-MONO2-16-ankle.dcm');
imtool(I,'DisplayRange',[])
info = dicominfo('CT-MONO2-16-ankle.dcm');
info.SeriesInstanceUID
max(I(:))
min(I(:))
Imodified = I;
Imodified(Imodified == 4080) = 32;
imshow(Imodified,[])
20
40
60
80
100
120
20
40
60
80
100
120
20
40
60
80
100
120
20
40
60
80
100
120
20
40
6080100120
20
40
6080
100
120
20406080100120
20
40
60
80
100
120
实验二 MRI 图像显示和读取
MRI 可获得人体横面、冠状面、矢状面及任何方向断面的图像，实现三维定位图像。 完整程序如下：
load mri ; D = squeeze(D); figure('Colormap',map) image_num = 8;
image(D(:,:,image_num)) axis image ; x = xlim; y = ylim
Image num=8
Image num=2 Image num=12 Imagenum=18
20
40
60
80
100
120
2-D Contour Slice
figure('Colormap',cm)
contourslice(D,[],[],image_num) axis ij xlim(x) ylim(y)
daspect([1,1,1])
Displaying 3-D Contour Slices
igure('Colormap',cm)
contourslice(D,[],[],[1,12,19,27],8); view(3); axis tight
Applying an Isosurface to the MRI Data
figure('Colormap',map) Ds = smooth3(D);
hiso = patch(isosurface(Ds,5),... FaceColor',[1,.75,.65],... EdgeColor','none'); isonormals(Ds,hiso)
hcap = patch(isocaps(D,5),... FaceColor','interp',... EdgeColor','none'); view(35,30) ; axis tight ;
20
40
60
80
100
120
daspect([1,1,.4]); lightangle(45,30);
set(gcf,'Renderer','zbuffer'); lighting phong set(hcap,'AmbientStrength',.6)
set(hiso,'SpecularColorReflectance',0,'SpecularExponent',50)
程序：
figure('Colormap',map)
Ds = smooth3(D);
hiso = patch(isosurface(Ds,5),... 'FaceColor',[1,.75,.65],... 'EdgeColor','none');
isonormals(Ds,hiso)
实验三平行束投影仿真
实验原理:sheep-logan 模型用来模拟头部断层图像，通过得到sheep-logan 头型，得到不同方向上投影数据，利用模拟X 光平行束重建图像。 程序如下：
clc; clear all ; close all ; N=256; I=phantom(N); figure;