專利名稱：基于彈性成像的圖像處理方法及系統(tǒng)的制作方法
基于彈性成像的圖像處理方法及系統(tǒng)技術(shù)領(lǐng)域：
本發(fā)明涉及圖像處理技術(shù)，特別是涉及一種基于彈性成像的圖像處理方法及系統(tǒng)。背景技術(shù)：
隨著科學(xué)技術(shù)的進(jìn)步，彈性成像技術(shù)在醫(yī)療衛(wèi)生領(lǐng)域得到了迅速的發(fā)展。彈性成像技術(shù)能夠得到組織機(jī)械力學(xué)的信息，豐富了模態(tài)成像內(nèi)容。彈性成像，是對組織施加一個外部激勵，在壓力-形變等物理規(guī)律作用下，組織將產(chǎn)生一個響應(yīng)，根據(jù)該響應(yīng)得到彈性成像圖。
彈性成像反映的是組織整體的彈性概況，目前尚沒有對于彈性成像的定量后處理，或?qū)椥猿上竦木植繄D像區(qū)域(或者是重點(diǎn)關(guān)注區(qū)域)處理的相關(guān)研究，在圖像數(shù)據(jù)的精確解析和利用方面尚有進(jìn)一步提升的空間。
發(fā)明內(nèi)容
基于此，有必要提供一種能夠?qū)椥猿上竦木植繀^(qū)域進(jìn)行處理且更加精確的基于彈性成像的處理方法。
另外，還有必要提供一種能夠?qū)椥猿上竦木植繀^(qū)域進(jìn)行處理且更加精確的基于彈性成像的處理系統(tǒng)。
一種基于彈性成像的圖像處理方法，包括以下步驟獲取彈性圖像；從所述彈性圖像中選取關(guān)注區(qū)域，并從所述關(guān)注區(qū)域提取灰度分布圖；對所述灰度分布圖的曲線進(jìn)行平滑處理；通過對所述灰度分布圖曲線上點(diǎn)的斜率劃分得到過渡區(qū)及對照區(qū)。
優(yōu)選地，所述對所述灰度分布圖的曲線進(jìn)行平滑處理的步驟為按照所述灰度分布圖曲線上的選取點(diǎn)分別提取相鄰點(diǎn)，所述相鄰點(diǎn)包括預(yù)設(shè)數(shù)量的上相鄰點(diǎn)和下相鄰點(diǎn)； 計(jì)算獲得上相鄰點(diǎn)和下相鄰點(diǎn)的灰度平均值，并分別與所述選取點(diǎn)的灰度值做差；判斷所述上相鄰點(diǎn)與所述選取點(diǎn)的差值的絕對值是否大于所述下相鄰點(diǎn)與所述選取點(diǎn)的差值的絕對值，是，則所述選取點(diǎn)的灰度值與下相鄰點(diǎn)的灰度平均值相同。
優(yōu)選地，所述通過對所述灰度分布圖曲線上點(diǎn)的斜率劃分得到過渡區(qū)及對照區(qū)的步驟為獲取所述灰度分布圖曲線上點(diǎn)的斜率值；將所述斜率值與預(yù)設(shè)的閾值范圍進(jìn)行比較，判斷所述斜率值是否超過所述閾值范圍，是，則所述斜率值對應(yīng)的所述灰度分布圖曲線的區(qū)域?yàn)閷φ諈^(qū)，否，則為過渡區(qū)。
優(yōu)選地，所述閾值范圍為0 0.5。
優(yōu)選地，所述通過對所述灰度分布圖曲線上點(diǎn)的斜率劃分得到過渡區(qū)及對照區(qū)的步驟之后還包括計(jì)算所述過渡區(qū)的面積和所述對照區(qū)的面積；所述過渡區(qū)的面積與所述對照區(qū)的面積的比值與預(yù)設(shè)的面積閾值進(jìn)行比較，判斷是否大于所述面積閾值，是，則所述過渡區(qū)為重點(diǎn)關(guān)注區(qū)。
一種基于彈性成像的圖像處理系統(tǒng)，包括彈性圖像模塊，用于加載彈性圖像；提取模塊，用于從所述彈性圖像中選取關(guān)注區(qū)域，并從所述關(guān)注區(qū)域提取灰度分布圖；平滑處理模塊，用于對所述灰度分布圖的曲線進(jìn)行平滑處理；劃分模塊，用于通過對所述灰度分布圖曲線上點(diǎn)的斜率劃分得到過渡區(qū)及對照區(qū)。
優(yōu)選地，所屬平滑處理模塊包括選取單元，用于按照所述灰度分布圖曲線上的選取點(diǎn)分別提取相鄰點(diǎn)，所述相鄰點(diǎn)包括預(yù)設(shè)數(shù)量的上相鄰點(diǎn)和下相鄰點(diǎn)；計(jì)算單元，用于計(jì)算獲得上相鄰點(diǎn)和下相鄰點(diǎn)的灰度平均值，并分別與所述選取點(diǎn)的灰度值做差；判斷單元， 用于判斷所述上相鄰點(diǎn)與所述選取點(diǎn)的差值的絕對值是否大于所述下相鄰點(diǎn)與所述選取點(diǎn)的差值的絕對值，是，則所述選取點(diǎn)的灰度值與下相鄰點(diǎn)的灰度平均值相同。
優(yōu)選地，所述劃分模塊包括斜率值獲取單元，用于獲取所述灰度分布圖曲線上點(diǎn)的斜率值；比較單元，用于將所述斜率值與預(yù)設(shè)的閾值范圍進(jìn)行比較，判斷所述斜率值是否超過所述閾值范圍，是，則所述斜率值對應(yīng)的所述灰度分布圖曲線的區(qū)域?yàn)閷φ諈^(qū)，否，則為過渡區(qū)。
優(yōu)選地，所述閾值范圍為0 0.5。
優(yōu)選地，還包括面積計(jì)算模塊，用于計(jì)算所述過渡區(qū)的面積和所述對照區(qū)的面積；面積判斷模塊，用于所述過渡區(qū)的面積與所述對照區(qū)的面積的比值與預(yù)設(shè)的面積閾值進(jìn)行比較，判斷是否大于所述面積閾值，是，則所述過渡區(qū)為重點(diǎn)關(guān)注區(qū)。
通過基于彈性成像的圖像處理方法及系統(tǒng)在彈性圖像中選取關(guān)注區(qū)域(局部區(qū)域)，并提取灰度分布圖，然后根據(jù)灰度分布圖曲線上點(diǎn)的斜率信息，獲得過渡區(qū)及對照區(qū)， 得到豐富的圖像定量數(shù)據(jù)，可以從繁多復(fù)雜的信息中有目的或有針對性的獲取信息，進(jìn)而獲得更加精確的信息。

