專利名稱：：醫學成像系統的制作方法
技術領域：
：本發明涉及一種醫學成像系統以及一種相應的方法。本發明尤其可應用于回波心動描記成像領域。
背景技術：
：在JournalofAmericanCollegeofCardiology1996Vol28的658—664頁中公開了一種醫學成像方法，其能夠計算表示器官的組織的良好健康狀態的心肌機能指標MPI，其也稱為TEI。在所公開的示例中，該器官是心臟。從血液流入以及流出左心室的速度的計時來計算MPI指標。心臟的組織處于良好健康狀態意味著，如果它們能夠有效地將血液射出左心室，則它們就是正確收縮的。MPI指標越小，則射血就越好。上述成像系統的一個缺點在于，MPI指標給出了對心臟組織的健康狀態的全局表示。如果具有較低的射血分數，并無法知道它來自組織的哪個個區域，并且不知道組織的哪個區域受到損傷以及哪個區域沒有受到損傷。
發明內容本發明的目的是提出一種系統，該系統能夠計算局部心肌機能指標，以便確定器官的哪個組織受到損傷或沒有受到損傷。為此，該系統包括用于控制以下操作的控制器-采集器官的圖像序列；-在所述圖像序列上確定組織區域；-確定表征所述組織區域的運動的參數；-基于所述參數，確定表征心動周期的階段組；-基于所述階段組，計算所述組織區域的局部心肌機能指標。因此，由于確定了局部心肌機能指標，就能夠定義在該器官中的組織區域是否有病癥。-根據第一實施例，所述參數是速度。其優點在于，在測量中實現了良好的時間分辨率。-根據第二實施例，所述參數是位移。其能夠獲得組織區域沿著兩個軸的運動的參數，這是更為準確的。-根據第三實施例，所述參數是速度與應變(strain)的組合。這就能夠更為準確地確定心動周期的一些階段。-根據一個非限定性實施例，所述控制器還用于控制以下操作在所述圖像序列上確定另一組織區域；以及確定相關的階段組。這能夠從參考組織區域自動確定其他組織區域。此外，這能夠確定其他組織區域的局部心肌機能指標，并且能夠具有以所述局部心肌機能指標為參數表示的所述心臟的整體顯示。-根據一個非限定性實施例，對另一組織區域的確定使用沿著所述器官的線。這能夠以更快的方式確定組織區域。-根據一個非限定性實施例，對另一組織區域的確定是前向和后向執行的。這能夠以更為準確的方式確定其他組織區域。-根據一個非限定性實施例，對其他階段組的確定是利用動態時間規整(dynamictimewarping)執行的。這能夠通過與參考階段組進行相關來自動確定其他階段組。本發明還涉及一種用于醫學成像的方法，該方法包括以下步驟-采集器官的圖像序列；-在所述圖像序列上確定組織區域；-確定表征所述組織區域的運動的參數；-基于所述參數，確定表征心動周期的階段組；-基于所述階段組，計算所述組織區域的局部心肌機能指標。本發明最后涉及一種包含用于執行上述方法的程序指令的計算機程序產品。參考以下所述的各個實施例，本發明的這些方面和其他方面將會變得清晰并且被闡明。現在將參考附圖，借助于多個非限定性示例來更為詳細地描述本發明，在附圖中-圖1是諸如心臟的器官的示意圖，使用根據本發明的系統從中采集圖像序列；-圖2是根據本發明的與探頭協同操作的系統的示意圖；-圖3是心臟的速度圖像序列的示意-圖4示出了從速度圖像序列，例如圖3中所示的速度圖像序列，所獲得的典型速度曲線；-圖5是速度圖像序列到位移圖像序列以及到應變速率(strainrate)圖像序列的變換的示意-圖6示出了從速度圖像序列，例如圖3中所示的速度圖像序列，所獲得的異常速度曲線；-圖7是在圖3的圖像序列中確定另一組織區域的一個非限定性實施例；-圖8是由用于將圖5的速度曲線與圖4的速度曲線相關的相關算法得到的矩陣；-圖9是示出了基于圖8的相關，在圖5的速度曲線上搜索一些階段的示意-圖IO是示出了與用圖7的確定操作所找到的一些組織區域相關的一些階段的示意-圖11顯示了根據本發明的用于醫學成像的方法的一個非限定性實施例的示意圖。具體實施例方式可以使用圖2中所示的系統SYS來采集任何器官(諸如心臟)的圖像。在以下描述中將會考慮心臟HRT的示例。要回憶的是，心臟HRT由左右心室LV和RV、主動脈AO、以及左右心房LA和RA組成，如圖1所示，動脈血從左心室LV流出到主動脈AO，同時右心室RV將從右心房RA接收的靜脈血輸出到肺動脈。主動脈AO可以由主動脈瓣AV閉合，左心房可以由二尖瓣MV閉合。圖2中描述了系統SYS。其與換能器陣列TAR及其相關電子電路協同操作，一起形成了探頭PRB。系統SYS包括用于控制以下操作的控制器CTRL:-采集器官HRT的圖像序列SQ(為此，控制器CTRL配置探頭PRB);-在所述圖像序列上確定組織區域PT;-確定表征所述組織區域PT的運動的參數P;-基于所述參數P，確定表征心動周期的階段組PH;-基于相關的階段組PH，計算所述組織區域PT的局部心肌機能指標LMPIo可任選地，系統SYS還包括心電圖觸發器ECG一T、用于顯示所采集的圖像序列的屏幕SCR，諸如LCD屏幕、以及用戶接口M一USER。要注意的是，控制器CTRL包括微處理器，能夠借助于指令對該微處理器進行預編程，或者能夠由系統SYS的用戶對該微處理器進行編程，例如經由接口M—USER。要注意的是，組織區域PT包含圖像序列中的至少一個點(例如像素或體素)或多個點(例如像素或體素)。在一個非限定性實施例中，控制器CTRL還用于控制以下操作-在所述圖像序列SQ上確定另一組織區域PT1;以及-確定表征所述另一組織區域PT1的心動周期的階段組PH。以下將詳細描述系統SYS的操作。1)采集圖像序列SQ。為了采集心臟HRT的圖像序列SQ，將探頭PRB應用到患者身體上，在一個實施例中是應用到心臟附近的頂點上。在第一實施例中，采集速度圖像序列SQVE。可以使用本領域技術人員公知的稱為TDI處理的彩色組織多普勒成像方法，來獲得組織區域的位移速度。結果得到了當心臟的組織區域移動時，例如朝向換能器陣列TAR移動時，心臟的組織區域的彩色編碼形式的速度圖像序列。通常，紅色與組織區域的收縮(朝向換能器的運動)相關，藍色與舒張(遠離換能器的運動)相關。因此，圖像越紅，則組織收縮并朝向換能器陣列TAR運動的越多。與之對照，圖像越藍，則組織舒張并遠離換能器陣列TAR運動的越多。該TDI處理具有的優點在于，在對組織區域的速度進行測量時實現了沿著換能器陣列TAR的軸的良好的時間分辨率。在第二實施例中，采集二維或三維的灰度級圖像序列SQGR。該第二實施例的優點在于-避免了如TDI處理一樣所必需的那種專用的協議過程，因為灰度級序列已經是臨床過程的一部分，不再需要額外的采集。-與僅在換能器方向上的運動可使用的TDI處理相比，沿著兩個軸的像素運動是可使用的，包括縱向軸和徑向軸。因此，使用該第二實施例，灰度級圖像序列的像素運動更為精確。為了使用戶能夠在這兩個實施例之間進行選擇，用戶接口M—USER包括用于在這兩個實施例之間進行選擇的模塊。2)確定組織區域PT。在第一實施例中，系統SYS的用戶能夠通過對如圖3中所示的圖像序列SQ的視覺評價來選擇組織區域PT。為此，用戶在該序列SQ上，在心壁內的與心臟的組織區域PT0相對應的某個位置上點擊光標CURS。在第二實施例中，能夠基于以下將會描述的相關參數自動進行選擇，并且用戶能夠選擇使該自動選擇有效或無效。為了使用戶能夠在這兩個實施例之間進行選擇，用戶接口M一USER包括用于在這兩個實施例之間進行選擇的模塊。3)確定表征所述組織區域PT運動的參數P。在第一實施例中，參數P是速度參數。在第二實施例中，參數P是位移參數。在第三實施例中，參數P是應變(也稱為變形)和速度參數的組合。在第四實施例中，參數P是應變速率(也稱為速度梯度)和速度參數的組合。根據該參數之一，定義了表示所選組織區域PT的運動的曲線C，并且可以將曲線C顯示在系統SYS的屏幕SCR上。對于速度曲線CVE，X軸表示時間，Y軸表示以厘米/秒計的組織區域PT的運動速度。這種曲線的一個示例在圖4中給出。在0下方，區域的運動是遠離傳感器陣列TAR運動，而在0上方，區域的運動是朝向傳感器陣列TAR運動。對于位移曲線CDI，X軸表示時間，Y軸表示以厘米計的區域PT的位移。其示出了從先前圖像到當前圖像的像素位移。對于應變和速度參數的組合，定義了速度曲線和應變曲線。對于應變/變形曲線CST，X軸表示時間，Y軸表示在所選區域與另一區域之間的以百分比計的變形。對于應變速率和速度參數的組合，定義了速度曲線和應變曲線。對于應變速率/速度梯度曲線CVG，X軸表示時間，Y軸被定義為秒—1。其表示在兩個區域之間的壓縮或膨脹的速度。為了使用戶能夠在這4個實施例之間進行選擇，用戶接口M—USER包括用于在這4個實施例之間進行選擇的模塊。要注意的是，這些曲線可以從步驟l)中采集的圖像序列直接地或非直接地推導出來，或者是從速度圖像序列SQVE或者是從灰度級圖像序列SQGR，一個圖像序列是與形成心臟的區域組PT相關的一組曲線C。速度曲線CVE是直接從速度圖像序列SQVE推導而來的。例如，使用本領域技術人員公知的斑點跟蹤(speckletracing)或紋理跟蹤，從灰度級圖像序列SQGR推導出位移曲線CDI。在另一個示例中，曲線C可以通過對一條曲線或多條曲線執行微分或積分(空間上或時間上的)來從另一曲線推導出。因此，位移曲線CDI可以通過時間積分而從速度曲線CVE推導出，速度曲線CVE可以通過時間微分而從位移曲線CDI推導出。應變/變形曲線CST可以通過空間微分而從兩條位移曲線CDI推導出。應變速率/速度梯度曲線CVG可以通過空間微分而從兩條速度曲線CVE推導出。并且最終，應變/變形曲線CST可以通過時間積分而從應變速率曲線CVG推導出。圖5示出了該變換。4)基于所述參數P，確定表征心動周期的階段組PH。在一個非限定性實施例中，表征一個心動周期CC的多個階段如圖4中所示。-在時刻to和tl之間的等容收縮階段IVC，-時刻tl和t2之間的射血階段EJC，-在時刻t2和t3之間的等容舒張階段IVR，以及-在時刻t3和t4之間的舒張階段RLX。在這4個階段期間，主動脈瓣AV和二尖瓣MV按如下所述地開啟或閉合。<table>tableseeoriginaldocumentpage10</column></row><table>因此，在等容收縮階段ivc期間，心臟在相同的血容積內收縮。左心室LV的內部壓力增加。當該內部壓力高于主動脈AO的外部壓力時，主動脈瓣AV開啟，導致血液從左心室LV射出到主動脈AO。在健康的心臟中，射出血液容積的2/3。這是射血階段EJC。在射血的同時，內部壓力降低。當內部壓力等于外部壓力時，主動脈瓣AV閉合。心臟膨脹。這是等容舒張階段IVR。當內部壓力低于左心房LA的外部壓力時，二尖瓣開啟，血液從左心房LA流到左心室LV。這是舒張階段RLX。最后，當左心室LV的內部壓力等于主動脈AO的外部壓力時，二尖瓣MV閉合。左心室LV充滿血液。兩個瓣AV和MV都閉合。等容階段的持續時間是肌肉強度的重要指標，因為較長的時間就暗示了肌肉不夠強壯，以致于不能快速地增加LV腔內的壓力來達到用于射血的主動脈AO中的壓力或者用于舒張的左心房LA中的壓力。要注意的是，該階段組是在一個心動周期CC過程中定義的。因此，為了獲得完整心動周期的速度曲線，可以使用例如患者的心電圖ECG，其示出了心臟的電子活動，更具體而言是收縮的發起。回憶到心臟收縮階段是心臟HRT進行收縮從而導致血液射入動脈的階段(心臟收縮階段包括等容收縮IVC階段和射血EJC階段)，心臟舒張階段是心臟HRT舒張的階段(心臟舒張階段包括等容舒張IVR階段和舒張RLX階段)。因此，圖像序列SQ的采集例如在所述心電圖ECG上進行同步。為了實現該同步，系統SYS必須使用ECG觸發器ECG—T。這種階段PH是從與所選的當前組織區域PT相關的參數P推導出的，就是說從相關的曲線C推導出的。在第一實施例中，通過對與所選的組織區域PT相關的曲線C的視覺評價來定義階段組PH。系統SYS的用戶通過自己的經驗來確定這些階段以及對應的時間t。在第二實施例中，使用算法ALG在與所選的組織區域相關的曲線C上定義階段組PH。在第一變化例中，其是基于示出了典型階段PHS的典型曲線CS的。典型曲線CS是示出了在來自健康心臟的一個心動周期CC期間的階段組PHS。圖4示出了典型的速度曲線CVE，圖6示出了有病的心臟的異常速度曲線CVE。然后通過這兩條曲線CS與C之間的相似度來從典型階段PHS推導出階段組PH。在第二個變化例中，算法ALG是基于所研究的曲線C的零交叉點和/或者最大值、最小值時刻以及諸如所述曲線C的斜率的最大值和最小值之類的參數。如上所述，階段組PH或者可以從速度曲線CVE推導出，或者可以從位移曲線CDI推導出，或者可以從速度曲線CVE和應變曲線CST的組合推導出、或者可以從速度曲線CVE和應變速率曲線CVG的組合推導出。因此，從速度曲線CVE可以直接推導出這4個階段IVC、EJC、IVR和RLX。從位移曲線CDI，可以直接在CDI上推導出這4個階段，或者通過對所述曲線CDI執行時間微分從而得到速度曲線CVE。從應變曲線CST，可以推導出兩個階段，等容收縮的開始(t0)和等容舒張的開始(t2)。然后，可以從速度曲線CVE推導出另兩個階段EJC和RLX。希望使用該組合，因為在應變曲線CST中更容易檢測兩個階段IVC和IVR，因此能夠通過在例如速度曲線CVE中簡單地查找峰值點而不是零交叉點來檢測它們。實際上，它們表示在一個心動周期過程中應變的最大值和最小值，其中，該最大值和最小值分別對應于心肌舒張的開始和心肌收縮的開始。這種應變峰值點的檢測在文檔WO2004/092766Al中以非限定性示例給出。從應變速率曲線CVG可以推導出等容收縮IVC的結束(tl)。然后從速度曲線CVE推導出其他階段。5)基于相關的階段組PH，計算所述組織區域PT的局部心肌機能指標LMPI。LMPI指標被估計為等容收縮時間與等容舒張時間之和除以射血時間。LMPI指標=(IVC+IVS)/EJC由于己經按照以上所述確定了所選組織區域PT的階段組PH，因此就能夠計算局部LMPI指標。根據所計算的局部LMPI指標的值，能夠推導出在所關注的組織區域PT中是否有任何病癥，諸如例如急性心肌梗塞。得到局部LMPI指標的好處在于顯示出在心臟中是否存在受損的組織區域PT。因此，希望的是確定其他組織區域PT的LMPI指標，這是如以下所述地執行的。6)在所述圖像序列上確定其他組織區域PT并且確定相關的階段組PH。6.1)確定另一區域PT在第一實施例中，使用沿著心肌層的線LX來確定另一組織區域，如圖7所示。在第一變化例中，在該線LX上以預定的間隔確定預定數量的組織區域，例如從參考組織區域PTR1開始，并且其在圖7上用點來表示。在第二變化例中，使用了本領域技術人員公知的前向傳播，并且基于至少一個參考組織區域。在圖7所示的示例中，已經在二尖瓣MV的兩個邊界上選擇了兩個參考區域PTR1和PTR2。該前向傳播沿著經過心肌層壁W的線LX從第一參考區域PTR1行進到第二參考區域PTR2。對于每個新區域PT，新的參考區域是沿著該線LX的上一個區域。區域PT1的參考區域是區域PTR1，區域PT2的參考區域是區域PT1，依此類推。在這兩個第一變化例的一個非限定性實施例中，選擇組織區域PT以使其在其長度方向上垂直于線LX，并且在所述區域PT內僅考慮一個像素PX來計算與所述區域PT相關的階段組PH。這一區域在圖7中示出為PT3。因此，唯一要計算的階段組將會代表所述區域PT的所有像素PX，因為這一區域的一個像素的運動類似于所述區域的其他像素的運動。因此，計算將會更為簡單和更為快速。在這兩個變化例的另一個非限定性實施例中，可以前向或后向地執行所述確定，如圖7中的箭頭所示。這就允許以下所述的對階段組的確定更為準確。在第二實施例中，可以使用前向傳播而不使用線LX。要注意的是，當定義了另一組織區域PT1、PT2等等時，按照以上步驟3)所述，自動確定與相關曲線C1、C2等等相對應的參數P。6.2)當選擇了另一組織區域PT1時，按照如下所述定義相關的階段組。在一個非限定性實施例中，基于本領域技術人員公知的動態編程，使用稱為DTW的動態時間規整來執行所述確定。DTW允許在兩條曲線之間執行相關。在第一步驟中，在例如通過前向傳播找到的當前組織區域PT1的曲線Cl與參考區域PTR的參考曲線CR之間執行相關。為了執行該相關，使用時間-時間矩陣來展示這兩條曲線CR和Cl之間的時間對準，如圖8所示，在圖8中參考曲線CR在側邊上向上繪出，曲線Cl在底部繪出。針對心動周期CC繪制了該時間矩陣。在該非限定性實施例中，已經使用了兩條速度曲線，并且使用了代表所述速度曲線的向量VR和VI。這些向量表征了參考區域PTR和當前區域PT1。在一個非限定性變化例中，這些向量由形成所述曲線的一些值組成。用相關路徑(correlationpath)PTH來表示所述時間對準。如果這兩條曲線CR和C1是相同的，則所述時間對準會由這樣的相關路徑PTHR來表示即該相關路徑是所述矩陣中斜率為45度的對角線，如圖8中所示。如果當前曲線C1與參考曲線CR相同，但是在時間上擴張，則所述時間對準會由這樣的相關路徑PTHR來表示即該相關路徑是所述矩陣中斜率小于45度的對角線(未示出)。在圖8的示例中，兩條曲線CR和C1并不相同。已經計算了相關路徑PTH1并且將其繪出以展示所述時間對準。在第二步驟中，在所述相關路徑PTH1的幫助下確定當前曲線C1的階段，如圖9所示。與參考曲線CR的第一階段等容收縮IVC相對應的間隔被投影到該相關路徑PTH1上，得到了間隔IVC1。對其他3個階段射血、等容收縮以及收縮執行相同的操作。因此，為與當前組織區域PT1相對應當前曲線C1找到了階段IVC1、EJC1、IVR1和RLX1。在所述DTW處理的第一變化例中，當使用沿著心肌層的線LX進行前向傳播時，可以如下所述地校正所計算的當前區域PT的階段組。圖10示出了為多個區域PT1到PTn-l以及用于前向傳播的開始的兩個參考區域(其是參考區域PTR1和PTR2)計算的階段組。這些階段組表示在水平軸上，而區域表示在垂直軸上。由于彼此接近的兩個區域的運動相似，因此它們的階段組應該差別不大。這就意味著它們各自的值應該彼此接近。因此，在前向傳播過程中不同階段組的進展應該相似。如果不是這種情況，就意味著在DTW的計算中可能存在錯誤，并且由此在對當前區域PT的階段定義中可能存在錯誤，并且可以再次執行該計算。在圖10中，對于不同的相鄰組織區域PT的每個階段IVC、EJC、IVC和RLX的發展分別用線L1、L2、L3、L4來表示。第一條線Ll以及最后兩條線L3和L4是均勻發展的。第二條線L2示出區域PT2的突變B。存在計算的錯誤。因此對該區域PT2再次計算DTW。在DTW的第二變化例中，可以計算一個區域PT的兩個階段組，并將平均值保留作為良好的階段組PH。這可以用于如上所述地以前向和后向行進的方式對區域的確定。第一和第二變化例可以合并在一起。在該DTW步驟的結束，對其他組織區域PT再一次執行步驟6和7，直到達到停止標準為止。在一個非限定性示例中，所述停止標準可以是一個組織區Z,該組織區Z包含已經由用戶劃定的多個組織區域PT。7)基于所述其他組織區域PT的相關階段組PH，計算所述其他組織區域PT的局部LMPI指標。因此，由于已經在圖像序列SQ上確定了不同組織區域PT的階段組PH，與這些區域相關的所有局部LMPI指標都是可用的，并且可以按照以上步驟5)中所述地計算得到。8)顯示以局部LMPI指標為參數表示的彩色圖像。系統SYS的用戶可以容易地看到可能受損或沒有受損的組織區域。例如，紅色可以與較大的LMPI值相關，藍色與較小的LMPI值相關，或者可能會有與LMPI值的平均值相關的彩色顯示。總之，圖ll示出了根據本發明的醫學成像方法，其中，能夠看到系統SYS所控制的不同操作。當然，一些操作可以并行執行。例如，步驟1和3可以并行執行。因此，已經描述的本發明的成像系統具有以下優點.--其考慮了心臟收縮功能和心臟舒張功能，而與心率和血壓無關。-心臟疾病具有多種病因，包括冠狀動脈疾病、初期的心肌和瓣膜心臟疾病，并且在心肌層面上的疾病可能是由于心臟收縮和心臟舒張機能不良的組合、單獨的心臟收縮機能不良或者單獨的心臟舒張機能不良引起的。因此，本系統比僅能夠測量左心室的心臟收縮功能的系統相比有著更大的優勢。-本系統使得能夠將患有臨床心臟病的病人與沒有心臟病但具有等效的心室機能不良的病人區分開。要注意的是，上述實施例是說明性的而不是限定本發明，本領域技術人員在不脫離附帶的權利要求所定義的本發明的范圍的情況下能夠設計很多可替換實施例。例如，能夠使用3D數據以3D方式計算LMPI。可以不僅為左心室LV計算LMPI，而且還為其他心臟腔，諸如右心室RV計算LMPI。還要注意，所述實施例并不限于任何成像形式。本發明可以應用于能夠采集速度信息或解剖圖像序列的任何系統。在一個非限定性實施例中，本系統SYS能夠應用于超聲成像。在此情況中，探頭PRB是超聲探頭，所采集的圖像序列SQ是超聲圖像序列。在權利要求中，括號中的任何參考標記都不應解釋為限制權利要求。在任何權利要求中以及整個說明書中，詞語"包含"和"包括"以及類似詞語并不排除有所列出的元件和步驟之外的其他元件和步驟的存在。元件的單數并不排除引用多個這種元件的情況，反之亦然。本發明可以借助于包含幾個元件的硬件來實現，并且可以借助于適當編程的計算機來實現。在枚舉出幾個模塊的裝置中，這些模塊之中的幾個模塊可以用同一個硬件來實現。在彼此不同的從屬權利要求中引用了特定手段這一事實并不表示不能使用這些手段的組合來獲得有利性。權利要求1、一種系統，包括用于控制以下操作的控制器(CTRL)-采集器官(HRT)的圖像序列(SQ)；-在所述圖像序列上確定組織區域(PT)；-確定表征所述組織區域(PT)的運動的參數(P)；-基于所述參數(P)，確定表征心動周期的階段組(PH)；-基于所述階段組(PH)，計算所述組織區域(PTR)的局部心肌機能指標(LMPI)。2、如權利要求l所述的系統，其中，所述參數是速度。3、如權利要求l所述的系統，其中，所述參數是位移。4、如權利要求l所述的系統，其中，所述參數是速度和應變的組合。5、如權利要求1所述的系統，其中，所述控制器還用于控制以下操作-在所述圖像序列(SQ)上確定另一組織區域(PT1);以及-確定相關的階段組(PH)。6、如權利要求5所述的系統，其中，所述確定另一組織區域的操作使用了沿著所述器官的線(LX)。7、如權利要求5所述的系統，其中，所述確定另一組織區域的操作是前向和后向執行的。8、如權利要求5所述的系統，其中，所述確定所述相關的階段組(PH)的操作是使用動態時間規整(DTW)執行的。9、一種用于醫學成像的方法，包括以下步驟-采集器官(HRT)的圖像序列(SQ);-在所述圖像序列上確定組織區域(PT);-確定表征所述組織區域(PT)的運動的參數(P);-基于所述參數(P)，確定表征心動周期的階段組(PH);-基于所述階段組(PH)，計算所述組織區域(PT)的局部心肌機能指標(LMPI)。10、如權利要求9所述的用于醫學成像的方法，還包括步驟通過視覺評價來選擇所述組織區域(PT)。11、一種計算機程序產品，其包含程序指令，所述程序指令用于當由處理器執行所述程序時執行如權利要求9所述的方法。全文摘要本系統涉及一種醫學成像系統。首先，采集器官(HRT)的圖像序列(SQ)。然后，在所述圖像序列上定義組織區域(PT)，并且定義表征所述組織區域(PT)的運動的參數(P)。然后，基于所述參數(P)，定義表征心動周期的階段組(PH)。最后，基于所述階段組(PH)，計算所述組織區域(PTR)的局部心肌機能指標(LMPI)。文檔編號G06T7/00GK101523435SQ200780036925公開日2009年9月2日申請日期2007年10月1日優先權日2006年10月4日發明者E·德尼,E·薩魯克斯,O·博納富,O·熱拉爾,P·阿蘭申請人:皇家飛利浦電子股份有限公司