舊題新作——基于MATLAB的切片三維重建

馬軍 高俊峰 劉天成

摘要：本文基于2001年數學建模國賽A題，通過圖片分割、圖像處理、骨架提取、數據擬合等方法得到了切片的三維重建。運用分水嶺分割算法提取輪廓，再用圖像骨架提取法對切片圖像進行骨架提取，求出骨架線上所有的點到輪廓邊緣上所有點間的距離，采用極大極小距離法求出切片最大內切圓圓心及其管道半徑。

關鍵詞：分水嶺分割算法;骨架提取;MATLAB圖像信息處理

一 問題重述（略）

二 模型假設

1.假設切片管道是由球心沿中軸線且半徑固定的球滾動包絡而成。

2.假設切片中軸線與每張切片有且只有一個交點，中軸線與切片表面不平行。

三 問題分析

問題一，管道中軸線與每張切片僅有一個交點，每張切片必然存在最大內切圓，且圓心為中軸線與每張切片的交點，內切圓半徑即為管道的半徑。

四 模型的建立與求解

4.1 問題一：管道的中軸線與半徑的求解

4.1.1 切片最大內切圓的證明

4.1.3 切片圖像骨架提取算法

（1）算法原理：骨架指切片的中軸線。骨架細化是指從原來的圖中去掉一些點，仍保持原圖的骨架的一種抽取技術。

（2）基本步驟：①骨架提取算法流程：開始-MATLAB讀圖-用graythreash函數計算閾值-用im2bw函數將圖像變為二值圖像-調用bwmorph去除邊界-骨架提取-結束。②運用MATLAB進行函數運算，得到切片圖像骨架提取結果。

五 模型的評價及改進

首先，分水嶺算法提取邊緣效果優于edge函數得到的結果。其次，最大內切圓我們給出了圖像以及證明。最后，利用多項式擬合曲線方程，能較為快捷的找到圓心以及圓的半徑，繼而較為快捷的找到中心軸的曲線方程。

參考文獻

[1].于萬波.《基于MATLAB的圖像處理》[M].北京：清華大學出版社，2008：67-72.

[2].張德豐.《詳解MATLAB數字圖像處理》[M].北京：電子工業出版社，2010：257-269.

（作者單位：陜西科技大學）