文档介绍：血管的三维重建三
断面可用于了解生物组织、器官等的形态。例如，将样本染色后切成厚约1m m的切片，在显微镜下观察该横断面的组织形态结构。如果用切片机连续不断地将样本切成数十、成百的平行切片， 可依次逐片观察。根据拍照并采样得到的平行切，该圆的半径等于滚动球的半径。
基于：
1）球的任意截面都是圆
2）经过球心的球截面是所有
截圆当中半径最大的圆
问题分析
9
方法2：利用patch()
for b=0:99
m(:,:,b+1)=[imread([int2str(b),'.bmp'],'bmp')]';
end
m1=m(1:4:512,1:4:512,:);
for i=1:size(m1,1)
for j=1:size(m1,2)
for k=1:size(m1,3)
if m1(i,j,k)==0
m1(i, j, k)=88;
else
m1(i,j,k)=0;
end, end, end,end
ms=smooth3(m1);
hiso=patch(isosurface(ms,5),'FaceColor',[1,.75,.65],'EdgeColor','none');
hcap=patch(isocaps(m1,5),'FaceColor','interp','EdgeColor','none');
colormap('default')
view(45,30), axis tight
grid
daspect([1,1,.4])
lightangle(45,30); lighting phong
isonormals(ms, hiso)
set(hcap,'AmbientStrength',.6)
set(hiso,'SpecularColorReflectance',0,'SpecularExponent',50)
rotate3d
15
5）如何由切片图叠合作出血管的立体图？
方法2：利用patch()
16
中轴线的建模：求中轴线与各横断面的交点和曲线拟合、逼近。参赛者使用的方法主要有：
1）枚举法
2）平行切线法
3）外推法
4）滚球法
5）投影法
6）变换法
7）细化法
建模方法思想
17
方法一：以黑色点为圆心，以r0为半径画圆，判断此圆周上的象素点是否为黑色．若全为黑色，则增大半径再判断，直到圆周上有白色象素点，记录该圆的圆心坐标和半径，然后，取下一个黑色象素点重复以上步骤，直至黑色区域中所有象素点都搜索完为止；所有记录点中对应于最大半径的圆心坐标，就是该切片的最大内切圆圆心坐标，此半径即是血管的半径
具体实现时出现的问题
18
方法一得出的中轴线投影：
具体实现时出现的问题
0
20
40
60
80
100
120
140
160
180
0
10
20
30
40
50
60
70
80
90
100
y
z
0
10
20
30
40
50
60
70
80
90
100
-200
-150
-100
-50
0
50
100
150
200
z
x
中轴线在ZX平面投影图 中轴线在YZ平面投影图
什么原因造成内切圆圆心的误差如此之大？
19
模型的检验
不检验只能说完成问题的一半。由于所给图象只是离散数据，以象素为单位。无论以何种方式建模，都是近似计算，效果如何，检验很必要。还可通过检验，发现模型的误差，修正模型，提高模型的正确性。阅卷中发现的检验方法主要有：
1）逐片比较
2）法平面法
3）滚动法
建模方法思想
20
模型检验时也发现内切圆圆心误差较大，特别是后面的切片。
建模方法思想
Z=30
Z=0
21
模型检验时也发现内切圆圆心误差较大，特别是后面的切片。
建模方法思想
Z=50
Z=40
22
模型检验时也发现内切圆圆心误差较大，特别是后面的切片。
建模方法思想
Z=60
Z=57
23
误差大的原因和改进途径
误差大的原因：
1）图象误差
实际图象边界上的点是连续的，在转换成bmp图象时，象素表示的图象边界是离散的，成锯齿状，与实际图象有误差（舍入误差）。
2）同一张切片上的最大内切圆不唯一
解决办法：
1）方法一：取平均
求出同一张切片上的所有最大内切圆的圆心，然后求平均值。
24
误差大的原因和改进途径
1）方法一：取平均
求出同一张切片上的所有最大内切圆的圆心，然后求平均值。