本申請是申請日為2013年5月22日、申請?zhí)枮?01310333790.0、名稱為“外科手術導航系統”的中國發(fā)明專利申請的分案申請。
本發(fā)明涉及使用外科手術導航系統以執(zhí)行外科手術。更具體地說,本發(fā)明涉及使用導航系統,以利用基準圖案實現治療方案。
背景技術:
電外科設備已經被廣泛地使用。電外科手術涉及在外科手術期間利用熱能和/或電能來對生物組織進行切割、解剖、消融、凝結、燒灼、密封或其他處理。電外科手術通常使用手持件進行,所述手持件包括適于在電外科手術期間傳輸能量至組織部位的外科設備(例如,末端執(zhí)行器或消融探頭)、可操作地輸出能量的遠程電外科發(fā)生器、以及將所述外科設備可操作地連接至遠程發(fā)生器的電纜組件。
某些疾病的治療需要破壞惡性組織(例如腫瘤)生長。在諸如癌癥的疾病治療中,已經發(fā)現某些類型的腫瘤細胞在與通常會對健康細胞有害的溫度相比稍低的高溫下變性。已知的治療方法,諸如高溫療法,通常包括將患病細胞加熱至高于41℃的溫度而保持鄰近健康細胞低于不可逆的細胞破壞發(fā)生的溫度。這些方法可以涉及施加電磁輻射以加熱、消融和/或凝結組織。有多種不同類型的電外科設備能夠用于執(zhí)行消融術。
可以使用二維(2d)術前計算斷層掃描(ct)圖像和“消融區(qū)圖”來執(zhí)行癌性腫瘤或良性腫瘤的腫瘤微創(chuàng)消融術,“消融區(qū)圖”通常描述消融針在實驗的、離體組織中在一系列輸入參數(功率、時間)上的特征。能量劑量(功率、時間)能夠與消融組織效應(體積、形狀)關聯以進行具體設計。能夠通過微波天線設計控制輸送至組織的能量劑量,例如,可利用天線扼流圈來提供從設備至組織的在已知位置的微波傳輸。在另一實例中,電介質緩沖使得能夠進行獨立于不同的或變化的組織特性的、從設備至組織的相對恒定的能量輸送。
在用戶確定應該使用哪種消融針來實現目標的處理后,用戶通過超聲引導來進行處理。通常,需要高水平技能來將外科設備置于由超聲識別的目標中。其中最重要的是,有選擇引導設備朝向超聲圖像表面(例如,正在成像的目標處)所需的角度和進入點的能力。
超聲引導介入涉及使用實時超聲成像(經腹、術中等)來精確地引導外科設備至它們的預定目標。這可以通過經皮應用和/或術中應用來進行。在每種情形中,超聲系統將包括換能器,其對患者組織成像,并用于識別目標以及預期和/或跟隨朝向目標的器械的路徑。
超聲引導介入目前廣泛用于穿刺活檢術以確定已檢測的可疑病變(乳房、肝臟、腎臟、以及其他軟組織)的惡性。另外,用來獲得頸靜脈通路(jugularaccess)以及允許輸送藥物的中心線布置是共同的。最后,新興用途包括腫瘤消融和器官(肝臟、肺臟、腎臟,等等)的手術切除。在腫瘤消融的情形中,在實現超聲引導定位后,可利用活檢類針來輸送能量(rf、微波、冰凍,等等)以殺死腫瘤。在器官切除的情形中,對解剖期間表層下解剖結構的密切了解、以及外科設備關于該解剖結構的顯示,對于獲得成功的手術切緣而避開致命結構來說是關鍵的。
在這些情形中的每種情形中,超聲引導通常提供二維圖像平面,其從施加到患者的換能器的遠端獲取。為了成功定位設備,對用戶而言最重要的是能夠可視化和表征目標、選擇器械角度和進入點以到達目標、以及觀察外科設備及其朝向目標的運動。目前,用戶對目標成像并使用高水平技能來選擇器械角度和進入點。用戶必須移動超聲換能器來觀察器械路徑(因此丟失目標部位)或假定路徑是準確的直至設備到達圖像平面。最重要的是,選擇引導設備朝向超聲圖像平面(例如,正在成像的目標處)所需的角度和進入點。
技術實現要素:
本發(fā)明可使用短語“在一個實施方式中”、“在各個實施方式中”、“在一些實施方式中”、或“在其他實施方式中”,其中每個指代根據本發(fā)明的一個或多個相同或不同的實施方式。出于描述目的,以形式“a/b”的短語意指a或b。出于描述目的,以形式“a和/或b”的短語意指“(a)、(b)、或(a和b)”。出于描述目的,以形式“a、b、或c中的至少一個”的短語意指“(a)、(b)、(c)、(a和b)、(a和c)、(b和c)、或(a、b和c)”。
如附圖所示的以及整個下文說明書所描述的,在指代外科設備上的相對定位時與傳統一樣,術語“近側”指代設備更接近用戶或發(fā)生器的端,而術語“遠側”指代設備遠離用戶或發(fā)生器的端。術語“用戶”指代任何醫(yī)務人員(即,醫(yī)生、護士、等),其執(zhí)行包括使用本文所述的本發(fā)明各個方面的醫(yī)療程序。
如本說明書中所使用的,術語“外科設備”通常指代外科手術工具,其施加電外科能量以治療組織。外科設備可包括但不限于針、探針、導管、內窺鏡器械、腹腔鏡器械、血管閉合設備、外科縫合器等。術語“電外科能量”通常指代任何形式的電磁、光、或聲能。
電磁(em)能通常通過頻率的增加或波長的降低分類為無線電波、微波、紅外線、可見光、紫外線、x射線和伽馬射線。如本文所使用的,術語“微波”通常指代處于300兆赫(mhz)(3×108循環(huán)/秒)至300千兆赫(ghz)(3×1011循環(huán)/秒)的頻率范圍的電磁波。如本文所使用的,術語“rf”通常指代具有比微波更低頻率的電磁波。如本文所使用的,術語“超聲”通常指代具有大于人類聽力上限的頻率的循環(huán)聲壓(cyclicsoundpressure)。
如本說明書中所使用的,術語“消融術”通常指代任何消融術,諸如微波消融、射頻(rf)消融或微波消融輔助切除。如本說明書所使用的,“能量施加器”通常指代能夠用于從諸如微波或rf電外科發(fā)生器的能量發(fā)生源傳輸能量至組織的任何設備。
如本說明書中所使用的,術語“能量源”和“電源”指代以適于使電子電路運轉的形式的任何電能量源(例如,電池)。如本說明書中所使用的,“傳輸線路”通常指代能夠用于使信號從一點至另一點傳播的任何傳輸介質。如本說明書所使用的,術語“開關”或“各開關”通常指代任意電動致動器,機械致動器,電子機械致動器(可旋轉致動器、可樞轉致動器、觸發(fā)器類致動器、按鈕等),光學致動器,或通常實現電子設備、或其組件、器械、設備、輸送線路或連接線及其上附件、或軟件的連接和斷開連接目的的任何合適設備。
如本說明書所使用的,“電子設備”通常指代利用在真空、氣體、或半導體中運動的電子或離子的特性的設備或對象。如本文所使用的,“電子線路”通常指代電子或離子運動的路徑,以及由設備或對象為電子或離子提供的方向。如本文所使用的,“電氣電路”或簡單來說“電路”通常指代多個電氣設備和導體的組合,它們在彼此連接時形成導電路徑以實現期望的功能。除了互連部分以外,電路的任何組成部件可以稱為“電路元件”,其可包括模擬和/或數字組件。
術語“發(fā)生器”可指代能夠提供能量的設備。這樣的設備可包括能量源以及電路,所述電路能夠將能量源輸出的能量改變?yōu)榫哂衅谕麖姸?、頻率、和/或波形的輸出能量。
如本說明書中所使用的,“用戶接口”通常指代任何視覺、圖形、觸覺、聽覺、感官、或其他機構以提供信息至用戶或其他實體,和/或從用戶或其他實體接收信息。本文所使用的術語“用戶接口”可指代人類用戶(或操作者)與一個或多個設備之間的接口,以使得能夠在用戶和設備(多個)之間通信。本發(fā)明的各個實施方式中可利用的用戶接口的實例包括但不限于:開關、電位計、按鈕、撥號盤、滑動器、鼠標、定點設備、鍵盤、小鍵盤、操縱桿、軌跡球、顯示屏、各種類型的圖形用戶接口(gui)、觸摸屏、麥克風、以及其他類型的可接收人類發(fā)出的某種形式的刺激并產生響應其的信號的傳感器或設備。如本文所使用的,“計算機”通常指代以目的方式傳輸信息的任何事物。
本文所述的系統還可以利用一個或多個控制器來接收各種信息并變換所接收的信息以產生輸出??刂破骺砂ㄈ魏晤愋偷挠嬎阍O備、計算電路、或能夠執(zhí)行存儲在存儲器中的一系列指令的任何類型的處理器或處理電路??刂破骺砂ǘ鄠€處理器和/或多芯中央處理單元(cpu),并可包括任何類型的處理器,諸如微處理器、數字信號處理器、微控制器等。控制器還可包括存儲器以存儲數據和/或算法以執(zhí)行一系列指令。
本文所述的任何方法、程序、算法或代碼可以轉化為編程語言或計算機程序、或以編程語言或計算機程序來表達。“編程語言”和“計算機程序”是用于向計算機指定指令的任何語言,并包括(但不限于)這些語言和它們的派生:assembler、basic、batchfiles、bcpl、c、c+、c++、delphi、fortran、java、javascript、machinecode、操作系統命令語言、pascal、perl、pl1、腳本語言、visualbasic、自身指定程序的元語言、以及所有的第一、第二、第三、第四及第五代計算機語言。還包括數據庫和其他數據模式、以及任意其他元語言。出于限定目的,在經解釋、編譯的語言,或使用編譯和解譯方法的語言之間沒有區(qū)別。出于限定目的,在程序的編譯版和源版本之間沒有區(qū)別。因此,對于程序,當編程語言能夠以多于一種狀態(tài)(諸如源、編譯、對象、或鏈接語言)存在時,可以指代任何的狀態(tài)以及所有狀態(tài)。該限定還包括有效指令和這些指令的含義。
本文所述的任何方法、程序、算法或代碼可以包含在一個或多個機器可讀介質或存儲器中。術語“存儲器”可包括提供(例如,存儲和/或傳輸)機器(諸如處理器、計算機、或數字處理設備)可讀形式的信息的機構。例如,存儲器可包括只讀存儲器(rom)、隨機存取存儲器(ram)、磁盤存儲介質、光存儲介質、閃存設備、或其他任何易失性或非易失性存儲器存儲設備。包含在它們之上的代碼或指令可通過載波信號、紅外信號、數字信號、和其他類似信號來表達。
如本說明書中所使用的,短語“治療計劃”指代選定的消融針、能量水平、和/或治療持續(xù)時間以實現目標治療。術語“目標”指代準備要治療的組織區(qū)域,并且可包括且不限于腫瘤、纖維瘤、和待消融的其他組織。短語“消融區(qū)”指代將消融的組織區(qū)域和/或組織體積。
如本說明書中所使用的,短語“計算機斷層掃描”(ct)或“計算機軸向斷層掃描”(cat)指代利用由計算機處理產生的斷層掃描的醫(yī)療成像方法。數字幾何處理用于從繞單個旋轉軸采集的大量二維x射線圖像來生成對象內部的三維圖像。
如本說明書中所使用的,術語磁共振成像(mri)、核磁共振成像(nmri)、或磁共振斷層掃描(mrt)指代在放射學中使用的用以可視化內部細節(jié)結構的醫(yī)療成像技術。mri利用核磁共振(nmr)特性來對體內的原子核成像。mri機器使用強力磁場來排列(align)體內的一些原子核的磁化,同時使用射頻場來系統地改變該磁化的排列。這引起核產生可由掃描器檢測的旋轉磁場,并且該信息被記錄以構建所掃描身體區(qū)域的圖像。
如本說明書中所使用的,術語“三維超聲”或“3d超聲”指代提供三維圖像的醫(yī)療超聲技術。
如本說明書所使用的,術語“醫(yī)學數字成像和通信”(dicom)指代用于處理、存儲、打印、和傳輸涉及醫(yī)療成像的信息的標準。它包括文件格式定義和網絡通信協議。通信協議是應用協議,其使用tcp/ip以在各系統間通信。dicom文件可以在能接收dicom格式的圖像和患者數據的兩個實體之間交換。
本文所述的任意系統和方法可通過有線網絡、無線網絡、點對點通信協議、dicom通信協議、傳輸線路、可移動存儲介質等在其間傳輸數據。
本文所述系統可利用被配置為檢測組織的一個或多個屬性和/或周圍環(huán)境的一個或多個傳感器。所述屬性包括但不限于:組織阻抗、組織類型、組織透明度、組織順應性、組織或鉗口構件(jawmember)的溫度、組織中的水含量、鉗開口角度、組織中的水移動性、能量輸送、以及鉗閉合壓。
在本發(fā)明的一個方面,提供了一種導航系統。所述導航系統包括:超聲設備,具有布置于其上的基準圖案,其被配置以獲得掃描平面中的超聲圖像;以及外科設備,其具有圖像獲取設備,被配置以獲取基準圖案的基準圖像??刂破鞅慌渲靡越邮粘晥D像和基準圖像,其中所述控制器基于基準圖像確定外科設備關于掃描平面的位置,并且顯示器被配置以基于外科設備關于掃描平面的位置來顯示外科設備的虛擬圖像和超聲圖像。
在導航系統中,基準圖案固定至超聲設備上的已知位置,并且圖像獲取設備固定至外科設備上的已知位置?;鶞蕡D案具有存在于已知拓撲結構中的具有相對位置和已知特征的多個標記。控制器將基準圖像對應至模型圖像,評估攝影機的位姿、并將外科設備轉換至模型坐標系??刂破鬟€針對透鏡畸變來校正基準圖像。另外,控制器能夠識別基準標記中的拓撲結構,其中拓撲結構參照兩個或更多個獨立唯一標識符,它們定位在標記上的單個圖案上的已知位置。
附圖說明
本發(fā)明的上述和其他方面、特征、和優(yōu)勢將在結合附圖時根據下面的詳細描述而變得更加明顯,其中:
圖1是根據本發(fā)明的一個實施方式的規(guī)劃和導航系統的系統框圖;
圖2a和2b是根據本發(fā)明的一個實施方式的消融針的示意圖;
圖3是圖2a和2b的消融針的輻射圖的示意圖;
圖4是根據本發(fā)明的一個實施方式的規(guī)劃系統的示意圖;
圖5是描繪根據本發(fā)明的一個實施方式的規(guī)劃系統的整體操作的流程圖;
圖6和7是用于根據本發(fā)明的一個實施方式的規(guī)劃系統中的圖形用戶界面的示意圖;
圖8是描繪根據本發(fā)明的一個實施方式的圖像分割和逆向規(guī)劃的算法的流程圖;
圖9是描繪根據本發(fā)明的一個實施方式的用于分割結節(jié)的算法的流程圖;
圖10a-圖10b是消融區(qū)和能量輸送之間的關系的圖形表示;
圖11a是根據本發(fā)明另一實施方式的血管和目標之間的關系的示意圖;
圖11b是根據本發(fā)明另一實施方式的替代劑量曲線的圖形表示;
圖12a-圖12c是根據本發(fā)明另一實施方式的規(guī)劃方法的示意圖;
圖13是根據本發(fā)明的一個實施方式的導航系統的示意圖;
圖14a和14b是用于圖13的導航系統的圖形用戶接口的示意圖;
圖15是描繪根據本發(fā)明的一個實施方式的基準跟蹤算法的流程圖;
圖16a和16b分別描繪了由攝影機獲得的圖像以及該圖像的校正版本;
圖17是描繪根據本發(fā)明的一個實施方式的用于找到白色圓圈的算法的流程圖;
圖18a-圖18c描繪了圖17中所描繪的算法的中間圖像結果;
圖19是描繪根據本發(fā)明的一個實施方式的用于找到黑色圓圈和黑色區(qū)域的算法的流程圖;
圖20a-圖20d描繪了圖19中所描繪算法的中間圖像結果;
圖21a是描繪根據本發(fā)明的一個實施方式的對應算法的流程圖;
圖21b是描繪根據本發(fā)明的一個實施方式的用于施加拓撲限制的算法的流程圖;
圖22a-圖22d是用于圖21a的算法的基準模型的示意圖;
圖23是根據本發(fā)明另一實施方式的集成的規(guī)劃和導航系統的示意圖;
圖24是根據本發(fā)明又一實施方式的集成的規(guī)劃和導航系統的示意圖;
圖25a和圖25b是適于與圖24的系統一起使用的導航系統的示意圖;以及
圖26-圖29是根據本發(fā)明各個實施方式的用于圖24的系統的圖形用戶界面的示意圖。
具體實施方式
本發(fā)明的特定實施方式將在下文中參照附圖進行描述;然而,應該理解的是,所公開的實施方式僅是本發(fā)明的實例并且可以各種形式體現。沒有詳細描述公知的功能或構造以避免以不必要的細節(jié)使本發(fā)明晦澀。因此,本文公開的具體結構和功能細節(jié)不能解釋為限制,而僅僅作為權利要求書的基礎,以及作為教導本領域技術人員以實際上任何合適的細節(jié)結構來不同地利用本發(fā)明的代表性基礎。在整個附圖描述中同樣的參考數字指代相似或相同的元件。
轉向附圖,圖1描繪了根據本發(fā)明各實施方式的規(guī)劃和導航系統的整體視圖。如圖1所示,患者“p”的術前圖像15通過圖像獲取設備10獲取。圖像獲取設備10可包括但不限于獲得二維(2d)或三維(3d)圖像的mri設備、cat設備、或超聲設備。圖像獲取設備10存儲傳輸至規(guī)劃系統100的術前圖像15。術前圖像15可以通過將圖像15上傳到網絡、通過無線通信裝置傳送圖像15至規(guī)劃系統100、和/或將圖像15存儲在插入規(guī)劃系統100的可移動存儲器中而傳輸至規(guī)劃系統100。在本發(fā)明的一個實施方式中,術前圖像15以dicom格式保存。在一些實施方式中,圖像獲取設備10和規(guī)劃系統100可結合為一個獨立單元。
規(guī)劃系統100(其將在下文更加詳細描述)接收術前圖像15并確定目標尺寸?;谀繕顺叽绾瓦x定的外科設備,規(guī)劃系統100確定設置,其包括實現目標治療的能量水平和治療持續(xù)時間。
導航系統200,其在下文將更加詳細描述,利用布置于醫(yī)療成像設備(例如,超聲成像設備)上的基準圖案來確定外科設備的體內位置。外科設備的體內位置在顯示設備上關于由醫(yī)療成像設備獲得的圖像來顯示。一旦外科設備定位在目標附近,用戶就基于由規(guī)劃系統確定的治療區(qū)設置來實現目標的治療。
在一些實施方式中,用戶使用規(guī)劃系統100來確定治療區(qū)設置,并在使用導航系統200實現治療時利用治療區(qū)設置。在其他實施方式中,當外科設備處于目標附近時,規(guī)劃系統100傳輸治療區(qū)設置至導航系統200以自動實現目標的治療。另外,在一些實施方式中,規(guī)劃系統100和導航系統200結合為單個獨立系統。例如,單個處理器和單個用戶接口可用于規(guī)劃系統100和導航系統200,單個處理器或多個用戶接口可用于規(guī)劃系統100和導航系統200,或者多個處理器和單個用戶接口可用于規(guī)劃系統100和導航系統200。
圖2a示出了根據本發(fā)明一個實施方式的外科設備的實例。具體而言,圖2a示出了具有扼流圈72的消融針60的變體的側視圖,以及圖2b示出了沿圖2a的2b-2b的剖面圖。消融針60示出了輻射部62,其通過饋線(或軸桿)64電連接至近側定位的耦合器66。輻射部62被示出為具有涂覆在節(jié)段62上的密封層68。扼流圈72被示出為部分地布置在饋線64的遠側節(jié)段上以形成扼流部70,其靠近輻射部62定位。
為了改善消融針60的能量聚焦,扼流圈72用于將場傳播或輻射圖案包含在消融針60的遠端。通常,扼流圈72布置在消融針60上靠近輻射段處。扼流圈72置于介電材料上,所述介電材料布置在消融針60上。扼流圈72是導電層,其可以通過管路或涂層覆蓋以迫使該導電層與底層消融針60相符合,進而形成更加遠離或接近輻射部62的電連接(或短路)。扼流圈72和底層消融針60之間的電連接還可以通過其他連接方法(諸如軟焊、焊接、銅焊、壓接、使用導電粘合劑等)實現。消融針60電耦接至發(fā)生器,其為消融針60提供電外科能量。
圖3是根據本發(fā)明的消融針60的一個實施方式的剖面圖,其示出以發(fā)射的輻射圖的圖形表示。
圖4-圖12c描述了根據本發(fā)明各個實施方式的規(guī)劃系統100的操作。轉向圖4,規(guī)劃系統100包括接收器102、存儲器104、控制器106、輸入設備108(例如,鼠標、鍵盤、觸摸板、觸屏等)、以及顯示器110。在規(guī)劃系統100的操作中,接收器102接收dicom格式的術前圖像15并將這些圖像存儲在存儲器104中??刂破?06隨后處理圖像15(這將在下文更詳細描述),并將經處理的圖像顯示在顯示器110上。使用輸入設備108,用戶可以導航通過圖像15,從圖像15中選擇一幅圖像,在所選擇圖像上選擇種子點,選擇消融針,調節(jié)能量水平,以及調節(jié)治療持續(xù)時間。由輸入設備108提供的輸入顯示在顯示器110上。
圖5描繪了由規(guī)劃系統100使用的用以確定治療方案的的算法的概圖。如圖5所示,在步驟120中,通過無線連接、網絡、或通過從可移動存儲介質下載圖像來獲取dicom格式的圖像,并將其存儲在存儲器104中??刂破?06隨后執(zhí)行圖像15的自動三維(3d)渲染,并在步驟122中顯示3d渲染的圖像(如圖6所示)。在步驟124,執(zhí)行圖像分割以對特定感興趣區(qū)域劃界并計算感興趣區(qū)域的體積。如下所述,分割可以是用戶驅動或自動的。在步驟126,控制器執(zhí)行逆向規(guī)劃操作,其還將在下文中更為詳細地描述,從而確定治療感興趣區(qū)域的治療算法。該治療算法可包括外科設備、能量水平、和/或治療持續(xù)時間的選擇。作為選擇,用戶可以選擇外科設備、能量水平、和/或治療持續(xù)時間以滿足治療醫(yī)生的意愿,其將包括治療目標的“邊緣值”以及周圍組織的邊緣。
圖6和7描繪了可在顯示器110上顯示的圖形用戶界面(gui)。如圖6和7所示,每個gui分割為多個區(qū)域(例如,區(qū)域132、134和136)以用于顯示經渲染的dicom圖像。例如,區(qū)域132顯示了患者“p”沿橫剖面的圖像,以及區(qū)域134顯示患者“p”沿冠狀剖面的圖像。區(qū)域136描繪了患者“p”的3d渲染。在其他實施方式中,矢狀面也可顯示在gui上。gui允許用戶在下拉菜單131中選擇不同的消融針。gui還允許用戶分別在區(qū)域133和135中調節(jié)功率和時間設置。此外,gui在區(qū)域137中具有多個附加工具,其包括但不限于:啟動種子點選擇的規(guī)劃工具、對比工具、縮放工具、拖動工具、用于滾動通過dicom圖像的滾動工具、以及用于顯示dicom數據集的體渲染的3d渲染工具。
圖8的流程圖描繪了用于執(zhí)行圖像分割步驟124和逆向規(guī)劃步驟126的基本算法。如圖8所示,用戶在步驟140中選擇種子點(參見圖6,其中,十字線位于區(qū)域132和134中目標“t”的中心)。在手動選擇種子點后,規(guī)劃系統100在步驟142中分割結節(jié)以對感興趣體積劃界。在其他實施方式中,可基于像素的強度值自動檢測種子點。
圖9描繪了用于分割結節(jié)的算法的流程圖。如圖9所示,一旦在步驟151中識別了種子點,該算法在步驟152中創(chuàng)建感興趣區(qū)域(roi)。例如,roi可包括4cm3的體積。在步驟153中,連接的閾值濾波器施加閾值并在存儲于存儲器104中的dicom圖像中尋找連接至種子點的所有像素。例如,在分割肺結節(jié)時,閾值可開始于-400houndsfield單位(hu)并終止于100hu。
在步驟154,控制器106應用幾何濾波器來計算對象的尺寸和形狀。幾何濾波器使得能夠測量標示體積中所有對象的幾何特征。該標示體積能夠例如表示分割為不同解剖結構的醫(yī)療圖像。這些對象的各個幾何特征的測量結果能夠提供對圖像的其他深入了解。
該算法在步驟155中確定是否檢測到預定形狀。如果未檢測到預定形狀,則算法行進至步驟156,其中閾值增加預定值。該算法重復步驟153至155直至檢測到預定的對象。
一旦檢測到預定對象,該算法在步驟157處終止,并且規(guī)劃系統100行進至步驟144以執(zhí)行體積分析。在體積分析期間,可通過控制器106計算球形對象的下述屬性:最小直徑、最大直徑、平均直徑、體積、球形度、最小密度、最大密度、以及平均密度。所計算的屬性可以在顯示器110上顯示,如圖7的區(qū)域139中所示。體積分析可使用幾何濾波器來確定最小直徑、最大直徑、體積、延伸率、表面積、和/或球形度。在步驟144,圖像強度統計濾波器還可與幾何濾波器結合使用。圖像強度統計濾波器計算最小密度、最大密度以及平均密度。
在步驟146中,針對劃界目標計算功率和時間設置。圖10描繪了對于給定時間段沉積至組織中的能量和所產生的消融區(qū)域之間的關系的各個曲線圖。該關系通過考慮目標組織(即,腫瘤、纖維瘤等)的維度和特征以及指定消融針的能量劑量/天線設計而允許逆向規(guī)劃。下表1示出了消融針的消融體積、功率、和時間之間的關系的示例。
表1
使用表1中的值,能夠從該表中導出線性等式,從而計算最佳功率和時間設置。例如,使用線性回歸分析,表1提供了下面等式:
體積=0.292381*功率+8.685714*時間-44.0762(1)
其能夠寫成
功率=(體積-8.685714*時間+44.0762)/0.292381(2)
能夠利用來自體積分析的最大直徑加上1厘米邊際來計算期望的體積,如下:
期望體積=4/3*pi*期望半徑^3(3)
其中,期望半徑可如下計算:
期望半徑=最大結節(jié)直徑/2+邊際(4)
將期望體積代入等式(1)或(2)留下兩個未知數:功率和時間。使用等式(2),控制器106能夠通過代入時間值解出功率??刂破?06選擇將功率保持低于70w或某個其他預定值的時間最小值,從而用戶可以盡可能快地執(zhí)行該過程并同時保持功率處于安全范圍。
一旦計算了功率和時間146,功率和時間顯示在顯示器110上,如圖7所示(參見133和135)。用戶可分別使用控件133和135調節(jié)計算的功率和/或時間,從而調節(jié)治療區(qū)138a和/或邊際138b。
存儲器104和/或控制器106可存儲對應于不同外科設備的多個等式。當用戶在下拉菜單131中選擇了不同的外科設備時,控制器106可執(zhí)行上面所述的相同分析來確定保持功率低于70w或某個其他預定值的最小時間值。
雖然上述過程描述了使用單個種子點來確定預定對象,一些目標可具有不規(guī)則形狀,其不能夠通過預定治療區(qū)來治療而不對其他組織造成損傷。在這樣的情形中,可使用多個種子點來創(chuàng)建不規(guī)則形狀的治療方案,其使用以多個位置重新定位的單個外科設備或可同時使用的多個外科設備來治療不規(guī)則形狀的區(qū)域。
在其他實施方式中,存儲器104和/或控制器106可存儲外科設備和治療區(qū)性能的目錄,其包括功率、時間、器械數量、以及到達體外或體內治療區(qū)所需的器械間隔?;趫D像分割和體積分析的結果,控制器可自動選擇設備類型、設備數量、多個設備的間隔、和/或針對每個設備的功率及時間設置,以對roi進行治療。作為選擇,用戶可以使用gui手動選擇治療roi的設備類型、設備數量、多個設備的間隔、針對每個設備的功率和/或時間設置,從而產生治療方案。
在根據本發(fā)明的另一實施方式中,規(guī)劃系統100還可對除了目標之外的器官和其他重要結構進行分割。器官和其他結構(如血管)的分割用于提供更加先進的治療方案。如上關于圖10所述的,治療區(qū)與能量輸送以規(guī)則方式相關聯。此外,公知的是大于三(3)毫米的血管可消極地影響治療區(qū)形成。血管的分割將允許血管和待評估目標之間的交互,包括血管直徑(d1)以及血管與提出的目標之間的距離(d2)(參見圖11a)。該交互可通過用戶手動地評估或通過控制器106自動地評估。使用血管直徑d1和距離d2,規(guī)劃系統100可自動地建議用于治療目的的替換劑量曲線,如圖11b所示。作為選擇,控制器106可通過顯示器110為用戶推薦以移動治療區(qū)。另外,不同的治療區(qū)投影可以顯示在顯示器110上。此外,在圖8的計算功率和時間設置步驟146中,控制器可基于血管直徑和到目標區(qū)域的距離來影響不同的曲線。
圖12a-圖12c描繪了使用器官分割的先進治療規(guī)劃。器官分割在規(guī)劃療程中具有至少兩個優(yōu)勢。在第一個情形中,通常選擇微創(chuàng)治療以進行器官保留。通過分割器官,控制器106可以計算器官體積160并減去確定的消融區(qū)162以確定保留的器官164的體積,如圖12a所示。如果控制器106確定保留的器官體積太小,則控制器106可以提醒用戶需要替換治療方案或者它可以建議替代治療方案。
圖12b和12c描繪了針對定位在器官表面上的目標“t”的治療方案。傳統上,通常避免靠近器官表面的治療,或需要另外的技術以在執(zhí)行治療前將該器官與其他器官隔離。在根據本發(fā)明的另一實施方式中,在分割器官后,還可以確定目標“t”的位置。如果治療方案中的治療區(qū)162凸出于器官表面外并且目標“t”位于表面上,則控制器106可以提醒用戶治療區(qū)162可影響在目標“t”附近的其他器官和/或結構,以及該治療方案需要變更。在另一實施方式中,控制器106可自動為用戶提出指示外科設備、能量水平、治療持續(xù)時間的建議。控制器106還可建議如圖12b中示出的較小治療區(qū)162,或者它可以建議如圖12c所示的移動治療區(qū)162。
在其他實施方式中,在分割目標、組織、器官、和其他結構后,已知的組織性質可歸因于這些結構。這些組織性質包括但不限于:跨頻率的電導率和介電常數、熱導率、熱對流系數,等等。圖8的規(guī)劃算法可使用歸因于分割腫瘤、組織、器官、和其他結構的組織性質來求解pennes生物傳熱方程,從而計算出消融選定目標所需的劑量。成功實現使用生物傳熱方程更全面求解的關鍵包括:利用處于穩(wěn)態(tài)的已知組織性質來預測初始空間溫度分布,利用溫度上升時的組織性質來根據溫度的升高調節(jié)空間性質,以及利用處于液-氣相變的組織性質。
轉向圖13,根據本發(fā)明一個實施方式的導航系統通常被示為200。一般而言,導航系統200結合固定至超聲換能器202的參考貼片或基準貼片204?;鶞寿N片204可印制在超聲換能器202上、通過粘結劑附接至超聲換能器202、或可移除地耦接至超聲換能器202。在一些實施方式中,基準貼片布置在支撐結構上,所述支撐結構被配置成可移除的固定,如“夾在”,超聲換能器的外殼上。超聲換能器202耦接至產生聲波的超聲發(fā)生器210。超聲換能器202和超聲發(fā)生器210可組合為獨立單元。超聲換能器202朝患者“p”發(fā)射聲波。所述聲波從患者“p”內的各結構反射離開并由超聲換能器202接收。超聲換能器202將所反射的聲波傳輸至超聲發(fā)生器210,其將所反射的聲波實時轉換為二維(2d)圖像。所述2d圖像傳輸至控制器212。控制器212處理2d圖像并在顯示器214上顯示該2d圖像作為包括目標220的圖像218。圖像218是可以包括目標“t”的掃描平面“s”的實時表示。
導航系統還結合有固定至外科設備206的攝影機208。攝影機208實時地獲取基準貼片204的圖像以確定外科設備206關于掃描平面“s”的位置。具體地講,基準貼片204具有針對掃描平面“s”的限定的空間關系。該限定的空間關系存儲在控制器212中。攝影機208還具有針對外科設備206的已知空間關系,其存儲在控制器212中。為了確定外科設備206和掃描平面“s”之間的空間關系,攝影機208獲取基準貼片204的圖像并將該圖像傳輸至控制器212。通過使用基準貼片204的圖像,控制器212可計算外科設備206和掃描平面“s”之間的空間關系。
在控制器212確定了外科設備206和掃描平面“s”之間的空間關系后,控制器212在顯示器214上顯示該關系。如圖13所示,顯示器214包括掃描平面“s”的圖像218,其包括目標“t”的目標圖像220。另外,控制器212關于圖像218疊加外科設備206的虛擬圖像206a以指示外科設備206關于掃描平面“s”的位置。基于消融針206的角度和位置,控制器212可計算外科設備206的軌跡,并顯示所計算的軌跡,通常以216示出。在一些實施方式中,十字線或目標可疊加在圖像218上以指示外科設備206將與掃描平面“s”相交的位置。在其他實施方式中,所計算的軌跡216可以紅色或綠色示出以指示導航狀態(tài)。例如,如果外科設備206處于將要與目標“t”相交的路徑上,所計算的軌跡216將以綠色示出。如果外科設備206未處于將要與交叉目標“t”相交的路徑上,所計算的軌跡216將以紅色示出。
控制器212還可以通過用戶輸入外科設備類型、能量水平、以及治療持續(xù)時間來控制。外科設備類型、能量水平、以及治療持續(xù)時間可以在顯示器214上顯示,如圖14a所示。在外科設備206與目標“t”相交時,虛擬消融區(qū)222投影在圖像218上,如圖14b所示。能量水平和治療持續(xù)時間則可由用戶調節(jié),并且控制器212將調節(jié)虛擬消融區(qū)222以反映能量水平和治療持續(xù)時間的變化。
下面將參照圖15-圖22描述基準跟蹤系統。在基準跟蹤系統中,控制器212接收來自攝影機208的基準圖像??刂破?12還包括先前存儲于其上的攝影機208的攝影機校準和畸變系數、基準系統模型、以及攝影機-天線校準數據。在其他實施方式中,攝影機208的攝影機校準和畸變系數、基準系統模型、以及攝影機-天線校準數據可以在導航過程期間輸入至控制器212中?;诨鶞蕡D像、攝影機208的攝影機校準和畸變系數、基準系統模型、以及攝影機-天線校準數據,控制器212可將消融針206的位置以及診斷幀速率、殘留誤差、和跟蹤狀態(tài)輸出至顯示器214。在一些實施方式中,攝影機208和基準貼片204之間的距離可以在約5厘米至約20厘米的范圍內。在一些實施方式中,攝影機208和基準貼片204之間的距離可以在約1厘米至約100厘米的范圍內。
圖15示出了控制器212所利用的基準追蹤算法的基本流程圖。如圖15所示,在步驟230中獲取圖像幀。在步驟231,控制器212使用攝影機校準和畸變系數來校正透鏡畸變。由攝影機208獲取的圖像可呈現透鏡畸變,如圖16a所示。因此,在圖像能夠用于進一步計算前,圖像需要對畸變進行校正。在導航過程中使用攝影機208前,攝影機208用于采集不同角度上棋盤狀圖案的多幅圖像。該多幅圖像和多個角度用于創(chuàng)建攝影機矩陣和畸變系數??刂破?12隨后利用攝影機矩陣和畸變系數來對透鏡畸變進行校正。
在步驟232中,控制器212在圖像幀中使用圖17的算法找到白色圓圈。如圖17所示,步驟241(圖18a)中所接收到的圖像幀在步驟243中使用動態(tài)閾值(參見圖18b)閾值化。當使用動態(tài)閾值時,在每個有效幀后,動態(tài)閾值算法針對下一幀使用在有效幀中發(fā)現的圓圈來計算新的閾值。使用在有效幀中發(fā)現的圓圈,控制器212基于下面公式(5)來計算新的閾值:
閾值=(黑圓圈強度平均+白圓圈強度平均)/2(5)
可使用預定閾值來獲取初始有效幀,其隨后用于計算新的閾值。
作為選擇,控制器212可通過測試閾值范圍來掃描尋找初始閾值,直至發(fā)現產生有效幀的閾值為止。一旦發(fā)現了初始閾值,控制器212將基于有效幀使用公式(5)以用于動態(tài)閾值。
在其他實施方式中,可使用固定閾值。固定閾值可以是存儲在控制器212中的預定數字,或者它可以這樣確定,即,測試閾值范圍直至發(fā)現導致有效幀的閾值為止。
在閾值和自動增益控制被施加至圖像后,在步驟244中執(zhí)行連通分量分析以發(fā)現閾值化圖像中的所有對象。在步驟245,將幾何濾波器應用至連通分量分析的結果以及圖像幀。幾何濾波器計算對象的尺寸和形狀,并僅保持圓形及差不多如圖18c中所示的合適尺寸的那些對象。對于所有圓形對象計算加權質心并將其存儲。
返回圖15,除了在步驟232中找到白色圓圈外,控制器212還在步驟233中使用圖19中描繪的算法來找到黑色圓圈。用于找到黑色圓圈的算法類似于圖17中所示的用于找到白色圓圈的算法。為了找到黑色圓圈,在步驟241中接收到圖像幀后(參見圖20a),控制器212在步驟242中反轉圖像幀的強度,如圖20b所示。然后,如上面關于圖17所述的,對圖像進行如圖20c所示的閾值化,并執(zhí)行連通分量分析,并應用幾何濾波器來獲得圖20d所示的圖像。在步驟248中對于所有的黑色圓圈計算加權質心并將其存儲。此外,在步驟245,控制器212應用幾何濾波器來確定除圖像幀內的黑色圓圈之外的黑色區(qū)域??刂破?12在步驟249中存儲已確定的黑色區(qū)域。
在圖15的步驟234,控制器212使用圖21a中所示的算法找到基準圖像和基準模型之間的對應關系。在圖21a的步驟251,控制器212使用拓撲限制來選擇四個白色圓圈,如圖21b所示。如圖21b所示,在步驟261,控制器212獲得圖19的步驟249中存儲的黑色區(qū)域并獲得圖17的步驟246中存儲的白色圓圈??刂破?12接著在步驟263中選擇第一黑色區(qū)域,并在步驟264中對第一黑色區(qū)域中的白色圓圈數量進行計數。在步驟265,控制器212確定位于選定黑色區(qū)域中的圓圈的數量是否匹配圓圈預定數量。如果圓圈數量不匹配圓圈預定數量,則算法行進至步驟266,在該步驟選擇下一個黑色區(qū)域,并且在步驟264中對下一黑色區(qū)域中的圓圈數量再次計數。該過程重復直至步驟264中計數的圓圈數目匹配圓圈預定數目為止。一旦步驟264中計數的圓圈數目匹配圓圈預定數目,則算法行進至步驟267,在該步驟完成拓撲限制算法。在其他實施方式中,控制器212通過選擇四個最圓的圓圈來選定四個白色圓圈。
在選擇了四個圓圈后,它們使用步驟252中的凸殼算法以順時針順序布置。針對實向量空間v中點集x的凸殼或凸狀包絡是包含x的最小凸集。如果所有點都位于一條直線上,則凸殼是連接最外側兩點的線段。在平面情形中,凸殼是凸多邊形,除非所有點位于相同直線上。類似地,在三維中,凸殼通常是最小凸多面體,其包含該集中的所有點。另外,模型中的四個匹配基準點也以順時針順序布置。
在步驟253,計算平面單應矩陣。在計算了平面單應矩陣后,使用該單應矩陣來將基準模型轉換為圖像坐標,其使用圖22中所示的四個對應基準模型以找到最近的匹配圖像基準點(步驟254和255)??刂破?12還在步驟256中計算殘留誤差。該算法使用產生的3d變換來將3d基準模型轉換為2d圖像。它隨后比較映射在2d圖像中的基準點與2d圖像中檢測到的基準點之間的距離。該殘留誤差是像素的平均距離。該誤差用于檢驗精度以及部分地確定紅色/綠色導航狀態(tài)??刂破?12則選擇具有最高匹配度和最小殘留誤差的模型。為了更加精確的結果,需要最小數量(例如,三個)的黑色基準點匹配。
在圖15的步驟235中,執(zhí)行攝影機的位姿評估。攝影機位姿評估包括通過迭代地將模型基準點轉換至基準圖像平面上以及最小化像素殘留誤差,來計算攝影機和選定模型之間的3d變換。目標是找到誤差函數的全局最小值??砂l(fā)生的一個問題是誤差函數中出現的顯著的局部最小值(例如,從左側成像的天線看上去類似于從右側成像的天線),這需要避免??刂破?12通過從多個起始點執(zhí)行最小化以及選擇具有最小誤差的結果來避免局部最小值。一旦計算了3d變換,控制器就可以使用3d變換來將外科設備206的坐標轉換至模型空間,并將外科設備206顯示為顯示器214中的虛擬外科設備206a。
由于對象邊界在不同照明條件下擴展和收縮,傳統的方角基準點位置可基于照明條件而變化?;鶞寿N片204使用黑色和白色圓圈,并因此不被該問題妨礙,因為圓圈的中心始終保持不變并持續(xù)良好地工作以計算加權質心。其他對比圖像或顏色也是能預料的。
在本發(fā)明的另一實施方式中,以及如圖23所示,提供了規(guī)劃和導航系統300。系統300包括規(guī)劃系統302和導航系統304,它們連接至控制器306??刂破?06連接至顯示器308,其可包括單個顯示屏或多個顯示屏(例如,兩個顯示屏)。規(guī)劃系統302類似于規(guī)劃系統100,并且導航系統304類似于導航系統200。在系統300中,顯示器308顯示上面所描述的規(guī)劃操作和導航操作。規(guī)劃操作和導航操作可以作為分區(qū)屏幕布置顯示在單個顯示屏上,規(guī)劃操作和導航操作可顯示在單獨的屏幕上,或者規(guī)劃操作和導航操作可顯示在相同屏幕上并且用戶可在各視圖之間切換??刂破?06可從規(guī)劃系統引入劑量設置,并在導航操作期間使用該劑量設置來顯示消融區(qū)維度。
在本發(fā)明的其他實施方式中,ct導航和軟件可以與規(guī)劃系統100集成。轉向圖24、25a和25b,規(guī)劃和導航系統通常以400示出。系統400包括圖像獲取設備402,其獲取患者“p”的ct圖像,患者“p”具有電磁參考點428和/或光學參考點438。ct圖像以dicom格式提供至規(guī)劃系統404,其與規(guī)劃系統100類似。規(guī)劃系統400用于確定如上所述的治療方案,并且該治療方案被提供至控制器408并在顯示器410上作為規(guī)劃屏412顯示,如圖26所示。
導航系統406可利用圖25a所示的電磁跟蹤系統、或圖25b所示的紅外線跟蹤系統或光學跟蹤系統。轉向圖25a,導航系統420包括電磁場發(fā)生器422、具有電磁換能器426的外科設備424、以及布置在患者上的電磁參考點428。場發(fā)生器422發(fā)射電磁波,該電磁波由外科設備424和電磁參考點428上的電磁傳感器(未明確示出)檢測,然后用于計算外科設備424和電磁參考點428之間的空間關系。該空間關系可由場發(fā)生器422計算,或場發(fā)生器422可提供數據至控制器408以計算消融針424和電磁參考點428之間的空間關系。
圖25b描繪了可替代的導航系統430,其類似于上面圖13所述的導航系統。在圖25b中,光學參考點或基準點438置于患者上。附接至外科設備424的攝影機436采集基準點438的圖像并傳輸該圖像至控制器408以確定消融針關于基準點438的位置。
在從導航系統406接收數據后,控制器408可將外科設備424的位置與ct圖像關聯,從而將外科設備424導航至目標“t”,如下文所述。在該情形中,患者參考點(任意類型)上也可具有不透射線的標記,以允許ct期間可視化。這允許控制器將患者ct圖像坐標系關聯至器械跟蹤坐標系。
控制器408和顯示器410彼此協作以將ct圖像顯示在導航屏440上,如圖27所示。如圖27所示,顯示屏440包括橫向視圖442、冠狀圖444、以及矢狀圖446。每個視圖包括目標“t”和消融區(qū)452(包括邊際)的視圖。橫向視圖442、冠狀圖444、以及矢狀圖446、消融區(qū)452全部從規(guī)劃系統404引入。另外,所有規(guī)劃元件(例如,設備選擇、能量水平、以及治療持續(xù)時間)都自動的傳輸至導航屏440。導航屏440也是圖形用戶界面,其允許用戶調節(jié)設備選擇、能量水平、以及治療持續(xù)時間。
導航引導屏448設置在顯示屏440上以輔助將消融針導航至目標“t”?;趶膶Ш较到y406接收的數據,控制器可確定外科設備424是否與目標“t”對準。如果外科設備424沒與目標“t”對準,則圓圈454將偏離外圓圈453的中心。用戶然后將調節(jié)外科設備424的進入角度直至圓形454的中心與外圓圈453的中心對準。在一些實施方式中,在圓圈454的中心與外圓圈453的中心不對準時,圓圈454可顯示為紅色圓圈,或者在圓圈454的中心與外圓圈453的中心對準時,圓圈454可顯示為綠色圓圈。另外,控制器408可計算目標“t”和外科設備424之間的距離。
在圖28描述的另一實施方式中,控制器408將虛擬外科設備424a疊加在3d渲染圖像上,并在屏幕462上顯示該組合圖像。類似于上面所述的方法,用戶可以將圓圈453的中心與圓圈454的中心對準以將外科設備424導航至目標“t”。作為選擇,用戶可以通過觀察屏幕462上的虛擬外科設備424a來確定外科設備424關于目標“t”的位置,從而將外科設備424導航至目標“t”。
圖29描繪了本發(fā)明的另一實施方式,類似于上面的屏幕462,在圖29的實施方式中,屏幕472描繪了與先前采集并渲染的ct圖像成空間關系的虛擬外科設備424a。ct圖像已經被體渲染以劃界目標“t”及另外的結構、血管、和器官。通過對目標“t”及另外的結構、血管、和器官進行體渲染,用戶能夠將外科設備424導航進入患者,而同時避開了另外的結構、血管、和器官以防止不必要的傷害。
應該理解的是,前面描述僅是本發(fā)明的示意性描述。本領域技術人員能夠做出各種替代和修改而不偏離本發(fā)明的范圍。因此,本發(fā)明旨在包含所有這些替代、修改和變體。參照附圖所述的實施方式僅呈現以示范本發(fā)明的某些實例。與上文和/或所附權利要求書中描述的那些非實質性不同的其他元件、步驟、方法和技術也旨在包括在本發(fā)明的范圍內。