[發明專利]基于B樣條的醫學圖像傾斜校正方法

在審

• 申請號: 201810153351.4
• 申請日: 2018-02-11
• 公開（公告）號: CN108460738A
• 公開（公告）日: 2018-08-28
• 發明（設計）人: 潘梅森
• 申請（專利權）人: 湖南文理學院
• 主分類號: G06T5/00
• 分類號: G06T5/00;G06T7/13
• 代理公司: 北京風雅頌專利代理有限公司 11403
• 代理人: 曾志鵬

一種基于B樣條的醫學圖像傾斜校正方法，包括：獲取待校正的醫學圖像，并對所述待校正的醫學圖像進行二值化處理以獲取二值化的醫學圖像；根據B樣條梯度算子對所述二值化的醫學圖像進行邊緣檢測，并獲得二值化的邊緣圖像；計算所述邊緣圖像的包圍盒，以得到所述邊緣圖像的實際大小，進而根據所述邊緣圖像的實際大小截取所述待校正的醫學圖像的子圖像；計算所述子圖像的質心和歸一化的中心矩，并得到旋轉角；將待校正的醫學圖像按照所述旋轉角反向旋轉，以獲取校正后的醫學圖像，能夠實現對灰度不均勻的傾斜圖像的校正。

技術領域

本發明涉及醫學圖像處理技術領域，特別涉及一種基于B樣條的醫學圖 像傾斜校正方法以及一種非臨時性計算機可讀存儲介質。

背景技術

在醫學圖像成像過程中，尤其在CT(Computed Tomography，電子計算 機斷層掃描)和MRI(Magnetic Resonance，磁共振成像)成像過程中，由于 存在相關設備欠穩定、患者隨意欲動或情緒緊張等問題，使生成的圖像常常 存在較為嚴重的傾斜現象，這給后續的圖像匹配和圖像融合帶來極大的困難。

為了解決圖像傾斜的問題，國內外專家學者進行了各種有益的探索，提 出一些較為實用的方法，例如投影法(projectionprofile)、霍夫變換法(Hough transform，霍夫變換)和近鄰法(component nearest neighbor clustering)。其中， 投影法是常用的傾斜校正方法之一，它是基于對投影圖形形狀進行分析的方 法，但是由于需要計算每個傾斜角的投影形狀，因此這種方法的計算量非常 大；Hough變換法是借助于Hough變換，用圖像空間目標像素的坐標去計算 參數空間中參考點的可能軌跡，該方法對于直線圖形效果較好，校正精度較 高；近鄰法是通過找出所有連通區中心點的K個最近鄰點，計算每對最近鄰 點的矢量方向并統計生成直方圖，直方圖的峰值對應于整個頁面的傾斜角， 由于在圖像中連通成分較多，因此也相當費時。

但是，上述方法存在的問題是，無法對灰度不均勻的傾斜圖像進行傾斜 面校正。

發明內容

本發明的目的是提供一種基于B樣條的醫學圖像傾斜校正方法，能夠實 現對灰度不均勻的傾斜圖像的校正。

為解決上述問題，本發明的第一方面提供了一種基于B樣條的醫學圖像 傾斜校正方法，包括以下步驟：獲取待校正的醫學圖像，并對所述待校正的 醫學圖像進行二值化處理以獲取二值化的醫學圖像；根據B樣條梯度算子對 所述二值化的醫學圖像進行邊緣檢測，并獲得二值化的邊緣圖像；計算所述 邊緣圖像的包圍盒，以得到所述邊緣圖像的實際大小，進而根據所述邊緣圖 像的實際大小截取所述待校正的醫學圖像的子圖像；計算所述子圖像的質心 和歸一化的中心矩，并得到旋轉角；將待校正的醫學圖像按照所述旋轉角反向旋轉，以獲取校正后的醫學圖像。

進一步地，所述根據B樣條梯度算子對所述二值化的醫學圖像進行邊緣 檢測采用如下公式：

其中，▽f為像素點(x,y)的梯度大小，B_x和B_y為B樣條梯度算子。

進一步地，所述計算所述二值化的邊緣圖像的包圍盒為計算所述二值化 的邊緣圖像的上、下、左、右的邊界。

進一步地，所述子圖像的質心通過所述子圖像的零階矩和一階矩計算。

進一步地，所述旋轉角采用如下公式獲得：

其中，M′_1,1、M′_2,0和M'_0,2為中心矩。

根據本發明的另一個方面提供一種非臨時性計算機可讀存儲介質，其上 存儲有計算機程序，該程序被處理器執行時實現所述的醫學圖像傾斜校正方 法。