版權(quán)通告本專利文檔的公開內(nèi)容的一部分包括受版權(quán)保護的材料。版權(quán)所有者不反對任何人照專利和商標辦公室專利文件或記錄原樣復(fù)制本專利文件或?qū)@_內(nèi)容,但除此之外版權(quán)所有者保留所有相關(guān)的版權(quán)。
背景技術(shù):
:1.
技術(shù)領(lǐng)域:
:本發(fā)明涉及身體生物電信號的檢測、測量或記錄。更具體地,本發(fā)明涉及生成與心律失常相關(guān)的電解剖標示圖。2.相關(guān)領(lǐng)域描述:表1-首字母縮略詞和縮寫cog重心cv傳導(dǎo)速度ecg心電圖egm電描記圖ic-egm心內(nèi)電描記圖lat局部激活時間mri磁共振成像內(nèi)部器官的三維圖像可用于許多基于導(dǎo)管的診斷和治療應(yīng)用中,并且實時成像被廣泛用于外科規(guī)程中。現(xiàn)在,常常使用包括用于標測心臟電活動的電生理傳感器的心臟導(dǎo)管來執(zhí)行對心臟中電勢的標測。通常,根據(jù)心臟內(nèi)的位置感測并記錄心內(nèi)膜中的時變電勢,然后將其用于標測局部電描記圖或局部激活時間。由于通過心肌傳導(dǎo)電脈沖所需的時間,心內(nèi)膜中各點的激活時間不同。該電傳導(dǎo)在心臟中任何點處的方向常規(guī)地由激活矢量(在本文中也被稱為傳導(dǎo)速度矢量)表示,該激活矢量垂直于等電激活波前,這兩者均可衍生自激活時間的標示圖。激活波前通過心內(nèi)膜中任意點的傳播速率可表示為傳導(dǎo)速度矢量。可通過觀察現(xiàn)象,諸如多個激活波前、激活矢量的反常集聚、或速度矢量的變化或矢量與正常值的偏差來識別心臟的激活信號的傳導(dǎo)中的局部缺陷。此類缺陷的示例包括內(nèi)曲區(qū)域,該內(nèi)曲區(qū)域可與稱為復(fù)雜碎裂電描記圖的信號圖案相關(guān)聯(lián)。一旦通過此類標測來定位缺陷,就可對其進行消融(如果其功能異常)或以其它方式治療以便盡可能恢復(fù)心臟的正常功能。文獻characterizationofleftventricularactivationinpatientswithheartfailureandleftbundle-branchblock,(auricchioetal.,circulation.2004;109:1133-1139)描述了心力衰竭患者的左心室激活順序以及左束支阻滯qrs形態(tài),其中在固有節(jié)律與非同步起搏期間同時應(yīng)用了三維接觸標測和三維非接觸標測。由于阻滯線位于前面、側(cè)面或下面,所以“u形”激活波前出現(xiàn)在大多數(shù)患者體內(nèi)。阻滯線的功能行為表現(xiàn)在:在不同部位且以不同周期長度進行非同步心室起搏期間,阻滯線的位置發(fā)生改變。技術(shù)實現(xiàn)要素:本發(fā)明的實施方案提供了用于通過限定區(qū)域內(nèi)的心內(nèi)電描記圖的同時測量值進行心房纖顫電解剖標測的方法和系統(tǒng)。激活時間基于單極電描記圖的完整負斜率而不是該斜率上的單個基準點,例如,最大斜率(最大-dv/dt)、中間振幅或時間。相反,確定表示局部激活的斜率中每個斜率的lat范圍。lat范圍由分別表示斜率的開始和結(jié)束的峰點和谷點劃分界限,限定了其中激活波經(jīng)過的時間窗口。本發(fā)明的另選實施方案還提供了斜率前或斜率后的另選激活周期的識別,從而提供短雙電位和長雙電位。根據(jù)本發(fā)明的實施方案,提供了一種方法,該方法通過以下步驟來執(zhí)行:將多電極探頭插入活體受檢者的心臟中;并發(fā)地記錄來自位于心臟中相應(yīng)位置處的電極的電描記圖;界定電描記圖中相應(yīng)的激活時間間隔;由激活時間間隔生成電傳播波標示圖;通過調(diào)整電描記圖的激活時間間隔內(nèi)的局部激活時間使波的相干性最大;以及報告所調(diào)整的局部激活時間。該方法的一個方面包括消融心臟中的組織以修改這些波。根據(jù)該方法的另外方面,激活時間間隔包括第一間隔和第二間隔,第一間隔由表示斜率的開始與結(jié)束的峰點和谷點界定,第二間隔包括關(guān)于第一間隔的窗口,其中調(diào)整局部激活時間在第二間隔內(nèi)執(zhí)行。根據(jù)該方法的另一方面,關(guān)于第一間隔的窗口為±40ms。根據(jù)該方法的一個方面,電極的相應(yīng)位置處的局部激活時間表示為從0至1變化的模糊電極隸屬函數(shù)(μe)。根據(jù)該方法的一個方面,生成電傳播波標示圖包括將電描記圖分段成在相應(yīng)時間的一系列幀,其中幀是電極讀數(shù)到值矩陣的相應(yīng)指配。該方法的另一方面包括:從值矩陣中為電極指配局部激活時間;將心臟的至少一部分建模成網(wǎng)格,其中網(wǎng)格中頂點的一部分對應(yīng)于電極的相應(yīng)位置;基于所指配的局部激活時間確定來自幀的頂點之間電傳播的傳導(dǎo)速度;以及使用傳導(dǎo)速度計算波的相干性。根據(jù)該方法的另外方面,使波的相干性最大包括:將傳導(dǎo)速度表示成傳導(dǎo)速度矢量;計算來自傳導(dǎo)速度矢量的相應(yīng)頂點處的速度偏差矢量;以及使速度偏差矢量的長度最小。該方法的又一方面包括:在網(wǎng)格上內(nèi)插與電極的相應(yīng)位置不對應(yīng)的頂點的激活時間。根據(jù)該方法的另外方面,網(wǎng)格的頂點處的激活時間表示為在0和1之間變化的模糊頂點隸屬函數(shù)(μu),每個頂點處的頂點隸屬函數(shù)包括其相鄰頂點的頂點隸屬函數(shù)的加權(quán)組合。根據(jù)該方法的另外方面,使波的相干性最大包括調(diào)整網(wǎng)格的頂點處的激活時間,以使網(wǎng)格中的頂點隸屬函數(shù)的平均值最大。該方法的又一方面包括:通過外推法將頂點隸屬函數(shù)指配至網(wǎng)格的非相鄰頂點。根據(jù)本發(fā)明的實施方案,還提供了一種設(shè)備,該設(shè)備包括多電極探頭和處理器,該多電極探頭適合于插入活體受檢者的心臟中,該處理器構(gòu)造成接收來自電極的電信號并執(zhí)行以下步驟:并發(fā)地記錄來自心臟中相應(yīng)位置處的電極的電描記圖;界定電描記圖中相應(yīng)的激活時間間隔;由激活時間間隔生成電傳播波標示圖;通過調(diào)整電描記圖的激活時間間隔內(nèi)的局部激活時間使波的相干性最大;以及報告所調(diào)整的局部激活時間。該設(shè)備可包括消融功率發(fā)生器,該消融功率發(fā)生器連接到探頭,用于消融心臟中的組織以修改波。附圖說明為更好地理解本發(fā)明,就本發(fā)明的詳細說明以舉例的方式做出參考,該詳細說明應(yīng)結(jié)合以下附圖來閱讀,其中類似的元件用類似的附圖標號來表示,并且其中:圖1為根據(jù)本發(fā)明的實施方案的用于評估活體受檢者的心臟中的電活動的系統(tǒng)的圖解示意圖;圖2是根據(jù)本發(fā)明的實施方案制備的心房纖顫的一組單極lat;圖3是根據(jù)本發(fā)明的實施方案的單極電極的lat優(yōu)化方法的流程圖;圖4是根據(jù)本發(fā)明的實施方案的示出了電極柵格的圖,該電極柵格示出了心房纖顫時的阻滯線;圖5是根據(jù)本發(fā)明的實施方案的波形幀序列;圖6是根據(jù)本發(fā)明的實施方案的示出了幀分段的流程圖;圖7是根據(jù)本發(fā)明的實施方案的示例性幀分段標示圖;圖8是根據(jù)本發(fā)明的實施方案的示出了幀分段矩陣和電解剖標示圖的復(fù)合圖;圖9是根據(jù)本發(fā)明的實施方案的示出了區(qū)生長的圖;圖10是根據(jù)本發(fā)明的實施方案生成的分段標示圖的一部分;圖11是根據(jù)本發(fā)明的實施方案的示出了激活波的表示框圖;圖12是根據(jù)本發(fā)明的實施方案的示出了網(wǎng)格上的電波的速度矢量計算的圖;圖13是根據(jù)本發(fā)明的實施方案的示出了三角形網(wǎng)格的頂點處的速度矢量的計算的圖;圖14是根據(jù)本發(fā)明的實施方案的示出了可用于波標測的某些變量的計算的圖;圖15是根據(jù)本發(fā)明的實施方案的用于標測心臟中相干傳播波的方法的流程圖;圖16是根據(jù)本發(fā)明的另選實施方案的用于標測心臟中相干傳播波的方法的流程圖;圖17是根據(jù)本發(fā)明的實施方案的可用于相對于圖16所述的方法中的數(shù)據(jù)的激活標示圖;圖18是根據(jù)本發(fā)明的實施方案的用于使電極網(wǎng)格的lat與在更精確網(wǎng)格上的心臟解剖相關(guān)的流程圖;圖19是示出了根據(jù)本發(fā)明的實施方案的圖18的方法中各個階段的一組圖;圖20是示出了根據(jù)本發(fā)明的實施方案的來自心內(nèi)電描記圖的模糊lat的起因的復(fù)合圖;圖21是根據(jù)本發(fā)明的實施方案的使模糊lat與心內(nèi)電描記圖的形態(tài)相關(guān)的復(fù)合圖;圖22是根據(jù)本發(fā)明的實施方案的距心內(nèi)膜不同距離處的電極的函數(shù)μe的系列圖;圖23是根據(jù)本發(fā)明的實施方案的從多電極導(dǎo)管獲取的函數(shù)μe的系列圖;圖24是根據(jù)本發(fā)明的實施方案的示出了將電極lat標測到網(wǎng)格的最頂點上的圖;圖25是根據(jù)本發(fā)明的實施方案的示出了網(wǎng)格中鄰域隸屬函數(shù)的發(fā)展的圖;圖26是根據(jù)本發(fā)明的實施方案的示出了網(wǎng)格中鄰域重疊的圖;圖27示出了根據(jù)本發(fā)明的實施方案的網(wǎng)格鄰域并呈現(xiàn)了兩幅曲線圖;圖28是根據(jù)本發(fā)明的實施方案的示出了模糊lat計算的復(fù)合圖;圖29是根據(jù)本發(fā)明的實施方案的具有外推頂點的網(wǎng)格;圖30是根據(jù)本發(fā)明的實施方案的模糊激活優(yōu)化方法的流程圖;圖31示出根據(jù)本發(fā)明的實施方案的優(yōu)化前和優(yōu)化后的網(wǎng)格的頂點中的函數(shù)μu的兩幅曲線圖;圖32是示出了根據(jù)本發(fā)明的實施方案的包括相鄰頂點的區(qū)域內(nèi)中心處的均勻傳導(dǎo)性優(yōu)化的圖;并且圖33是根據(jù)本發(fā)明的實施方案的示出了包括相鄰頂點的區(qū)域內(nèi)中心處的均勻傳導(dǎo)性優(yōu)化的圖。具體實施方式在下文的具體實施方式中,示出了許多具體細節(jié),以便提供對本發(fā)明的各種原理的全面理解。然而,對于本領(lǐng)域的技術(shù)人員而言將顯而易見的是,并非所有這些細節(jié)都是實施本發(fā)明所必需的。在這種情況下,未詳細示出熟知的電路、控制邏輯、以及用于常規(guī)算法和過程的計算機程序指令的細節(jié),以免不必要地使一般概念模糊不清。以引用方式并入本文的文獻將被視作本申請的整體部分,不同的是,就任何術(shù)語在這些并入文件中以與本說明書中明確或隱含地作出的定義矛盾的方式定義而言,應(yīng)僅考慮本說明書中的定義。定義“注釋”或“注釋點”是指被認為表示感興趣事件的在電描記圖上的點或候選點。在本公開中,事件通常是電極所感測的電波的傳播的局部激活時間。電描記圖中的“活動”在本文中用于表示電描記圖信號中的不同突發(fā)性區(qū)或波動變化。此類區(qū)可被識別為在基線信號的區(qū)之間是突出的。在本公開中,“活動”更多地是指穿過心臟的一個或多個電傳播波在電描記圖上的表現(xiàn)?!安ā笔侵感呐K中的連續(xù)電傳播。心臟中電傳播波的“相干性”是指波在不同空間點或時間點上相位恒定以及頻率相等的度量值?!白铚木€”是指對心臟中的電傳播的阻礙或阻滯。此類線可對波劃分界限。波本身可包含稱為“波內(nèi)阻滯”的阻滯線。電描記圖的“一次斜率”是與在電極下通過的激活波的局部激活時間相關(guān)的斜率?!岸涡甭省笔桥c未在電極下通過、即來自遠端激活波諸如遠場活動的波相關(guān)的斜率。當一個斜率和另一個斜率在限定的時間窗口內(nèi)一致發(fā)生時,該斜率被“聯(lián)接”到另一個斜率。“阻滯點”是具有小于使用者限定的值(通常為0.2m/s)的傳導(dǎo)速度的點。除此之外或另選地,阻滯點是位于兩個電極之間的點,其中離開第一電極的激活波到達第二電極處,結(jié)果發(fā)現(xiàn)第二電極先前已在使用者限定的時間間隔例如100ms內(nèi)激活,該時間間隔就在到達之前并且在第二電極的不應(yīng)期開始之后?!白铚木€”或“阻滯線”是阻滯點的集合。“幀”是網(wǎng)格電極讀數(shù)中各個并發(fā)讀數(shù)至值矩陣的指配。網(wǎng)格表面上的電極或頂點的“相鄰者”是指以該電極或頂點為中心的3×3柵格中的表面上的其它電極或頂點。綜述現(xiàn)在轉(zhuǎn)到附圖,首先參見圖1,其為根據(jù)本發(fā)明的公開實施方案構(gòu)造和操作的用于在活體受檢者的心臟12上執(zhí)行消融規(guī)程的系統(tǒng)10的圖解示意圖。該系統(tǒng)包括導(dǎo)管14,由操作者16將導(dǎo)管14經(jīng)由皮膚穿過患者的血管系統(tǒng)插入心臟12的心室或血管結(jié)構(gòu)中。操作者16,通常為醫(yī)師,將導(dǎo)管的遠側(cè)末端18例如在消融靶部位處與心臟壁接觸??筛鶕?jù)公開于美國專利6,226,542和6,301,496中和公開于共同轉(zhuǎn)讓的美國專利6,892,091中的方法來制備電激活標示圖,這些專利的公開內(nèi)容以引用方式并入本文。系統(tǒng)10可包括用合適的軟件編程以用于實施下文所述功能的通用或嵌入式計算機處理器。因此,盡管本文中的其它附圖所示的系統(tǒng)10的部分示出為包括多個單獨的功能框,但這些單獨的功能框未必為單獨的物理實體,而是可代表例如在可由處理器訪問的存儲器中存儲的不同計算任務(wù)或數(shù)據(jù)對象。這些任務(wù)可在運行于單個處理器上或運行于多個處理器上的軟件中實施。該軟件可在有形非暫態(tài)介質(zhì)諸如cd-rom或非易失性存儲器上提供給所述一個處理器或多個處理器。另選地或除此之外,系統(tǒng)10可包括數(shù)字信號處理器或硬連線邏輯。一種采用系統(tǒng)10的元件的商品可以3系統(tǒng)購自biosensewebster,inc.(3333diamondcanyonroad,diamondbar,ca91765)。該系統(tǒng)可由本領(lǐng)域的技術(shù)人員進行修改以體現(xiàn)本文所述的本發(fā)明的原理??梢酝ㄟ^施加熱能對例如通過電激活標示圖評估確定為異常的區(qū)域進行消融,例如,通過將射頻電流通過導(dǎo)管中的線傳導(dǎo)至遠側(cè)末端18處的一個或多個電極,這些電極將射頻能量施加到心肌。能量在組織中被吸收,從而將組織加熱到一定溫度(通常為約50℃),在該溫度下組織永久性地失去其電興奮性。在規(guī)程成功后,此規(guī)程在心臟組織中形成非傳導(dǎo)性消融灶,該非傳導(dǎo)性消融灶中斷導(dǎo)致心律失常的異常電通路。本發(fā)明的原理可應(yīng)用于不同的心臟腔室,以診斷并治療多種不同的心律失常。導(dǎo)管14通常包括柄部20,柄部20具有在該柄部上的合適的控制器,以使操作者16能夠根據(jù)消融手術(shù)的需要對導(dǎo)管的遠側(cè)端部進行導(dǎo)向、定位和取向。為了幫助操作者16,導(dǎo)管14的遠側(cè)部分容納向位于控制臺24中的處理器22提供信號的位置傳感器(未示出)。處理器22可以履行如下所述的若干處理功能。導(dǎo)管14為多電極導(dǎo)管,該導(dǎo)管可以是如球囊37的右部分中所示的籃形導(dǎo)管或如左部分中所示的樣條導(dǎo)管。在存在多個電極32的任何情況下,這些電極用作感測電極,并在籃形或樣條電極上具有已知位置,并且彼此關(guān)系已知。因此,一旦導(dǎo)管例如通過構(gòu)建當前位置標示圖定位于心臟中,則心臟中每個電極32的位置為已知的。一種用于生成當前位置標示圖的方法描述于授予bartal等人的共同轉(zhuǎn)讓美國專利8,478,383中,該專利以引用方式并入本文??梢允闺娦盘柦?jīng)由纜線34從位于導(dǎo)管14的遠側(cè)末端18處或附近的電極32,在心臟12和控制臺24之間來回傳送??梢酝ㄟ^纜線34和電極32將起搏信號和其它控制信號從控制臺24傳送至心臟12。線連接件35將控制臺24與體表電極30和用于測量導(dǎo)管14的位置和取向坐標的定位子系統(tǒng)的其它部件連接在一起。處理器22或另一個處理器(未示出)可以是定位子系統(tǒng)的元件。電極32和體表電極30可用于按照以引用方式并入本文的授予govari等人的美國專利7,536,218中所教導(dǎo)的在消融部位處測量組織阻抗。溫度傳感器(未示出),通常為熱電偶或熱敏電阻器,可安裝在導(dǎo)管14的遠側(cè)末端18附近。控制臺24通常包括一個或多個消融功率發(fā)生器25。導(dǎo)管14可適于利用任何已知的消融技術(shù)將消融能量(例如,射頻能量、超聲能量和激光產(chǎn)生的光能)傳導(dǎo)至心臟。共同轉(zhuǎn)讓的美國專利6,814,733、6,997,924和7,156,816中公開了此類方法,這些專利以引用方式并入本文。在一個實施方案中,定位子系統(tǒng)包括磁定位跟蹤構(gòu)造,該磁定位跟蹤構(gòu)造利用生成磁場的線圈28,通過以預(yù)定工作空間生成磁場并感測導(dǎo)管處的這些磁場來確定導(dǎo)管14的位置和取向。合適的定位子系統(tǒng)在以引用方式并入本文的美國專利7,756,576以及上述美國專利7,536,218中有所描述。如上所述,導(dǎo)管14聯(lián)接到控制臺24,這使得操作者16能夠觀察并調(diào)控導(dǎo)管14的功能??刂婆_24包括處理器,優(yōu)選為具有適當信號處理電路的計算機。處理器被耦合以驅(qū)動監(jiān)視器29。信號處理電路通常接收、放大、過濾并數(shù)字化來自導(dǎo)管14的信號,這些信號包括由上述傳感器和位于導(dǎo)管14遠側(cè)的多個位置感測電極(未示出)生成的信號??刂婆_24和定位系統(tǒng)接收并使用數(shù)字化信號,以計算導(dǎo)管14的位置和取向,并如下文進一步詳細所述分析來自電極的電信號。通常,系統(tǒng)10包括為簡明起見而未示出于附圖中的其它元件。例如,系統(tǒng)10可包括心電圖(ecg)監(jiān)視器,其被聯(lián)接以從一個或多個體表電極接收信號,以便為控制臺24提供ecg同步信號。如上文提及,系統(tǒng)10通常還包括基準位置傳感器,其位于附接于受檢者身體外部的外加基準補片上,或者位于插入心臟12中并相對于心臟12保持在固定位置的內(nèi)置導(dǎo)管上。系統(tǒng)10可從外部成像模態(tài)諸如mri單元等接收圖像數(shù)據(jù),并且系統(tǒng)10包括圖像處理器,該圖像處理器可結(jié)合在處理器22中或由處理器22調(diào)用以用于生成并顯示如下文描述的圖像。心房纖顫的標測:心房纖顫的特征在于不具有周期性或重復(fù)圖案的復(fù)雜傳播圖案。可存在多條阻滯線,從而將各種形式的分離波分開。將心房激活時間標測至心房電極網(wǎng)格的嘗試導(dǎo)致測量誤差?;趤碜詷藴y導(dǎo)管的電極讀數(shù)的空間分辨率不足以用于評估復(fù)雜的心房纖顫激活圖案。生成重構(gòu)心房纖顫時所標測區(qū)域的活動的可信電解剖標示圖需要采取若干步驟。這些步驟通常包括:對所獲取的電描記圖進行預(yù)處理,檢測局部激活時間,然后將與標測陣列上各個電極下的組織的活動有關(guān)的lat組合成可解釋性標示圖或影片。一個波標測方法描述于標題為lineofblockdetection的共同轉(zhuǎn)讓的美國專利公布2016/0045123中,該專利公布以引用方式并入本文。其中所描述的規(guī)程對幀,即幀分段情況下的心房波進行檢測與標測。這些規(guī)程對于標測阻滯線所描繪的波特別有用。根據(jù)本發(fā)明的實施方案,激活時間基于單極電描記圖的完整負斜率而不是該斜率上的單個基準點,例如,最大斜率(最大-dv/dt)、中間振幅或時間。相反,確定表示局部激活的斜率中每個斜率的lat范圍。該lat范圍由表示斜率的開始和結(jié)束的峰點和谷點劃分界限,限定了其中激活波經(jīng)過的時間窗口?,F(xiàn)在參見圖2,該圖是根據(jù)本發(fā)明的實施方案制備的心房纖顫的一組單極lat。根據(jù)本發(fā)明的實施方案,這些值源自于由具有64(8×8)個電極的心內(nèi)膜標測導(dǎo)管所產(chǎn)生的讀數(shù)。在此描述中,每個正方形通過右上角處的位置編號39(從1至64變化)進行識別。正方形41示出粗體的單極lat,該lat通常從最大負斜率(dv/dt)獲得并且是可信的。這些正方形是固定的,并且不能改變。正方形43示出落入lat時間窗口內(nèi)的普通字體lat,該時間窗口具有上限和下限,上限和下限分別在正方形的上方和下方示出。這些限值源自于對斜率劃分界限的峰點和谷點的時間關(guān)聯(lián)性,lat在該斜率上進行限定(最大-dv/dt)。這些限值在lat之前與之后,并且在執(zhí)行下述算法前進行確定。在這些限值內(nèi),可通過下述算法更改正方形45的值。具有下劃線字體數(shù)字的正方形45不可信。沒有正方形的柵格上的位置也表示讀數(shù)不可靠。在一些情況下,例如,正方形47,lat時間窗口的另選上限和下限均可用,并且分別在正方形右側(cè)與左側(cè)處示出。另選限值通?;诙屉p電位或長雙電位,即,另選限值源自于之前和之后的與在其上限定lat的斜率(最大-dv/dt)有關(guān)的斜率。這些限值在執(zhí)行算法前進行確定。本文所公開的本發(fā)明原理適用于具有其它數(shù)目的電極的導(dǎo)管。lat的優(yōu)化:如上所述,一些電極的lat可指配在給定時間窗口的邊界內(nèi)。另選地,可以從心內(nèi)電描記圖中的另選斜率指配lat。用于指配的算法被稱作波傳播相干性函數(shù),該函數(shù)對于每個電極而言使傳播矢量的最小速度偏差角與該電極的lat相關(guān)?,F(xiàn)在參見圖3,該圖是根據(jù)本發(fā)明的實施方案的單極電極的lat優(yōu)化方法的高級流程圖。為了呈現(xiàn)清楚,在圖3和本文其它流程圖中以特定線性序列示出了過程步驟。然而,將顯而易見的是,這些步驟中的多個可并行地、異步地或以不同的順序執(zhí)行。本領(lǐng)域的技術(shù)人員還應(yīng)當理解,另選地,過程可被表示為多個相關(guān)的狀態(tài)或事件,例如,在狀態(tài)圖表中。此外,可能不需要所有示出的過程步驟來實施該方法。在初始步驟49中,心臟傳統(tǒng)上用多電極標測導(dǎo)管進行導(dǎo)管插入。導(dǎo)管諸如可得自biosensewebster的nav或?qū)Ч苓m用于初始步驟49中。將導(dǎo)管的電極放置成與心房之一中的相應(yīng)位置電流接觸。接著,在步驟51中,可以將電極及其位置的標示圖呈現(xiàn)為電極網(wǎng)格。該激活標示圖可顯示阻滯線。現(xiàn)在參見圖4,該圖是根據(jù)本發(fā)明的實施方案的示出了電極柵格的示例性圖,該電極柵格示出了心房纖顫時的阻滯線。在正方形53中指示局部激活時間。電極位置55,一些被正方形53遮蓋,形成柵格。低質(zhì)量的電極信號由星號57指示。這些信號對于lat確定是不可靠的。阻滯點59具有位于感興趣時間窗口之外的局部激活時間。阻滯線由阻滯點集合61指示。接著,在步驟63中,通過導(dǎo)管的多個電極,同時記錄心房電活動,每個電極具有相應(yīng)的位置,該位置可使用系統(tǒng)10(圖1)的位置追蹤能力確定。在可能的情況下,使用傳統(tǒng)方法,諸如最大-dv/dt偏轉(zhuǎn)法,最初檢測到心房去極化??墒褂霉餐D(zhuǎn)讓的專利申請序列號14/693,042中所教導(dǎo)的方法排除心室遠場活動,該專利申請以引用方式并入本文。將電描記圖分段成相應(yīng)時間下的一系列幀。接著,在步驟65中,將激活標示圖轉(zhuǎn)變?yōu)閭鞑ゲ耸緢D。這通過下述的幀分段過程來實現(xiàn)。該過程可涉及將心房建模成(通常但不是一定的)三維三角形網(wǎng)格。用于生成網(wǎng)格的方法包括樞轉(zhuǎn)球算法(ball-pivotingalgorithm)。另選地,網(wǎng)格可生成為狄羅尼(delaunay)三角形劃分,包括多個三角形。結(jié)合網(wǎng)格的幀分段實現(xiàn)了該轉(zhuǎn)變。下文在圖6的討論中對幀分段進行了描述?,F(xiàn)在參見圖5,該圖為根據(jù)本發(fā)明的實施方案的范圍為237-399ms的10個波幀的序列。如附圖左側(cè)處兩個幀代表性所示,區(qū)域67顯示沒有活動波。區(qū)域69表示由阻滯線73分開的波71。幀的左下角的數(shù)字表示以ms為單位的時間。回到圖3,在步驟75中,優(yōu)化電極的激活時間。這是迭代規(guī)程,該規(guī)程結(jié)合圖15進行描述??刂魄斑M至最終步驟77。對于其中獲得有效讀數(shù)的每個電極而言,報告固定lat或產(chǎn)生最小速度偏差角的lat。幀分段:現(xiàn)在參見圖6,該圖是示出了如步驟65所示(圖3)的幀分段的流程圖。該圖描述根據(jù)本發(fā)明的實施方案的幀分段。幀填充有大部分邏輯上彼此相關(guān)的lats。在這個過程內(nèi),計算傳導(dǎo)速度。參考后續(xù)幀中的對應(yīng)lat來確定傳導(dǎo)阻滯。在執(zhí)行算法的過程中保持源列表。源列表包含電極編號和相關(guān)聯(lián)的lat,以相對于相鄰lat對該lat進行檢驗以用于阻滯或傳導(dǎo)。被發(fā)現(xiàn)具有傳導(dǎo)性的電極編號被添加到源列表并且在下一輪算法中進行檢驗。以此方式,算法使屬于同一個波的電極編號的區(qū)生長。在框79中確定源幀。對這個框的輸入是幀結(jié)構(gòu)、距離矩陣和從源電極獲得的lat???9的輸出是源電極的幀編號。幀編號的指配基于源電極的lat處的幀空缺。對于所有空缺的幀來說,使用如表2中所示的matlab例程來計算以下特性,以便支持指配決定。表2基于每個空缺幀的特性,以列表1中的偽碼給定判定規(guī)則。列表1如果沒有可用的空缺幀充當minmap將源{ele/lat}指配到下一幀否則如果帶有一個或多個傳導(dǎo)性第一階相鄰激活的幀是可用的將源{ele/lat}指配到lat最接近的幀否則如果沒有可用的帶有傳導(dǎo)性相鄰激活的幀如果最近的lat<100ms(af循環(huán))將源{ele/lat}指配到lat最接近的幀否則將源{ele/lat}指配到下一幀結(jié)束結(jié)束。現(xiàn)在參見圖7,該圖是根據(jù)本發(fā)明的實施方案的由使用matlab例程的上述算法產(chǎn)生的示例性幀分段標示圖。該幀中的空白區(qū)域可歸因于缺失電極,或歸因于由于當在幀分段過程中檢測到波間阻滯時進行值的再指配而駐留在下一幀中的波?,F(xiàn)在參見圖8,該圖是根據(jù)本發(fā)明的實施方案的示出由上述算法和matlab例程產(chǎn)生的示例性幀分段矩陣81和電解剖標示圖的復(fù)合圖。阻滯線由點,例如矩陣81上的點83和標示圖87、89上的點85指示。區(qū)生長:通過電極網(wǎng)格進行的心房纖顫波檢測涉及區(qū)生長算法以及幀生成和分段算法?,F(xiàn)在參見圖9,該圖是根據(jù)本發(fā)明的實施方案的解釋區(qū)生長過程的圖。該過程是迭代的。例如通過識別碼來識別圍繞3×3柵格中心的相鄰電極是方便的。出于區(qū)生長的目的,使用3x3柵格中的中心電極和八個相鄰電極的lat來計算歸一化傳導(dǎo)速度。在框93中識別電極的3x3正方形柵格91,示出為由虛線描繪的正方形。接著,在框95中,在階段97處評估正方形柵格91中的傳導(dǎo)。這個過程需要:(1)計算正方形柵格91中的中心電極99和相鄰電極之間的3維距離;(2)確定中心電極99和相鄰電極之間的局部激活時間間隔。另外的信息可用于區(qū)的擴展:(1)lat時間窗口。這些時間窗口提供lat不準確度的指示。(2)3×3柵格內(nèi)的四個2×2正方形的傳導(dǎo)速度矢量。(3)相鄰ic-ecg的一次注釋和ff斜率(二次注釋)。(4)ic-ecg的質(zhì)量和lat質(zhì)量?,F(xiàn)在可基于cvnorm=d(lat)/d(loc)來確定傳導(dǎo)完整性或傳導(dǎo)阻滯,其中l(wèi)oc是指心臟內(nèi)電極的位置cvnorm≥cv。當cvnorm≤cvnorm_min時,指示阻滯,在這種情況下:cv≤cvnorm_min。另選的傳導(dǎo)檢測策略包括僅針對高質(zhì)量的ic-ecg和lat來確定傳導(dǎo)速度矢量的量值。這種方法具有對lat不準確度的敏感性問題。另一個另選的傳導(dǎo)檢測策略涉及使用標準方法在lat上擬合雙二次曲面的3x3適合性。這產(chǎn)生超定解決方案,但更能抵抗lat不準確度。關(guān)于幀分段和區(qū)生長的另外細節(jié)見于上述美國專利公布2016/0045123中。心房激活的標測:現(xiàn)在參見圖10,該圖是根據(jù)本發(fā)明的實施方案的使用上述規(guī)程而產(chǎn)生的常見幀分段標示圖的一部分。單個幀可包括多于一條波,并且給定波可擴展多個幀。在此類情況下,可通過影線圖案將波彼此區(qū)分開來。本發(fā)明的實施方案中采用內(nèi)插法,以在更精確的導(dǎo)管網(wǎng)格上產(chǎn)生激活波,如上所述。內(nèi)插法的優(yōu)點包括:當存在相對較大的信號振幅和負斜率時,更好的時空限定激活時間。當情況并非如此時,可顯示另選(之前或之后)的激活時間。這些另選激活時間通過lat附近的較大負斜率表征。由于對于顯示目的而言電極網(wǎng)格為紋路相對粗糙的,三維標測系統(tǒng),諸如上述carto3系統(tǒng),通過使用加權(quán)平均值具有內(nèi)插的取樣點,并且在偽著色標示圖中已經(jīng)顯示了結(jié)果,該加權(quán)平均值配置成與表面上的點之間的測地距離成反比。其它內(nèi)插法技術(shù)包括拉普拉斯內(nèi)插法。然而,應(yīng)當注意,傳統(tǒng)的表面拉普拉斯內(nèi)插法可不提供平滑的lat/波圖案,這是因為解剖學(xué)網(wǎng)格不規(guī)則,例如頂點之間的距離不同。這些問題可通過應(yīng)用名稱為highdefinitioncoloringofheartchambers的共同轉(zhuǎn)讓的專利申請序列號15/009,285的教導(dǎo)內(nèi)容來解決,該專利申請以引用方式并入本文。本發(fā)明的實施方案中所采用的內(nèi)插法技術(shù)有利于以下情況下的治療:其中信號缺少作為“模糊lat”的振幅突變。這些技術(shù)描述如下。現(xiàn)在參見圖11,該圖為根據(jù)本發(fā)明的實施方案的示出了激活波的表示的框圖。在框101中,同時發(fā)生如上所述的心內(nèi)電描記圖的獲取和幀分段???03以電極網(wǎng)格的分辨率生成激活波。框103接收來自電極網(wǎng)格的電描記圖,該電描記圖以信號105(elemesh)從框101中輸出。以下激活標示圖產(chǎn)生自框101并輸出為信號107:actmap每個電極最可能的lat,actfixmap如果斜率和振幅大于閾值,則lat固定,actinvalid沒有可用于該電極的有效lat。另外地,以下激活標示圖可用于其中l(wèi)at不固定的電極,但是可根據(jù)本文所公開的算法而變化:actearliermap另選的之前的lat,actlatermap另選的之后的lat,actminmaplat下限,actmaxmaplat上限。另外的波標示圖和傳導(dǎo)速度(cv)標示圖作為信號111從框101輸出???13將框103中產(chǎn)生的波標示圖轉(zhuǎn)變成分辨率更高的形式,通常使用更詳細的網(wǎng)格和上述規(guī)程。框113中執(zhí)行的內(nèi)插法規(guī)程產(chǎn)生相干的激活波,這能夠涉及,當明確定義的lat不可用時,改寫激活時間范圍、或之前或之后的激活時間。標示圖內(nèi)插法:如上所述,創(chuàng)建相干波涉及標示圖內(nèi)插法。假設(shè)該規(guī)程相對于三角形網(wǎng)格進行描述。可將其它網(wǎng)格構(gòu)造分解成三角形?,F(xiàn)在參見圖12,該圖是根據(jù)本發(fā)明的實施方案的示出了網(wǎng)格上的電波的速度矢量計算的圖。三角形115具有邊,包括邊117、119。每條邊都存在速度矢量。例如,邊117的速度矢量由以下公式給出:其中d12是三角形115的頂點121和頂點123之間的距離,并且lat1和lat2是頂點121、頂點123處的激活時間。用類似的方式計算邊119的速度矢量。整個三角形115的速度是沿著邊117、邊119的速度之和:現(xiàn)在參見圖13,該圖是示出了根據(jù)本發(fā)明的實施方案的三角形網(wǎng)格127的一部分的頂點125處的速度矢量的計算。頂點125由連接的三角形129、三角形131、三角形133形成,這些三角形具有相應(yīng)的根據(jù)圖12的公式計算的速度矢量。示出了三角形129、三角形131的速度矢量135、137。將連接到頂點125的所有三角形的速度矢量,包括速度矢量135、137,投影到平面139上,該平面由垂直于頂點125的法線141限定。例如,矢量143、145分別是速度矢量135、137在平面139上的投影。速度矢量147是連接到頂點125的三角形的所有投影速度矢量的平均值?,F(xiàn)在參見圖14,該圖是根據(jù)本發(fā)明的實施方案的示出了可用于波標測的某些變量的計算的圖。所投影的速度矢量,包括矢量149、矢量151及其差值矢量153,相對于頂點155定位。兩個所投影的速度矢量(例如,矢量149、矢量151)之間的角157由αk表示。頂點處的速度偏差矢量vd由下式給出:其中tr1和tr2是頂點處連接的三角形(未示出),并且是它們的速度矢量。根據(jù)本發(fā)明的實施方案,為了從心內(nèi)電描記圖提取傳播波,使速度偏差矢量最小。在另一應(yīng)用中,當需要限定非傳播波(即,阻滯波)時,可使速度偏差矢量最大。一般來講,當傳導(dǎo)速度低于最小值,通常為0.2m/s時,假定存在阻滯。繼續(xù)參見圖14,頂點155處的速度偏差角通過以下公式給出:其中n為形成共頂點155的連接三角形的數(shù)目。一般來講,當速度偏差角很小時,波為相干波。為了獲得傳播波,使速度偏差角最小。相干標測算法:回到圖11,假設(shè)信號107、信號111中所包含的標示圖可用,并且已經(jīng)將電極表示成三角形網(wǎng)格。以下規(guī)程是步驟75(圖3)的一個實施方案的詳細流程圖,這可在框103中執(zhí)行。現(xiàn)在參見圖15,該圖是根據(jù)本發(fā)明的實施方案的用于標測心臟中相干傳播波的方法的流程圖。在初始步驟159中,選擇從導(dǎo)管的電極中的一者獲得的信息,包括lat。接著,在決定步驟161中,確定從當前電極獲得的lat是否需要進行優(yōu)化。如果決定步驟161中的確定為否定的,即,讀數(shù)無效,或者lat通過明確定義的負斜率(最大–dv/dt)進行計算并且因此是固定的,則不需要采取進一步行動。控制直接前進至下述決定步驟167。如果決定步驟161中的確定為肯定的,則控制前進至步驟163,其中對應(yīng)于當前電極的網(wǎng)格上的頂點的速度偏差角最小。該優(yōu)化通過調(diào)整用于計算連接的三角形的速度矢量的lat(公式1)來實現(xiàn)。當前電極的調(diào)整限值從上述標示圖actminmap和actmaxmap獲得。接著,在決定步驟165中,確定在步驟163的當前迭代中計算的速度偏差角相對于前一次迭代中計算的值的改變是否已經(jīng)超過了預(yù)定公差。如果該確定為肯定的,則繼續(xù)進行該算法。控制返回到步驟163,以進行迭代循環(huán)。如果決定步驟165中的確定為否定的,則確定速度偏差角的值已經(jīng)集中,并且不需要繼續(xù)進行迭代步驟163。控制前進至決定步驟167,此時確定是否仍然需要處理來自更多電極的數(shù)據(jù)。如果確定為肯定的,則控制回到初始步驟159,并且選擇另一電極。如果決定步驟167中的確定為否定的,則控制前進至最終步驟169。報告優(yōu)化的lat。第一替代實施方案:在該實施方案中,lat斜率前和/或后的另選激活時間經(jīng)識別,提供短雙電位和長雙電位。這些另選lat可涉及最小速度偏差角的優(yōu)化?,F(xiàn)在參見圖16,該圖是根據(jù)本發(fā)明的另選實施方案的用于標測心臟中相干傳播波的方法的流程圖,在該實施方案中尋找替代激活時間。對每個電極重復(fù)該實施方案的規(guī)程。圖16中所示的規(guī)程可以用相對于圖15所述的方法進行協(xié)調(diào),例如,當包括步驟163和決定步驟165的迭代循環(huán)沒有集中或另外地沒有導(dǎo)致對特定電極的lat的可靠優(yōu)化時。當信號的振幅較低或質(zhì)量較低,或者受到無法完全消除的遠場電位的干擾時,發(fā)生此類情況。在決定步驟171中,確定特定電極的lat是否固定或讀數(shù)無效。如果決定步驟171中的確定為肯定的,則忽略算法的其余部分,并且控制直接前進至最終步驟177。在一些實施方案中,為了節(jié)省計算機資源,不自動生成與另選lat相關(guān)的功能標示圖(actearliermap,actlatermap),并且另選激活可能還不可用。如果決定步驟171中的確定為否定的,則在決定步驟173中確定當前電極的另選lat是否已經(jīng)存在。如果決定步驟173中的確定為肯定的,則控制前進至優(yōu)化步驟109,其中使用另選可能中對應(yīng)于最小速度偏差角的另選lat來執(zhí)行步驟163和決定步驟165中的上述迭代循環(huán)。當之前和之后的lat均可用時,針對它們中的兩者執(zhí)行優(yōu)化,并且從結(jié)果中選擇優(yōu)化的對應(yīng)于速度偏差角的更小值的lat。如果決定步驟173中的確定為否定的,則在步驟175中計算之前和之后的另選lat。這些lat根據(jù)二次斜率確定,二次斜率見于關(guān)于一次斜率的±40ms窗口內(nèi)。此后,基于如上所述所計算的另選lat執(zhí)行優(yōu)化步驟109。執(zhí)行優(yōu)化步驟109之后,在最終步驟177中,報告電極網(wǎng)格的相干lat?,F(xiàn)在參見圖17,該圖是根據(jù)本發(fā)明的實施方案的可用于相對于圖16所述的方法中的數(shù)據(jù)的示例性激活標示圖。電極網(wǎng)格的激活時間在附圖的左上角示出。呈現(xiàn)電極的電壓描記線,其中指出了常見的一次斜率和二次斜率??煽啃匀狈χ甘酒饕约安缓细裰甘拘畔⒎謩e作為圍繞描記線的橢圓和正方形示出。精確網(wǎng)格和解剖學(xué)標示圖的內(nèi)插:上述規(guī)程在相對粗糙的電極網(wǎng)格上產(chǎn)生了相干lat?,F(xiàn)在參見圖18,該圖是根據(jù)本發(fā)明的實施方案的用于使電極網(wǎng)格的lat與更精確三角形網(wǎng)格上表示的心臟解剖相關(guān)的流程圖。在初始步驟179中,使用上述lat優(yōu)化規(guī)程產(chǎn)生電極網(wǎng)格的相干lat。接著,在步驟181中,為了產(chǎn)生具有更高分辨率的標示圖,對于每個2×2電極正方形,使用雙線性內(nèi)插算法來內(nèi)插電極網(wǎng)格上的相干lat,這考慮了波間阻滯與波內(nèi)阻滯,這可使用上述美國專利公布2016/0045123的教導(dǎo)內(nèi)容進行識別。該步驟導(dǎo)致粒度相對細的內(nèi)插網(wǎng)格。然后在最終步驟183中,將內(nèi)插網(wǎng)格標示到心臟解剖學(xué)標志中。將電極的位置標示到解剖學(xué)模型中,例如,解剖學(xué)網(wǎng)格,僅僅當電極位置與解剖學(xué)網(wǎng)格點之間的距離低于20mm時如此標示。當前的位置感測設(shè)備,諸如上述carto系統(tǒng),能確定1mm-2mm內(nèi)的導(dǎo)管電極的三維位置。然而,該標測可導(dǎo)致解剖學(xué)網(wǎng)格點(平均)過采樣或欠采樣。當導(dǎo)管電極網(wǎng)格點的數(shù)目超過相近的解剖學(xué)網(wǎng)格點的數(shù)目(間距<20mm)時,發(fā)生過采樣。如果相近的解剖學(xué)網(wǎng)格點的數(shù)目超過導(dǎo)管電極網(wǎng)格點的數(shù)目,發(fā)生欠采樣。一旦所有靠近解剖學(xué)網(wǎng)格點的導(dǎo)管電極網(wǎng)格點已經(jīng)被標測,則將其余的導(dǎo)管電極網(wǎng)格點標記為無效。應(yīng)當優(yōu)選地進行未分配的解剖學(xué)網(wǎng)格點(<20mm)的內(nèi)插,因為解剖學(xué)網(wǎng)格相對于電極網(wǎng)格(n=64)欠采樣。導(dǎo)管能使心臟變形,使得如果離線而非治療期間制備解剖學(xué)模型,解剖學(xué)模型將不準確。失真隨導(dǎo)管模型而變化。建議基于心內(nèi)電描記圖的斜率分析反復(fù)識別以及定位靠近心內(nèi)膜表面或與心內(nèi)膜表面接觸的導(dǎo)管電極,并且根據(jù)治療期間的需要,通過調(diào)整解剖學(xué)網(wǎng)格來對失真進行補償。這樣,在整個治療周期,能準確地使傳播波與心臟解剖結(jié)構(gòu)相關(guān)聯(lián)?,F(xiàn)在參見圖19,該圖是根據(jù)本發(fā)明的實施方案的示出了圖18的方法中各個階段的一組圖。附圖右側(cè)處的標號185所指的圖表示編號后導(dǎo)管電極網(wǎng)格點187(為了清楚起見,省略了一些點標識符)的變平滑后的標示圖。該標示圖示出了標號185所指的圖的下部分中的傳播波189。標號191所指的圖示出了適配于心房的解剖學(xué)模型的導(dǎo)管電極。導(dǎo)管電極的分布作為不規(guī)則線性標記193示于模型上。中心區(qū)域195表示標測導(dǎo)管的包殼,該包殼可以與心內(nèi)膜表面接觸或不可與心內(nèi)膜表面接觸。標號197所指的圖為標示圖,該標示圖表示標測到解剖學(xué)模型的導(dǎo)管電極網(wǎng)格點的位置,該模型如上所述為粒度較細的三角形網(wǎng)格。示出了若干傳播波199。第二另選實施方案:前述實施方案中所述的范圍內(nèi)的lat調(diào)整是指定激活時間的更廣泛技術(shù)的特殊情況,這在本文稱為“模糊激活時間”或“模糊lat”。模糊lat是指根據(jù)模糊邏輯、模糊數(shù)學(xué)規(guī)則和模糊集的一般理論的模糊數(shù)值,這些理論最初由loftizadeh以及隨后由bartkosko等人描述。在此概念下,為變量指配真實值,該值可以是任何介于0和1之間的實數(shù)。相比之下,在布爾邏輯中,變量的真實值不是0就是1。在該實施方案中,將真實值指配給以下lat:當偏離于所檢測到的lat時下降的lat。此外,可使用模糊數(shù)學(xué)規(guī)則對模糊lat進行組合。使用模糊lat而產(chǎn)生的相干標示圖被稱作“模糊相干標示圖”。該概念已經(jīng)縮小成一系列l(wèi)at、概率,這些lat、概率包括用以創(chuàng)建相干標示圖的之前和之后的另選lat。用于電極網(wǎng)格的注釋:在標測陣列中,ne是陣列中的電極數(shù)目:late是電極e的lat(ms),1≤e≤ne,μe是用于電極e處的激活的模糊隸屬函數(shù)。對于網(wǎng)格表面,μe=1。模糊隸屬函數(shù)限定模糊lat。與最大–dv/dt相關(guān)聯(lián)的lat被指配了值1。三角形隸屬函數(shù)從峰點(之前)和谷點(之后)處逐漸變至零,峰點和谷點與負斜率以及l(fā)at相關(guān)聯(lián)。隸屬函數(shù)下的區(qū)域歸一化為1。因此,相對陡峭的負斜率提供明顯的高振幅隸屬函數(shù),該函數(shù)具有與對負斜率劃分界限的峰點和谷點之間的時間相等的窄基底。相比之下,較緩和的負斜率提供基底較寬的低振幅隸屬函數(shù)。對于一組標測在解剖學(xué)網(wǎng)格上的局部激活時間:nu是網(wǎng)格中的頂點數(shù)目。latu是網(wǎng)格頂點u(ms)處的lat,1≤u≤nu。μu是用于該頂點處的激活的模糊隸屬函數(shù),并且相當于網(wǎng)格表面上的函數(shù)μe,μu=1。然而,模糊隸屬函數(shù)μu是加權(quán)(去缺口)的三角形隸屬函數(shù)的復(fù)合。相鄰頂點的隸屬函數(shù)根據(jù)如圖25所示的在最頂點與相鄰頂點之間傳播的波的距離和傳導(dǎo)速度進行去缺口。為了獲得網(wǎng)格頂點的清晰lat,需要例如使用下述重心(cog)方法進行“去模糊化”?,F(xiàn)在參見圖20,該圖是根據(jù)本發(fā)明的實施方案的示出了來自心內(nèi)電描記圖的模糊lat的起因的復(fù)合圖。標號201所指的電描記圖具有峰值時間203(pe)、谷值時間205和中間時間207。標號209所指的隸屬函數(shù)μe的圖為三角形,該三角形分別在時間late處具有最大值,并在峰值時間pe和谷值時間ve處具有零值?,F(xiàn)在參見圖21,該圖是根據(jù)本發(fā)明實施方案的使模糊lat與心內(nèi)電描記圖的形態(tài)相關(guān)的復(fù)合圖。心內(nèi)電描記圖211的負斜率提供了感興趣的時間窗口的參考。在心內(nèi)電描記圖213中,較寬的傳播纖顫波生成了3ms的斜率,該斜率由豎直虛線界定,分別對應(yīng)于峰值時間pe和谷值時間ve。在圖215中,最大隸屬值max(μ)應(yīng)當為1。曲線217下的面積然后變成1.5??梢允购瘮?shù)μe的曲線圖的形態(tài)與電極和心內(nèi)膜之間的接觸狀態(tài)相關(guān)。現(xiàn)在參見圖22,該圖是根據(jù)本發(fā)明的實施方案的距心內(nèi)膜不同距離處的電極的函數(shù)μe的系列圖。在標號219所指的圖中,電極與心內(nèi)膜接觸。在標號221所指的圖、標號223所指的圖中,電極位于距心內(nèi)膜越來越遠處。曲線越來越寬,并且函數(shù)μe的最大值隨著與心內(nèi)膜距離的增加而減小?,F(xiàn)在參見圖23,該圖是根據(jù)本發(fā)明的實施方案的從多電極導(dǎo)管得到的函數(shù)μe的系列曲線圖。具有尖峰的三角形,例如三角形225,對應(yīng)于與心內(nèi)膜接觸的電極。其它三角形不接觸。一般來講,三角形的基底寬度越大,該函數(shù)的最大值越小,電極與心內(nèi)膜之間的距離越大。用于模糊相干標測目的的網(wǎng)格可以是導(dǎo)管電極網(wǎng)格或?qū)π呐K解剖結(jié)構(gòu)進行建模的網(wǎng)格。針對本公開的目的,“最頂點”是指網(wǎng)格中最接近地對應(yīng)于導(dǎo)管電極中一個導(dǎo)管電極的位置的頂點。模糊相干標測的注釋:3d空間中電極的位置:pe電極位置。spe標測陣列上的電極位置(ne)集。網(wǎng)格定義:tui網(wǎng)格三角形,1≤i≤ntu,ntu網(wǎng)格三角形的數(shù)目。vuj網(wǎng)格頂點,1≤j≤nvu,nvu網(wǎng)格頂點的數(shù)目。mutu和vu集,對于網(wǎng)格頂點:pv節(jié)點最頂點節(jié)點;nb節(jié)點鄰域節(jié)點ex節(jié)點外推節(jié)點(既不是最頂點節(jié)點,也不是鄰域節(jié)點);每個頂點均屬于以下集中的一個集:集a無lat(ex節(jié)點),集b正好一個模糊lat(pv或nb節(jié)點),集c多于一個模糊lat(pv或nb節(jié)點)。計算電極與解剖學(xué)網(wǎng)格頂點之間的距離矩陣:deu距離矩陣(nu×ne)。標測到最頂點上:給定deu和μe,將每個電極激活late的貢獻標測到網(wǎng)格上具有最大真實值μupv的最近側(cè)頂點(最頂點;pv)上,如圖24所示。late→latupv;μe→μupv使用隸屬函數(shù)限定最頂點周圍的近鄰者。相鄰latμu通過latupv周圍的量dunb,pv/cv整體確定。然而,還不知道激活的方向。模糊隸屬函數(shù)是寬度為2×dunb,pv/cv整體的平頂。該運算示于圖25中。鄰域可重疊,在網(wǎng)格頂點上創(chuàng)建多于一個模糊lat。現(xiàn)在參見圖26,該圖為根據(jù)本發(fā)明的實施方案的示出了網(wǎng)格中鄰域重疊的圖。附圖左側(cè)處是三個最頂點的組227,這些最頂點的鄰域通過相應(yīng)的圓圈來表示,在所述圓圈中以最頂點為中心。附圖右側(cè)處是組229,其中最頂點的相鄰頂點已經(jīng)如上所述進行定義。一些相鄰頂點,例如頂點231、頂點233,處于這些鄰域中的多于一個鄰域中?,F(xiàn)在參見圖27,該圖示出了根據(jù)本發(fā)明的實施方案的網(wǎng)格鄰域并呈現(xiàn)了兩幅圖。該附圖示出,網(wǎng)格頂點上多個模糊lat可具有不同的真實值。傳播波通過箭頭235指示。模糊lat237、239、241處于三個獨立鄰域的組243中,并且具有如標號251所指的圖上所示的相應(yīng)的真實值245、247、249。標號253所指的圖示出了隸屬函數(shù),該函數(shù)是加權(quán)、去缺口三角形隸屬函數(shù)的復(fù)合,使用“重心”(cog)方法計算組243的模糊lat,其中隸屬函數(shù)被分成兩部分,使得曲線下的面積相等?,F(xiàn)在參見圖28,根據(jù)本發(fā)明的實施方案,該圖是示出根據(jù)一組最頂點的cog進行模糊lat計算的復(fù)合圖。以圖形方式作為曲線圖示出三個最頂點255、257、259的函數(shù)μu,并且在圖像顯示器261上,該函數(shù)被頂點周圍的區(qū)域離散化。cog263具有頂點255、257、259的函數(shù)μu的復(fù)合曲線圖265。許多頂點(ex節(jié)點)可能沒有指配函數(shù)μu。ex節(jié)點的函數(shù)μu的值通過評估從最頂點(pv節(jié)點)下面經(jīng)過的外推激活波來指配?,F(xiàn)在參見圖29,該圖示出根據(jù)本發(fā)明的實施方案的具有外推頂點的網(wǎng)格。如下計算頂點的函數(shù)μ:1.計算ex節(jié)點與網(wǎng)格上所有最頂點之間的最短距離。2.計算函數(shù)μuex=f(μupv,dpv,ex)的值。該計算與上述相鄰節(jié)點的初始化相同,不同的是僅考慮具有五個最高真實值的函數(shù)μuex的情況,如曲線圖267所示。3.使用cog方法通過步驟2中的組合函數(shù)μuex來計算模糊lat。模糊相干標測中的激活優(yōu)化:現(xiàn)在參見圖30,該圖是根據(jù)本發(fā)明的實施方案的模糊激活優(yōu)化方法的流程圖。在初始步驟269中,將函數(shù)μu的初始值指配至最頂點以及網(wǎng)格激活幀的鄰域節(jié)點。這些規(guī)程在上文已有所描述,并且在此不再重復(fù)。接著,在步驟271中,評估成本函數(shù)。成本函數(shù)的值是網(wǎng)格頂點的所有最大模糊真實值μu的平均值。通過降低近側(cè)頂點的隸屬函數(shù)的加權(quán)來增加該結(jié)果。接著,在決定步驟273中,確定當前迭代中的成本函數(shù)是否小于前一次迭代的成本函數(shù)(在預(yù)定義公差內(nèi))。如果確定為否定的,則控制前進至最終步驟275,并且規(guī)程終止。如果在決定步驟273,確定為肯定的,則控制前進至步驟277。基于包括使當前頂點的速度偏差矢量最小的優(yōu)化規(guī)程來更新函數(shù)μu的值。該規(guī)程與結(jié)合圖14所述的規(guī)程類似。為簡潔起見,此處不再重復(fù)細節(jié)。影響是使激活波的傳導(dǎo)性的一致性最大,該激活波以至少0.2m/s的速度行進。此類波應(yīng)當通常以連續(xù)方式傳播。完成步驟277之后,存在具有相應(yīng)函數(shù)μu的模糊lat集。隨著成本函數(shù)在逐次的迭代中進行了優(yōu)化,函數(shù)μu的值增加,并且其伸展變緩,如圖31所示,該圖是根據(jù)本發(fā)明實施方案的優(yōu)化前和優(yōu)化后網(wǎng)格頂點中函數(shù)μu的兩幅曲線圖?;氐綀D30,接著,在步驟279中,通過函數(shù)μu的當前值計算值latu。然后該規(guī)程在步驟271迭代。如果相鄰頂點中的至少一個頂點的μu>μumin,標記這些頂點用于修正。在每次迭代結(jié)束的同時,更新其中函數(shù)μu經(jīng)標記用于修正的所有模糊lat。在每次迭代步驟中,調(diào)整頂點處的模糊lat,以使周圍三角形的速度矢量與這些矢量的平均值之間的角度最小,如上文在公式1-4的討論中所述。為了消除阻滯線,為速度設(shè)置下限?,F(xiàn)在參見圖32,該圖是根據(jù)本發(fā)明的實施方案的示出了包括相鄰頂點的區(qū)域內(nèi)中心處的均勻傳導(dǎo)性的優(yōu)化的圖。該圖示出了模糊lat為28ms的頂點281。傳導(dǎo)速度在箭頭旁邊示出。具有模糊lat(lat2–lat6)的五個相鄰頂點圍繞具有如所示的相應(yīng)傳導(dǎo)矢量的頂點281。一般來講,當平均cv為正時,波以比其靠近時更快的速度遠離某一點,并且cog處的lat(lat1)必定下降。由于三角形之間的時差,lat的變化導(dǎo)致速度矢量的旋轉(zhuǎn)。為了處理該問題,速度偏差角通過使用逐步增加和逐步縮短頂點281處的中心lat的方式來計算。選擇產(chǎn)生較小差值角的方向?,F(xiàn)在參見圖33,該圖為根據(jù)本發(fā)明的實施方案的類似于圖32的圖。在該示例中,存在傳導(dǎo)阻滯283,并且函數(shù)μu沒有更新。當傳導(dǎo)速度低于預(yù)定最小值時,移除頂點之間的連接,之后進行算法迭代。本領(lǐng)域的技術(shù)人員應(yīng)當理解,本發(fā)明并不限于上文中具體示出和描述的內(nèi)容。相反,本發(fā)明的范圍包括上文所述各種特征的組合與子組合兩者,以及本領(lǐng)域的技術(shù)人員在閱讀上述說明時可能想到的未在現(xiàn)有技術(shù)范圍內(nèi)的變型和修改。當前第1頁12