TW202496B - - Google Patents
Download PDFInfo
- Publication number
- TW202496B TW202496B TW080100256A TW80100256A TW202496B TW 202496 B TW202496 B TW 202496B TW 080100256 A TW080100256 A TW 080100256A TW 80100256 A TW80100256 A TW 80100256A TW 202496 B TW202496 B TW 202496B
- Authority
- TW
- Taiwan
- Prior art keywords
- data
- signal
- amplitude
- phase
- noise
- Prior art date
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52077—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging with means for elimination of unwanted signals, e.g. noise or interference
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
- G01S7/52025—Details of receivers for pulse systems
- G01S7/52026—Extracting wanted echo signals
- G01S7/52028—Extracting wanted echo signals using digital techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52046—Techniques for image enhancement involving transmitter or receiver
- G01S7/52047—Techniques for image enhancement involving transmitter or receiver for elimination of side lobes or of grating lobes; for increasing resolving power
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52046—Techniques for image enhancement involving transmitter or receiver
- G01S7/52049—Techniques for image enhancement involving transmitter or receiver using correction of medium-induced phase aberration
Description
^02496 A6 B6 經濟邾中央樣準局負工消費合作社印製 五、發明説明(1 ) 技術領域 本發明大體上關於超音成像裝置與方法。更明確地說 ,本發明闊於回波圖形裝置與方法,產生固定及移動靶目 標之影像,諸如通常在生物醫學應用中被掃描者。 發明背景 超音成像機常用於各種診斷和治療醫學以及臨床程序 ,例如心臓血管疾病,婦科及産科應用,腫瘤研究及肺部 疾病。初期妊娠診斷、胎兒成長及出生為一種最廣泛使用 超音成像設備之情形。醫事人員靠産生之影像來觀察並研 究早期懷孕和胎兒後續的生長、成長與移動。此種超音監 控亦可偵測出異常。婦産科學院現都推薦婦女接受例行産 前超音波評估,因此,超音成像變為産前撿査之主要工作 〇 然而,超音成像機也有其限制,並非全無危險。目前 之超音波裝置需要直接傳送高頻、高功率、音波信號,通 常在3—7MHz的範圍。産科中,超音波直接進入母體 腹部子宮内,或以探針插入陰道,以産生具有充分解析度 與清楚之影像以便適當評估。然而,對曝露於超音波之組 織的實驗研究指出,長久曝露於此種高能波會傷害胎兒或 母體組織。 標準的習知超音回波圖形成像組合顯示於圖1。探針 1 00載有多數値成已知幾何陣列的轉換器1 0 1。波/ 脈衝産生器102驅動轉換器10 1,同時將一超音倍號 .張尺度適用中國國家標準(CNS)甲4規格(2〗0 X 297公釐) 81.9.25,000 -3 " (請先閲讀背面之注意事項再填寫本頁) -装_ 訂· .绵. 03488 A6 B6 經濟部中央標準局員工消費合作社印製 五、發明説明(2 ) 傳向要成像的靶(未圖示)。靶的回波信號由轉換器 1 Ο 1所接收,並經由類比對數位變換器1 〇 4而通信到 諸如電腦1 0 3的信號處理裝置。傳輸的信號亦直接從波 /脈衝産生器102並經由變換器104,或通常從轉換 器1 ◦ 1而通信到信號處理裝置。處理傳輸的超音信號和 接收的回波信號通常牽涉到各轉換器之傳輸和接收信號的 自動相關,如自動相關器盒1 ◦ 6所示,並可包含各自動 相關之資料的交互相關(資料用於彼此的自動相關),如 交互相關器盒1 0 7所示。然後,自動相關和交互相關的 資料都儲存於資料矩陣1 〇 8。利用此種陣列,使三維影 像成為可行。回波信號到達的時間包含兩値分量,亦即範 圍與方位位置。 自動和交互相關所産生的資料在無線電天文學中稱為 複數可見度函數,由''可見度振幅〃資料,A i j'及 ''可 見度相位〃資料,炎i j τ組成。參見例如1 9 8 6年John Wiley & Sons 出販之Thompson 等人所著的'、Interferometry and Synthesis in Radio Astronomy "第 1 4 頁 。雖 然這些資料在時域,但傳統上是以頻域中的等效變數來描 述。本文以一般習慣來描述頻域,時域的傅立葉轉換為頻 域的傅立葉轉換。如超音成像領域中所公知,可見度振幅 度資料遠比可見度相位資料可靠,後者受雜訊(諸如複數 旁瓣、条統相位雜訊,及校準問題)所駸重破壞。校準問 題,尤其是因超音波在非均勻介質内傳播之相差所引起者 ,導致所諝的*相位像差",通常視為是良好回波圖形成 .張尺度適用中國國家橒準(CNS)甲4规格(210 X 297公犛) 81.9.25,000 —^ — (請先閲讀背面之注意事項再填寫本頁) 丨裝. 訂. 線- A6 B6 024S8 五、發明説明(3 ) (請先閲讀背面之注意事項再填寫本頁) 像最大的阻礙。 對相位像差問題提出兩種措施。一是不用或摒除相位 資料,這是因為可造成很大雜訊,其使用無法提高信號對 雜訊比或産生之影像的解析度及動態範圍。另一措施是提 供試圖校正相位像差之条統,例如利用濾波器、時間延遲 或其他近似法。 當摒除相位資料時,在習知超音波裝置中藉由測繪可 見度振幅資料來成像,如測繪器盒1 1 1所示。測繪器 1 1 1將資料與有關轉換器位置的位置資訊相结合,以合 成對應於探針轉換器陣列的孔。一旦測繪了影像,則靶的 影像便可顯示於輸出裝置1 1 2,最好是視頻顯示端。 不幸地,可見度振幅資料仍含有可觀的雜訊。設備校 準和旁瓣效應造成重大雜訊,所産生的影像遠低於最佳之 解析度和動態範圍。 經濟部中央標準局貝工消費合作社印製 當相位資料用於超音成像時,最常用來減少相位像差 的其中一種措施是使用根據交互相關技術之相位像差的適 應降低。給Has si er等人的美國專利48 1 76 14號以 及給O’Donnell的4 8 3 5 6 8 9號都利用適應降低。多 値轉換器陣列的信號交互相關用來使時間延遲校正插在每 個轉換器,以降低相位像差。此措施令人聯想到光學中之 ''橡皮鏡"措施。然而,此種設計僅是企圖得到真實的無 雜訊可見度相位。其結果可隨著掃描、時間、器官、患者 不同而變化,因此重復性及臨床效果有限。亦請參見,例 如,O’Donnell等人的相關技術論文'' Phase Aberration ,張尺度適用中國國家標準(CNS)甲4规格(2]0Xf297公龙) 81.9.25.000 -5 - 1;〇2490 A6 B6 五、發明説明(4 ) (請先閲讀背面之注意事嘹再場寫本頁)
Measurements in Medical Ultrasound 、、Human Studies ",Ultrasonic Imaging,第1 0册1 一 1 1頁( 1 9 8 8) ; O’Donnell 等人的 '、Aberration Correct ion without the Need for a Beacon Signal ” IEEE U-Itrasonics Symposium* 833-837頁(1988); O’Donnell等 人、、Phase-Aberration Correction Using Signals from Point Ref lectors and Diffuse Scatterers: Mea s urem-ent 〃, IEEE Transact ions on Ultrasonics, F reeoe 1 -* ectrics and Frequency Control, 3 5册 6期7 6 8 — 77 4 頁(1 9 8 8)。也參見 Hayakawa 的 '、Multifrequency Echo scopy for Quantitative Acoustical Characterization of Living Tissues *·、J· Acoustical So- cietyofAmerica, 69 册(6) , 1838— 184〇 頁(1 9 8 1 ),其中開發了人體組織中之衰減僳數的近 似。 經濟部中央標準局員工消費合作社印製 在 Somer 等人的 '、Real-Time Improvement of both Lateral and Range Resolution by Optical Signal Processing. Ultrasonics Symposium Proceedings, 1〇〇2—1〇〇5頁(1977),描述了超音影像增 強處理,其中在交互相關處理中使用相干光學濾波來增進 橫向和軸向解析度。這是達到近似相位像差校正的光學措 施。 頒給Luthra等人的美國專利號碼4 6 0 4 6 9 7號和 4553437號中,由在多數値頻率之轉換器陣列之振 81.9.25,000 本紙張尺度適用中西因豕揉準(CNS)甲4規格(2】0 X 297公楚) -6 - A6 B6 ^02490 五、發明説明(5 ) 幅和相位的向量相加産生了混合影像。整値影像由部份影 像相加而成。Mat sinnoto的美國專利4 5 8 6 1 3 5號中 ,利用相位資料來應用旁瓣由合成孔徑技術來提供重建的 全相資料組。 自動與交互相關亦用於Galbraith的美國專利 4 3 9 7 ◦ 0 6號,以決定地震學中之雜訊降低的數位時 域濾波參數。 加強信號對雜訊比的另一常用措施在翳學應用中使用 相當局的功率及高頻超音波。增加功率傾向於增加相對於 雜訊的倍號強度。至於增加頻率以降低雜訊,超音波束的 擴張與波束頻率成反比。因此,在生物醫學領域,常需高 頻超音成像裝置,試圖減少波束擴張並增大影像解析度。 然而,高頻波束穿過身髏組織時大為衰減,此亦減少信號 對雜訊比。事實上,在身體組織深度超過10公分以上之 超音影像的解析度很差,具有有限的臨床價值。在較大深 度,低頻波束衰減較少。但當頻率下降時,波束擴張且解 析度減少。因此,常使用頻率折衷以使特定應用的解析度 最佳,但一般使用皆較最佳頻率為高,犧牲成像之深度及 清晰。 因此,在所得影像的解析度品質,目前的超音回波圖 形或成像裝置受限。,影像通常非常模糊且有陰影,需要 準確影像以便偵側之組織細部常無法觀察到。因此,使用 目前的超音波設備需要相當的經驗與技巧,即使有此種經 驗和技巧,從回波圖形影像可獲得之資訊仍非常有限,需 &張尺度適用中S困家樣準(CNS)甲4規格(210 X 297公釐) 81.9.25,000 -7 - (請先閲讀背面之注意事項再填寫本頁) 丨裝· 訂· -—線· 經濟部中央標準局8工消費合作社印« 02498 A6 B6 經濟部中央標準局员工消費合作社印製 五、發明説明(6 ) 要主觀判斷。 結果,即使是高頻率、高功率,目前最常用於醫療的 超音波裝置在患者身髏與轉換器組只相距幾公分的徑向距 離,也僅能産生模糊的組織靶的影像。 另做其它嚐試來加強超音波影像的清晰和解析度,但 只有有限的成功。SaSak i的美國專利4 4 7 8 0 8 5號中 ,在陣列中改變超音轉換器之厚度以減少波束擴張。O’D-〇nne- 11的美國專利4 4 7〇3 ◦ 5號利用環狀陣列的超 音轉換器及時間延遲的脈衝模擬在近場具有尖銳聚焦之號 角叭轉換器,利用增進聚焦的此種糸統可達到高達2 0公 分,但在3MHz達成成像。使用可變的接收信號增益, 試圖降低影像中的旁瓣雜訊,在0 ’ Donne 1 1專利中達成深 度之聚焦的改良。 Coursant的美國專利467798 1號中,利用超音 轉換器的極化待性,設法改善超音回波圖形影像聚焦。此 措施的缺點是缺少起始極化資訊及總強度。對超音波影像 的動態範圍與解析度改善不大。 Brisken等人的美國專利44427 1 5號亦利用可 變頻率超音波掃描,音調變化用於Yano的美國專利 4 6 6 4 1 2 2號的裝置。都卜勒偏移亦用來偵測掃描靶 之蓮動,如Seo的美國專利4 5 0 9 5 2 5號所示。 最後,在論文的摘要,我建議在水下聲音成像中,線 性和非線性的成像技術可用來幫助恢復相位可觀察量,以 增加動態範圍和影像準確度。我也建議可應用諸如光學及 張尺茂適用中國國家櫺準(CNS)甲4規格(210 X 297公釐) 81.9.25,000 -8 - (請先閲讀背面之注意事項再蜞寫本頁) |裝· 訂· 線· A6 1,02486 _ B6 五、發明説明(7 ) 無線電天文學之成像方法的技術。C〇hen的'' Phase Recovery and Calibration with Underwater Acoustic Arrays ^ , J . Acoustical Society of America,補 1, 82卷574 — 575頁(1 987)。如何應用及其醫 療成像之合適性的技術並未在摘要中述及。 雖然以習知影像加強技術達成適度改善,但未獲致産 生具有良好品質之影像的超音回波圖形裝置。結合此種信 號的高位準雜訊使持缠、準確的成像不可能。所以,需要 能産生具有高動態範圍之高解析度成像的安全回波圖形裝 置。 「發明概要」 所以,本發明之主要目的是提供一種回波圖形裝置和 方法,在不均勻介質中使用超音波信號來産生固定和移動 靶目標物之高解析度影像。 本發明另一目的是提供一種回波圖形裝置和方法,能 産生具有高解析度和較目前回波圖形裝置所産生者高二至 十倍之靶位置準確度的影像。 本發明另一目的是提供一種超音回波圖形裝置及方法 ,具有較大動態範圍和較大信號對雜訊比。 本發明另一目的是提供一種用於臨床與診斷研究而産 生身髏組織之高解析度影像的回波圖形裝置與方法。 本發明再一目的是提供一種診斷回波圖形裝置和方法 ,能在患者體内較深處産生高解析度影像。 L張尺度適用中國國家櫺準(CNS)甲4規格(2】0 X 297公釐) -9 - (請先閲讀背面之注意事項再塡寫本頁) 丨裝. 訂· 丨線. 經濟部t央標準局员工消費合作社印製 81.9.25,000 L03496 A6 B6 烴濟部中央標準局S工消費合作社印製 五、發明説明(8 ) 本發明又一目的是提供一種裝置和方法,以産生回波 圖形影像,充分接近及時顯示以用於手術中之影像。 本發明又一目的是提供一種裝置和方法,以産生在較 低頻率具有較高解析度之回波圖形影像。 本發明再一目的是提供一種回波圖形成像裝置與方法 ,加強移動靶之影像解析度。 本發明另一目的是提供一種超音成像裝置和方法,可 用於目前的超音轉換器陣列及相關器,以加強所産生之影 像的解析度和動態範圍。 為達到根據本發明之目的的上述及其它目標,提供産 生固定和移動靶目標之高解析度影像之回波圖形裝置和方 法。該裝置和方法包括一種組合,將多數値超音信號傳入 組織之選擇區,並接收多數個超音回波信號。信號以相關 器裝置處理,將傳送和接收之超音信號自動及交互相關。 本發明的裝置和處理中,不僅測繪從交互相關所造成之可 見度振幅,也測繪相位資料。因此,可見度振幅和可見度 相位資料用來産生微分相位、閉合相位及閉合振幅度之資 料組,並將這些資料組區分以進行非線性雜訊降低處理。 雜訊降低最好包括三重可見度相位資料之相加及四値可見 度振幅資料之組的比率,及至少一非線性雜訊降低處理。 本發明的一觀點中,由反覆旁瓣減法技術而從資料降低旁 瓣雜訊,而本發明的另一觀點中,由反覆測繪處理來降低 雜訊。也可结合這些雜訊降低處理,進一步增加動態範圍 。此外,一旦依據本發明的處理將資料區分及雜訊降低, .張尺度適用中因國家樣準(⑽甲4規格,1〇297_公笼) 81.9.25,000 (請先閲讀背面之注意事項再堉寫本頁) |裝_ 訂· .線. 02493 A6 B6 經濟部t央標準局KT工消費合作社印製 五、發明説明(9 ) 則隨後收集之資料可被處理然後交互相關,以顯示靶蓮動 。較佳實施例中,該裝置和方法做為臨床診斷工具,以産 生無外傷、低功率之身體組織的高解析度成像。 本發明其它目的和特點可配合圖式而從以下的詳細說 明及申請專利範圍而更明顯,其中: 圖1為習知回波圖形裝置之方塊圖。 圖2為根據本發明一實施例之利用旁瓣雜訊降低處理 以産生高解析度,高動態範圍影像之回波圖形裝置的方塊 圖。 圖3為圖2之条統中的一對轉換器之交互相關輸出的 三維域圖。 圖4為顯示所有交互相關轉換器之振幅度分佈之頻域 圔。 圖5為在時域之圖1之裝置所産生的合成圖。 圖6為轉換器陣列之孔徑之時域中的波束旁瓣型樣圖 〇 圔7A為圖5之圖上的一點減去圖6之波束旁瓣型樣 後之圖5的圖的部分清除式樣。 圖7 B為圖5之圖上的第二點減去波束旁瓣型樣後之 圖5之的圖的進一步清除式樣。 圖8為根據本發明第二實施例之利用混合測繪處理來 産生高解析度、高動態範圍之回波圖形裝置之方塊圖。 圖9為根據本發明所構成之回波圖形裝置的方塊圖, 用於産生顯示靶目標蓮動的影像。 (請先閲讀背面之注意事項再填窝本頁〕 •丨裝· 訂· •線_ 本紙張尺度適用中a S家標準(CNS)甲4规格(2】0 X 297公釐) -11 - 〇024δο Α6 Β6 經濟部中央標準局S工消费合作社印製 五、發明説明(10 ) 「較佳實施例的說明」 如同用於傳統超音回波圖形裝置和方法,本發明之高 解析度回波圖形糸統採用多數値超音信號,從轉換器陣列 在選擇的靶目標(諸如選擇之身體·組織的區域)之方向傳 送。一組N個轉換器接收自靶組織彈回的反射回波信號。 這些回波信號經處理而産生靶區的三維影像圖。本發明中 ,超音信號之處理産生大為增高之解析度和改善的動態範 圍。 依據本發明之圖2的第一實施例中,利用旁瓣雜訊減 法處理。另一實施例中,如圖8所示,利用混合測繪成像 處理。最後,圖9的處理利用順序成像技術來偵測靶移動 Ο 現在參照圖2,顯示依據本發明一實施例之適於施行 資料區隔、設備雜訊降低及旁瓣減法處理之裝置的縱觀。 如圔2所示,為簡化而未顯示圖1的波/脈波産生器及類 比對數比變換器,這是由於此種元件為習知,且不形成本 發明的新穎部分。 「旁瓣減法縱覽」 超音回波圖形裝置10包括傳送及接收超音轉換器列 或組12,含轉換器T;至T 〃。本發明的處理同時自 T ~値轉換器獲得波束響應。在傳輸信號輪廓與接收信 號輪廓之時間的自動相關及轉換器的時間後,宜將轉換器 (請先閲讀背面之注意事項再填寫本頁) .装· 訂. -線_ 本纸張尺度適用中國S家標準(CNS)甲4規格(210 X 297公穿) -12 - 81.9.25,000 -02486 A6 B6 經濟部中央標準局員工消費合作社印製 五、發明説明(11 ) 對之輸出在各種延遲而交互相關且傅立葉轉換,産生各對 轉換器之''可見度振幅〃和、、可見度相位"。可見度振 幅〃與''可見度相位"的表示式定義如下。 雖然廣泛知道使用轉換器陣列,但本發明之裝置和方 法的首要恃點在使用轉換器陣列,這是因為其可産生較使 用單一之大轉換器之轉換器陣列所可捕捉的資料還多。此 外,捕捉之資料可用來對照以消除雜訊,此點單一轉換器 無法辦到。 可見度振幅與相位資料係以傳統方式自轉換器陣列取 得,標示為1 4 ,將各轉換器之傳輸信號的回波信號自動 相關。圖2的較佳形式中,交互相關器16使自動相闊之 信號交互相關,産生各對轉換器之可見度振幅和可見度相 位。. 然而,本發明的雜訊降低裝置和方法適於不含交互相 關器16之現存超音波裝置。因此,若輸入轉換器對的各 U、V位置,則將自動相關器1 4所産生的時域影像傅立 葉反轉,也可得到可見度振幅和可見度相位。可從自動相 開的輸出,合成可見度振幅度與可見度相位。此措施的缺 點在於不容許微分相位資料組的産生,但由形成、x閉合相 位〃、''閉合振幅〃、旁瓣減法混合測繪而達到雜訊降低 Ο 因此,本文所用之 ''相關"須包括單獨之自動相關及 先自動相關再交互相關之含意。 由於增加用於測繪的資料有利,且交互相闊器容許微 .張尺度適用中國S家標準(CNS)甲4規格(210 >$497公51· ) 81.9.25,000 -13 - (請先閱讀背面之注意事項再填寫本頁) .装. 訂· .線. 20248b A6 B6 五、發明説明(i2 ) (請先閲讀背面之注意事項再填寫本頁) 分相位資料組的産生,所以本發明的較佳裝置與方法包括 交互相關器1 6 ,使自動相關轉換器信號交互相關。時域 交互相關資料表示成各自的空間頻率,然後與空間頻率一 起儲存於資料矩陣2 2 ,否則描述為頻域中的U、V點對 應於轉換器對之相對位置。 在本發明之改良處理和裝置中,可見度振幅度與可見 度相位資料分別在矩陣2 6 a與2 6 b進一步區分,以容 許雜訊降低並産生測繪的更多有用資料。第一,在矩陣 2 6 c得到在不同延遲之可見度相位之差異比較的''微分 相位〃。第二,進一步區分之完成俗相加三個可見度相位 之時域組及四個可見度振幅之比值組,産生稱為a閉合振 幅〃(在矩陣2 6 d )及'''閉合相位"(在矩陣2 6 e ) 的新資料組。也産生閉合相位和閉合振幅的對應頻域資料 。此種相加與比值形成以閉合相位與閉合振幅資料大為降 低設備校準雜訊,藉以提供兩組具有降低之雜訊位準之新 資料。 經濟部中央櫺準局貝工消費合作杜印製 在圖2之處理中,所有五組頻域資料(26a-26e)用來産生在時域之影像圖(圖5)。可見度振幅 度資料一如傳統有些雜訊,但通常相當準確。可見度相位 在習知糸統中的用途有限,但使用對資料庫引起重大雜訊 。徹分相位資料一樣有可觀的雜訊,以前並未用於超音成 像。然而,若可降低雜訊,則可見度和微分相位資料含有 效資料點。閉合振幅與閉合相位資 '料由相加及比值來降低 雜訊,使資料庫增加更多有效資料點。 81.9.25,000 木纸張尺度通用中國國家標準(CNS)甲4规格(2]0X ^97公穿) 一 14 - 202496 A6 B6 經濟部中央標準局员工消費合作社印製 五、發明説明(13 ) 然而,由測繪此五組資料所産生的圖(圖5 )仍有很 多雜訊。儘管如此,仍比習知處理中所發現者具有更多有 效資料,亦即閉合資料及在雜訊可見度相位資料中之有效 資料。這些資料組的傅立葉反轉所産生的圔在無線電天文 學傳統上稱為 ''髒圖"。為了從資料組所得之髒圖除去雜 訊,圖2顯示使用反覆旁瓣減法處理,從資料(主要是相 位資料)清除雜訊的旁瓣分量,以大大加強成像。旁瓣減 法係基於先前用來從無線電天文學和雷達信號去除雜訊的 超解析光點識別。參見例如Hogborn的、、CLEAN as a Pattern Recognition Procedure ” , Indirect Imaging, C-a m b r i d g e U n i v e r s i t y P r e s s,p p . 2 4 7 - 2 5 4 (19 8 3 );以 及 '、Aperture Synthesis with a Non-Regular Distribution of Interferometer Baselines,、、 A s t r onomy A~ strophys ics , Supp 1 . 15,pp.417-426(1974)〇 因此,光點(圔5之60)在髒圖上被識別。(藉由 識別器3 6 ),假設光點顯示靶構造之實際部份。由於已 知任何用來合成孔徑之轉換器陣列12因不完整和/或非 連續孔徑之故,而在有效資料點附近産生複旁瓣雜訊,所 以産生用來産生時域資料之轉換器T之特別陣列12的特 徵旁瓣雜訊型樣(圖6)。波束旁瓣孔徑型樣的合成係由 測繪五資料矩陣(26a_26e)之U — V資料所達成 ,亦卽在所有立體角作傅立葉或類似傅立葉反轉,以産生 時域X、Y圖。時域中之孔徑代表的傅立葉反轉有效産生 具有旁瓣的波束型樣。 .張尺度適用中國國家標準(CNS)甲4规格(2】0 X 297公犛) 81.9.25,000 一 15 _ (靖先閲讀背面之注意事項再填寫本頁) |裝- 訂· 線· L024G6 A6 B6 經濟部中央標準局MC工消費合作社印製 五、發明説明(14 ) 然後在減法盒34,從髒圖(圖5)上之選擇光點 6 0旁的資料減去旁瓣孔徑資料圖。此減法宜在頻域進行 ,消除在時域圖上之資料點6 ◦附近之旁瓣雜訊,藉以清 除圖。注意時域減法亦可。 然後比較器4 1査看部份雜訊降低之,在大部分的情 形,發現仍有太大雜訊。然後對另一光點(例如6 0 a) 掃描清潔圖自第二光點6 0 a附近的資料減去相同旁瓣雜 訊型樣。重複此處理,直至比較器41發現圖(圖7B) 符合比較器的去除雜訊之標準為止,一般為接近熱雜訊限 制,在該點,清潔圖顯示於輸出裝置4 2。 利用閉合振幅度與閉合相位之有效資料點的相加,加 上可見度和微分相位的保持以及相位資料和可見度振幅之 旁瓣雜訊之清潔,導致顯示的影像有較多資料點、較高解 析度及較低雜訊、較大動態範圍。解析度之改善與先前技 藝比較,估計為2至1 0倍,而動態範圍之改良估計為 1 0至1 ◦ ◦倍。 産生回波圖形影像所需的功率與靶和轉換器間之距離 的四次方成正比。因此,若動態範圍、信號對雜訊比可利 用本處理來增加1 ◦倍,則此增加可用來增加成像深度或 減少所需功率,或兩者的組合。因此,若相同動態範圍可 接受,則成像距離可由動態範圍增加之四次方根來增加, 亦即動態範圍增加10倍,則距離增加70至80%,動 態範圍增加1 0 0倍,則距離增加3 %。可減少之功 率與動態範圍的增加成正比。因此,10倍動態範圍的增 ,張尺度適用中國國家櫺準(CNS)甲4规格(210 X 297公蹵) 81.9.25,000 -16 - (請先閲讀背面之注意事項再頊寫本頁) 裝. 訂. .線· 02493 Α6 Β6 絰濟部中央標準局8工消費合作社印製 五、發明説明(15 ) 加致使相同影像之功率的1 〇倍降低。同樣地,1 〇 〇倍 動範圍的增加致使産生相同影像所需之功率的1 〇 〇倍降 低。 「信號相開」 在圖2中,N個轉換器T的陣列或組1 2在靶目標( 未示出),諸如人體組織,之方向傳送信號。轉換器的數 目N越大,則得到之幾何孔徑影像愈大。轉換器1 2不需 排成圖1之線性陣列,通常轉換器在共同框架或探針上排 成二雒陣列,在預定固定之相對位置。 在本發明的較佳形式中,轉換器Τι至各傳送一 寛帶脈衝或都卜勒偏移連續波信號,也接收回到轉換器之 回波或反射信號。所傳送的信號最好是已知的振幅度與相 位,自驅動轉換器之産生器(未圖示)輸入至電腦60, 如輸入6 1所示。自各轉換器傳輸的超音波信號經自由空 間而傳播,直至接觸到組織之選擇區並彈回為止,反射波 回頭朝向轉換器源而成為回波信號,並由同一轉換器來接 收。然而,回波信號也會被相鄰轉換器接收,但所收到之 信號不同。因此,各轉換器的不同投射距離、不同性質及 雜訊持性,及插入不均勻介質之差異層所造成之反射性分 佈差異,都造成接收信號之不同。 須知傳送器和接收器之轉換器的分離組亦可使用於本 發明之方法和裝置中,只要各相位和幅度資料能相關。 値別轉換器在對應於靶之位置、距離、輪廓結構以及 (請先閲讀背面之注意事項再填寫本頁) t r % 本紙張尺度適用中國國家標準(CNS)甲4規格(2】0 Χ^97公釐) 81.9.25,000 -17 - OC4C0 A6 B6 經濟部中央標準局8工消費合作社印« 五、發明説明(16 ) 介質之折射性的時間延遲範圍上接收反射信號。延遲與傳 播介質之音速有關。各轉換器所接收之回波信號之相位和 振幅度在各瞬間來量測。在傳送和接收信號間之振幅和相 位差的相關用來構成選擇組織之三維圖,如圖3所示。 自動相關器,以方塊14表示,用於使各轉換器之傳 送和接收信號在傳播徑向時間延遲(r )上自動相關。傳 播時間延遲由下式而得: r = 2 D X C s 其中為了適當單位,r是以秒表示的時間,D是以公 分表示的距離,Cs是介質傳播音速。自動相關器執行下 列公知的演算(或其等效),以施行各轉換器之傳送和接 收信號間之自動相關。 C 1 ] A 7 (r) = /R (t) R (t-r) dt 其中, R ( t )為在時間t所傳送之傳輸信號, R ( t — r )為回波信號, t為目前之初時, r為時間偏差,定義為在轉換器傳輸和接收間的信號 傳播期間, A 7為N値轉換器之組中的轉換器1至N之自動相關 張尺度適用中國國家標準(CNS)甲4規格(210 X 297公嫠) -18 - 81-9.25,000 (請先閲讀背面之注意事項再塡寫本頁) 訂. .線. 024S; A6 B6 烴濟部中央標準局8工消費合作社印« 五、發明説明(17 ) 信號。 自動相關之結果決定各轉換器之時間(ί)的偏差期 間中之傳送和接收之脈波信號間其振幅度頻率共同之處。 交互相關器步驟,如圖2的方塊1 6所示,用於使成 對之自動相關信號交互相關。交互相關器1 6執行下列演 算(或其等效),以施行兩自動相關信號A ,和A 2之間 的交互相開: [2 ] A i A 2 ( r ) = / A i ( t ) A 2 ( t - r ') d t 其中, r ’為轉換器對之間之傳播的幾何時間延遲。所以, r ·由靶源横向位置所造成之到達時間差異所定義,由下 式而得: r ' = ( D i 2/C s ) ( co s Θ icos0 i+cos0 i s i η φ i -cos θ 2 cos<f> 2~cos0 2sin0 2) 其中,D u等於轉換器T ;與T 2間之距離,C s為 傳播音速,(θ< , )對應於描繪靶位置的二維方向 餘弦角,分別從轉換器T ;和T 2所見。 每一指定的r ·交互相關器施行在N値自動相關信號 間之可能不同之交互相關的總數。例如,若N等於5,則 施行下列交互相闊: (請先閲讀背面之注意事項再填寫本頁) i τ
H t E3a i.«^(CNS) Ψ4*5.^(21_0 X 297_α« ) 81.9.25,000 A6 B6 r0:4S〇 五、發明説明(18 ) 1 X 2 2 X 3 3X4 4X5 1 X 3 2 X 4 3X5 1 X 4 2 X 5 1 X 5 各數字(1 一 5)在N = 5的例子中代表來自轉換器 組1 2的自動相關信號。 交互相關之目的是決定振幅和相位以及由轉換器組 1 2所界定之整値幾何孔徑的密度分布。交互相關1 6決 定關於靶之交互相關轉換器對之間的相對物理位置差所引 起之相位偏差,該差異因轉換器組12的形態而固定。 參照圖3,顯示交互相關轉換器之選擇對的三維時域 圔。三維圔中,其時間偏差(t)對應於X軸,r對應於 Y軸,振幅峰值5 2繪於Z軸。各校正線5 0對應於指定 延遲偏差之對之間的交互相關的輸出。在特定座標(X, Υ, Z)之峰值52的存在對應於轉換器對之靶目標内之 構造的識別。幅峰52亦伴有無可避免之旁瓣雜訊峰值5 4 ,若顯示成為合成孔徑的部分圖,則其存在破壞了圖3 之解析度。當顯示成為光域掃描型之顯示時,則對於本文 所討論的各種圖,也可顯示灰度或假色圖素為基礎的顯示 0 一種似傅立葉的轉換,定義為拉普拉斯轉換,將函數 轉變為頻域(例如,線性函數之傅立葉轉換或傅立葉轉換 本纸張尺度適用中國因家樣準(CNS)甲4規格<210 X 297公釐) -20 - (請先閲讀背面之注意事項再填寫本頁) t 丁 經濟部中央樣準扃8工消费合作社印¾ 81.9.25,000 0:4的 A6 B6 經濟部中央櫺準局貝工消費合作社印« 五、發明説明(19 ) 之非線性函數之頻域轉換),如方塊1 8所示,將圖轉變 成頻率,來産生圖3之時域圖的功率密度函數(PDF) 。此函數定義為: AJ2-(U,V)-/AiA2(r ')ei2Wudz ' 變換器18傅立葉轉換各單獨之交互相關信號,將整 値陣列的轉換器之時域圖變換為頻域。變換器1 8執行傅 立葉轉換以産生空間頻率資料(U - V),用來合成波束 與旁瓣m應。 波束與旁瓣遒應傜轉換器陣列12之特定形態所固有 的。各可見度振幅和相位分配一位置(U , V ),對應於 陣列1 2所産生之幾何孔徑中的位置。 波束與旁瓣繼應之決定像以分配二維位置予各轉換器 ,並執行此延遲陣列之傅立葉轉換,達到二維空間頻譜, 因而各轉換器對對應至一(U, V)點。此種二維空間頻 率資料的集合在本文稱為'"u_v 〃平面。u_v平面之 傅立葉反轉産生波束與旁瓣響應。 各(u, v)點被分配其對應的轉換器對交互相關振 幅度與相位。在傅立葉轉換頻域中之u - v平面的代表産 生資料,然後再轉換而産生 ''髒波束〃。交互相關料資料 的傅立葉轉換頻域産生資訊,可再經傅立葉反轉而産生、、 勝圖〃 〇 為方便起見,下表綜述以標準孔徑合成方法形成影像 (請先閲讀背面之注意事項再填寫本頁) -—装· 訂· 線· 本纸張尺度適用中國國家標準(CNS)甲4規格(210 ;^f97公笼) 81.9.25,000 ^,02498 A6 B6 經濟部中央標準局8工消費合作社印褽 五、發明説明(20 ) 之資料的主要定義域。 資料 階段1 階段2 階段3 轉換器位置位置差異 U-V平面 髂波 轉換器對 交互相關 可見度振幅,相位 髒圖 之輸出 階段2為頻域量。 圖4中,顯示所有交互相關之轉換器對之交互相關振 幅的頻域分布圖。可見度振幅度資料以各別空間頻率繪於 U-V平面上,成為平滑正弦狀之頻率波形5 8。如箭頭 55所示的振幅顯示構造的存在,而圖3之峰值52出現 在圖4中成為其中一條線5 8上的一點。 可見度振幅與相位在各種空間頻率儲存於圖2之矩陣 2 4及U — V資料矩陣2 5。可見度振幅和相位資料及相 關空間頻率傜以複數符號(i , j )表示,為轉換器陣列 元件之函數。可見度振幅A 及可見度相位0 i“資料 是資料的頻域表示。 注意,本干涉量測方法像作用於近場資料,所以使用 近場校正於相闋器,為公知之干涉量測學。 圖1中,可見度振幅資料A 通常是用於習知回波 圖形成像測繪的分量。 在本發明之改良方法中,不僅使用可見度相位資料 張尺度適用中國國家樣準(CNS)甲4規格(210 X 297公笼) -22 — 81.9.25,000 (請先閲讀背面之注意事項再瑱寫本頁) 丨装· 訂. 丨線. 02498 A6 B6 經濟部中央標準局貝工消費合作社印製 五、發明説明(21 ) 炎,而且産生並使用另外三値資料組。 最初,可見度相位的差異比較,如先前交互相關在不 同增量時間延遲r所測定者,會産生'''微分相位〃資料 0 ί ί"。這些資料係産生自吵雜的可見度相位,所以同樣 相當吵雜訊,但包含有效資料點。微分相位儲存於矩陣 2 6 C 〇 在無線電天文學的領域中,中曾証明轉換器陣列之任 何三値可見度相位的向量相加會抵銷設備校準雜訊。而且 很重要的,由傳播介質中之不均勻所引起的任何相差(> 相位像差〃)將被抵消。非代表靶結構和/或位置之相位 項所污染的可見度相位,由新型的相位所替換,該新型相 位係基於三個轉換器而非兩個。 新相位資料由三重可見度相位組成,稱為N'閉合相位 "資料。可由實際可見度相位之各分童來表示: Φ w i j m Φ i j +( θ i ~ Θ j ) + ( 7 i~ Ί j ) + ( β i~ β J ) 其中,<6 〇 =靶的真實可見度相位, θ „ =轉換器Τ η的内部相位, 7 η =對應於轉換器之位置不確定之相位項, 冷„=轉換器Τη量測之傳播介質中之不均勻所引起 的相位偏差。 因此,閉合相位為: 張尺度適用中國國家櫺準(CNS)甲4规格(210 Χ^97公楚) -23 - 81.9.25,000 (請先閲讀背面之注意事項再填寫本頁) -裝- 訂· 線- 202498 A6 B6 經 濟 部 中 央 標 準 局 Ά 工 消 費 合 作 社 印 « 五、發明説明(22 ) Φ w i j k = Φ i j + φ j k + φ hi + ( θ i-e.j) + (e j-e *) + ( θ h-e i) + (7 i-7 j ) + ( 7 j - y h) + ( y k- Ί i ) + ( β t-β j ) + ( β j-β 〇 + (/3 ,-/3 n =Φ i J + φ j h + <f> k i 很清楚地看到污染項(特別是産生相位像差之/3 „項 )消失在中。因此,利用閉合相位是一種確實消除相位像 的方法。因此,陣列中之三個可見度位的所有可能組合以 向量相加,並將和儲存於矩陣2 6 e作為、、閉合¥位" Φ W i j h 〇 同樣地,已知的無線電天文學中,若作四個可見度振 幅度的比,則消除可見度振幅資料中之設備校準誤差。此 可由轉換器k、 1、m、η之閉合振幅度Α υμ表示,其 中A <,為轉換器對丨、>]之量測之可見度振。六〇下式 而得: A i y = G i G J ^ i j 其中,Gi與G,為具有不確定之轉換器i和j的增 益因數,而「〇為轉換器對i j之真實可見度振幅度。因 此, (請先閲讀背面之注意事項再填寫本頁) 裝- 訂. 丨線· 本纸張夂度通用中國國家標準(CNS)甲4规格<2巧1^97_公梦) 81.9.25.000 A,0249b A6 B6 绶濟部中央標準局貝工消費合作社印製 五、發明説明(23 ) GfcGlGmGn^fct ^mn A h l m η ~ "' G * G m G £ G π P k m l n 此式表示可見度振幅度之比,其中消除載有雜訊之增 益因數。簡言之,閉合振幅A 有轉換器增益中之不 確定的雜訊。因此,陣列中之四個可見度振幅的所有可能 組合儲存於資料矩陣2 6 d ,作為 ''閉合振幅"A " t* t ο 在閉合相位中及非單位比值之閉合振幅的任何非零值 代表指示結構存在於靶之資料。因此,這些資料可用於影 像測繪,當作以前超音波成像未用過之新資料點。此外, 由於向量相加及比值消除了設備校準雜訊,所以振幅度資 料相當可靠,且相位資料戯刻性改善,尤其是與可度見和 微分相位資料比較。 相位相加及振幅比值中,利用本方法可消除一轉換器 相對於其他之雜訊。閉合振幅和閉合相位為新一組的資料 可觀測值,與可見度振幅度和可見度相位資料分離。因此 ,以本方法所包括之閉合振幅和閉合相位資料,仍有其他 資料可用來處理信號,以增加信號對雜訊比。 因此,在相關後,本發明之處理的下一步驟是將可見 度振幅和可見相位予以向量相加、比值形成及分隔可見度 振幅資料可儲存於矩陣2 6 a ,可見度相位資料存於矩陣 26b,微分相位資料在於矩陣26c,閉合振幅資料存 (請先閲讀背面之注意事項再填寫本頁) 丨裝. 訂. •線. 本紙張尺度適用中國困家標準(⑽甲4規格,xf_公楚) 81.9,25,000 ^02496 M A6 B6 五、發明説明(24 ) 於矩陣2 6 d ,閉合相位資料存於矩陣2 6 e。 表示可見度和閉合振幅及相位和微分相位之波束與旁 瓣鎏應的U_V儲存於對應之子矩陣中。 因此,本發明之處理中之資料的向量相加産生額外有 用之資料。區分步驟本身為新穎步驟,提供具有顯著之雜 訊和相位像差之降低的資料組。 資料區分亦提供一資料組,獨無受相位像差效果之影 響,因為在中間介質之折射效應所造成之相位的相對差異 會在閉合相位相加中加以消除,所以造成相位像差之分量 不會出現於閉合相位。因此,利用測繪中的閉合相位,産 生不受相位像差干擾之影像。
「藉由旁瓣減法的雜訊降低J 雖然在傳統方法中大量增加矩陣2 6 a — 2 6 e的資 料,但仍含有極大雜訊,尤其是在相位資料。然而,可僅 測繪資料而無需另外處理,可獲得改善之解析度與動能範 圍。若要測繪可見度振幅度、閉合振幅及閉合相位,則只 要將閉合資料點加至傳統所用的可見度振幅即可達成較大 影像解析度。 然而,但在本發明中,使用另外之雜訊降低處理,以 加強所有資料,尤其是吵雜的相位資料。如圖2所示,此 處理為合成孔徑中固有之複旁瓣雜訊減法的形式。 已知當多數轉換器用於轉換器陣列時,合成轉換器的 波束或孔徑比填充(單一轉換器)孔徑吵雜。因此,當使 (請先閲讀背面之注意事項再填寫本頁) τ % 經濟部中央橾準局R工消費合作社印Μ衣 本紙張尺度適用中國困家標準(CNS)甲4规格(2】_0备^9=公穿) 81.9.25,000 Α6 Β6 經濟部中央標準局負工消費合作社印製 五、發明説明(25 ) 用合成轉換器波束(通常在超音波成像中)時,在此種轉 換器成為合成孔徑中的固定點時,來自各轉換器之相對相 位與振幅資料中總是有一定程度之不確定。雖然合成超音 波束或孔徑具有較填充孔徑更大的波束持性適用,但波束 雜訊對信號比的增加便減弱了此種增強之適用。 因此,在信號處理之下値步驟是將合成孔徑區域之a 髒"圖予以合成。髒圖傜定義為可見度振幅和相位的傅立 葉反轉,也代表閉合振幅和閉合相位之反轉與模式。例如 ,在僅傅立葉反轉可見度振幅與相位之簡單情形下,髒圖 I (X,y)由下式而得(對於各轉換器丁〃) ·· K I(x,y)= Σ V(U^,VA)exp{-i2^ (U^X+V^Y)} N=1 其中V (U〃,V〃)係指定之(u, v)對的可見 度函數(可見度振幅和可見度相位)。 圔2之測繪器2 7所示的測繪,使用子矩陣2 6 a — 26e中之頻域資料,標準髒圖6 1顯示於圖5。將矩陣 26a_26e中之資料映到合成孔徑的圖61傜公知技 藝,不在此複述。 圔6 1含多重峰值60、6〇a,與靶目標之光亮結 構對應。注意此時圖6 1仍含高位準雜訊,如峰值6 2與 6 2 a所示,由孔徑形態固定之旁瓣鬱應所産生。 (請先閲讀背面之注意事項再填寫本頁) —裝· 訂. 線· 衣紙張尺度適用中國國家櫺準(CNS)甲4規格(2】0 X 297公梦) -27 - 81.9.25,000 ^〇24δ'〇 A6 B6 經濟部中央櫺準局霣工消費合作社印« 五、發明説明(26 ) 旁瓣減處理使用矩陣26a—26e中之U, V資料 ,模擬轉換器陣列對單點源之合成波束髒波束)辔應 。髒波束是陣列之合成波束或點饗應。否則陣列的轉移函 數,由下式而得: K Ρ(Χ,Υ)= Σ βχΡ{-ΐ2π (U^X+V^Y)} N-l 注意I (X, Y)代表髒波束與可見度函數之褶積。 波束合成器40將整値U,V資料轉換或測繪成在頻 域的波束功率型樣。若靶是點結構或源,則可自U、V資 料測繪的波束功率型樣有效代表固有産生的旁瓣。因此, 在探針上特定轉換器陣列12之特有的幾何形狀總會在點 源各側産生相同之旁瓣雜訊。若重新建構轉換器陣列12 ,則會改變旁瓣雜訊,但對任何固定陣列,在一點附近僅 有一旁瓣雜訊的波束功率型樣。 若測繪此種U、V資料,則圖6表示典型的波束饗應 圖65,有一中央峰值64代表點源和旁瓣68在孔徑幾 何特性之位置和振幅。 為從圖5的圖61去除旁瓣62雜訊本處理用反覆減 法技術。參照圖7 ,顯示在執行反覆旁瓣減法處理後之圔 5之圖的部分清潔式樣。反覆處理包括下列步驟: 1.圖5之圖6 1上的亮度峰值6 ◦以識別器3 6識 末紙張尺廑適用中國國家標準(CNS)甲4规格(210 X 297公釐) 81.9.25,000 (請先閲讀背面之注意事項再填寫本頁) 装- 訂. .線. ^0:49β Α6 Β6 經濟部中央標準局W工消费合作社印製 五、發明説明(27 ) 別。此由尋找大峰值而識別數學性施行於圖6 1之資料。 2. 圖5之圖以似傅立葉反轉轉換至頻域,如轉換器 方塊32所示,以準備於旁瓣減法步驟來除去頻域中的雜 訊。 3. 在旁瓣減法步驟中,點源的U、 V波束饗應資料 自包圍圖5之亮度峰值60的資料減去,如減法器方塊3 4所示。此步驟從圖5的圔去除點6◦附近之旁瓣62, 産生除去旁瓣繼應的較清潔之圖6 7,如圖7A所示。 4 .然後比較器4 1掃描或比較圖7 A之較清潔的圖 67,標準是雜訊位準。輸入裝置43能指定圖67之動 態範圍、信號對雜訊位準的最大範圍必須超過一定量,或 繼缠清潔,直到至少是最亮點60之強度一半之圔5之圖 6 1上的所有光點已清除為止。圖7A的圖67中,從光 點6 0週圍清除旁瓣雜訊6 2對第二高光點6 ◦ a附近之 旁瓣雜訊6 2 a並無什麼效果。 5 .所以,若解析度未達比較器4 1標準,則圖 6 7再反轉至時域,並以識別器3 6識別第二光點6 0 a Ο 6.較清潔之圖67再反轉至頻域,減法器34從包 圍下一最亮點6 0 a的資料減去相同旁瓣型樣,産生圔7 B之圖6 9。 7 .新清理的圖6 9再以比較器比較,若未符合標準 ,則圖回至另一光點選擇及波束旁瓣減法。 8.若新清理的圖之比較産生可接受之結果,則圖7 張尺度適用中困國家標準(CNS)甲4規格(210 X 297公«· ) 81.9.25,000 -29 - (請先閲讀背面之注意事項再填窝本頁) —裝. 訂· 線·
A6 B6 經濟部中央標準局IW工消費合作社印¾ 五、發明説明(28 ) B之清潔圖以似傅立葉轉換器39自頻域轉換成時域,並 顙示於顯示輸出裝置4 2。 利用反覆旁瓣減法,去除相位資料中的大量雜訊,以 及可見度振幅資料和閉合振幅資料中的一些雜訊。結果可 拾取未用於超音波成像之五値源的增加之可用資料,並改 善所有資料之信號對雜訊比。 利用本方法和裝置所顯示的所得影像較習知回波圖形 影像有很大的改善。明確地說,本發明所産生之影像增加 了 1 0到1 0 0倍的動態範圍及2到1 0倍的解析度。 〔藉由混合測繪的雜訊降低」 本發明另一觀點中,以混合測繪方法來減少雜訊。如 圖8所示,方法亦包括微分相位、閉合振幅與閉合相位, 額外有用資料組的産生。然而,圖8之方法與裝置採用混 合測繪,另一使用於無線電天文學的信號處理技術,代替 由旁瓣減法來清潔所有資料。 圖8中,顯示依據本發明第二實施例之從混合測繪處 理來産生高解析度、高動態範圍影像之第二回波圖形裝置 的方塊圖。混合測繪處理器提供了信號處理從轉換器排1 2所得之資料的第二裝置。當可見度相位受糸統雜訊破壞 時,此第二處理特別有用。 除去影像旁瓣雜訊,影像解析度和動態範圍經由混合 映測繪或似自動校準方式加以改善。混合或模式資料用來 代替不要之實際相位資料,並用來産生圖。以此方式,去 (請先閲讀背面之注意事項再蜞窝本頁) .裝· 訂. .線· 衣紙張尺度適用中國國家樣準(CNS)甲4規格(210 X 297公釐) -30 一 81.9.25,000 A6 B6 經濟t中央標準局β工消費合作社印製 五、發明説明<29 ) 除製造雜訊之壊相位資料點,嵌入模式資料以改善所得影 像之解析度與動態範圍。 轉換器組12、自動相關器14、交互相闊器16、 傅立葉變換器1 8、儲存可見度相位和振幅資料及U — V 資料之資料儲存裝置2 2、以及將此種資料分成五値子矩 陣26a — 26e之操作與圖2的第一實施例相同。在此 階段,圖6 1 (圖5 )含高位準雜訊,通常為相位雜訊, 導致不良之解析度和動態範圍。 然而,混合測繪中,基本上假設振幅資料較正確,而 相位資料被破壞。此與傳統上用於圖1之習知超音成像處 理的假設相同。 所以,藉由混合測繪在減少雜訊之第一步驟最好以測 繪器2 7 a在頻域中僅測繪可見度振幅和閉合振幅資料, 産生類似於圖5之圖61的圖,但因無相位資料存在,故 測繪資料較少。 識別器3 6 a再度識別光點,等於圖5之6 0的點。 然後此最亮點用來産生模式相位資料。指定最亮點為零相 位點,並假設孔徑中僅有一亮點,如似傅立葉反轉方塊 5 2所示,將最亮點資料反轉成頻域,産生恃定孔徑中的 一點構造之模式振幅和模式相位資料組,這些模式資料點 得自傳統上用於無線電天文學之亮度分布演算法。參考H-ogbom 的 '、CLEAN as a Pattern Recognition Procedure, "Indirect Imaging, Cambridge Univ· Press , pp 247 -254(1983)0 (請先閲讀背面之注意事項再填寫本頁) -裝. 訂. 線. 衣紙張尺度適用中囷國家標準(CNS)甲4規格(210 X 297公笼) -31 - 81.9.25,000 0:4ϋδ Α6 Β6 經濟部中央樣準局S工消費合作杜印« 五、發明説明(30 ) 模式振幅與相位資料由比較器裝置4 1 a與矩陣 26a — 26e中之頻域資料比較,如括弧53所示。模 式振幅資料與實際振幅資料通常極為吻合,但模式相位資 料則與實際相位資料差異很大。然後資料替換完成,如混 合資料輸入方塊5 1所示,模式相位資料取代在矩陣 26b、 26c與26e的實際相位資料。以比較器 4 1 a作適配之比較,也可將一些模式振幅度資料點替代 矩陣26 a與26 d中之振幅資料。 一當替代模式資料,則産生新振幅為基礎的圖。若孔 徑中僅有一靶點(不會發生於生物醫學應用上),則模式 資料替換唯一亮點的實際資料會結束資料處理。然而,在 幾乎所有的情況下,識別到第二亮點60 a。現在使用第 二亮點60 a,配合第一亮點6 ◦,産生更複雜的兩點模 式資料組。這是由利用換亮度分布演算法之52的似傅立 葉轉換並假設有兩亮點來逹成。再一次比較模式資料並替 換入矩陣26a — 26e。重複此處理並識別第三亮點, 進行三亮點反轉以産生更複雜之楔式資料組。 當模式資料替換時,雜訊在圖上所造成之亮點會減弱 並消失,留下對應於結構之亮點。因此,若每一迴圈加入 一亮點,則由於有些可能是會消失的雜訊,所以不能由識 別1 0値亮點來開始。 資料替換之反覆處理繼續,直到振幅度資料與模式振 幅度資料高度配合(通常介於10至100次重複之間) 為止。適配可由K平方(CHI-squared )餘數的最小化或 (請先閲讀背面之注意事項再填寫本頁) 丨裝. 訂· 線. 本紙張尺度適用中S國家橒準(CNS)甲4規格(210 X 297公釐) -32 - 81.9.25,000 ,〇:498 A6 B6 經濟部中央標準局貝工消費合作社印製 五、發明説明(31 ) 其他標準來界定。當高度配合發出時,大部份的相位資料 和一些振幅資料會由混合或模式資料來替換,相位資料中 之旁瓣雜訊大量減少。 比較器使矩陣2 6 a — 2 6 e中之所有資料被測繪, 如括弧5 0所示,並顯示於輸出裝置42 a,結果一樣是 以混合測繪經由增加資料及降低雜訊來增強解析度。 「組合處理」 可利用旁瓣減法或混合測繪來減少超音信號雜訊。也 可使用兩種處理。較佳形式中,先以旁瓣減法來降低雜訊 ,再以混合測繪來處理清潔之資料。此順序較反序處理為 佳。 「蓮動偵測」 本發明的裝置與方法可用來決定並顯示在靶區之結構 的蓮動。圖9顯示回波圖形信號處理条統,其中可準識別 在回波圖内之任何特性的蓮動或構造變化。本發明的此觀 點中,連绩産生之資料用來偵測蓮動。 圖9的處理中,使傳送與回波信號相關並區分,如配 合圖2所述。其後,資料以旁瓣減法(圖2 )來清潔,或 混合測繪(圖8 )到雜訊減少及資料增加之圖完成準備顯 示的點。 不用圖2或圖8的處理來顯示圖,可見度振幅、可見 度相位、閉合振幅,閉合相位及相差資料以自動相關比較 ,張尺度適用中固西家樣準(CNS)甲4規格(210 X 297公釐) -33 - 81.9.25,〇〇〇 (請先閲讀背面之注意事!?再«寫本頁} -·裝. 訂- .線. ^02490 A6 B6 經濟部中央棣準局W工消費合作社印製 五、發明説明(32 ) 來進一步處理,在第一組資料後之順序時間捕捉相同資料 。因此,圖9顯示在時間=1之矩陣70a,、7〇bi 、70c;、70d;、70e;中之雜訊降低的資料。 除了已從資料除去雜訊外,餘皆與圖2的矩陣26a— 2 6 e全等。 時間=2時,相同的雜訊降低資料存於矩陣7 0 a 2 —70e2,程序繼缅至時間為止。因此可看出圖9 之資料組分為兩種,亦即同時時間(不同延遲)及順序時 間(不同延遲)。因此,有不同亮持性、構造之資料組作 為相位參考,且稍後同型之資料組也具有此種持點之資料 Ο 各同時時間之資料組以配合演算法反轉,如方塊 7 1 7 1 〃示,並用相位差找出準確位置及各資料組 時間的亮持性之間的構造差異。位置資訊以特性間相差编 碼。各時間之圖上之亮特性的相對位置資料由相差決定) 與特特性構和亮度一起存於矩陣7 22 _7 2w。 下個步驟使儲存於矩陣72: -72〃之順序相鄰資 料交互相關,以識別固定特性。,若時間順序資料之交互 相關僳零偏移,則恃性不無移動且固定。不動特性的資料 與其他資料區隔並存於矩陣73: — 73*^ ,以降低電 腦處理時間。 不在零偏移之交互相關的資料代表從某一時間到下一 時間之構造之蓮動。在非零偏移資料中通常可發現兩種蓮 動,亦即真實且要研究的運動,和設備感應之雜訊。雜訊 (請先閲讀背面之注意事項再填寫本頁) .装. 訂· .緯· 木纸張尺度適用中因國家標準(⑽甲4規格⑵f ~ 81.9.25,000 ^02496 A6 B6 經濟部中央標準局β工消費合作社印製 五、發明説明(33 ) 蓮動使影像外形移動或模糊。整値孔徑看似稍有位移。 對照下,構造之運動限於框架内的物體。所以,在標 示''擴展/偏置〃方塊7 3 , — 7 3 〃 -』的上半部,找到 相對於定義之參考資料組之所有順序資料組之擴展偏移函 數。此函數定義各順序之模糊或震動及位置不準確。 一旦識別擴展/偏移後,則使資料再次交互相關,以 改正資料並消除震動。因此,方塊75i — 75〃之資料 係對應於矩陣7 2 7 - 7 2 〃中之資料的蓮動資料,但已 改正而消除震動。 下値步驟使矩陣75 / — 75〃中無震動的資料順序 交互相關,産生方塊7 6 /之時間=1至時間=2之結構 待性的動資料,方塊762之時間=2至時間=3等等, 直至方塊76〃-』之時間=N — 1至時間=N為止。然後 將此蓮動直接顯示於輸出裝置7 8,以濾波器77; — 7 7 #過濾,或如同矩陣73i—73 (''固定特性 ")内所含之固定目標的資料所顯示者。可根據各種標準 來過濾,例如,若固定資料亦通過震動之消除,則其可由 矩陣73』一 73#-/之輸入來過濾掉。 圖9的方法主要採用圖的雜訊降低順序圖像、作比較 、消除設備雜訊、然後依序顯示蓮動和構造變化。 雖然本發明參照一些特定實施例加以說明,尤其關於 生物醫學應用,但這是說明本發明而非限制本發明熟悉技 藝人士就本發明附列之申請專利範圍所定義者作各種之修 改,並不脱離本發明之精神與範畴。 (請先閲讀背面之注意事項再填寫本頁) 裝. 訂· .線. 紙張尺度適用中國园家標準咖甲4规格公货) 81,9.25,000
Claims (1)
- ,0:4δδ 卜年il 補充 B: c· D- 六、申諳專利苑ffl 1. 一種用以産生靶目檫之回波圖形影像的方法,包 含以下步驟:將多數個超音信號傳向該靶目標,·從該靶目 標接收對應之多數個回波信號;使該超音信號及該回波倍 號相關,以産生包含振幅資料的相關資料,足以測繪該_靶 目標及相位資料的影像;測繪該振幅資料,其改良包括以 下步驟: 將該相關資料分成多數個資料組,包含可見度振輻資 料組和可見度相位資料組; 由下列至少其中之一來降低存在於該相闊資料之雜訊 :(i )旁瓣減法處理,(i i )混合測繪處理,(i i i ) 三個可見度相位資料值之可能組合的向量相加,以産生閉 合相位資料組,及(iv)形成四値可見度振幅資料值之可 能組合的比,以産生閉合振幅資料組;及 在該測繪步驟中,在該降低雜訊步驟後測结該振幅資 料及該相位資料。 2. 如申請專利範圍第1項之方法,其中, 經濟部中央標準局貝工消費合作社印製 該降低步驟由三個可見度相位資料值之所有可能組合 的向量相加來達成,.以産生無相位像差之閉合相位資料組 ,及 在該測繪步驟中,测繪該閉合相#資料組。 3. 如申請專利範圍第2項之方法,其中, 該降低步驟另外包含形成四値可見度振幅資料值之所 有可能組合之比的步驟,以産生閉合振幅資料组,及 在該測繪步驟中,測繪'該閉合相位資料組及該閉合振 81. 2. 2,500(2) (請先閱讀背面之注意事項再填窝本頁) 本紙張尺度进用中Η B家撑半(CNS) Μ規格(210x297公釐) -36 - 經濟部中央標準局貝工消费合作社印製 81. 2 . 2.500(H) ^02403 g __D7_ 六、申請專利範面 幅資料组。 4. 如申請專利範圍第3項之方法,另包括以下步驟 • . 在該測繪步驟之前,産生微分相位資料組,及 在該測繪步驟中,測繪該徼分相位資料。 5. 如申請專利範圍第2項之方法,其中 在該向量相加步驟後,由使用該閉合相位資料組來旁 瓣減法處理來降低雜訊。 6. 如申請專利範圍第2項之方法,其中 在該向量相加步驟後,由使用該閉合相位資料组的混 合測繪處理來降低雜訊。 7 .如申請專利範圍第4之方法,其中 在向童相加、形成比值和産生的該步驟後,由使用該 閉合相位資料組、該閉合振幅資料組、該徼分相位資料组 、該可見度振幅資料組及該可見度相位資料組的旁瓣減法 處理來降低雜訊。 8. 如申請專利範圍第4項之方法,其中 在向置相加、形成比值和産生微分相位資料組的該 步驟後,由使利用該可見度振幅資料組、該閉合振幅資料 組、該閉合相位資料组的混合測繪處理來降低雜訊。 9. 如申請專利範圍第8項之方法,其中 在該混合測繪處理時,利用該可見度相位資料組及該 微分相位資料組。 10. 如申請專利範圍第1項之方法,其中, 本紙張尺度进用中® R家標平(CNS) VM規格(2〗0x297·^· j, (請先閱讀背面之注意事項再填宵本百) k. •訂· .缘· ,02496 B7 C7 D7 々'申請專利苑S1 該傳輸步驟傜自超音轉換器陣列傳送超音信號來逹成 ,以合成孔徑以將該靼目標成像;及 (請先閲讀背面之注意事項再填寫本百) 該降低步驟傜以旁瓣減法處理來逹成,其中從該資料 組的選擇區減去該超音轉換器陣列所合成之該孔徑中固有 的旁瓣響應c 1 1.如申請專利範圍第1 0項之方法,其中, 該旁瓣減法處理包含以下步驟: (i )將該多數資料組之振幅及相位資料測繪成時域 之髒圖; (ii)在該髒圖中識別一亮點; (i i i )産生頻域中之該孔徑的旁瓣饗應型樣; (^)將該髒圖轉換成頻域; (v )自頻域中之識別亮點周圍的資料點減去該旁瓣 鎏應型樣,産生清潔之頻率資料; (V i )將該清潔頻率資料轉換到時域中之清潔圔;及 在該旁瓣減法處理後,顯示該清潔圖。 12. 如申請專利範圍第11項之方法,另包括以下 步驟: 濟 部 + 央 標 準 局 Ά X 消 费 合 作 杜 印 製 比較該清潔圖與該髒圖以降低雜訊,及 對該髒映圖中的其他亮點重複該琴瓣減法處理。 13. 如申諳專利範圍第1項之方法,其中, 該降低步驟由混合測繪處理來完成,其中指定資料點 的模式資料取代在區分資料組中的實際資料。 14. 如申請專利範圍第13項之方法,其中, 81. 2. 2,500(H) 本紙Jt尺度適ffl中田困家標iMCNS)下4規格(210x297公" -38 一 ^〇24δδ AT Β7 Cl D7_ 六 '申請專利苑SI 該混合測繪處理由下列步驟來達成·· (i )測繪在時域中的振幅資料,産生第一振幅圔; (i i )識別該第一振幅圖中之亮資料點; (iii )使用亮度分布演算法以反轉換該亮資料點至 頻域,産生該亮資料點之模式振幅及相位資料; (iv)以該模式振幅和相位資料替換在該區分資料组 中的對應之實際振幅與相位資料;及 然後測繪該區分資料組。 15.如申請專利範圍第14項之方法,另包括下列 步驟: (i )在該替換步驟後,自振幅資料組測繪第二振幅 圖; (i i)識別在該第二振幅圖中與該第一亮資料點不同 之第二亮資料點; (iii )使用亮點分布演算法以反轉換該第二亮資料 點至頻域,産生該第二亮資料點之模式振幅及相位資料; 及 (iv)以該模式振幅和相位資料替換在該區分資料組 中的對應之實際振幅與相位資料。 1 6 .如申請專利範圍第1 4項之_方法,另包括下列 步驟: 在該區分步驟中産生可見度振幅,和可見度相位之資 料組,在該降低步驟中産生閉合振幅和閉合相位之資料組 本紙張尺度適扣中田Η家標苹(CNS)TM规格(210x2974^—)% 81. 2 . 2.500(H) -39 - (請先閱讀背面之注意事項再填寫本頁} •故· 經濟部中央樣準局貝工消费合作社印製> ^0245b 7 7 7 7 Λ B c D 濟 部 中 央 標 準 貝 工 消 费 合 作 社 印 製 六、申請專利範a 17.如申諳專利範圍第is項之方法,其中, 測繪振幅資料的該步驟由測繪可見度振幅資料和閉合 振幅資料來完成。 1 8 .如申請專利範圍第1 6項之方法,其.中, 在該區分步驟中産生微分相位的資料組。 19·如申請專利範圍第15項之方法,另包括以下 步驟: 以適配比較來比較模式振幅資料與實際振幅資料; 重複振幅測繪,亮資料點識別、亮資料點的反轉及模 式資料替換步驟,直至該比較步驟産生特定配合為止;及 當該比較步驟找到該振幅資料匹配特定配合時, 測繪所有資料组。 20. 如申請專利範圍第19項之方法,另每括下列 步驟: 顯示所有資料组之圖。 21. 如申請專利範圍第1項之方法,另包括下列步 m * 該超音信號以多數信號組在多數順序時間期間傳輸; 各組信號及對應組的回波信號以該區分和降低步驟來 處理; 顯示測繪振幅資料和相位資料;其中, 該顯示步驟由下列步驟來達成: (i )測繪所有該資料组; (ii)識別各該圖上的多數亮待性,並儲存關於各該 本紙張尺度適;I) t 田困家標iMCNS)〒4規格(210x297/i^_)4Q _ 81. 2_ 2.500(H) (請先閲讀背面之注意窣項再填寫本百) •装. •訂. ’線· 2〇24δδ c? _____D7_ 六、申請專利氣SI 亮特性之位置與亮度的待性資料; (請先閱讀背面之注意事項再填寫本頁) (i i i )以下一組特性資料的該特性資料使毎一特性 資料組的待性資料交互相關; (iv)識別特性資料組之間的資料變化以顯示蓮動; (V )改正資料,以從各特性資料組除去震動; (v i)使順序相鄰改正之資料交互相關,以識別順序 蓮動;及 (vi i )順序顯示改正之交互相關的資料組。 2 2 .如申請專利範圍第2 1項之方法,另包括下列 步驟: 分離第一交互相關步驟所造成之不改愛之資料。 2 3 .如申請專利範圍第2 2項之方法,另包活下列 步驟: 顯示該不改變之資料。 24. 如申請專利範圍第22項之方法,另包括下列 步驟: 利用該不改變之資料來過濾該剩餘資料。 經濟部中央標準局貝工消费合作社印製 W 25. 如申請專利範圍第1 1項之方法,其中, 該超音信號以多數信號組在多數順序時間期間傳輸; 各組信號及對應組的回波信號以寧區分及降低步驟來 處理;其中, 該顯示步驟由下列步驟來達成: (i )測繪所有該資料組; (ii)識別各該圖上之多數亮特性,並儲存關於各該 81. 2 . 2.500(H) 本紙張尺度適JD中B Η家標堆(CNS) TM規格(210χ29Ρϋ ' 41 ^ AT ‘‘0248b B7 C7 _ D7 六'申請專利範® 亮特性之位置和亮度的特性資料; (請先閱锖背面之注意事項再滇寫本頁) (iii )以下一組特性資料的該特性資料使每一特性 資料的特性資料交互相關; (iv) 識別特性資料組之間的資料變化以顯示蓮動; (v) 改正資料,以從各特性資料組除去震動; (vi) 使順序相鄰改正之資料交互相關,以識別順序 蓮動;及 (vi i )順序顯示改正之交互相關的資料組。 26. 如申請專利範圍第20項之方法,其中, 該超音信號以多數信組在多數順序時間期間傳輸; 各组信號及對應組的回波信號以該區分和降低步驟來 處理;及其中, 該顯示步驟由下列步驟來達成: (i)測繪所有該資料組; (Π)識別各該圖上的多數亮特性,並儲存關於各該 亮待性之位置與亮度的特性資料; ,經濟部中央標準局貝工消背合作社印製 (iii )以下一組特性資料的該特性資料使每一特性 資料的特性資料交互相關; (iv)識別特性資料組之間的資料變化以顯示蓮動; (V)改正資料,以從各特性資料組除去震動; (vi)使順序相鄰改正之資料交互相關,以識別順序 蓮動;及 (vii )順序顯示改正之交互相開的資料組。 27. 如申請專利範圍第1項之方法,其中, 本纸張疋度適用中國國家樣準(CNS>甲4規格(210X297公釐> -42 - A7 ” B7 ^〇249b c? _ D7__ 六、申說專利範ffl 該降低步驟包括該旁瓣減法處理和該混合測繪處理兩 者。 28. 如申請專利範圍第27項之方法,其中, 先執行該旁瓣減法處理,隨後由該混合測繪處理來處 理該旁瓣減法處理的輸出。 29. 如申請專利範圍第28項之方法,其中, 該超音信號以多數信號組在多數順序時間期間傳輸; 各組信號及對應組的回波信號由該區分步驟及由該旁 瓣減法處理和該混合測繪處理來處理;及下列步驟: (i )测繪所有該資料組; (ii)識別各該圖上的多數亮特性,並儲存顔於各該 亮特性之位置與亮度的特性資料; (iii )以下一組特性資料的該特性資料使每一特性 資料的特性資料交互相關; (i v)識別待性資料組之間的資料變化以顯示蓮動; (V)改正資料,以從各特性資料組除去震動; 經濟部中央搮準局β工消费合作社印製 (請先閱讀背面之注意事項再填寫本頁} (Vi)使順序相鄰改正之資料交互相關,以識別順序 蓮動;及 (Vi i )順序顯示改正之交互相關的資料組。 30. —種用以産生回波圖形影像之裝置,包括超音 轉換器,聯結以驅動該轉換器産生器裝置,聯結至該轉換 器以接收其回波信號資料且聯結以接收傳送之超音信號資 料的信號處理裝置,該信號處理裝罝相蹰並測繪該超音信 號資料和該回波信號資料之振幅資料,以及聯結至該信號 本紙張尺度適用中Β國家標準(CNS)T4規格(210X297公釐) -43 - (〇24SS A 7 B7 C7 D7 經濟部屮央標準局貝工消費合作社印製 、申請專利範圊 處理裝置以從該信號處理裝置接收測繪資料並顯示該測繪 資料而成為影像的顯示裝置,其改良包含: 該裝置包括多數超音波轉換器裝成陣列並聯結至該産 生器裝置和該信號處理裝置;及 該信號處理裝置包括處理相位資料的裝置,該處理相 位資料的該裝置包括降低該相位資料内的雜訊之裝置,降 低雜訊之該裝置由下列其中之一提供:(i)自該相位資 料減法旁瓣雜訊之裝置,及(Π)替換實際相位和振幅資 料為混合圖資料之裝置;及 在降低雜訊後,該處理裝置進一步測繪幾乎所有之該 振幅資料和該相位資料。 31. 如申請專利範圍第30項之裝置,其中, 該倍號處理裝置包括減去旁瓣雜訊之裝置及替換混合 圖資料之裝置。 32. 如申請專利範圔第3 1項之裝置,其中, 減去混合資料之該裝置接收減去旁瓣雜訊的該裝置之 輸出。 33. 如申請專利範圍第30項之裝置,其中, 該信號處理裝置包括處理多數順序組超音信號之裝置 ,處理各組以降低雜訊,及 該信號處理裝置包括使各組信號之特性資料交互相關 並順序顯示該交互相閬特性資料以顯示靶目標運動的裝置 3 4 . 種用以産生回波圖形影像之裝置,包括裝成 本纸張尺度適用中國國家《準(CNS)甲4規格(210X297公釐) (請先聞讀背面之注意事項再填萁本页) •装. •打. •線 -44 - ,02498 A7 B7 C7 D7 經濟部中央標準局貝工消费合作社印製 六'申請專利範園 陣列的多數超音轉換器,聯結以驅動該轉換器的信號産生 器裝置,聯結以接收該轉換器之回波信號並聯結以接收回 波信號的信號處理裝置,該信號處理裝置回應以使該回波 信號相對該傳輸信號而相關,産生適於影像測繪之振幅信 號,以及聯結至該信號處理裝置以接收該幅度信號及其響 應以産生影像的顯示裝置,該裝置之改良包含: 該倍號處理裝置回應該傳輸信號及該回波信號,使該 回波信號相對於該傳輸信號而相關,産生振幅信號和相位 信號兩者; 該信號處理裝置更包括雜訊降低裝置,回應該相位信 號和該振幅信號,以降低該相位信號和該振幅信號中之雜 訊; 該雜訊降低裝置包括下列至少一種:(i)自該振Hi 信號和該相位信號減去旁瓣雜訊之裝置,(π)替換該振 幅信號和該相位信號以混合圖像信號之裝置,(i i i )向 量相加三相位信號組以産生閉合相位信號之裝置,及(iv )形成四振幅信號組之比以産生閉合振幅信號之裝置; 在進行雜訊降低後,該顯示裝置聯結至該信號處理處 置,以接收在一位置的信號。 35. 如申請專利範圍第34項之裝置,其中, 該雜訊降低裝置包括裝置回應相位信號和該振幅倍號 的裝置,以通信至具有降低之旁瓣雜訊之顯示裝置。 36. 如申請專利範圍第34項之裝置,其中, 該雜訊降低裝置包括從該振幅和該相位信號減去旁瓣 <請先閱讀背面之注意事項再堪荈本頁) ·«.· •訂· .線 本紙張尺度適《中國B家樣準(CNS)<f 4規格(210x297公釐) -45 - 經濟部中央搮準局貝工消费合作社印製 _____D7_ 六、申锖專利範® 雜訊之裝置。 37. 如申請專利範圍第34項之装置,其中, 該雜訊降低裝置包括將該振幅和該相位信號換以混合 圖信號的裝置。 38. 如申請專利範圍第34項之裝置,其中, 該雜訊降低裝置包括形成四振幅信號組之比以産生閉 合振幅信號的裝置。 39. 如申請專利範圍第34項之裝置其中, 該雜訊降低裝置包括向童相加三相位信號組以産生閉 合相位信號的裝置。 40. 如申請專利範圍第39項之裝置,其中, 該信號處理裝置包括將該振幅信號、該相位信號及該 閉合相位信號分成分離信號組的裝置,及 該信號處理裝置係形成來測繪該影像,以回應該分離 信號組。 4 1.如申請專利範圍第40項之裝置,其中, 該雜訊降低裝置亦包括形成四振幅倍號組之比以産生 閉合振幅信號的裝置,及 該區分裝置區分該振幅信號與其餘信號。 42. 如申請專利範圍第4 1項之裝置,其中, 該倍號處理裝置傜以回應該相位信號,産生一組撖分 相位倍號,及 該區分裝置區分該徹分相位信號與其餘信號。 43. 如申請專利範圍第42項之裝置,其中, f請先«I讀背面之注意事项再琪其本页} ·«.· .打. •線 糸紙張尺度適用中B Η家搮準(CNS) T4規格(210X297公釐) -46 - B7 C7 D7 03496 六、申請專利範80 該雜訊降低裝置包括替換混合信^的裝置,聯結於該 區分裝置以接收其信號及饗應,以混合測繪信號替換該分 離信號組中之信號。 44.如申請專利範圍第42項之裝置,其中, 該雜訊降低裝置包括減去旁瓣雜訊的裝置,聯結於該 區分裝置以接收其信號,並回應以降低該分離信號組中之 信號的旁瓣雜訊。 (請先閱讀背面之注意事項再填寫本頁) 經濟部屮央標準局負工消費合作社印製 夂紙張尺度適用中面β家揉準(CNS>T4規格(210X297公釐) -47 -
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US34104789A | 1989-04-20 | 1989-04-20 | |
US07/508,131 US5111823A (en) | 1989-04-20 | 1990-04-16 | Apparatus and method for generating echographic images |
Publications (1)
Publication Number | Publication Date |
---|---|
TW202496B true TW202496B (zh) | 1993-03-21 |
Family
ID=26992356
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
TW080100256A TW202496B (zh) | 1989-04-20 | 1991-01-09 |
Country Status (8)
Country | Link |
---|---|
US (1) | US5111823A (zh) |
EP (1) | EP0469076A4 (zh) |
JP (1) | JP2801396B2 (zh) |
AU (1) | AU646805B2 (zh) |
CA (1) | CA2053249A1 (zh) |
NZ (1) | NZ233395A (zh) |
TW (1) | TW202496B (zh) |
WO (1) | WO1990012540A1 (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10271821B2 (en) | 2014-12-23 | 2019-04-30 | Industrial Technology Research Institute | Method of ultrasound imaging and ultrasound scanner |
Families Citing this family (75)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5383457A (en) * | 1987-04-20 | 1995-01-24 | National Fertility Institute | Method and apparatus for processing images |
US5299577A (en) * | 1989-04-20 | 1994-04-05 | National Fertility Institute | Apparatus and method for image processing including one-dimensional clean approximation |
JP3000233B2 (ja) * | 1990-11-19 | 2000-01-17 | 株式会社日立メデイコ | 超音波診断装置 |
JP2950632B2 (ja) * | 1991-03-26 | 1999-09-20 | 株式会社東芝 | 超音波診断装置 |
JP3109749B2 (ja) * | 1991-04-17 | 2000-11-20 | 株式会社東芝 | 超音波映像化装置 |
US5278757A (en) * | 1991-11-15 | 1994-01-11 | The Trustees Of The University Of Pennsylvania | Synthetic aperture ultrasonic imaging system using a minimum or reduced redundancy phased array |
US5482044A (en) * | 1992-01-14 | 1996-01-09 | Diasonics Ultrasound, Inc. | Direct demodulation in ultrasound instruments |
US5271404A (en) * | 1992-06-25 | 1993-12-21 | Cardiometrics, Inc. | Method and apparatus for processing signal data to form an envelope on line |
US5357964A (en) * | 1993-02-08 | 1994-10-25 | Spivey Brett A | Doppler imaging device |
US5515852A (en) * | 1994-06-06 | 1996-05-14 | Hewlett-Packard Company | Method and apparatus for a detection strength spatial filter in an ultrasound imaging system |
US5551433A (en) * | 1994-08-05 | 1996-09-03 | Acuson Corporation | Method and apparatus for a geometric aberration transform in an adaptive focusing ultrasound beamformer system |
US5862269A (en) * | 1994-11-23 | 1999-01-19 | Trustees Of Boston University | Apparatus and method for rapidly convergent parallel processed deconvolution |
US6005827A (en) * | 1995-03-02 | 1999-12-21 | Acuson Corporation | Ultrasonic harmonic imaging system and method |
US6009046A (en) * | 1995-03-02 | 1999-12-28 | Acuson Corporation | Ultrasonic harmonic imaging system and method |
US6104670A (en) * | 1995-03-02 | 2000-08-15 | Acuson Corporation | Ultrasonic harmonic imaging system and method |
US5608690A (en) * | 1995-03-02 | 1997-03-04 | Acuson Corporation | Transmit beamformer with frequency dependent focus |
US6027448A (en) * | 1995-03-02 | 2000-02-22 | Acuson Corporation | Ultrasonic transducer and method for harmonic imaging |
JP4237256B2 (ja) * | 1996-02-29 | 2009-03-11 | シーメンス メディカル ソリューションズ ユーエスエイ インコーポレイテッド | 超音波トランスジューサ |
US5846202A (en) * | 1996-07-30 | 1998-12-08 | Acuson Corporation | Ultrasound method and system for imaging |
US7104956B1 (en) | 1996-11-08 | 2006-09-12 | Research Corporation Technologies, Inc. | Finite amplitude distortion-based inhomogeneous pulse echo ultrasonic imaging |
US6030344A (en) * | 1996-12-04 | 2000-02-29 | Acuson Corporation | Methods and apparatus for ultrasound image quantification |
US5961460A (en) * | 1997-04-11 | 1999-10-05 | Acuson Corporation | Ultrasound imaging enhancement methods and systems |
US5882306A (en) * | 1997-04-11 | 1999-03-16 | Acuson Corporation | Ultrasound imaging methods and systems |
US6110120A (en) | 1997-04-11 | 2000-08-29 | Acuson Corporation | Gated ultrasound imaging apparatus and method |
US6050944A (en) * | 1997-06-17 | 2000-04-18 | Acuson Corporation | Method and apparatus for frequency control of an ultrasound system |
US6193659B1 (en) | 1997-07-15 | 2001-02-27 | Acuson Corporation | Medical ultrasonic diagnostic imaging method and apparatus |
US5913823A (en) * | 1997-07-15 | 1999-06-22 | Acuson Corporation | Ultrasound imaging method and system for transmit signal generation for an ultrasonic imaging system capable of harmonic imaging |
US5833614A (en) * | 1997-07-15 | 1998-11-10 | Acuson Corporation | Ultrasonic imaging method and apparatus for generating pulse width modulated waveforms with reduced harmonic response |
US6132374A (en) * | 1997-08-01 | 2000-10-17 | Acuson Corporation | Ultrasonic imaging method and system |
US6023977A (en) * | 1997-08-01 | 2000-02-15 | Acuson Corporation | Ultrasonic imaging aberration correction system and method |
US6312379B1 (en) | 1997-08-15 | 2001-11-06 | Acuson Corporation | Ultrasonic harmonic imaging system and method using waveform pre-distortion |
US5944666A (en) * | 1997-08-21 | 1999-08-31 | Acuson Corporation | Ultrasonic method for imaging blood flow including disruption or activation of contrast agent |
US5928151A (en) * | 1997-08-22 | 1999-07-27 | Acuson Corporation | Ultrasonic system and method for harmonic imaging in three dimensions |
US5873830A (en) * | 1997-08-22 | 1999-02-23 | Acuson Corporation | Ultrasound imaging system and method for improving resolution and operation |
US6106465A (en) * | 1997-08-22 | 2000-08-22 | Acuson Corporation | Ultrasonic method and system for boundary detection of an object of interest in an ultrasound image |
US5860931A (en) * | 1997-10-10 | 1999-01-19 | Acuson Corporation | Ultrasound method and system for measuring perfusion |
US5935069A (en) * | 1997-10-10 | 1999-08-10 | Acuson Corporation | Ultrasound system and method for variable transmission of ultrasonic signals |
US5897500A (en) * | 1997-12-18 | 1999-04-27 | Acuson Corporation | Ultrasonic imaging system and method for displaying composite fundamental and harmonic images |
US6066095A (en) * | 1998-05-13 | 2000-05-23 | Duke University | Ultrasound methods, systems, and computer program products for determining movement of biological tissues |
DE19822957C1 (de) * | 1998-05-22 | 2000-05-25 | Deutsch Zentr Luft & Raumfahrt | Verfahren zur Detektion und Unterdrückung von Störsignalen in SAR-Daten und Einrichtung zur Durchführung des Verfahrens |
US5957852A (en) * | 1998-06-02 | 1999-09-28 | Acuson Corporation | Ultrasonic harmonic imaging system and method |
US6511426B1 (en) | 1998-06-02 | 2003-01-28 | Acuson Corporation | Medical diagnostic ultrasound system and method for versatile processing |
US6116244A (en) * | 1998-06-02 | 2000-09-12 | Acuson Corporation | Ultrasonic system and method for three-dimensional imaging with opacity control |
US6048316A (en) * | 1998-10-16 | 2000-04-11 | Acuson Corporation | Medical diagnostic ultrasonic imaging system and method for displaying composite fundamental and harmonic images |
DE19915583A1 (de) * | 1999-04-07 | 2000-10-12 | Miroslaw Wrobel | Vorrichtung und Verfahren für die 3D-Echtzeitsonographie |
US6775400B1 (en) | 1999-10-29 | 2004-08-10 | Acuson Corporation | Medical diagnostic ultrasonic imaging method and apparatus for suppressing electronic noise |
DE10027828A1 (de) * | 2000-06-05 | 2001-12-06 | Sonem Gmbh | Aktives Ultraschall-Sichtgerät |
US6508764B1 (en) | 2000-10-31 | 2003-01-21 | Koninklijke Philips Electronics N.V. | Aberration correction apparatus and methods |
DE10112034A1 (de) * | 2001-03-14 | 2002-10-02 | Sonem Gmbh | Anordnung zur Bildwiedergabe für Computertomographen mit Ultraschall |
IL144110A0 (en) * | 2001-07-02 | 2002-05-23 | Sharony Reuven | Methods and apparatus for objective fetal diagnosis |
US20030050648A1 (en) | 2001-09-11 | 2003-03-13 | Spiration, Inc. | Removable lung reduction devices, systems, and methods |
US6592594B2 (en) | 2001-10-25 | 2003-07-15 | Spiration, Inc. | Bronchial obstruction device deployment system and method |
US20030181922A1 (en) | 2002-03-20 | 2003-09-25 | Spiration, Inc. | Removable anchored lung volume reduction devices and methods |
US6852108B2 (en) * | 2002-05-14 | 2005-02-08 | Spiration, Inc. | Apparatus and method for resecting and removing selected body tissue from a site inside a patient |
AU2003285747A1 (en) | 2002-12-15 | 2004-07-09 | Ultramove Ltd. | System and method for determination of fetal movement |
US7100616B2 (en) | 2003-04-08 | 2006-09-05 | Spiration, Inc. | Bronchoscopic lung volume reduction method |
US7998073B2 (en) * | 2003-08-04 | 2011-08-16 | Imacor Inc. | Ultrasound imaging with reduced noise |
US7533671B2 (en) | 2003-08-08 | 2009-05-19 | Spiration, Inc. | Bronchoscopic repair of air leaks in a lung |
CN100566664C (zh) * | 2003-11-26 | 2009-12-09 | 普瑞斯玛医药技术有限责任公司 | 使用窄探针的经食管的超声波 |
JP2007513726A (ja) * | 2003-12-16 | 2007-05-31 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 浸透、解像度及びフレームレートの自動制御を有する超音波画像診断システム |
US8234923B2 (en) * | 2004-09-20 | 2012-08-07 | Innervision Medical Technologies Inc. | Systems and methods for ultrasound imaging |
US20060106309A1 (en) * | 2004-11-16 | 2006-05-18 | Siemens Medical Solutions Usa, Inc. | Aberration correction beam patterns |
KR101316796B1 (ko) * | 2005-04-14 | 2013-10-11 | 베라소닉스, 인코포레이티드 | 화소 배향 프로세싱을 갖는 초음파 이미지화 시스템 |
US20070047743A1 (en) * | 2005-08-26 | 2007-03-01 | Step Communications Corporation, A Nevada Corporation | Method and apparatus for improving noise discrimination using enhanced phase difference value |
US7472041B2 (en) | 2005-08-26 | 2008-12-30 | Step Communications Corporation | Method and apparatus for accommodating device and/or signal mismatch in a sensor array |
US7619563B2 (en) * | 2005-08-26 | 2009-11-17 | Step Communications Corporation | Beam former using phase difference enhancement |
US7415372B2 (en) * | 2005-08-26 | 2008-08-19 | Step Communications Corporation | Method and apparatus for improving noise discrimination in multiple sensor pairs |
US7691151B2 (en) | 2006-03-31 | 2010-04-06 | Spiration, Inc. | Articulable Anchor |
US8193967B2 (en) * | 2008-12-10 | 2012-06-05 | The United States Of America As Represented By The Secretary Of The Army | Method and system for forming very low noise imagery using pixel classification |
US7796829B2 (en) * | 2008-12-10 | 2010-09-14 | The United States Of America As Represented By The Secretary Of The Army | Method and system for forming an image with enhanced contrast and/or reduced noise |
US9250323B2 (en) | 2008-12-10 | 2016-02-02 | The United States Of America As Represented By The Secretary Of The Army | Target detection utilizing image array comparison |
US8665132B2 (en) * | 2008-12-10 | 2014-03-04 | The United States Of America As Represented By The Secretary Of The Army | System and method for iterative fourier side lobe reduction |
KR101123005B1 (ko) | 2010-06-14 | 2012-03-12 | 알피니언메디칼시스템 주식회사 | 초음파 진단장치, 거기에 이용되는 그래픽 환경 제어장치 및 그 제어방법 |
US8795241B2 (en) | 2011-05-13 | 2014-08-05 | Spiration, Inc. | Deployment catheter |
JP6340086B2 (ja) * | 2014-12-26 | 2018-06-06 | 古野電気株式会社 | 超音波体組織検出装置、超音波体組織検出方法、および、超音波体組織検出プログラム |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4442715A (en) * | 1980-10-23 | 1984-04-17 | General Electric Company | Variable frequency ultrasonic system |
US4397006A (en) * | 1980-12-31 | 1983-08-02 | Mobil Oil Corporation | Cross trace coherent noise filtering for seismograms |
JPS5829455A (ja) * | 1981-08-18 | 1983-02-21 | 株式会社東芝 | 超音波診断装置 |
JPS5897347A (ja) * | 1981-12-03 | 1983-06-09 | 株式会社東芝 | 超音波診断装置 |
US4470305A (en) * | 1982-09-27 | 1984-09-11 | General Electric Company | Annular array used as a horn transducer |
JPS5994092A (ja) * | 1982-11-20 | 1984-05-30 | Toshiba Corp | 開口合成法に基づく信号処理装置 |
US4604697A (en) * | 1983-08-05 | 1986-08-05 | Interspec, Inc. | Body imaging using vectorial addition of acoustic reflection to achieve effect of scanning beam continuously focused in range |
US4553437A (en) * | 1984-01-30 | 1985-11-19 | Imaging Associates | Hybrid non-invasive ultrasonic imaging system |
JPS6131135A (ja) * | 1984-07-25 | 1986-02-13 | 株式会社東芝 | 超音波診断装置 |
FR2570838B1 (fr) * | 1984-09-25 | 1986-11-21 | Labo Electronique Physique | Appareil d'examen de milieux par echographie ultrasonore avec focalisation angulaire |
EP0256481B1 (de) * | 1986-08-20 | 1994-01-12 | Siemens Aktiengesellschaft | Verfahren und Einrichtung zur adaptiven Fokussierung bei einem medizinischen Ultraschall-Bildgabegerät |
US4835689A (en) * | 1987-09-21 | 1989-05-30 | General Electric Company | Adaptive coherent energy beam formation using phase conjugation |
US4989143A (en) * | 1987-12-11 | 1991-01-29 | General Electric Company | Adaptive coherent energy beam formation using iterative phase conjugation |
-
1990
- 1990-04-16 US US07/508,131 patent/US5111823A/en not_active Expired - Fee Related
- 1990-04-19 EP EP19900907895 patent/EP0469076A4/en not_active Withdrawn
- 1990-04-19 AU AU55696/90A patent/AU646805B2/en not_active Ceased
- 1990-04-19 CA CA002053249A patent/CA2053249A1/en not_active Abandoned
- 1990-04-19 JP JP2507235A patent/JP2801396B2/ja not_active Expired - Lifetime
- 1990-04-19 WO PCT/US1990/002021 patent/WO1990012540A1/en not_active Application Discontinuation
- 1990-04-20 NZ NZ233395A patent/NZ233395A/en unknown
-
1991
- 1991-01-09 TW TW080100256A patent/TW202496B/zh active
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10271821B2 (en) | 2014-12-23 | 2019-04-30 | Industrial Technology Research Institute | Method of ultrasound imaging and ultrasound scanner |
Also Published As
Publication number | Publication date |
---|---|
AU646805B2 (en) | 1994-03-10 |
JP2801396B2 (ja) | 1998-09-21 |
US5111823A (en) | 1992-05-12 |
AU5569690A (en) | 1990-11-16 |
JPH04505114A (ja) | 1992-09-10 |
NZ233395A (en) | 1992-04-28 |
EP0469076A4 (en) | 1992-08-19 |
CA2053249A1 (en) | 1990-10-21 |
EP0469076A1 (en) | 1992-02-05 |
WO1990012540A1 (en) | 1990-11-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
TW202496B (zh) | ||
AU682816B2 (en) | Method and apparatus for processing images | |
Bercoff | Ultrafast ultrasound imaging | |
US8840554B2 (en) | Ultrasonic 3-dimensional image reconstruction method and ultrasonic wave system thereof | |
Qi et al. | Image quality enhancement using a deep neural network for plane wave medical ultrasound imaging | |
US20220003721A1 (en) | Methods and systems for non-invasively characterising a heterogeneous medium using ultrasound | |
EP2763585A1 (en) | Nonlinear imaging with dual band pulse complexes | |
US11776526B2 (en) | Method and system for ultrasonic characterization of a medium | |
KR20150112057A (ko) | 공간 일관성 기초 초음파 신호 처리 모듈 및 그에 의한 초음파 신호 처리 방법 | |
CA3127822A1 (en) | Method and system for ultrasonic characterization of a medium | |
CA3127924A1 (en) | Method and system for ultrasonic characterization of a medium | |
Lu | Limited diffraction array beams | |
US20240036004A1 (en) | Method and system for ultrasonic characterization of a medium | |
CA3127873A1 (en) | Method and system for ultrasonic characterization of a medium | |
Løvstakken | Signal processing in diagnostic ultrasound: Algorithms for real-time estimation and visualization of blood flow velocity | |
Robert et al. | The time-reversal operator with virtual transducers: Application to far-field aberration correction | |
EP3603526A1 (en) | Ultrasonic diagnostic device and method for controlling ultrasonic diagnostic device | |
Haun et al. | Efficient three-dimensional imaging from a small cylindrical aperture | |
CN108303699A (zh) | 超声相控阵远场超分辨率成像方法、装置、存储介质和系统 | |
Rodríguez et al. | The signal projection method: A general approach to active imaging with arrays | |
Deshmukh | REAL-TIME ELASTOGRAPHY SYSTEMS | |
Svet | Acoustical imaging in heterogeneous environments | |
Ibáñez et al. | The Signal Projection Method: A General Approach to Active Imaging with Arrays. | |
Otawa et al. | Visualization of Untrained Target Shape and Distance Using a Multi-element Ultrasound Transducer and Neural Network | |
Mehra | Frequency modulated pulse for ultrasonic b-scan imaging in attenuating medium |