JP2024511870A - 発散レンズアレイ - Google Patents
発散レンズアレイ Download PDFInfo
- Publication number
- JP2024511870A JP2024511870A JP2023560710A JP2023560710A JP2024511870A JP 2024511870 A JP2024511870 A JP 2024511870A JP 2023560710 A JP2023560710 A JP 2023560710A JP 2023560710 A JP2023560710 A JP 2023560710A JP 2024511870 A JP2024511870 A JP 2024511870A
- Authority
- JP
- Japan
- Prior art keywords
- ultrasound
- transducer
- imaging
- curvature
- transducers
- Prior art date
- Legal status (The legal status 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 status listed.)
- Pending
Links
- 238000002604 ultrasonography Methods 0.000 claims abstract description 148
- 238000003384 imaging method Methods 0.000 claims abstract description 75
- 238000012285 ultrasound imaging Methods 0.000 claims abstract description 34
- 239000000523 sample Substances 0.000 claims description 50
- 238000012937 correction Methods 0.000 claims description 21
- 239000011159 matrix material Substances 0.000 claims description 18
- 238000000034 method Methods 0.000 claims description 17
- 230000005540 biological transmission Effects 0.000 claims description 14
- 239000007787 solid Substances 0.000 claims description 7
- 239000000463 material Substances 0.000 claims description 6
- 230000004044 response Effects 0.000 claims description 6
- 239000002131 composite material Substances 0.000 claims description 5
- 230000001427 coherent effect Effects 0.000 claims description 3
- 238000013329 compounding Methods 0.000 claims description 3
- 230000004807 localization Effects 0.000 claims description 3
- 230000000644 propagated effect Effects 0.000 claims 3
- 210000004556 brain Anatomy 0.000 abstract description 10
- 230000035945 sensitivity Effects 0.000 abstract description 10
- 238000003491 array Methods 0.000 abstract description 4
- 210000005166 vasculature Anatomy 0.000 description 12
- 210000000056 organ Anatomy 0.000 description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 8
- 230000015654 memory Effects 0.000 description 7
- 201000010099 disease Diseases 0.000 description 5
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 5
- 210000000988 bone and bone Anatomy 0.000 description 4
- 239000007788 liquid Substances 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000001902 propagating effect Effects 0.000 description 4
- 210000004204 blood vessel Anatomy 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 239000004205 dimethyl polysiloxane Substances 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 229920000435 poly(dimethylsiloxane) Polymers 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 230000000747 cardiac effect Effects 0.000 description 2
- -1 for example Substances 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000035515 penetration Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 210000001519 tissue Anatomy 0.000 description 2
- 208000006386 Bone Resorption Diseases 0.000 description 1
- 239000004593 Epoxy Substances 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 230000024279 bone resorption Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 239000000499 gel Substances 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000004089 microcirculation Effects 0.000 description 1
- 238000000386 microscopy Methods 0.000 description 1
- 210000004088 microvessel Anatomy 0.000 description 1
- 238000002610 neuroimaging Methods 0.000 description 1
- 238000010895 photoacoustic effect Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 230000005236 sound signal Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4483—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/483—Diagnostic techniques involving the acquisition of a 3D volume of data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5269—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B06—GENERATING OR TRANSMITTING MECHANICAL VIBRATIONS IN GENERAL
- B06B—METHODS OR APPARATUS FOR GENERATING OR TRANSMITTING MECHANICAL VIBRATIONS OF INFRASONIC, SONIC, OR ULTRASONIC FREQUENCY, e.g. FOR PERFORMING MECHANICAL WORK IN GENERAL
- B06B1/00—Methods or apparatus for generating mechanical vibrations of infrasonic, sonic, or ultrasonic frequency
- B06B1/02—Methods or apparatus for generating mechanical vibrations of infrasonic, sonic, or ultrasonic frequency making use of electrical energy
- B06B1/06—Methods or apparatus for generating mechanical vibrations of infrasonic, sonic, or ultrasonic frequency making use of electrical energy operating with piezoelectric effect or with electrostriction
- B06B1/0607—Methods or apparatus for generating mechanical vibrations of infrasonic, sonic, or ultrasonic frequency making use of electrical energy operating with piezoelectric effect or with electrostriction using multiple elements
- B06B1/0622—Methods or apparatus for generating mechanical vibrations of infrasonic, sonic, or ultrasonic frequency making use of electrical energy operating with piezoelectric effect or with electrostriction using multiple elements on one surface
- B06B1/0629—Square array
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B06—GENERATING OR TRANSMITTING MECHANICAL VIBRATIONS IN GENERAL
- B06B—METHODS OR APPARATUS FOR GENERATING OR TRANSMITTING MECHANICAL VIBRATIONS OF INFRASONIC, SONIC, OR ULTRASONIC FREQUENCY, e.g. FOR PERFORMING MECHANICAL WORK IN GENERAL
- B06B1/00—Methods or apparatus for generating mechanical vibrations of infrasonic, sonic, or ultrasonic frequency
- B06B1/02—Methods or apparatus for generating mechanical vibrations of infrasonic, sonic, or ultrasonic frequency making use of electrical energy
- B06B1/06—Methods or apparatus for generating mechanical vibrations of infrasonic, sonic, or ultrasonic frequency making use of electrical energy operating with piezoelectric effect or with electrostriction
- B06B1/0644—Methods or apparatus for generating mechanical vibrations of infrasonic, sonic, or ultrasonic frequency making use of electrical energy operating with piezoelectric effect or with electrostriction using a single piezoelectric element
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B06—GENERATING OR TRANSMITTING MECHANICAL VIBRATIONS IN GENERAL
- B06B—METHODS OR APPARATUS FOR GENERATING OR TRANSMITTING MECHANICAL VIBRATIONS OF INFRASONIC, SONIC, OR ULTRASONIC FREQUENCY, e.g. FOR PERFORMING MECHANICAL WORK IN GENERAL
- B06B1/00—Methods or apparatus for generating mechanical vibrations of infrasonic, sonic, or ultrasonic frequency
- B06B1/02—Methods or apparatus for generating mechanical vibrations of infrasonic, sonic, or ultrasonic frequency making use of electrical energy
- B06B1/06—Methods or apparatus for generating mechanical vibrations of infrasonic, sonic, or ultrasonic frequency making use of electrical energy operating with piezoelectric effect or with electrostriction
- B06B1/0644—Methods or apparatus for generating mechanical vibrations of infrasonic, sonic, or ultrasonic frequency making use of electrical energy operating with piezoelectric effect or with electrostriction using a single piezoelectric element
- B06B1/0651—Methods or apparatus for generating mechanical vibrations of infrasonic, sonic, or ultrasonic frequency making use of electrical energy operating with piezoelectric effect or with electrostriction using a single piezoelectric element of circular shape
-
- 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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8909—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
- G01S15/8915—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array
- G01S15/8925—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array the array being a two-dimensional transducer configuration, i.e. matrix or orthogonal linear arrays
-
- 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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8993—Three dimensional imaging systems
-
- 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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8997—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using synthetic aperture 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
-
- 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/52079—Constructional features
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K11/00—Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/18—Methods or devices for transmitting, conducting or directing sound
- G10K11/26—Sound-focusing or directing, e.g. scanning
- G10K11/32—Sound-focusing or directing, e.g. scanning characterised by the shape of the source
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K11/00—Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/18—Methods or devices for transmitting, conducting or directing sound
- G10K11/26—Sound-focusing or directing, e.g. scanning
- G10K11/34—Sound-focusing or directing, e.g. scanning using electrical steering of transducer arrays, e.g. beam steering
- G10K11/341—Circuits therefor
- G10K11/346—Circuits therefor using phase variation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0808—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of the brain
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0883—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of the heart
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0891—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of blood vessels
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/42—Details of probe positioning or probe attachment to the patient
- A61B8/4209—Details of probe positioning or probe attachment to the patient by using holders, e.g. positioning frames
- A61B8/4218—Details of probe positioning or probe attachment to the patient by using holders, e.g. positioning frames characterised by articulated arms
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/42—Details of probe positioning or probe attachment to the patient
- A61B8/4209—Details of probe positioning or probe attachment to the patient by using holders, e.g. positioning frames
- A61B8/4227—Details of probe positioning or probe attachment to the patient by using holders, e.g. positioning frames characterised by straps, belts, cuffs or braces
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4483—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer
- A61B8/4488—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer the transducer being a phased array
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/481—Diagnostic techniques involving the use of contrast agent, e.g. microbubbles introduced into the bloodstream
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5238—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
- A61B8/5246—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from the same or different imaging techniques, e.g. color Doppler and B-mode
- A61B8/5253—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from the same or different imaging techniques, e.g. color Doppler and B-mode combining overlapping images, e.g. spatial compounding
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/54—Control of the diagnostic device
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B06—GENERATING OR TRANSMITTING MECHANICAL VIBRATIONS IN GENERAL
- B06B—METHODS OR APPARATUS FOR GENERATING OR TRANSMITTING MECHANICAL VIBRATIONS OF INFRASONIC, SONIC, OR ULTRASONIC FREQUENCY, e.g. FOR PERFORMING MECHANICAL WORK IN GENERAL
- B06B2201/00—Indexing scheme associated with B06B1/0207 for details covered by B06B1/0207 but not provided for in any of its subgroups
- B06B2201/70—Specific application
- B06B2201/76—Medical, dental
-
- 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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8979—Combined Doppler and pulse-echo imaging systems
- G01S15/8981—Discriminating between fixed and moving objects or between objects moving at different speeds, e.g. wall clutter filter
-
- 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/52036—Details of receivers using analysis of echo signal for target characterisation
- G01S7/52038—Details of receivers using analysis of echo signal for target characterisation involving non-linear properties of the propagation medium or of the reflective target
- G01S7/52039—Details of receivers using analysis of echo signal for target characterisation involving non-linear properties of the propagation medium or of the reflective target exploiting the non-linear response of a contrast enhancer, e.g. a contrast agent
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Pathology (AREA)
- Animal Behavior & Ethology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biophysics (AREA)
- Heart & Thoracic Surgery (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Mechanical Engineering (AREA)
- Gynecology & Obstetrics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Nonlinear Science (AREA)
- Multimedia (AREA)
- Hematology (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Transducers For Ultrasonic Waves (AREA)
Abstract
医療分野で3D超音波撮像が強力なツールになっている一方で、主な欠点は、主に、変換器の2Dアレイの寸法に関する、大きな3D立体を撮像する難しさである。空間解像度を損なわないため、変換器のアレイを使用する必要があり、変換器のサイズは、超音波の波長を超えない。そのような要件によって、現在の技術では達成可能でない、または費用が高くなり過ぎる大きな3D体積を撮像するためのアレイの寸法になる。本開示は、より大きな変換器を使用することによって、上の技術的制約を克服し、各変換器が曲がった形状を有する受信面を有し、または、音響レンズに嵌合される。変換器のそのような構成は、脳または心臓のような大きな3D体積を、高解像度および高感度で撮像するのに好適な変換器の2Dアレイをもたらす。
Description
本発明は、超音波撮像技法に関し、具体的には、大きな視野での3D立体撮像のための技法に関する。
今日では、医療分野で超音波の使用に対して関心があることは、もはや証明する必要がない。
一般的に、超音波は、個別に超音波を生成することが可能な複数の変換器のアレイを備えるプローブから生成することができる。これらの超音波は、応答して後方散乱信号を生成することができる媒体に向けることができ、後方散乱信号は、順に同じ変換器のアレイまたは異なるアレイによって記録することができる。後方散乱信号から、媒体の2Dまたは3D画像を再構築するのを可能とすることができる。超音波撮像分野での最近の進歩によって、アニメーション化した3D超音波撮像を得ることを可能にする4D超音波撮像を行う可能性がもたらされている。
医療分野で3Dまたは4D超音波撮像が強力なツールになっている一方で、それらの主な欠点は、大きな3D立体を撮像する難しさである。実際に、大きな3D立体を撮像するための容量が、変換器のアレイの寸法に直接関係する。空間解像度を損なわないため、変換器のアレイを使用する必要があり、ここで変換器のサイズは、超音波の波長を超えない。この要件によって、現在の技術では達成可能でない、または費用が高くなり過ぎるアレイの構成になる。たとえば、3D撮像での良好な解像度および大きな視野を保つ目的で、10cm×10cmの開口(すなわち、10cm×10cmのアレイ)および(6MHzの中心周波数に対応する)約250マイクロメートルの変換器のサイズを有するプローブによって、160000個の変換器のアレイとなる。
大きな視野で3D立体撮像を実施するために、従来技術では、より少ない変換器を有するアレイを使用するためにある粗なマトリックスを使用する。粗なマトリックスが3D立体撮像を実施するのを可能にするが、より少ない数の変換器によって、3D撮像の品質に影響がおよぶ。
さらに、古典的アレイまたは粗なアレイの両方は、通常では、骨の背後にあってもよい臓器全体の3Dまたは4D超音波立体撮像を、撮像された臓器における特定の場所(たとえば、疾患)を決定するのに十分な解像度および感度で実施することが可能ではない。実際には、骨が強く吸収するために、従来技術は、通常では、骨における(人工的または自然な)開口を使用して超音波3D撮像を実施し、大きな視野の撮像および高い横方向解像度で特定の場所を決定する可能性を難しくさせる。
J.A. Jensen、S.I. Nikolov、K.L. Gammelmark、M.H. Pedersen著、「Synthetic Aperture Ultrasound Imaging」、Ultrasonics、Elsevier、The Netherlands、vol. 44、no.1、e5~e16頁, 2006年
したがって、画像品質が損なうことなく、大きな視野で3D超音波立体撮像を実施する必要があり、特に、ある種の場合には骨の背後に配置される場合がある臓器の3D超音波立体撮像を実施する必要がある。
このために、本開示は、変換器を有するアレイを備えることができる超音波3D撮像プローブを提案し、各変換器が、少なくとも、媒体中の撮像する領域から前記媒体における超音波放出に続く後方散乱超音波を取得するように構成することができ、前記後方散乱超音波が前記媒体中で伝播するときのある波長λに対応するある中心周波数を前記後方散乱超音波が有することができ、前記アレイが前記媒体と接触して配置されるように適合されるアレイ表面を形成し、各変換器が少なくとも前記後方散乱超音波を受け取るための受信面を有することができ、
各変換器が、2λよりも大きな直径を有する内接円を規定するアレイ表面に形状を有することができ、各変換器が
- 前記受信面に垂直な少なくとも1つの平面中に、湾曲を有する受信面の曲がった形状を有することができ、
- および/または、前記受信面に垂直な少なくとも1つの平面中に、湾曲を有する曲がった形状を有する変換器をエミュレートするように設計された音響レンズと嵌合することができ、
- および/または、独立に制御されるように適合された変換器要素のマトリックスによって形成される。
各変換器が、2λよりも大きな直径を有する内接円を規定するアレイ表面に形状を有することができ、各変換器が
- 前記受信面に垂直な少なくとも1つの平面中に、湾曲を有する受信面の曲がった形状を有することができ、
- および/または、前記受信面に垂直な少なくとも1つの平面中に、湾曲を有する曲がった形状を有する変換器をエミュレートするように設計された音響レンズと嵌合することができ、
- および/または、独立に制御されるように適合された変換器要素のマトリックスによって形成される。
したがって、有利なことに、2λよりも大きな直径を有する内接円を規定するアレイ表面に形状を有する変換器を使用すると、超音波撮像で大きな視野を実施するのに必要な大きな表面を有するプローブの使用を妨げる技術的制約を克服することが可能である。さらに、変換器用により大きな受信面を使用することと、曲がった形状または/および受信面用に曲がった変換器をエミュレートするように設計された音響レンズ(たとえば、発散レンズ)を使用することとを組み合わせることによって、通常では、2λよりも大きな直径を有する内接円を規定するアレイ表面に形状を有する変換器の場合であるはずがない、横方向および時間解像度を向上させることが可能となる。
(たとえば、変換器からの)超音波の放出、または超音波(すなわち、媒体からの後方散乱超音波)の受信は、焦点をぼかすことができる。曲がった形状を有する、もしくは曲がった変換器をエミュレートするように設計された音響レンズを有する、または、独立に制御されるように適合された変換器要素のマトリックスによって形成される各変換器は、変換器の開口角を広げる(すなわち、指向性を低下させる)ように構成することができ、すなわち、発散レンズとして働くことができ、したがって、受信または/および放出において、超音波の焦点をぼかすように構成することができる。
2λよりも大きな直径を有する内接円を規定するアレイ表面に形状を有すると共に、曲がった形状を有する、もしくは曲がった変換器をエミュレートするように設計された音響レンズを有する、または、独立に制御されるように適合された変換器要素のマトリックスによって形成される変換器を有する超音波3D撮像プローブは、2λよりも大きな直径に起因して大きな振幅を有する超音波を送信および/もしくは受信することが可能であってもよく、媒体中で発散する(または、発散)超音波を送信することによって低い指向性を維持することが可能であってもよく、このことによって、視野、解像度、および感度の項目で画像/立体品質を向上させることが可能となることができる。
こうして、そのような性質によって、大きな視野で3Dまたは4D超音波撮像を実施する実現性がもたらされる。
曲がった形状を有する変換器によって、曲がった送信または/および受信面を有する変換器が理解されるべきである。
媒体は、たとえば水のような媒体または/および柔らかい組織であってもよく、それらを含んでもよい。
1つまたは複数の実施形態では、各変換器の受信面が曲がった形状を有してもよく、湾曲は、λと40λの間に含まれる少なくとも1つの曲率半径を有してもよい。
1つまたは複数の実施形態では、各変換器の受信面が曲がった形状を有してもよく、曲がった形状は凸であってもよい。
1つまたは複数の実施形態では、各変換器の受信面が、曲がった形状を有する変換器をエミュレートするように設計された音響レンズと嵌合することができ、音響レンズが、変換器の前記受信面の反対側に外面を有することができ、前記外面は凸または凹であってもよい。
1つまたは複数の実施形態では、外面は凸であってもよく、前記外面に垂直な少なくとも1つの平面に湾曲を有してもよく、前記湾曲が、λと40λの間に含まれる少なくとも1つの曲率半径を有してもよい。
1つまたは複数の実施形態では、外面は凹であってもよく、前記外面に垂直な少なくとも1つの平面に湾曲を有してもよく、前記湾曲が、λと40λの間に含まれる少なくとも1つの曲率半径を有してもよい。
1つまたは複数の実施形態では、音響レンズは、それぞれがいくつかの音速に関連するいくつかの材料を含む複合音響レンズであってもよい。
1つまたは複数の実施形態では、アレイは、1000λ2よりも大きくてもよい表面積を有することができる。
1つまたは複数の実施形態では、直径は30マイクロメートルよりも大きくてもよい。
1つまたは複数の実施形態では、直径は100マイクロメートルと2000マイクロメートルの間であってもよい。
1つまたは複数の実施形態では、中心周波数は、0.5MHzと100MHzの間、好ましくは、1MHzと10MHzの間に含まれてもよい。
本開示は、超音波3D撮像デバイスにやはり関し、超音波3D撮像デバイスは、少なくとも、
- 本開示にしたがった超音波3D撮像プローブと、
- 各送信された超音波に応答して撮像する領域から変換器によって受信された前記後方散乱超音波から信号を取得するために構成される制御ユニットであって、取得した信号に基づいて3D画像を生成するためさらに構成することができる、制御ユニットと
を備えることができる。
- 本開示にしたがった超音波3D撮像プローブと、
- 各送信された超音波に応答して撮像する領域から変換器によって受信された前記後方散乱超音波から信号を取得するために構成される制御ユニットであって、取得した信号に基づいて3D画像を生成するためさらに構成することができる、制御ユニットと
を備えることができる。
1つまたは複数の実施形態では、制御ユニットは、前記後方散乱超音波を生成するため撮像する領域中に前記変換器に超音波を送信させるようにさらに構成することができる。
したがって、有利なことに、従来型変換器より大きな寸法の変換器を使用することによって、媒体中により多いエネルギーを送信し、受信においてより良好な感度を得ることが可能になる。
1つまたは複数の実施形態では、制御ユニットは、各送信超音波の送信後に得られた取得信号に補正遅延を適用するようにさらに構成することができ、前記補正遅延は、送信面の曲がった形状を考慮に入れる。
1つまたは複数の実施形態では、超音波3D撮像プローブの各変換器は、独立に制御されるように適合された変換器要素のマトリックスによって形成することができ、制御ユニットは、変換器規模で曲がった形状を有する変換器をエミュレートするため変換器のマトリックスの変換器要素に時間遅延を適用するように構成することができる。
本開示は、撮像する領域中を伝播する後方散乱超音波を取得するための本開示にしたがった超音波3D撮像プローブを使用するステップを含む超音波撮像方法にやはり関し、後方散乱超音波は波長λを有することができる。
本開示は、各送信超音波に応答して撮像する領域中を伝播する後方散乱超音波からの信号を取得するため、および、取得した信号に基づいて3D画像を生成するための本開示にしたがった超音波3D撮像デバイスを使用するステップを含む超音波撮像方法にやはり関し、後方散乱超音波は波長λを有することができる。
1つまたは複数の実施形態では、超音波は、合成開口法にしたがって、撮像する領域中に送信される。
1つまたは複数の実施形態では、超音波3D撮像デバイスを使用するステップを含む超音波撮像方法は、
- 補正した取得信号を得るために、各送信超音波の後に得られる取得信号に補正遅延を適用するステップであって、前記補正遅延が、送信面の曲がった形状を考慮に入れるステップと、
- 動的合焦および3D画像の生成を実施するために補正した取得信号にビーム形成を適用するステップであって、前記ビーム形成が、各送信波の受信において実施されて、コヒーレントコンパウンディングまたは合成開口撮像を実施することができる、ステップと
をさらに含むことができる。
- 補正した取得信号を得るために、各送信超音波の後に得られる取得信号に補正遅延を適用するステップであって、前記補正遅延が、送信面の曲がった形状を考慮に入れるステップと、
- 動的合焦および3D画像の生成を実施するために補正した取得信号にビーム形成を適用するステップであって、前記ビーム形成が、各送信波の受信において実施されて、コヒーレントコンパウンディングまたは合成開口撮像を実施することができる、ステップと
をさらに含むことができる。
1つまたは複数の実施形態では、1つの3D撮像を実施するために、ある数の送信超音波を使用することができ、前記数は、100未満であってもよい。
1つまたは複数の実施形態では、1つの3D撮像を実施するため使用される送信超音波の数は、30未満であってもよい。
1つまたは複数の実施形態では、撮像する媒体中にマイクロバブルが存在する場合があり、スーパーローカリゼーションアルゴリズムを使用して3D画像中で位置特定および/または追跡することができる。
1つまたは複数の実施形態では、スーパーローカリゼーションアルゴリズムは、
(a)以下すなわち、
- 組織から信号を除去してマイクロバブルだけを保つために、3D Bモード立体にSVDクラッタフィルタを適用するステップと、
- 各3D Bモード立体中の最大値を検出するステップと、
- 偽陽性マイクロバブルを除去するためにマイクロバブルの平均エネルギーレベルに基づいたエネルギーレベルしきい値を適用するステップと、
- 偽陽性マイクロバブルを除去するために3D点広がり関数(PSF)相関に基づいた相関値のしきい値を適用するステップと
を含む、マイクロバブルを検出するステップと、
(b)以下すなわち、
- 前記マイクロバブルの最大値周辺のボクセルに3D放物面補間を使用することによりマイクロバブルの正確な中心を検出するステップと、
- 正確な中心の座標を使用して、フレーム毎にバブルを追跡するステップと、
- 追跡の位置および速度を回復するステップと
によってマイクロバブルを追跡するステップと、
(c)以下すなわち、
- グリッド上に超解像立体を再構築するステップであって、超解像立体の密度が、グリッドの各ボクセル中のマイクロバブルの数であり、超解像立体の速度が、グリッドの各ボクセルにおける平均速度である、ステップ
によって、超解像立体を生成するステップと、
(d)超解像立体に基づいて3D超解像画像を生成するステップと
を含むことができる。
(a)以下すなわち、
- 組織から信号を除去してマイクロバブルだけを保つために、3D Bモード立体にSVDクラッタフィルタを適用するステップと、
- 各3D Bモード立体中の最大値を検出するステップと、
- 偽陽性マイクロバブルを除去するためにマイクロバブルの平均エネルギーレベルに基づいたエネルギーレベルしきい値を適用するステップと、
- 偽陽性マイクロバブルを除去するために3D点広がり関数(PSF)相関に基づいた相関値のしきい値を適用するステップと
を含む、マイクロバブルを検出するステップと、
(b)以下すなわち、
- 前記マイクロバブルの最大値周辺のボクセルに3D放物面補間を使用することによりマイクロバブルの正確な中心を検出するステップと、
- 正確な中心の座標を使用して、フレーム毎にバブルを追跡するステップと、
- 追跡の位置および速度を回復するステップと
によってマイクロバブルを追跡するステップと、
(c)以下すなわち、
- グリッド上に超解像立体を再構築するステップであって、超解像立体の密度が、グリッドの各ボクセル中のマイクロバブルの数であり、超解像立体の速度が、グリッドの各ボクセルにおける平均速度である、ステップ
によって、超解像立体を生成するステップと、
(d)超解像立体に基づいて3D超解像画像を生成するステップと
を含むことができる。
SVDによって、特異値分解を理解されたい。
1つまたは複数の実施形態では、変換器規模で曲がった形状を有する変換器をエミュレートするため変換器のマトリックスの変換器要素に時間遅延を適用することができる。
他の特徴、詳細、および利点は、以下の詳細な説明および図面に示されることになる。
図では、同じ参照記号は、同一または同様の要素を示す。
図1および2は、超音波撮像のための装置の例を図示する。
図1および図2に示される装置は、たとえば、媒体の領域を撮像するため、超高速3Dまたは4D超音波撮像に適合することができる。媒体は、たとえば、人間の患者といった、生物体の体であってもよい。撮像する領域は、生物体の脳もしくは脳の部分であってもよく、生物体の心臓もしくは心臓の部分であってもよい。
装置は、たとえば、少なくとも2Dアレイ超音波プローブ2および制御システムを含むことができる。装置は、生物体の脳の一部だけ、または脳全体をカバーすることができる。
2Dアレイ超音波プローブ2は、たとえば、1mmよりも狭いピッチで、数百から数千の数N個の変換器Tijを有することができる。2Dアレイ超音波プローブの変換器は、XY平面において異なる形状を呈する場合がある。一実施形態では、変換器の形状は、正方形、長方形、または円形であってもよい。変換器の厚さは、圧電材料に依存して、10マイクロメートルと1000マイクロメートルの間に含まれてもよい。
2Dアレイ超音波プローブ2は、XY平面に垂直な軸Zに沿って超音波を送信する、2つの垂直な軸X、Yに沿ったマトリックスとして配設されるn*nの変換器を有することができる。1つの具体的な例では、2Dアレイ超音波プローブ2は、6mmピッチで、256個の変換器Tij(16×16)を有することができる。
さらに、2Dアレイ超音波プローブは、受信器、または放出器、または放出器および受信器として働くように構成することができる。
超音波の周波数(または、中心周波数)は、たとえば0.5MHzと100MHzの間、たとえば1MHzと10MHzの間に含まれてもよい。
さらに、制御システムは、超音波が毎秒100超音波以上、たとえば毎秒数百から数千の超音波の割合で送信されるようにプログラム(または構成)することができる。
制御システムは、たとえば、制御ユニット3およびコンピュータ4を含むことができる。この例では、制御ユニット3は、2Dアレイ超音波プローブ2を制御し、それから信号を取得するために使用される一方で、コンピュータ4は、制御ユニット3を制御し、制御ユニット3によって取得された信号から3D画像シーケンスを生成し、それから数量化パラメータを決定するために使用される。変形形態では、単一の電子デバイスが、制御ユニット3およびコンピュータ4の全部の機能性を満足させることができる。
図2に示されるように、制御ユニット3は、たとえば以下を含むことができる。
- 2Dアレイ超音波プローブ2のn個の変換器Tijに個別に接続される、n*nのアナログ/デジタル変換器5(ADij)
- n*nのアナログ/デジタル変換器5にそれぞれ接続される、n*nのバッファメモリ6(Bij)
- バッファメモリ6およびコンピュータ4と通信する中央処理ユニット7(CPU)
- 中央処理ユニット7に接続されるメモリ8(MEM)
- 中央処理ユニット7に接続されるデジタル信号プロセッサ9(DSP)
- 2Dアレイ超音波プローブ2のn個の変換器Tijに個別に接続される、n*nのアナログ/デジタル変換器5(ADij)
- n*nのアナログ/デジタル変換器5にそれぞれ接続される、n*nのバッファメモリ6(Bij)
- バッファメモリ6およびコンピュータ4と通信する中央処理ユニット7(CPU)
- 中央処理ユニット7に接続されるメモリ8(MEM)
- 中央処理ユニット7に接続されるデジタル信号プロセッサ9(DSP)
図3は、合成開口手法における取得および撮像の例を示す。
超音波2Dアレイプローブ2は、撮像する領域10(たとえば、脳または脳の部分)に近い媒体の表面1(たとえば、生物体の頭の外面)に接触することができる。制御ユニット3によって制御される超音波2Dアレイは、変換器に、媒体中に少なくとも1つの超音波を送信させることができる(この例では、1度に1つの変換器)。超音波の各々の放出後に、N個の変換器が、媒体中の撮像する領域10で生成された後方散乱信号を受信するように構成することができる。たとえば、送信される超音波は、球形波であってもよく、撮像する全領域10をカバーすることが可能であってもよい。したがって、時刻t0において、変換器301が媒体10の中に球形波303を送信することができ、発生した後方散乱信号は、N個すべての変換器320によって受信することができる。t1(t0+Δt)において、変換器305が媒体中に球形波307を送信することができ、発生した後方散乱信号は、N個すべての変換器320によってtk(t0+k.Δt)に受信することができ、変換器N309が媒体10中に球形波311を送信することができ、発生した後方散乱信号は、N個すべての変換器320によって受信することができる、などである。
一実施形態によれば、超音波2Dアレイは、放出用に、K個の変換器毎にただ1つの変換器を使用する場合がある。1つの好ましい実施形態では、Kは、2以上であってもよい。たとえば、超音波2Dアレイは、放出において、2つの変換器毎(K=2)または4つの変換器毎(K=4)に使用することができる。1つまたは複数の実施形態では、超音波2Dアレイは、放出において、P個の変換器の少なくとも1つのクラスタが使用することができ、Pは、64と10000の間の変換器に含まれてもよい。したがって、この実施形態では、媒体中に超音波を送信するため各変換器を順に使用する従来型合成開口手法とは反対に、K個の変換器毎に発信することによって、合成開口手法でさえ、速い撮像率を保つことが可能になる。
変換器による各々の送信された信号について、送信された信号に応答して媒体で後方散乱信号が発生する場合がある。発生した後方散乱信号は、変換器301、305、309による各々の放出後の、撮像する領域10中の散乱についての情報を含む。発生した後方散乱信号は、N個すべての変換器320によって受信することができ、次いで、アナログデジタル変換器5によってそれぞれデジタル化されて、n個のバッファメモリ6中に記憶することができる。超音波(たとえば、球形波)の各々の発信後にn個のバッファメモリに記憶されるデータは、以降では、生の信号データと呼ばれる場合がある。
次いで、変換器による1回の放出後に得られる生の信号データは、1回の放出後の低解像度画像を形成するため、すべての画像点においてステアリングおよび合焦(すなわち、ビーム形成)することができる。したがって、たとえば、低解像度画像L-img_01は、変換器301によって送信された超音波信号、および、すべての変換器320によって受信される発生した後方散乱信号から作ることができる。
変換器301、305、309などによって順に、複数の低解像度画像L1、L2、Ljを得ることが可能であり、これらは、完全な高解像度画像
(たとえば、3D画像)を再構築するため、350で加算することができる。
合焦は、変換器によって受信された信号をコヒーレントに加算することによって、画像中の点毎に行うことができる。
超音波2Dアレイ中の放出の場所がわかっているため、放出から受信の正確な時間を計算して、受信における合焦で使用することができる。球形波の発生源から画像中の点へ、および受信変換器へと戻る距離は、以下であってもよい。
この式で、
が画像中の点に対応することができ、
が送信変換器の中間の場所であってもよく、
が受信変換器の中間の場所に対応することができる。したがって、単一の放出についての受信における合焦は、同相ですべての変換器からの受信信号g(t,i)を加算することによって実施することができ、
は、i番目の受信変換器の位置であってもよい。したがって、合焦した信号Lは、次式となってもよい。
Cが音の速度であり、iが受信変換器番号であり、Nrが受信している変換器の数である(Nr≦N)。これは、低解像度画像をもたらすために、結果として得られる画像中のあらゆる点
について行うことができる。あらゆる点
についての合焦は、前方遅延または後方遅延を考慮して実施することができる。前方遅延は、距離
を媒体中の音の速度Cで割ったものに対応することができ、後方遅延は、距離
を媒体中の音の速度Cで割ったものに対応することができる。こうして、各低解像度画像は、それぞれの低解像度画像のすべての点に適用される式(2)から得ることができる。
次いで、高解像度画像は、低解像度画像の加算によって得ることができる。したがって、高解像度画像
の1つの点
について、次式にしたがった、すべての放出(または、発信)についての信号を加えることによって、合成合焦を得ることができる。
ここで、Ntは、送信イベントの数であり、
は、高解像度画像が示される。
合成開口手法は、以下の文書にさらに詳細に記載される。
- [1] J.A. Jensen、S.I. Nikolov著、米国特許第6689063号、「Method and apparatus for acquiring images by recursive ultrasound images」
- [2] J.A. Jensen、S.I. Nikolov、K.L. Gammelmark、M.H. Pedersen著、「Synthetic Aperture Ultrasound Imaging」、Ultrasonics、Elsevier、The Netherlands、vol. 44、no.1、e5~e16頁, 2006年
- [1] J.A. Jensen、S.I. Nikolov著、米国特許第6689063号、「Method and apparatus for acquiring images by recursive ultrasound images」
- [2] J.A. Jensen、S.I. Nikolov、K.L. Gammelmark、M.H. Pedersen著、「Synthetic Aperture Ultrasound Imaging」、Ultrasonics、Elsevier、The Netherlands、vol. 44、no.1、e5~e16頁, 2006年
従来型合成開口撮像における基本的な問題は、画像における不十分な信号対雑音比である。というのは、単一の変換器が放出用に使用されるためである。これによって、従来の撮像で完全な開口を使用することと比較して、はるかに低い放出エネルギーがもたらされ、したがって、透過の深さが制限される。
合成開口撮像におけるそのような制限は、本開示に提案され下で記載されるような変換器の特定の寸法の使用によって克服することができる。
より正確には、各変換器は、2λよりも大きな直径を有する内接円を規定するアレイ表面に形状を有し、波長λは、超音波の周波数(たとえば、中心周波数)に対応する。
水中で(または、人体中で)伝播する1MHzの周波数を有する超音波では、波長は、1.54mmに対応する。
1つまたは複数の実施形態では、内接円の直径は、たとえば、2λと10λの間に含まれてもよい。
1つまたは複数の実施形態では、内接円の直径は、一般的に30マイクロメートルよりも大きく、たとえば、100マイクロメートルと2000マイクロメートルの間であってもよい。
1つまたは複数の実施形態では、内接円の直径よりもむしろ、表面アレイ中の変換器の対角の寸法を使用することが可能であり、対角は、変換器の中心点を通過する。そのような例では、変換器の対角の長さは、波長の2倍よりも大きくてもよい。
別の例によれば、変換器が円形形状を呈する場合に、アレイ表面の円形形状の直径が、内接円の直径に対応する。
こうして、より大きな変換器の(すなわち、より大きな表面積の)使用のおかげで、通常では技術的な視点から費用がかかり困難である、変換器の数を増やすことなく、大きな視野での3Dまたは4D超音波立体撮像を実現するのに必要な大きな表面(たとえば、10cm×10cmよりも大きい)を備えるプローブで、そのような変換器を覆う(または、カバーする)ことを可能にすることができる。加えて、より大きな変換器の使用によって、従来の変換器の寸法(たとえば、0.5×0.5λ、λ=250μm)と比較して、変換器の感度を向上させることが可能となり、上述の合成開口の制限を克服することがやはり可能となる。より大きな変換器の使用によって、より高いエネルギーを送信すること、したがって、媒体中の超音波の透過の深さを増やすことがやはり可能になる。
しかし、波長の2倍よりも大きな寸法の使用によって、各変換器の強い指向性がもたらされる可能性がある。たとえば、たとえば水中で1.54mmの波長λに対応する1MHzの周波数で、4λ×4λの寸法を有する正方形変換器、すなわち、4λの直径を有する内接円に対応するものについての指向性は、0度と45度の間に含まれてもよい。
したがって、オフアクシス後方散乱信号を検出することがより困難になる可能性があり、このことによって得られる3D画像の空間解像度を低下させる場合がある。
したがって、波長の2倍よりも大きな変換器の寸法の使用を補正するため、各変換器は、曲がった形状を有する(たとえば、凸)送信/受信面を呈する場合がある、または/および、曲がった形状を有する変換器をエミュレートするように設計され、(変換器の)送信/受信面上に配置される音響レンズを備える場合がある、または/および、変換器規模で曲がった変換器表面をエミュレートするために、独立に制御されるように適合される、変換器要素のマトリックスの特別配置に細分化する場合がある。
図4aから図4dは、変換器の送信/受信面用の曲がった形状の一例を概略的に図示する。
図4aは、曲がった形状401を有する受信(または、送信/受信)面を備える2Dアレイ2の正方形変換器400の概略斜視図に対応することができる。図4bは、図4aの上面図に対応することができ、2λよりも大きな直径を有する内接円を規定する変換器400の形状を図示する。図4cおよび図4dは、図4aで提示された変換器400の、平面XZおよびYZにおける概略断面図を表すことができる。
図4aを参照して、変換器400は、曲がった形状を有する受信面401を備えることができる。たとえば、この曲がった形状は、凸形状または凹形状であってもよい。
図4bを参照して、以前に説明したように、表面アレイ中の正方形変換器の内接円407の直径405は2λよりも大きな。正方形変換器の場合に、正方形変換器の寸法420、425は、内接円407の直径と等しくてもよい。
図4cを参照して、曲がった形状は、XZ平面における第1の湾曲410によって規定することができる。一実施形態では、第1の湾曲410の曲率半径は、λと40λの間に含まれてもよい。好ましい実施形態では、第1の湾曲の曲率半径は、2λと20λの間に含まれてもよい。
この第1の湾曲に加えて、図4dを参照して、送信/受信面の曲がった形状は、平面ZY中に第2の湾曲415を備えることができる。一実施形態では、第2の湾曲415の曲率半径は、λと40λの間に含まれてもよい。1つまたは複数の実施形態では、第2の湾曲の曲率半径は、2λと20λの間に含まれてもよい。1つまたは複数の実施形態では、第2の湾曲の曲率半径は、第1の湾曲の曲率半径と等しくてもよい。
変換器の受信(または、受信/送信)面のそのような性質によって、スネルの法則およびそれぞれの(変換器中および水などの液体中の)超音波速度にしたがって、発散レンズとして働くことにより、変換器の指向性を下げることを可能にすることができる。変換器の焦点、虚焦点を変換器の後に配置することができ、曲がった形状の曲率半径を選択することによって、調整することができる。
そのような性質によって、受信または放出において、焦点をぼかした超音波を得ることを可能にすることができる。
一例によれば、(たとえば、水中で)1.54mmの波長λに対応する1MHzの周波数で、4λの直径を有する内接円、2λの2つの曲率半径を(X軸の周りおよびY軸の周りに)有する凸に曲がった形状を規定する表面アレイ中の正方形変換器では、150度と180度の間に含まれる変換器の指向性を得るのが可能になることができる。
別の例によれば、水中で1.54mmの波長に対応する1MHzの周波数で、2λ以上の直径、4λの2つの曲率半径によって規定される凸に曲がった形状を有する受信面を有する円形変換器は、60度と120度の間に含まれる変換器の指向性を得るのが可能になることができる。
図5a、図5b、および図5cは、1つまたは複数の実施形態において、変換器の受信面上に配置され、曲がった形状を有する変換器をエミュレートするように設計される音響レンズの一例を概略的に図示する。
音響レンズは、たとえば、PDMS(ポリジメチルシロキサン)またはエポキシまたは音響レンズとして使用されるのに好適な任意の材料のような、少なくとも1つの材料で作ることができる。1つまたは複数の実施形態では、音響レンズは、複数の層でできてもよく、そのために各層は、それぞれの音速を有するそれぞれの材料でできている。
図5aは、発散レンズとして働く平凸音響レンズ500を有する2dアレイ2の円形変換器の概略斜視図に対応することができる。以前に説明したように、内接円の直径に対応するアレイ表面における円形変換器501の直径525は、2λよりも大きくてもよい。図5bおよび図5cは、図5aで提示された変換器および音響用平凸音響レンズの平面XZおよびYZにおける概略断面図を表すことができる。
図5bを参照すると、平凸音響レンズは、平面502および凸面503によって規定することができる。平凸レンズの平面502は、変換器501の受信(または送信/受信)面と接触してもよい。音響レンズは、変換器の受信面と完全にまたは部分的に一致する(または嵌合する)ように構成することができる。
一実施形態では、平凸音響レンズの平面は、変換器の形状と同じ形状を有することができる。たとえば、平凸レンズの平面は、長方形(または正方形または円形)形状を有する変換器について、長方形(または正方形または円形)であってもよい。
図5bおよび図5cを参照して、曲がった形状を有する受信面について前に記載したように、凸面503は、1つまたは複数の湾曲510、515によって規定することができる。
図5bを参照して、凸面503は、平面XZにおける第1の湾曲510によって規定することができる。1つまたは複数の実施形態では、第1の湾曲の曲率半径は、λと40λの間に含まれてもよい。
この第1の湾曲に加えて、図5cを参照して、凸面503は、平面ZYにおける第2の湾曲515を備えることができる。1つまたは複数の実施形態では、第2の湾曲の曲率半径は、λと40λの間に含まれてもよい。1つまたは複数の実施形態では、凸面503の各湾曲のそれぞれの半径が等しくてもよい。
1つまたは複数の実施形態では、凸面503の第1の湾曲または/および第2の湾曲の曲率半径は、2λと40λの間に含まれてもよい。
音響レンズ(たとえば、平凸音響レンズ)のそのような構成によって、2λと40λの間に含まれる変換器の後に配置される虚焦点を得ることを可能にすること、および、60度と150度の間に含まれる変換器の指向性を得ることを可能にすることができる。
一例によれば、液体(たとえば水)中で1.54mmの波長λに対応する1MHzの周波数で、2λよりも大きなアレイ表面における直径を有する円形変換器、および、3λの2つの曲率半径を有するPDMS中の音響平凸レンズによって、60度と150度の間に含まれる変換器の指向性を得るのを可能にすることができる。
1つまたは複数の実施形態では、音響レンズは、短い焦点を有する平凹レンズであってもよい。平凹レンズの焦点は、0.1mmと1cmの間に含まれてもよく、または/および、平凹レンズの少なくとも1つの曲率半径は、λと10λの間に含まれてもよい。短い焦点の使用によって、発散音響レンズのように働く音響レンズを得ることが可能である。というのは、平凹(すなわち、収束レンズ)によって合焦される超音波は、平凹レンズの焦点と交差した後に発散するためである。
そのような性質によって、受信または放出において、焦点をぼかした超音波を得ることを可能にすることができる。
1つまたは複数の実施形態では、独立に制御されるように適合された変換器要素のマトリックスを使用することによって、曲がった形状を有する変換器をエミュレートすることも可能にすることができる。変換器要素のマトリックスの各々の下位の変換器は、2Dアレイの変換器と同様の方法で、音信号(たとえば、超音波および/または後方散乱信号)を送信/受信するように構成することができる。変換器規模で曲がった形状を有する変換器、たとえば、曲面(たとえば、受信/送信面)を有する変換器をエミュレートするため、変換器のマトリックスの各変換器要素に、それぞれの電子的な時間遅延を適用することができる。1つまたは複数の実施形態では、変換器要素の表面は、0.1λ2とλ2の間に含まれてもよい。
それによって、図4aから図4dおよび図5aから図5cに提示された実施形態は、変換器の開口角を広げること(すなわち、指向性を低下させること)、したがって、より多くのオフアクシス後方散乱信号を集めることによって、横方向解像度および時間解像度を向上させることを可能にすることができる。より多くのオフアクシス後方散乱信号を集めるということによって、高解像度画像
(たとえば、高3D解像度画像)の各点について実施される合成合焦を改善することができる。
しかし、曲がった形状を有する変換器の使用または音響レンズ(たとえば、発散レンズ)の使用によって、送信において超音波信号が歪み、受信において合焦したビーム形成を構築するのに困難がもたらされる可能性がある。したがって、各送信後にN個の変換器によって受信される各後方散乱信号にそれぞれ遅延を加えることによって、曲がった形状の送信/受信面または音響レンズの使用を補償することが必要な場合がある。
図6は、1つまたは複数の実施形態中で、音響レンズについて、遅延を補正するための原理を図示する。
無補正遅延60は、撮像する領域と変換器601、602の中心との間の時間距離に対応してもよい。1つまたは複数の実施形態では、音響レンズについて補正した遅延は、送信変換器601から撮像する領域へ、および、撮像する領域から受信変換器602への異なる移動時間を考慮に入れることに対応することができる。
音響レンズを使用するときの補正遅延は、上で述べた式(1)を以下の式で置き換えることによって考慮に入れることができる。
この式(4)において、
が画像中の点に対応することができ、
が送信変換器曲率中心の位置であってもよく、
が送信ビーム603と音響レンズ601が交差する点Aの標準投影の送信変換器上の位置であってもよく、
が受信変換器曲率中心の位置であってもよく、
が受信ビーム604と音響レンズ602が交差する点Bの標準投影の受信変換器上の位置であってもよい。tltが、
からAへビーム(すなわち、超音波)が音響レンズ601中を移動する時間に対応することができ、ttが、Aから
へビームが媒体中を移動する時間に対応することができ、trが、
からBへビームが媒体中を移動する時間に対応することができ、tlrが、Bから
へビームが音響レンズ602中を移動する時間に対応することができる。
こうして、単一の発信についての合焦は、すべての変換器からの受信信号g(t,i)を同相で加算することによって実施することができ、
がi番目の受信変換器の位置であってもよい。合焦した信号Lは、次式となってもよい。
合成合焦は、次いで、上で述べた式(3)にしたがって得ることができる。
1つまたは複数の実施形態では、水中聴音器で遅延を測定することによって、実際の遅延を実験的に評価することもできる。
1つまたは複数の実施形態では、曲がった形状(たとえば、凸形状)を有する変換器について、tltおよびtlrは、上の式中で、ゼロと考えることができる。
図7は、3D超音波撮像のための受信における、2Dアレイの変換器の異なる構成のシミュレーションによる比較を図示する。
各画像703、705、707は、正方形形状の変換器についての異なる構成を有する表面アレイ(5×5cm)の、169個の正方形形状の変換器からの500個の散乱のイメージングシミュレーションの結果であってもよい(3Dファントムイメージングとも呼ばれる)。たとえば合成開口による、送信源は、1.03mmの波長λに対応する1.5MHzに設定することができ、9個の変換器によって実施することができる。これらの9個の変換器は、受信のために使用される2Dアレイからの変換器であってもよく、または、別の2Dアレイからの変換器であってもよい。
画像703は、正方形形状の変換器の2D超音波アレイを使用することによって得ることができ、正方形形状の変換器は、波長の4倍に等しい直径を有する内接円を規定することができる。次いで、変換器の寸法は、表面アレイ中で、4λ×4λに対応する。
画像705は、変換器の寸法が表面アレイ中で0.5λ×0.5λであってもよい、従来型変換器の2D超音波アレイを使用することによって、得ることができる。
画像707は、本開示の1つまたは複数の実施形態にしたがって、曲がった形状(凸形状)を有する受信面(または、送信/受信面)を備える正方形形状の変換器の2D超音波アレイを使用することによって得ることができる。より具体的には、正方形形状の変換器が、波長の4倍に等しい直径を有する内接円を規定し、受信面(または、送信/受信面)の曲率半径が3λである。
最終的に、グラフ711は、各々の前に述べた構成について、x=0およびz=3.8センチメートルに配置された散乱から来る後方散乱信号の、X軸にしたがった振幅を提示する。したがって、曲線720が画像707を与える変換器の構成に対応し、曲線730が画像703を与える変換器の構成に対応し、曲線740が画像705を与える変換器の構成に対応する。
グラフ711から、湾曲あり720または湾曲なし730で、変換器の寸法がより大きなものを使用すると、変換器740の従来型の寸法と比較して、感度を改善すること(すなわち、信号対雑音比を向上させること)を可能にすることができることに気づくことができる。
同様に、曲がった形状707を有する受信面を備える2Dアレイの変換器は、従来型の変換器705と比較して、または、同じ寸法だが曲がった形状の受信面を有さない変換器703と比較して、感度ならびにコントラストを改善すること(画像703、705、707上の白矢印)も可能にすることができる。
図8は、2Dアレイの変換器の受信面が曲がった形状を備えるときの、補正遅延の効果を図示する。
各画像803および805は、2λよりも大きな直径を有する、すなわち、2λよりも大きな内接円の直径に対応する、169個の円形変換器を備える表面アレイ(5×5cm)の同じ構成についてからの、800個の散乱のイメージングシミュレーションの結果であってもよい(3Dファントムイメージングとも呼ばれる)。合成開口による、送信源は、水中で1.04mmの波長λに対応する1.5MHzに設定することができ、9個の凸形円形変換器によって実施することができる。
画像803は、変換器で受信した後方散乱信号からもたらされる生の信号データに補正遅延を適用しないことによって得ることができる。画像805は、生の信号データに補正遅延を適用することによって得ることができる。一例によれば、曲がった形状(この場合、凸)を考慮に入れる補正遅延は、1マイクロ秒と50マイクロ秒の間に含まれてもよい。
最後に、グラフィック807は、画像803、805の両方について、x=0、z=3.8センチメートルに配置された散乱から来る後方散乱信号の中心振幅を提示する。
グラフ807から、生の信号データに補正遅延810を適用することによって、補正遅延815を適用しない同様の構成の凸形変換器との比較において、感度を改善すること(すなわち、信号対雑音比を向上させること)を可能にできることに気づくことができる。
同様に、画像803および805で気づくように、変換器の送信/受信面の凸形状を考慮に入れるため補正遅延を適用することによって、補正遅延を適用しない同様の構成の凸形変換器との比較において、感度ならびにコントラストを改善することを可能にできる(画像803および805上の白矢印)。
用途の例
上で提示された超音波3D撮像デバイスを有するプローブのそのような性質によって、3D(または4D)超音波撮像を実施するのを可能にすることができる。具体的には、たとえば、脳または心臓のように、骨の背後に配置される、または、皮膚表面から深くに配置される場合がある臓器全体の、高解像度での3D超音波撮像を実施することを可能にすることができる。実際に、それぞれ曲がった形状と組み合わされた大きな変換器の使用によって、高いエネルギー送信ならびに高い受信感度を保証することができる一方で、大きな開口、したがって、撮像する臓器の高解像度および大きな視野を保つことができる。
上で提示された超音波3D撮像デバイスを有するプローブのそのような性質によって、3D(または4D)超音波撮像を実施するのを可能にすることができる。具体的には、たとえば、脳または心臓のように、骨の背後に配置される、または、皮膚表面から深くに配置される場合がある臓器全体の、高解像度での3D超音波撮像を実施することを可能にすることができる。実際に、それぞれ曲がった形状と組み合わされた大きな変換器の使用によって、高いエネルギー送信ならびに高い受信感度を保証することができる一方で、大きな開口、したがって、撮像する臓器の高解像度および大きな視野を保つことができる。
たとえば、図9は、1つまたは複数の実施形態における3D超音波撮像デバイスの用途の第1の例を図示する。
本開示にしたがった3D超音波撮像デバイスのプローブ901は、患者の頭部902の表面に位置決めすることができる。1つまたは複数の実施形態では、プローブは、手動で配置すること、または、操作者によって制御されるロボットアームによって配置することができる。放出の周波数(または中心周波数)は、0.5MHzと3MHzの間に含まれてもよい。頭部の表面とプローブの間での、超音波の送信および受信を有利にするため、液体/ゲルを使用することができる。
そのような性質によって、大きな視野で、患者の脳の3Dまたは4D超音波撮像(または、超音波超高速撮像)を実施することを可能にすることができる。
さらに、プローブの場所を案内するロボットアームを使用して、脳の完全な3D画像を迅速で効果的な方法で再構築するために、頭部のいくつかの場所で撮像を実施することができる。
一実施形態によれば、プローブは、図8に提示されるように平面であってもよく、または、患者の頭部に配置されるヘルメット中に配置することができる。
一実施形態によれば、プローブは、平面または曲がった形状であってもよい。アレイが曲がった形状である場合、アレイの平面に対する上の言及は、撮像する領域と接触して配置されるように適合されるアレイ表面と理解するべきである。
1つまたは複数の実施形態によれば、超音波2Dアレイの表面は、視野の拡大を可能にする1000λ2よりも大きくてもよい。
1つまたは複数の実施形態では、たとえば、患者の皮膚の表面に、脳の撮像のため頭部上に、または心臓の撮像のため胸部上に位置決めすることができる。
1つまたは複数の実施形態では、プローブは、単に受信器として使用する場合がある。そのような実施形態では、第2のプローブを使用して、媒体中に超音波を発生させることができる。1つまたは複数の実施形態では、超音波は、光音響効果によって生成することができ、たとえば、媒体中で生成される超音波音響波のためレーザを使用することができる。そのような場合、プローブ901の変換器の送信/受信面は、後方散乱超音波を受診するためだけに使用することができる。
1つまたは複数の実施形態では、超音波3D撮像方法は、以下のように実行することができる。
- 各々の送信された超音波に対応して撮像する領域中を伝播する後方散乱超音波から信号を取得するステップ
- 取得信号に基づいて3D画像を生成するステップであって、前記後方散乱超音波が前記波長λを有する、ステップ
- 補正した取得信号を得るために、各送信された超音波の後に得られる取得信号に補正遅延を適用するステップであって、前記補正遅延が、送信面の曲がった形状を考慮に入れるステップ
- 動的合焦および3D画像の生成を実施するために補正した取得信号にビーム形成を適用するステップであって、ビーム形成が、各送信された波の受信において実施されて、コヒーレントコンパウンディングまたは合成開口撮像を実施することができる、ステップ
- 各々の送信された超音波に対応して撮像する領域中を伝播する後方散乱超音波から信号を取得するステップ
- 取得信号に基づいて3D画像を生成するステップであって、前記後方散乱超音波が前記波長λを有する、ステップ
- 補正した取得信号を得るために、各送信された超音波の後に得られる取得信号に補正遅延を適用するステップであって、前記補正遅延が、送信面の曲がった形状を考慮に入れるステップ
- 動的合焦および3D画像の生成を実施するために補正した取得信号にビーム形成を適用するステップであって、ビーム形成が、各送信された波の受信において実施されて、コヒーレントコンパウンディングまたは合成開口撮像を実施することができる、ステップ
超音波は、前に記載された合成開口法にしたがって、撮像する領域中に送信することができる。
生成される3D超音波画像は、たとえば、臓器全体の血管系の3D画像であってもよい。臓器全体の血管系の3D画像は、スーパーローカリゼーションアルゴリズム(または超音波ローカリゼーションアルゴリズム)を使用して、たとえば、撮像する臓器の血管系といった、撮像する領域の血流中に注入されたマイクロバブルを位置特定すること、または/および、追跡することによって得ることができる。マイクロバブルでスーパーローカリゼーションアルゴリズムをそのように使用することによって、高い空間解像度を有する血管系マップを生成することができる。実際に、超音波回折限界よりも1桁以上小さな10マイクロメートルほどの小ささ、および伝統的な周波数で制限される撮像深度よりはるかに大きな深さで、臓器の血管系の解像度を得ることを可能にすることができる。スーパーローカリゼーションアルゴリズムを使用して生成した3D超音波画像は、3D超解像超音波画像と呼ぶことができる。
そのようなスーパーローカリゼーションアルゴリズムは、以下のように実行することができる。
(a)以下すなわち、
- 組織から信号を除去してマイクロバブルだけを保つために、3D Bモード立体にSVDクラッタフィルタを適用するステップと、
- 各3D Bモード立体中の最大値を検出するステップと、
- 偽陽性マイクロバブルを除去するためにマイクロバブルの平均エネルギーレベルに基づいたエネルギーレベルしきい値を適用するステップと、
- 偽陽性マイクロバブルを除去するために3D点広がり関数(PSF)相関に基づいた相関値のしきい値を適用するステップと
を含む、マイクロバブルを検出するステップ
(b)以下すなわち、
- 前記マイクロバブルの最大値周辺のボクセルに3D放物面補間を使用することによりマイクロバブルの正確な中心を検出するステップと、
- 正確な中心の座標を使用して、フレーム毎にバブルを追跡するステップと、
- 追跡の位置および速度を回復するステップと
によってマイクロバブルを追跡するステップ
(c)以下すなわち、
- グリッド上に超解像立体を再構築するステップであって、超解像立体の密度が、グリッドの各ボクセル中のマイクロバブルの数であり、超解像立体の速度が、グリッドの各ボクセルにおける平均速度である、ステップ
によって、超解像立体を生成するステップ
(d)超解像立体に基づいて3D超解像画像を生成するステップ
(a)以下すなわち、
- 組織から信号を除去してマイクロバブルだけを保つために、3D Bモード立体にSVDクラッタフィルタを適用するステップと、
- 各3D Bモード立体中の最大値を検出するステップと、
- 偽陽性マイクロバブルを除去するためにマイクロバブルの平均エネルギーレベルに基づいたエネルギーレベルしきい値を適用するステップと、
- 偽陽性マイクロバブルを除去するために3D点広がり関数(PSF)相関に基づいた相関値のしきい値を適用するステップと
を含む、マイクロバブルを検出するステップ
(b)以下すなわち、
- 前記マイクロバブルの最大値周辺のボクセルに3D放物面補間を使用することによりマイクロバブルの正確な中心を検出するステップと、
- 正確な中心の座標を使用して、フレーム毎にバブルを追跡するステップと、
- 追跡の位置および速度を回復するステップと
によってマイクロバブルを追跡するステップ
(c)以下すなわち、
- グリッド上に超解像立体を再構築するステップであって、超解像立体の密度が、グリッドの各ボクセル中のマイクロバブルの数であり、超解像立体の速度が、グリッドの各ボクセルにおける平均速度である、ステップ
によって、超解像立体を生成するステップ
(d)超解像立体に基づいて3D超解像画像を生成するステップ
図10は、1つまたは複数の実施形態において、スーパーローカリゼーションアルゴリズムを使用することにより得られる3D超解像超音波撮像の例を図示する。
図10は、0.3λ×0.3λの各変換器用のサイズを有する従来型3D超音波撮像プローブ(10cm×10cm)変換器を使用することと、マイクロバブルおよびスーパーローカリゼーションアルゴリズムの使用を組み合わせることによって得られた、分離されて脈動するラットの心臓についての、3D超音波冠動脈微小循環撮像に対応する。
図10上で3Dを実施するため、高い体積/率で体積全体にわたって3D画像を再構築するために、平面波が高い繰返し率でマトリックス変換器によって放出された。
図10上で見ることができるように、心臓の血流中へのマイクロバブルの注入は、次いで位置特定および/または追跡されるが、10μmよりも小さな解像度で心臓の血管系の高解像度撮像を得ることを可能にすることができる。高い繰返し率のために、心臓の血管系中の血液循環を追従することを可能にすることができる。
さらに、10μmよりも小さな解像度のおかげで、血管系中で高い精度で潜在的な疾患を位置特定することを可能にすることができるが、伝統的な3D撮像技法で識別されなかった早期の疾患を識別することも可能にすることができる。
図11は、媒体中に超音波を送信する発散レンズを有する変換器の実験的な例を図示する。
図11では、たとえば、発散レンズといった、曲がった形状と組み合わされたときの、変換器の指向性への影響が図示される。
図11の画像Aは、変換器1110が発散レンズ(すなわち、音響レンズ)とは組み合わされずに媒体中に超音波を送信する場合を図示する。画像Aは、媒体中で、X軸が横の次元に対応し、Z軸が深さに対応する、XZウィンドウにしたがった媒体における切断面を提示する。
図11の画像Bは、変換器1120が発散レンズ(すなわち、音響レンズ)1120aと組み合わされて媒体中に超音波を送信する場合を図示する。画像Bは、媒体中で、X軸が横の次元に対応し、Z軸が深さに対応する、XZウィンドウにしたがった媒体における切断面を提示する。
図11からの各々の画像AおよびBは、媒体中で送信された超音波に応答して、各空間点で水中聴音器によって記録された、送信された最大圧力信号に基づいて得ることができる。超音波の周波数は1MHzである。任意単位(u.a.)のスケールは、たとえば、デシベルでの超音波の強度を表すことができる。
画像AおよびBから、発散レンズ1120の使用によって、変換器の開口角を広げること(すなわち、指向性を低下させること)、すなわち、横方向寸法を広げることを可能にすることができ、したがって、受信または/および放出において、超音波の焦点をぼかすことを可能にすることができることに気づくことができる。
図12は、発散レンズを有する変換器を有する超音波プローブ、および、チューブの中のバブルを検出および追跡するために超音波ローカリゼーション顕微鏡検査アルゴリズムを使用してそのような超音波プローブを使用することによって得られるチューブの中のマイクロバブルの画像を図示する。
図12の画像A中に提示される超音波プローブ1210は、図1および図2で提示された超音波撮像のための装置で使用するように構成および配置することができる。
超音波プローブ1210は、たとえば、16個の変換器(または、超音波変換器)といった複数の変換器を有するマトリックスアレイを備えることができる。マトリックスアレイの各変換器は、発散レンズ1220と組み合わせることができる。超音波プローブは、図1および図2で提示された超音波撮像のための装置に接続することができる。
そのような超音波プローブと組み合わせた超音波撮像のための装置は、1.5mmである音の波長よりも小さなチューブの中に含まれる液体(たとえば、水)の中のマイクロバブルを撮像および追跡するために使用することができる。マイクロバブル追跡後のチューブの画像が、図12の画像Bに提示される。チューブおよびその直径は、血管系の血管を表すことができる。画像Bから、音の波長よりも小さなチューブ中のマイクロバブル1250を追跡および撮像して高解像度画像を得ることが可能であり、したがって、血管系(または、脈管構造)の血管中で高解像度(典型的には、波長/10、たとえば、150μm未満)を得ることを可能にすることができることを示し、したがって、血管系中で高い精度で潜在的な疾患を位置特定することを可能にすることができるが、伝統的な3D撮像技法で識別されなかった早期の疾患を識別することもできることを理解することができる。方法は、血管中の流速を評価し、微小血管網の機能情報を得るために使用することもできる。
2 2Dアレイ超音波プローブ、超音波3D撮像プローブ
3 制御ユニット
4 コンピュータ
5 アナログ/デジタル変換器
6 バッファメモリ
7 中央処理ユニット
8 メモリ
9 デジタル信号プロセッサ
10 撮像する領域、媒体
301 変換器
303 球形波
305 変換器
307 球形波
309 変換器
311 球形波
320 変換器
400 正方形変換器
401 曲がった形状、受信面
405 直径
407 内接円
410 第1の湾曲
415 第2の湾曲
420 寸法
425 寸法
500 平凸音響レンズ
501 円形変換器
502 平面
503 凸面、曲がった形状
510 第1の湾曲
515 第2の湾曲
525 直径
600 無補正遅延
601 送信変換器、音響レンズ
602 受信変換器、音響レンズ
603 送信ビーム
604 受信ビーム
703 画像、変換器
705 画像、従来型の変換器
707 画像、曲がった形状
711 グラフ
720 曲線、湾曲あり
730 曲線、湾曲なし
740 曲線、変換器
803 画像
805 画像
807 グラフィック
810 補正遅延
815 補正遅延
901 プローブ
902 頭部
1110 変換器
1120 変換器
1120a 発散レンズ、音響レンズ
1210 超音波プローブ
1220 発散レンズ
1250 マイクロバブル
3 制御ユニット
4 コンピュータ
5 アナログ/デジタル変換器
6 バッファメモリ
7 中央処理ユニット
8 メモリ
9 デジタル信号プロセッサ
10 撮像する領域、媒体
301 変換器
303 球形波
305 変換器
307 球形波
309 変換器
311 球形波
320 変換器
400 正方形変換器
401 曲がった形状、受信面
405 直径
407 内接円
410 第1の湾曲
415 第2の湾曲
420 寸法
425 寸法
500 平凸音響レンズ
501 円形変換器
502 平面
503 凸面、曲がった形状
510 第1の湾曲
515 第2の湾曲
525 直径
600 無補正遅延
601 送信変換器、音響レンズ
602 受信変換器、音響レンズ
603 送信ビーム
604 受信ビーム
703 画像、変換器
705 画像、従来型の変換器
707 画像、曲がった形状
711 グラフ
720 曲線、湾曲あり
730 曲線、湾曲なし
740 曲線、変換器
803 画像
805 画像
807 グラフィック
810 補正遅延
815 補正遅延
901 プローブ
902 頭部
1110 変換器
1120 変換器
1120a 発散レンズ、音響レンズ
1210 超音波プローブ
1220 発散レンズ
1250 マイクロバブル
Claims (15)
- 変換器を有するアレイを備える超音波3D撮像プローブ(2)であって、各変換器が、少なくとも、媒体中の撮像する領域(10)から前記媒体における超音波放出に続く後方散乱超音波を取得するように構成され、前記後方散乱超音波が前記媒体中で伝播したときのある波長λに対応するある中心周波数を前記後方散乱超音波が有し、前記アレイが前記媒体と接触して配置されるように適合されるアレイ表面を形成し、各変換器が少なくとも前記後方散乱超音波を受け取るための受信面を有し、
各変換器が、2λよりも大きな直径を有する内接円を規定する、前記アレイ表面における形状を有し、各変換器が
- 前記受信面に垂直な少なくとも1つの平面中に湾曲を有する前記受信面の曲がった形状(401)を有し、
- および/または、前記受信面に垂直な少なくとも1つの平面中に、湾曲を有する曲がった形状(503)を有する変換器をエミュレートするように設計された音響レンズと嵌合され、
- および/または、独立に制御されるように適合された変換器要素のマトリックスによって形成される、超音波3D撮像プローブ(2)。 - 各変換器の前記受信面が曲がった形状を有し、前記湾曲が、λと40λの間に含まれる少なくとも1つの曲率半径を有し、
または、各変換器の前記受信面が曲がった形状を有し、前記曲がった形状が凸である、請求項1に記載の超音波3D撮像プローブ。 - 各変換器の前記受信面が曲がった形状を有する変換器をエミュレートするように設計された音響レンズと嵌合され、前記音響レンズが、前記変換器の前記受信面の反対側に外面を有し、前記外面が凸または凹である、請求項1に記載の超音波3D撮像プローブ。
- 前記外面が凸であり、前記外面に垂直な少なくとも1つの平面に湾曲を有し、前記湾曲が、λと40λの間に含まれる少なくとも1つの曲率半径を有し、
または、前記外面が凹であり、前記外面に垂直な少なくとも1つの平面に湾曲を有し、前記湾曲が、λと40λの間に含まれる少なくとも1つの曲率半径を有する、請求項3に記載の超音波3D撮像プローブ。 - 前記音響レンズが、それぞれがいくつかの音速に関連するいくつかの材料を含む複合音響レンズである、請求項1に記載の超音波3D撮像プローブ。
- 前記直径が30マイクロメートルよりも大きい、請求項1に記載の超音波3D撮像プローブ。
- 少なくとも、
- 請求項1から6のいずれか一項に記載の超音波3D撮像プローブ(2)と、
- 各送信された超音波に応答して前記撮像する領域(10)から前記変換器によって受信された前記後方散乱超音波から信号を取得するために構成される制御ユニット(3)であって、取得信号に基づいて3D画像を生成するためさらに構成される、制御ユニット(3)と
を備える、超音波3D撮像デバイス。 - 前記制御ユニットが、前記後方散乱超音波を生成するため前記撮像する領域中に前記変換器に超音波を送信させるようにさらに構成される、請求項7に記載の超音波3D撮像デバイス。
- 前記制御ユニットが、各送信された超音波の送信後に得られる前記取得信号に補正遅延(603,604)を適用するようにさらに構成され、前記補正遅延が、送信面の前記曲がった形状を考慮に入れる、請求項8に記載の超音波3D撮像デバイス。
- 前記超音波3D撮像プローブの各変換器が、独立に制御されるように適合された変換器要素のマトリックスによって形成され、前記制御ユニットが、変換器規模で曲がった形状を有する変換器をエミュレートするため変換器の前記マトリックスの前記変換器要素に時間遅延を適用するように構成される、請求項7に記載の超音波3D撮像デバイス。
- 撮像する領域中を伝播した後方散乱超音波を取得するための請求項1から6のいずれか一項に記載の超音波3D撮像プローブを使用するステップを含む超音波撮像方法であって、前記後方散乱超音波が前記波長λを有する、超音波撮像方法。
- 各送信された超音波に応答して撮像する領域中を伝播した後方散乱超音波からの信号を取得するため、および、前記取得信号に基づいて3D画像を生成するための請求項7に記載の超音波3D撮像デバイスを使用するステップを含む超音波撮像方法であって、前記後方散乱超音波が前記波長λを有する、超音波撮像方法。
- 前記超音波が合成開口法にしたがって前記撮像する領域中に送信される、請求項12に記載の、超音波3D撮像デバイスを使用するステップを含む、超音波撮像方法。
- - 補正した取得信号を得るために、各送信された超音波の後に得られる前記取得信号に補正遅延を適用するステップであって、前記補正遅延が、前記送信面の前記曲がった形状を考慮に入れる、ステップと、
- 動的合焦および前記3D画像の生成を実施するために前記補正した取得信号にビーム形成を適用するステップであって、前記ビーム形成が各送信された波の受信で実施され、コヒーレントコンパウンディングまたは合成開口撮像を実施する、ステップと
をさらに含む、請求項12に記載の超音波3D撮像デバイスを使用するステップを含む、超音波撮像方法。 - マイクロバブルが、前記撮像する媒体中に存在し、スーパーローカリゼーションアルゴリズムを使用して前記3D画像中で位置特定および/または追跡され、
前記スーパーローカリゼーションアルゴリズムが、
(a)
- 組織から信号を除去してマイクロバブルだけを保つために、3D Bモード立体にSVDクラッタフィルタを適用するステップ、
- 各3D Bモード立体中の最大値を検出するステップ、
- 偽陽性マイクロバブルを除去するためにマイクロバブルの平均エネルギーレベルに基づいたエネルギーレベルしきい値を適用するステップ、および
- 偽陽性マイクロバブルを除去するために3D点広がり関数(PSF)相関に基づいた相関値のしきい値を適用するステップ
を含む、前記マイクロバブルを検出するステップと、
(b)
- 前記マイクロバブルの最大値周辺のボクセルに3D放物面補間を使用することによりマイクロバブルの正確な中心を検出するステップ、
- 正確な中心の座標を使用して、フレーム毎にバブルを追跡するステップ、および
- 追跡の位置および速度を回復するステップ
によって前記マイクロバブルを追跡するステップと、
(c)
- グリッド上に超解像立体を再構築するステップであって、前記超解像立体の密度が、前記グリッドの各ボクセル中のマイクロバブルの数であり、前記超解像立体の速度が、前記グリッドの各ボクセルにおける平均速度である、ステップ
によって、超解像立体を生成するステップと、
(d)前記超解像立体に基づいて3D超解像画像を生成するステップと、
を含む、請求項11に記載の超音波撮像方法。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP21305434 | 2021-04-02 | ||
EP21305434.9 | 2021-04-02 | ||
PCT/EP2022/058808 WO2022207928A1 (en) | 2021-04-02 | 2022-04-01 | Divergent lens array |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2024511870A true JP2024511870A (ja) | 2024-03-15 |
Family
ID=75588164
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2023560710A Pending JP2024511870A (ja) | 2021-04-02 | 2022-04-01 | 発散レンズアレイ |
Country Status (7)
Country | Link |
---|---|
US (1) | US20240180524A1 (ja) |
EP (1) | EP4312794A1 (ja) |
JP (1) | JP2024511870A (ja) |
KR (1) | KR20230164169A (ja) |
CN (1) | CN117320635A (ja) |
CA (1) | CA3214124A1 (ja) |
WO (1) | WO2022207928A1 (ja) |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ATE245301T1 (de) | 1999-05-10 | 2003-08-15 | B K Medical As | Rekursive ultraschallabbildung |
US7285094B2 (en) * | 2002-01-30 | 2007-10-23 | Nohara Timothy J | 3D ultrasonic imaging apparatus and method |
US7263888B2 (en) * | 2003-10-16 | 2007-09-04 | General Electric Company | Two dimensional phased arrays for volumetric ultrasonic inspection and methods of use |
EP2701607B1 (en) * | 2011-04-26 | 2020-09-16 | University Of Virginia Patent Foundation | Bone surface image reconstruction using ultrasound |
US9883848B2 (en) * | 2013-09-13 | 2018-02-06 | Maui Imaging, Inc. | Ultrasound imaging using apparent point-source transmit transducer |
FR3051693B1 (fr) * | 2016-05-31 | 2018-05-11 | Imasonic | Reseau d'elements transducteurs ultrasonores |
-
2022
- 2022-04-01 KR KR1020237037830A patent/KR20230164169A/ko unknown
- 2022-04-01 WO PCT/EP2022/058808 patent/WO2022207928A1/en active Application Filing
- 2022-04-01 CA CA3214124A patent/CA3214124A1/en active Pending
- 2022-04-01 CN CN202280036059.9A patent/CN117320635A/zh active Pending
- 2022-04-01 JP JP2023560710A patent/JP2024511870A/ja active Pending
- 2022-04-01 US US18/553,591 patent/US20240180524A1/en active Pending
- 2022-04-01 EP EP22720417.9A patent/EP4312794A1/en active Pending
Also Published As
Publication number | Publication date |
---|---|
KR20230164169A (ko) | 2023-12-01 |
CA3214124A1 (en) | 2022-10-06 |
CN117320635A (zh) | 2023-12-29 |
WO2022207928A1 (en) | 2022-10-06 |
US20240180524A1 (en) | 2024-06-06 |
EP4312794A1 (en) | 2024-02-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10653392B2 (en) | Ultrasound imaging using apparent point-source transmit transducer | |
CN108778530B (zh) | 具有稀疏阵列探测器的超声成像 | |
OˈReilly et al. | A super‐resolution ultrasound method for brain vascular mapping | |
Montaldo et al. | Coherent plane-wave compounding for very high frame rate ultrasonography and transient elastography | |
US4183249A (en) | Lens system for acoustical imaging | |
CA2658063C (en) | Method and apparatus for ultrasound synthetic imaging | |
JP6489797B2 (ja) | 被検体情報取得装置 | |
Jensen | Ultrasound imaging and its modeling | |
JP3868724B2 (ja) | 超音波血管内視鏡システム | |
Diamantis et al. | Resolving ultrasound contrast microbubbles using minimum variance beamforming | |
KR20180013956A (ko) | 단일 추적 위치 전단파 탄성 이미징을 위한 방법, 시스템 및 컴퓨터 프로그램 제품 | |
CN107205720B (zh) | 一种超声适应性波束形成方法及其对经颅成像的应用 | |
JP2007513672A (ja) | 2次元アレイトランスデューサを用いる立体的超音波画像化システム | |
JP2017003587A (ja) | ハイブリッド光学音響断層撮影及び超音波検査に関するデバイスおよび方法 | |
Nabavizadeh et al. | Multi-source and multi-directional shear wave generation with intersecting steered ultrasound push beams | |
JP2024511870A (ja) | 発散レンズアレイ | |
KR20150118732A (ko) | 초음파 장치 및 그 제어 방법 | |
CN117158911B (zh) | 一种多声速自适应光声层析图像重建方法 | |
Caudoux et al. | Curved toroidal row column addressed transducer for 3D ultrafast ultrasound imaging | |
EP4394441A1 (en) | A beamforming method and system | |
Aubry et al. | Skull surface detection algorithm to optimize time reversal focusing through a human skull | |
Gopalakrishnan et al. | Extending the Capability of Linear Array Ultrasound Probe to Concave Array Using Low-Cost Acoustic Lens for High Frame Rate Focused Imaging | |
Dalgarno et al. | Development of Super-resolution Sharpness-based Axial Localization for Ultrasound Imaging | |
Park | Integrated system for ultrasonic, elasticity and photoacoustic imaging | |
Omidvar et al. | Shape Estimation of Ultrasound Arrays Using Spatial Coherence: A Preliminary Study |