圖1為一實(shí)施例中基于彈性成像的圖像處理方法的流程圖2為一實(shí)施例中步驟S30的具體流程圖3為舉例說明示意圖4為一實(shí)施例中步驟S40的具體流程圖5為另一實(shí)施例的方法流程圖6為測試組織的實(shí)物圖7為測試組織的沿指示線的灰度分布度；
圖8為測試組織的磁共振Tl加權(quán)的圖9為測試組織的磁共振Tl加權(quán)圖沿指示線的灰度分布圖10為測試組織的磁共振彈性圖像；
圖11為測試組織的磁共振彈性圖像的灰度分布圖12為另一組測試組織的實(shí)物圖13為另一組測試組織的沿指示線的灰度分布度；
圖14為另一組測試組織的磁共振Tl加權(quán)的圖15為另一組測試組織的磁共振Tl加權(quán)圖沿指示線的灰度分布圖16為另一組測試組織的磁共振彈性圖像；
圖17為另一組測試組織的磁共振彈性圖像的灰度分布圖18為附圖17進(jìn)行平滑處理的圖19為基于彈性成像的圖像處理系統(tǒng)的邏輯圖20為一實(shí)施例中基于彈性成像的圖像處理系統(tǒng)的邏輯圖。
具體實(shí)施方式
為了解決目前還沒有對彈性成像定量的后處理，或?qū)椥猿上竦木植繄D像區(qū)域 (或者是重點(diǎn)關(guān)注區(qū)域)處理的相關(guān)研究，在圖像數(shù)據(jù)的解析方面的缺陷，提供了一種基于彈性成像的圖像處理方法，詳細(xì)方案如下
結(jié)合附圖1，一種基于彈性成像的圖像處理方法，具體步驟如下
步驟SlO 獲取彈性圖像。彈性成像技術(shù)是一種較為成熟的成像技術(shù)，根據(jù)該技術(shù)可以獲得彈性圖像，該彈性圖像可以基于MRI (磁共振)或超聲彈性成像獲取的圖像。
步驟S20 從彈性圖像中選取關(guān)注區(qū)域，并從關(guān)注區(qū)域提取灰度分布圖。具體地， 在彈性圖像中在關(guān)注區(qū)域上，選取一“關(guān)注線”，沿該“關(guān)注線”從彈性圖像中提取灰度分布圖。
步驟S30 對灰度分布圖的曲線進(jìn)行平滑處理。對獲得的灰度分布圖多數(shù)為鋸齒狀曲線，不易于從該灰度分布圖獲取數(shù)據(jù)。通過對灰度分布圖的曲線進(jìn)行平滑處理后，易于對該關(guān)注區(qū)域的曲線進(jìn)行處理，得到豐富的圖像數(shù)據(jù)。
在一實(shí)施例中，結(jié)合附圖2，步驟S30具體為
步驟S31 按照灰度分布圖曲線上的選取點(diǎn)分別提取相鄰點(diǎn)，相鄰點(diǎn)包括預(yù)設(shè)數(shù)量的上相鄰點(diǎn)和下相鄰點(diǎn)。具體地，在灰度分布曲線上任意處獲得選取點(diǎn)，同時在該選取點(diǎn)相鄰點(diǎn)處選取預(yù)設(shè)數(shù)量的上相鄰點(diǎn)和下相鄰點(diǎn)，該預(yù)設(shè)數(shù)量可以為2個或3個。
步驟S32 計(jì)算獲得上相鄰點(diǎn)和下相鄰點(diǎn)的灰度平均值，并分別與選取點(diǎn)的灰度值做差。具體地，在灰度分布圖曲線上獲得預(yù)設(shè)數(shù)量的上相鄰點(diǎn)和下相鄰點(diǎn)的灰度值，同時根據(jù)算數(shù)平均法計(jì)算出上相鄰點(diǎn)和下相鄰點(diǎn)的灰度平均值。上相鄰點(diǎn)和下相鄰點(diǎn)的灰度值與選取點(diǎn)的灰度值做差，獲得差值。
步驟S33 判斷上相鄰點(diǎn)與選取點(diǎn)的差值的絕對值是否大于下相鄰點(diǎn)與選取點(diǎn)的差值的絕對值，是，則選取點(diǎn)的灰度值與下相鄰點(diǎn)的灰度平均值相同。具體地，獲取上相鄰點(diǎn)的差值與選取點(diǎn)的差值的絕對值，以及下相鄰點(diǎn)的差值與選取點(diǎn)的差值的絕對值，并判斷絕對值大小。是，則選取點(diǎn)的灰度值與下相鄰點(diǎn)的灰度平均值相同；否，則選取點(diǎn)的灰度值與上相鄰點(diǎn)的灰度平均值相同。
結(jié)合具體例子說明，選取點(diǎn)的相鄰點(diǎn)的預(yù)設(shè)數(shù)量為2個，結(jié)合附圖3，具體如下
在灰度分布圖曲線上隨機(jī)的獲得選取點(diǎn)，具體的選取點(diǎn)為A，其2個上相鄰點(diǎn)為 A1和A2, 2個下相鄰點(diǎn)為L1和A_2 ；所對應(yīng)的上相鄰點(diǎn)的灰度值為F1和F2, 2個下相鄰點(diǎn)的灰度值為F—和F_2，A的灰度值為F。上相鄰點(diǎn)對應(yīng)的灰度平均值為Fs = (FJF2)/2 ；下相鄰點(diǎn)對應(yīng)的灰度平均值為I7x = (F_i+F_2)/2。若F與Fs的差值的絕對值大于F與!^x的差值的絕對值，則選取點(diǎn)A所對應(yīng)的F變換為Fx，即與下相鄰點(diǎn)的!^x值相同。
步驟S40 通過對灰度分布圖曲線上點(diǎn)的斜率劃分得到過渡區(qū)及對照區(qū)。具體地， 獲取灰度分布圖曲線上點(diǎn)的斜率值，根據(jù)斜率值劃分區(qū)域數(shù)據(jù)分析區(qū)域，即過渡區(qū)及對照區(qū)，進(jìn)一步的從彈性圖像中的關(guān)注區(qū)域獲取豐富的圖像數(shù)據(jù)。
在一實(shí)施例中，結(jié)合附圖4，步驟S40具體為
步驟S41 獲取灰度分布圖曲線上點(diǎn)的斜率值。具體地，在灰度分布曲線的點(diǎn)上獲取對應(yīng)的點(diǎn)的斜率。
步驟S42 將斜率值與預(yù)設(shè)的閾值范圍進(jìn)行比較，判斷斜率值是否超過閾值范圍， 是，則斜率值對應(yīng)的灰度分布圖曲線的區(qū)域?yàn)閷φ諈^(qū)，否，則為過渡區(qū)。具體地，將斜率值與預(yù)設(shè)的閾值范圍進(jìn)行比較，在其它實(shí)施例中，該閾值范圍為0 0. 5，若超過該閾值范圍，則該斜率值對應(yīng)的灰度分布圖曲線與橫坐標(biāo)(距離)所包圍的區(qū)域?yàn)閷φ諈^(qū)；否，則為過渡區(qū)。
通過本方案在彈性圖像中選取關(guān)注區(qū)域，并提取灰度分布圖，然后根據(jù)灰度分布圖曲線上點(diǎn)的斜率信息，獲得過渡區(qū)域及對照區(qū)域?？梢?，對彈性圖像的關(guān)注區(qū)域進(jìn)行圖像處理，得到了豐富的圖像數(shù)據(jù)，而不僅僅是簡單的提供彈性圖像。
在另一實(shí)施例中，結(jié)合附圖5，還包括如下步驟
步驟S50 計(jì)算過渡區(qū)的面積和對照區(qū)的面積。具體地，采用微積分方法、填補(bǔ)發(fā)或估算法分別得到過渡區(qū)的面積和對照區(qū)的面積。
步驟S60 過渡區(qū)的面積與對照區(qū)的面積的比值與預(yù)設(shè)的面積閾值進(jìn)行比較，判斷是否大于面積閾值，是，則過渡區(qū)為重點(diǎn)關(guān)注區(qū)。具體地，根據(jù)渡區(qū)的面積與對照區(qū)的面積的比值，若超過面積閾值，例如為1/3，則確認(rèn)該過渡區(qū)為重點(diǎn)關(guān)注區(qū)域，為進(jìn)一步獲取信息提供準(zhǔn)確的指引；若沒有超過面積閾值，則為非重點(diǎn)關(guān)注區(qū)。
結(jié)合具體附圖6 11的實(shí)驗(yàn)數(shù)據(jù)，說明本申請的特點(diǎn)，具體地
附圖6為測試組織的實(shí)物圖，附圖7為該測試組織的沿指示線的灰度分布度；
附圖8為測試組織的磁共振Tl加權(quán)的圖，附圖9為測試組織的磁共振Tl加權(quán)圖沿指示線的灰度分布附圖10為測試組織的磁共振彈性圖像，附圖11為測試組織的磁共振彈性圖像的灰度分布圖。
清楚的看到實(shí)物圖的灰度分布圖為鋸齒狀，相應(yīng)的磁共振圖的灰度分布圖為小幅度波動狀，趨于均一，進(jìn)一步提取組織特征信息的難度較大；而磁共振彈性圖像的灰度分布圖為有一定斜率趨勢的曲線，信息豐富。
同時，實(shí)物圖的灰度分布圖、磁共振的灰度分布圖以及磁共振彈性圖像的灰度分布圖都不相一致，而實(shí)物圖和磁共振的灰度分布圖僅是從圖像本身反映對實(shí)物或圖像的認(rèn)知，而無法獲取豐富的數(shù)據(jù)。磁共振彈性圖像的灰度分布圖從另一個角度進(jìn)行分析，通過對該灰度分布圖的處理獲取豐富的數(shù)據(jù)。
結(jié)合另一組附圖12 17的實(shí)驗(yàn)數(shù)據(jù)，進(jìn)一步說明本申請的優(yōu)勢，具體地
附圖12為另一組測試組織的實(shí)物圖，附圖13為另一組測試組織的沿指示線的灰度分布度；
附圖14為另一組測試組織的磁共振Tl加權(quán)的圖，附圖15為另一組測試組織的磁共振Tl加權(quán)圖沿指示線的灰度分布附圖16為另一組測試組織的磁共振彈性圖像，附圖17為另一組測試組織的磁共振彈性圖像的灰度分布圖。
實(shí)物圖的灰度分布圖和磁共振常規(guī)圖像的灰度分布圖與磁共振彈性圖像的灰度分布相差較大。本申請方案從彈性成像中所獲取的信息較為豐富準(zhǔn)確，可對關(guān)注區(qū)域進(jìn)一步分析。
現(xiàn)例舉一具體方案，結(jié)合附圖18，劃分過渡區(qū)和對照區(qū)，以及辨別重點(diǎn)關(guān)注區(qū)和非重點(diǎn)關(guān)注區(qū)。
附圖18為附圖17進(jìn)行平滑處理的圖，由此可見，能獲得較為光滑的曲線圖。根據(jù)斜率值在0 0. 5范圍內(nèi)所對應(yīng)的灰度分布圖曲線的區(qū)域?yàn)閷φ諈^(qū)，不在該范圍內(nèi)的為過渡區(qū)。
同時，每個區(qū)域所對應(yīng)的面積為Si、S2和S3?？梢?，S1/S2和S3/S2大于1/3，則屬于重點(diǎn)關(guān)注區(qū)。
基于上述基于彈性成像的圖像處理方法，還有必要建立基于彈性成像的圖像處理系統(tǒng)，具體如下
結(jié)合附圖19，一種基于彈性成像的圖像處理系統(tǒng)，具體步驟如下
彈性圖像模塊10，用于加載彈性圖像。彈性成像技術(shù)是一種較為成熟的成像技術(shù)，根據(jù)該技術(shù)可以獲得彈性圖像并加載，該彈性圖像可以是磁共振彈性成像或超聲彈性圖像。
提取模塊20，用于從彈性圖像中選取關(guān)注區(qū)域，并從關(guān)注區(qū)域提取灰度分布圖。具體地，在彈性圖像中在關(guān)注區(qū)域上，選取一“關(guān)注線”，沿該“關(guān)注線”從彈性圖像中提取灰度分布圖。
平滑處理模塊30，用于對灰度分布圖的曲線進(jìn)行平滑處理。對獲得的灰度分布圖多數(shù)為鋸齒狀曲線，不易于從該灰度分布圖獲取數(shù)據(jù)。通過對灰度分布圖的曲線進(jìn)行平滑處理后，易于對該關(guān)注區(qū)域的曲線進(jìn)行處理，得到豐富的圖像數(shù)據(jù)。
在一實(shí)施例中，平滑處理模塊30包括
選取單元31，用于按照灰度分布圖曲線上的選取點(diǎn)分別提取相鄰點(diǎn)，相鄰點(diǎn)包括預(yù)設(shè)數(shù)量的上相鄰點(diǎn)和下相鄰點(diǎn)。具體地，在灰度分布曲線上任意處獲得選取點(diǎn)，同時在該選取點(diǎn)相鄰點(diǎn)處選取預(yù)設(shè)數(shù)量的上相鄰點(diǎn)和下相鄰點(diǎn)，該預(yù)設(shè)數(shù)量可以為2個或3個。
計(jì)算單元32，用于計(jì)算獲得上相鄰點(diǎn)和下相鄰點(diǎn)的灰度平均值，并分別與選取點(diǎn)的灰度值做差。具體地，在灰度分布圖曲線上獲得預(yù)設(shè)數(shù)量的上相鄰點(diǎn)和下相鄰點(diǎn)的灰度值，同時根據(jù)算數(shù)平均法計(jì)算出上相鄰點(diǎn)和下相鄰點(diǎn)的灰度平均值。上相鄰點(diǎn)和下相鄰點(diǎn)的灰度值與選取點(diǎn)的灰度值做差，獲得差值。
判斷單元33，用于判斷上相鄰點(diǎn)與選取點(diǎn)的差值的絕對值是否大于下相鄰點(diǎn)與選取點(diǎn)的差值的絕對值，是，則選取點(diǎn)的灰度值與下相鄰點(diǎn)的灰度平均值相同。具體地，獲取上相鄰點(diǎn)的差值與選取點(diǎn)的差值的絕對值，以及下相鄰點(diǎn)的差值與選取點(diǎn)的差值的絕對值，并判斷絕對值大小。是，則選取點(diǎn)的灰度值與下相鄰點(diǎn)的灰度平均值相同；否，則選取點(diǎn)的灰度值與上相鄰點(diǎn)的灰度平均值相同。
結(jié)合具體例子說明，選取點(diǎn)的相鄰點(diǎn)的預(yù)設(shè)數(shù)量為2個，具體如下
在灰度分布圖曲線上隨機(jī)的獲得選取點(diǎn)，具體的選取點(diǎn)為A，其2個上相鄰點(diǎn)為 A1和A2, 2個下相鄰點(diǎn)為L1和A_2 ；所對應(yīng)的上相鄰點(diǎn)的灰度值為F1和F2, 2個下相鄰點(diǎn)的灰度值為F—和F_2，A的灰度值為F。上相鄰點(diǎn)對應(yīng)的灰度平均值為Fs = (FJF2)/2 ；下相鄰點(diǎn)對應(yīng)的灰度平均值為I7x = (F_i+F_2)/2。若F與Fs的差值的絕對值大于F與!^x的差值的絕對值，則選取點(diǎn)A所對應(yīng)的F變換為Fx，即與下相鄰點(diǎn)的!^x值相同。
劃分模塊40，用于通過對灰度分布圖曲線上點(diǎn)的斜率劃分得到過渡區(qū)及對照區(qū)。 具體地，獲取灰度分布圖曲線上點(diǎn)的斜率值，根據(jù)斜率值劃分區(qū)域數(shù)據(jù)分析區(qū)域，即過渡區(qū)及對照區(qū)，進(jìn)一步的從彈性圖像中的關(guān)注區(qū)域獲取豐富的圖像數(shù)據(jù)。
在一實(shí)施例中，劃分模塊40包括
斜率值獲取單元41，獲取灰度分布圖曲線上點(diǎn)的斜率值。具體地，在灰度分布曲線的點(diǎn)上獲取對應(yīng)的點(diǎn)的斜率。
比較單元42，用于將斜率值與預(yù)設(shè)的閾值范圍進(jìn)行比較，判斷斜率值是否超過閾值范圍，是，則斜率值對應(yīng)的灰度分布圖曲線的區(qū)域?yàn)閷φ諈^(qū)，否，則為過渡區(qū)。具體地，將斜率值與預(yù)設(shè)的閾值范圍進(jìn)行比較，在其它實(shí)施例中，該閾值范圍為O 0. 5，若超過該閾值范圍，則該斜率值對應(yīng)的灰度分布圖曲線與橫坐標(biāo)(距離)所包圍的區(qū)域?yàn)閷φ諈^(qū)；否， 則為過渡區(qū)。
通過本方案在彈性圖像中選取關(guān)注區(qū)域，并提取灰度分布圖，然后根據(jù)灰度分布圖曲線上點(diǎn)的斜率信息，獲得過渡區(qū)域及對照區(qū)域?？梢姡瑢椥詧D像的關(guān)注區(qū)域進(jìn)行圖像處理，得到了豐富的圖像數(shù)據(jù)，而不僅僅是簡單的提供彈性圖像。
在另一實(shí)施例中，結(jié)合附圖20，還包括
面積計(jì)算模塊，用于計(jì)算過渡區(qū)的面積和對照區(qū)的面積。具體地，采用微積分方法、填補(bǔ)發(fā)或估算法分別得到過渡區(qū)的面積和對照區(qū)的面積。
面積判斷模塊，用于過渡區(qū)的面積與對照區(qū)的面積的比值與預(yù)設(shè)的面積閾值進(jìn)行比較，判斷是否大于面積閾值，是，則關(guān)注區(qū)域?yàn)橹攸c(diǎn)關(guān)注區(qū)。具體地，根據(jù)渡區(qū)的面積與對照區(qū)的面積的比值，若超過面積閾值，例如為1/3，則確認(rèn)該關(guān)注區(qū)域?yàn)橹攸c(diǎn)關(guān)注區(qū)域，為進(jìn)一步獲取信息提供準(zhǔn)確的指引；若沒有超過面積閾值，則為非重點(diǎn)關(guān)注區(qū)。
以上所述實(shí)施例僅表達(dá)了本發(fā)明的幾種實(shí)施方式，其描述較為具體和詳細(xì)，但并不能因此而理解為對本發(fā)明專利范圍的限制。應(yīng)當(dāng)指出的是，對于本領(lǐng)域的普通技術(shù)人員來說，在不脫離本發(fā)明構(gòu)思的前提下，還可以做出若干變形和改進(jìn)，這些都屬于本發(fā)明的保護(hù)范圍。因此，本發(fā)明專利的保護(hù)范圍應(yīng)以所附權(quán)利要求為準(zhǔn)。
權(quán)利要求
1.一種基于彈性成像的圖像處理方法，其特征在于，包括以下步驟獲取彈性圖像；從所述彈性圖像中選取關(guān)注區(qū)域，并從所述關(guān)注區(qū)域提取灰度分布圖；對所述灰度分布圖的曲線進(jìn)行平滑處理；通過對所述灰度分布圖曲線上點(diǎn)的斜率劃分得到過渡區(qū)及對照區(qū)。
2.根據(jù)權(quán)利要求1所述的基于彈性成像的圖像處理方法，其特征在于，所述對所述灰度分布圖的曲線進(jìn)行平滑處理的步驟為按照所述灰度分布圖曲線上的選取點(diǎn)分別提取相鄰點(diǎn)，所述相鄰點(diǎn)包括預(yù)設(shè)數(shù)量的上相鄰點(diǎn)和下相鄰點(diǎn)；計(jì)算獲得上相鄰點(diǎn)和下相鄰點(diǎn)的灰度平均值，并分別與所述選取點(diǎn)的灰度值做差；判斷所述上相鄰點(diǎn)與所述選取點(diǎn)的差值的絕對值是否大于所述下相鄰點(diǎn)與所述選取點(diǎn)的差值的絕對值，是，則所述選取點(diǎn)的灰度值與下相鄰點(diǎn)的灰度平均值相同。
3.根據(jù)權(quán)利要求2所述的基于彈性成像的圖像處理方法，其特征在于，所述通過對所述灰度分布圖曲線上點(diǎn)的斜率劃分得到過渡區(qū)及對照區(qū)的步驟為獲取所述灰度分布圖曲線上點(diǎn)的斜率值；將所述斜率值與預(yù)設(shè)的閾值范圍進(jìn)行比較，判斷所述斜率值是否超過所述閾值范圍， 是，則所述斜率值對應(yīng)的所述灰度分布圖曲線的區(qū)域?yàn)閷φ諈^(qū)，否，則為過渡區(qū)。
4.根據(jù)權(quán)利要求3所述的基于彈性成像的圖像處理方法，其特征在于，所述閾值范圍為 0. 5。
5.根據(jù)權(quán)利要求3所述的基于彈性成像的圖像處理方法，其特征在于，所述通過對所述灰度分布圖曲線上點(diǎn)的斜率劃分得到過渡區(qū)及對照區(qū)的步驟之后還包括計(jì)算所述過渡區(qū)的面積和所述對照區(qū)的面積；所述過渡區(qū)的面積與所述對照區(qū)的面積的比值與預(yù)設(shè)的面積閾值進(jìn)行比較，判斷是否大于所述面積閾值，是，則所述過渡區(qū)為重點(diǎn)關(guān)注區(qū)。
6.一種基于彈性成像的圖像處理系統(tǒng)，其特征在于，包括彈性圖像模塊，用于加載彈性圖像；提取模塊，用于從所述彈性圖像中選取關(guān)注區(qū)域，并從所述關(guān)注區(qū)域提取灰度分布圖；平滑處理模塊，用于對所述灰度分布圖的曲線進(jìn)行平滑處理；劃分模塊，用于通過對所述灰度分布圖曲線上點(diǎn)的斜率劃分得到過渡區(qū)及對照區(qū)。
7.根據(jù)權(quán)利要求6所述的基于彈性成像的圖像處理系統(tǒng)，其特征在于，所屬平滑處理模塊包括選取單元，用于按照所述灰度分布圖曲線上的選取點(diǎn)分別提取相鄰點(diǎn)，所述相鄰點(diǎn)包括預(yù)設(shè)數(shù)量的上相鄰點(diǎn)和下相鄰點(diǎn)；計(jì)算單元，用于計(jì)算獲得上相鄰點(diǎn)和下相鄰點(diǎn)的灰度平均值，并分別與所述選取點(diǎn)的灰度值做差；判斷單元，用于判斷所述上相鄰點(diǎn)與所述選取點(diǎn)的差值的絕對值是否大于所述下相鄰點(diǎn)與所述選取點(diǎn)的差值的絕對值，是，則所述選取點(diǎn)的灰度值與下相鄰點(diǎn)的灰度平均值相同。
8.根據(jù)權(quán)利要求7所述的基于彈性成像的圖像處理系統(tǒng)，其特征在于，所述劃分模塊包括斜率值獲取單元，用于獲取所述灰度分布圖曲線上點(diǎn)的斜率值； 比較單元，用于將所述斜率值與預(yù)設(shè)的閾值范圍進(jìn)行比較，判斷所述斜率值是否超過所述閾值范圍，是，則所述斜率值對應(yīng)的所述灰度分布圖曲線的區(qū)域?yàn)閷φ諈^(qū)，否，則為過渡區(qū)。
9.根據(jù)權(quán)利要求8所述的基于彈性成像的圖像處理系統(tǒng)，其特征在于，所述閾值范圍為 0. 5。
10.根據(jù)權(quán)利要求7所述的基于彈性成像的圖像處理系統(tǒng)，其特征在于，還包括 面積計(jì)算模塊，用于計(jì)算所述過渡區(qū)的面積和所述對照區(qū)的面積；面積判斷模塊，用于所述過渡區(qū)的面積與所述對照區(qū)的面積的比值與預(yù)設(shè)的面積閾值進(jìn)行比較，判斷是否大于所述面積閾值，是，則所述過渡區(qū)為重點(diǎn)關(guān)注區(qū)。
全文摘要
一種基于彈性成像的圖像處理方法及系統(tǒng)，包括以下步驟獲取彈性圖像；從所述彈性圖像中選取關(guān)注區(qū)域，并從所述關(guān)注區(qū)域提取灰度分布圖；對所述灰度分布圖的曲線進(jìn)行平滑處理；通過對所述灰度分布圖曲線上點(diǎn)的斜率劃分得到過渡區(qū)及辨別區(qū)。通過對彈性圖像的關(guān)注區(qū)域進(jìn)行過渡區(qū)提取，得到圖像數(shù)據(jù)的定量信息，可以從繁多復(fù)雜的信息中有目的或有針對性的獲取特征。同時，也可以根據(jù)豐富的數(shù)據(jù)進(jìn)行更加精確的數(shù)據(jù)分析。
文檔編號G06T7/00GK102521823SQ20111035038
公開日2012年6月27日 申請日期2011年11月8日 優(yōu)先權(quán)日2011年11月8日
發(fā)明者劉新, 張麗娟, 鄭海榮, 鐘耀祖, 黃文慧 申請人:中國科學(xué)院深圳先進(jìn)技術(shù)研究院