TW202496B - - Google Patents

Download PDF

Info

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
Application number
TW080100256A
Other languages
English (en)
Original Assignee
Nat Fertility Inst
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nat Fertility Inst filed Critical Nat Fertility Inst
Application granted granted Critical
Publication of TW202496B publication Critical patent/TW202496B/zh

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52077Details 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52023Details of receivers
    • G01S7/52025Details of receivers for pulse systems
    • G01S7/52026Extracting wanted echo signals
    • G01S7/52028Extracting wanted echo signals using digital techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52046Techniques for image enhancement involving transmitter or receiver
    • G01S7/52047Techniques for image enhancement involving transmitter or receiver for elimination of side lobes or of grating lobes; for increasing resolving power
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52046Techniques for image enhancement involving transmitter or receiver
    • G01S7/52049Techniques 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)

  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 -
TW080100256A 1989-04-20 1991-01-09 TW202496B (zh)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Cited By (1)

* Cited by examiner, † Cited by third party
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