血管的三维重建

据血管序列切片的二值图像特征,利用MATLAB丰富的矩阵运算和图像处理命令,将血管三维重建过程分成半径搜索、交点定位和轴线拟合三个主要步骤,并编制通用M程序包实现从数据采集到模型渲染的全程自动计算.最后应用给出的方法完成了对100张序列切片图像的计算机三维重建.

%附录7：求拟合图与原切片图的重合度的matlab程序代码（baifenbi1.m）
function baifenbi=baifenbi1（pnjjpn）%输出拟合图与原切片图的重合度
%填充新图
juzheng=pnjj; %pnjj为通过dian.m得到的轮廓边界二值矩阵
%先填充左边界的右半部分
for i=1:511
     for j=1:511
         if（pnjj（ij）==0&pnjj（ij+1）~=0）
        pnjj（ij+1）=pnjj（ij）;        
    end
  end
end
you=pnjj; 
%再填充右边界的左半部分 
for i=1:512   
     for j=512:-1:2          
         if（juzheng（ij）==0&juzheng（ij-1）~=0） 
             juzheng（ij-1）=juzheng（ij）;  
         end
     end
end
zuo=juzheng; 
shijijuzheng=you|zuo; %通过矩阵的或运算得到填充后的新图 
imshow（you|zuo） 
%原图的黑点的个数 
biaozhunjuzheng=imread（strcat（int2str（pn）‘.bmp‘））;
nbiao=0; 
for i=1:512     
    for j=1:512          
        if（biaozhunjuzheng（ij）==0）            
            nbiao=nbiao+1;         
        end
    end
end
%新图与原图重合部分黑点的个数 
chonghegs=0;
for i=1:512     
    for j=1:512
        if（biaozhunjuzheng（ij）==0&shijijuzheng（ij）==0）
           chonghegs=chonghegs+1;
        end
    end
end
%求百分比
baifenbi=chonghegs/nbiao;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1139  2015-04-15 09:56  三维重建\baifenbi1.m
     文件        1472  2015-04-15 09:49  三维重建\dian.m
     文件        1112  2015-04-15 15:50  三维重建\ff.m
     文件         280  2015-04-14 18:32  三维重建\pczx.m
     文件         452  2015-04-15 15:58  三维重建\xueguansanweitu.m
     文件         702  2015-04-15 09:58  三维重建\zhongzhouxiantouying.m
     文件         231  2015-04-15 09:58  三维重建\zhuanhua.m
     目录           0  2015-04-15 16:03  三维重建\