WO2010017508A1 - Device, system, and method to measure abdominal aortic aneurysm diameter - Google Patents

Device, system, and method to measure abdominal aortic aneurysm diameter Download PDF

Info

Publication number
WO2010017508A1
WO2010017508A1 PCT/US2009/053201 US2009053201W WO2010017508A1 WO 2010017508 A1 WO2010017508 A1 WO 2010017508A1 US 2009053201 W US2009053201 W US 2009053201W WO 2010017508 A1 WO2010017508 A1 WO 2010017508A1
Authority
WO
WIPO (PCT)
Prior art keywords
ultrasound
scan
transceiver
availability
abdominal aorta
Prior art date
Application number
PCT/US2009/053201
Other languages
French (fr)
Inventor
Jongtae Yuk
Gerald Mcmorrow
Original Assignee
Verathon Inc.
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 Verathon Inc. filed Critical Verathon Inc.
Priority to CA2732997A priority Critical patent/CA2732997C/en
Priority to EP09805626.0A priority patent/EP2323559A4/en
Priority to JP2011522292A priority patent/JP5658151B2/en
Publication of WO2010017508A1 publication Critical patent/WO2010017508A1/en

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/52023Details of receivers
    • G01S7/52036Details of receivers using analysis of echo signal for target characterisation
    • G01S7/52038Details of receivers using analysis of echo signal for target characterisation involving non-linear properties of the propagation medium or of the reflective target
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/02007Evaluating blood vessel condition, e.g. elasticity, compliance
    • A61B5/02014Determining aneurysm
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0858Detecting organic movements or changes, e.g. tumours, cysts, swellings involving measuring tissue layers, e.g. skin, interfaces
    • 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/52079Constructional features
    • G01S7/52082Constructional features involving a modular construction, e.g. a computer with short range imaging equipment
    • 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/52098Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging related to workflow protocols
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/12Edge-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • G06T7/62Analysis of geometric attributes of area, perimeter, diameter or volume
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/107Measuring physical dimensions, e.g. size of the entire body or parts thereof
    • A61B5/1075Measuring physical dimensions, e.g. size of the entire body or parts thereof for measuring dimensions by non-invasive methods, e.g. for determining thickness of tissue layer
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4444Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to the probe
    • A61B8/4455Features of the external shape of the probe, e.g. ergonomic aspects
    • 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10132Ultrasound image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30101Blood vessel; Artery; Vein; Vascular

Definitions

  • the aorta artery in the abdomen carries blood from the heart to the abdominal region.
  • One disorder of the abdominal aorta is known as an abdominal aortic aneurysm, which is a permanent localized dilation of the arterial wall of the abdominal aorta.
  • an abdominal aortic aneurysm which is a permanent localized dilation of the arterial wall of the abdominal aorta.
  • dilation of the arterial wall is greater than 1.5 times the typical, i.e. nominal, diameter, it is referred to as an aneurysm.
  • FIGURE 1-1 shows a normal abdominal aorta at.
  • FIGURE 1-2 shows a typical aortic aneurysm at.
  • An aortic aneurysm is usually located below the renal arteries and the kidney arteries and above the aorta-iliac bifurcation. Below the aortic-iliac bifurcation are additional arteries. Abdominal aortic aneurysms are a fairly common disorder, occurring in approximately 5- 7% of the population over age 60. Since effective screening programs are not yet established, the diagnosis of AAA is frequently made at the time of rupture or impending rupture, which leads to a dramatic increase of post-operative morbidity. (Daly, et al., 2004)
  • the accuracy is within three millimeters of the actual size of the aneurysm, using conventional CT or MRI processing.
  • CT or MRI processing are very expensive, both to purchase/lease and to maintain.
  • a trained sonographer is necessary to interpret the results of the scans. As a result, many aneurysms go undetected and/or are not consistently monitored after discovery, until rupture resulting in death of the patient.
  • FIGURE 1-1 shows a normal abdominal aorta
  • FIGURE 1-2 shows a typical abdominal aortic aneurysm
  • FIGURES 2-1 thru 2-4 depict a partial schematic and a partial isometric view of a transceiver, a scan cone comprising a rotational array of scan planes, and a scan plane of the array of an ultrasound harmonic imaging system;
  • FIGURE 3-1 is a side elevation view of an ultrasound transceiver 1OA that includes an inertial reference unit
  • FIGURE 3-2 illustrates a side and partial isometric that schematically depicts an ultrasound scanner employing C-mode and B-mode ultrasound modalities
  • FIGURE 3-3 illustrates a partial isometric and schematic view of an ultrasound scanner system
  • FIGURE 4 is a schematic illustration of a server-accessed local area network and/or internet in communication with a plurality of ultrasound imaging systems
  • FIGURE 5 shows exemplary probe locations along the midline of a patient's torso.
  • FIGURE 6-1 is a B-mode image of an AAA is in the same field as a shadow region ;
  • FIGURES 6-2 and 6-3 are histograms of the A-mode lines 65 (FIGURE 6-2) and 65 (FIGURE 6-3);
  • FIGURES 7-1 thru 7-12 are twelve B-mode images showing shadows and non-availability
  • FIGURE 8 schematically illustrates an availability plot in C-mode (view from the top of the scan cone);
  • FIGURE 9-1 and 9-2 illustrates an aiming direction scheme; and a decision tree for the arrow feedback from the indicator 22 of FIGURE 3-3;
  • FIGURE 10 illustrates a screen shot of an exemplary aiming guide
  • FIGURE 11 depicts a flowchart of a aortic aneurysm detection algorithm employing fundamental ultrasound energies
  • FIGURES 12-1 and 12-2 are exemplary B-mode images of longitudinal section and cross section of AAA phantom with their typical segmentations (red dotted lines) for volume estimation;
  • FIGURES 12-3 is a schematic diagram illustrating an exemplary embodiment of the limited segmentation
  • FIGURE 13 is a graphic diagram for the conversion from volume to diameter assuming an object shape is sphere.
  • FIGURE 14 depicts a flowchart of the Find Initial Walls sub-algorithm of FIGURE 11.
  • Embodiments of the invention may be described in the general context of computer-executable instructions, such as program modules, executed by one or more computer processors or other devices, computer-readable media on which such
  • program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types.
  • functionality of the program modules may be combined or distributed as desired in various embodiments.
  • Embodiments of the invention may include or otherwise utilize at least some form of computer readable media, which may be associated with one or more processors and/or memory devices.
  • Computer readable media can be any available media that can be accessed by one or more components of such operating environment.
  • Computer readable media may comprise computer storage media and communication media.
  • Computer storage media includes volatile and nonvolatile, removable and non-removable media implemented in any method or technology for storage of information such as computer readable instructions, data structures, program modules or other data.
  • Computer storage media includes, but is not limited to, RAM, ROM, EEPROM, flash memory or other memory technology, CD- ROM, digital versatile disks (DVD) or other optical storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other medium which can be used to store the desired information and which can be accessed by one or more components of such operating environment.
  • Communication media typically embodies computer readable instructions, data structures, program modules or other data in a modulated data signal such as a carrier wave or other transport mechanism and includes any information delivery media.
  • modulated data signal means a signal that has one or more of its characteristics set or changed in such a manner as to encode information in the signal.
  • communication media includes wired media such as a wired network or direct-wired connection, and wireless media such as acoustic, RF, infrared and other wireless media. Combinations of the any of the above should also be included within the scope of computer readable media.
  • wired media such as a wired network or direct-wired connection
  • wireless media such as acoustic, RF, infrared and other wireless media. Combinations of the any of the above should also be included within the scope of computer readable media.
  • Disclosure below includes systems and methods to detect and measure an AAA involving transmitting ultrasound energy having at least one of a fundamental frequency to the AAA, collecting ultrasound echoes returning from the AAA and generating signals from the ultrasound echoes, and identifying within the ultrasound signals those attributable to fundamental ultrasound frequencies. Thereafter, the fundamental-frequency-derived signals undergo signal processing via computer executable program instructions to present an image of the AAA on a display, and calculate the volume and diameter of the AAA.
  • the signal processing applied to the transceiver echoic fundamental ultrasound signals include an algorithm having computer readable instructions for ascertaining the certainty that a given scan line traverses a AAA, a shadow region, or both a AAA and a shadow region using a grading algorithm for predicting the scan line's AAA or shadow classification.
  • the ultrasound transceivers or distal collection devices are capable of collecting in vivo three-dimensional (3-D) cone-shaped ultrasound images of a patient.
  • a radio frequency pulsed ultrasound field is transmitted into the body, and the back-scattered "echoes" may be transducer-detected and presented as a one-dimensional (1-D) voltage trace, which may be also referred to as a radio frequency (RF) line.
  • RF radio frequency
  • RF radio frequency
  • the system includes an ultrasound transceiver positioned to deliver ultrasound energy and receive echoes of the ultrasound energy across a plurality of scan planes; an algorithm configured to signal process the received echoes and characterize detected signals across a plurality of scan planes, wherein a "percentage of availability" measurement may be made. Percentage of availability (POA) may be based on whether relevant information is contained within the scanplane as compared with shadows or other regions not of interest.
  • the system includes a display for presenting a visual depiction of availability contained in a plurality of scan planes; and display for indicating positioning information of the ultrasound transceiver based on the visual depiction of availability.
  • FIGURES 2-1- 2-4 depict a partial schematic and a partial isometric view of a transceiver, a scan cone comprising a rotational array of scan planes, and a scan plane of the array of various ultrasound systems capable of collecting RF line analysis.
  • FIGURE 2-1 is a side elevation view of an ultrasound transceiver 1OA that includes an inertial reference unit 38, according to an embodiment of the invention.
  • the transceiver 1OA includes a transceiver housing 35 having an outwardly extending handle 33 suitably configured to allow a user to manipulate the transceiver 1OA relative to a patient.
  • Ultrasound transducers operating within the transceiver 1OA can be equipped to collect and ready signals for ultrasound fundamental and/or harmonic frequency analysis.
  • the handle 33 includes a trigger 37 that allows the user to initiate an ultrasound scan of a selected anatomical portion.
  • the transceiver 1OA also includes a transceiver dome 34 that contacts a surface portion of the patient when the selected anatomical portion is scanned (See FIGURE 5).
  • the dome 34 generally provides an appropriate acoustical impedance match to the anatomical portion and/or permits ultrasound energy to be properly focused as it is projected into the anatomical portion.
  • the transceiver 1OA further includes one transducer element, or preferably an array of separately excitable ultrasound transducer elements, positioned within or otherwise adjacent to the housing 35.
  • the transducer elements can be suitably positioned within the housing 35 or otherwise to project ultrasound energy outwardly from the dome 34, and to permit reception of acoustic reflections generated by internal structures within the anatomical portion.
  • the one or more array of ultrasound elements can include a one- dimensional, or a two-dimensional array of piezoelectric elements that can be moved within the housing 35 by a motor.
  • the array can be stationary with respect to the housing 35 so that the selected anatomical region can be scanned by selectively energizing the elements in the array.
  • a directional indicator panel or aiming guide panel 32 includes a plurality of arrows that can be illuminated for initial targeting and guiding a user to access the targeting of an organ or structure within a region of interest (ROI).
  • ROI region of interest
  • the directional indicator panel 32 of transceiver 1OA has a virtual equivalent in the form of a targeting icon screenshot 77B that appears on an display 76 on the console 74; both the indicator panel 32 and display 76 (displaying targeting icon 77B) can function to guide a transceiver user to place the transceiver to obtain a "good scan" of the abdominal aortic region of interest (e.g., for the detection of AAA).
  • the directional arrows will be not illuminated. If the AAA is off-center or a shadow appears in the field of interest, an arrow or set of arrows can be illuminated to direct the user to reposition the transceiver 1OA, 1OC acoustically at a second or subsequent dermal location of the subject.
  • the acoustic coupling can be achieved by liquid sonic gel applied to the skin of the patient or by sonic gel pads against which the transceiver dome 34 is placed.
  • the directional indicator panel 32 can also be presented on the display 54 of computer 52 in imaging subsystems described in FIGURES 4 below, and/or presented on the transceiver display 36.
  • transceiver 1OA can, in one embodiment, include an inertial reference unit that includes an accelerometer and/or gyroscope 38 positioned preferably within or adjacent to housing 35.
  • the accelerometer 38 can be operable to sense an acceleration of the transceiver 1OA, preferably relative to a coordinate system, while the gyroscope can be operable to sense an angular velocity of the transceiver 1OA relative to the same or another coordinate system.
  • the gyroscope can be of conventional configuration that employs dynamic elements, or it can be an optoelectronic device, such as the known optical ring gyroscope.
  • the accelerometer and the gyroscope can include a commonly packaged and/or solid-state device.
  • One suitable commonly packaged device can be the MT6 miniature inertial measurement unit, available from Omni Instruments, Incorporated, although other suitable alternatives exist.
  • the accelerometer and/or the gyroscope can include commonly packaged micro-electromechanical system (MEMS) devices, which are commercially available from MEMSense, Incorporated.
  • MEMS micro-electromechanical system
  • the accelerometer and the gyroscope cooperatively permit the determination of positional and/or angular changes relative to a known position that is proximate to an anatomical region of interest in the patient.
  • the transceiver 1OA shown in FIGURE 2-1 includes a display 36 and/or is capable of being in signal communication with a device with a display, FIGURES 3-3 and 4) operable to allow viewing of processed results from an ultrasound scan, and/or to allow an operational interaction between the user and the transceiver 1OA.
  • the display 36 of the handheld transceiver 1OA can be configured to display alphanumeric data that indicates a proper and/or an optimal position of the transceiver 1OA relative to the selected anatomical portion.
  • Display 36 can be used to view two- or three-dimensional images of the selected anatomical region.
  • the display 36 can be a liquid crystal display (LCD), a light emitting diode (LED) display, a cathode ray tube (CRT) display, or other suitable display devices operable to present alphanumeric data and/or graphical images to a user.
  • LCD liquid crystal display
  • LED light emitting diode
  • CRT cathode ray tube
  • the transceiver dome 34 of the transceiver 1OA can be positioned against a surface portion of a patient that is proximal to the anatomical portion to be scanned. See for example FIGURE 5 for exemplary positioning for AAA detection.
  • the user actuates the transceiver 1OA by depressing a trigger 37.
  • the transceiver 1OA transmits ultrasound signals into the body, and receives corresponding return echo signals that can be at least partially processed by the transceiver 1OA to generate an ultrasound image of the selected anatomical portion.
  • the transceiver 1 OA transmits ultrasound signals in a range that extends from approximately about two megahertz (MHz) to approximately about ten MHz. Ultrasound energies beyond 10 MHz can be utilized.
  • the transceiver 1OA can be operably coupled to an ultrasound system that can be configured to generate ultrasound energy at a predetermined frequency and/or pulse repetition rate and to transfer the ultrasound energy to the transceiver 1OA.
  • the system also includes a processor that can be configured to process reflected ultrasound energy that is received by the transceiver 1OA to produce an image of the scanned anatomical region.
  • the system generally includes a viewing device, such as a cathode ray tube (CRT), a liquid crystal display (LCD), a plasma display device, or other similar display device, that can be used to view the generated image.
  • a viewing device such as a cathode ray tube (CRT), a liquid crystal display (LCD), a plasma display device, or other similar display device, that can be used to view the generated image.
  • the system can also include one or more peripheral devices that cooperatively assist the processor to control the operation of the transceiver 1OA, such a keyboard, a pointing device, or other similar devices.
  • the transceiver 1OA can be a self-contained device that includes a microprocessor positioned within the housing 35 and software associated with the microprocessor to operably control the transceiver 1OA, and to process the reflected ultrasound energy to generate the ultrasound image.
  • the display 36 can be used to display the generated image and/or to view other information associated with the operation of the transceiver 1OA.
  • the information can include alphanumeric data that indicates a preferred position of the transceiver 1OA prior to performing a series of scans.
  • the transceiver 1OA can be operably coupled to a general-purpose computer (see FIGURE 4), such as a laptop or a desktop computer that includes software that at least partially controls the operation of the transceiver 1OA, and also includes software to process information transferred from the transceiver 1OA, so that an image of the scanned anatomical region can be generated.
  • the transceiver 1OA can also be optionally equipped with electrical contacts to make communication with receiving cradles 50 as illustrated in FIGURE 4 below.
  • transceiver 1OA of FIGURE 2-1 can be used in any of the foregoing embodiments, other transceivers can also be used.
  • the transceiver can lack one or more features of the transceiver 1OA.
  • a suitable transceiver need not be a manually portable device, and/or need not have a top-mounted display, and/or can selectively lack other features or exhibit further differences.
  • FIGURE 2-2 is a graphical representation of a plurality of scan planes that form a three-dimensional (3D) array having a substantially conical shape.
  • An ultrasound scan cone 40 formed by a rotational array of two-dimensional scan planes 42 projects outwardly from the dome 34 of the transceivers 1OA.
  • Other transceiver embodiments of transceiver 1OA can also be configured to develop a scan cone 40 formed by a rotational array of two-dimensional scan planes 42.
  • the pluralities of scan planes 40 can be oriented about an axis 43 extending through the transceivers 1 OA.
  • each of the scan planes 42 can be positioned about the axis 43, preferably, but not necessarily at a predetermined angular position ⁇ .
  • the scan planes 42 can be mutually spaced apart by angles ⁇ ⁇ and ⁇ 9 2 .
  • the scan lines within each of the scan planes 42 can be spaced apart by angles ⁇ j and ⁇ 2 .
  • angles ⁇ ⁇ and 6* 2 are depicted as approximately equal, it is understood that the angles ⁇ ⁇ and ⁇ 2 can have different values.
  • the angles ⁇ ⁇ and ⁇ 2 are shown as approximately equal, the angles ⁇ ⁇ and ⁇ 2 can also have different values.
  • Other scan cone configurations are possible. For example, a wedge-shaped scan cone, or other similar shapes can be generated by the transceiver 1OA.
  • FIGURE 2-3 is a graphical representation of a scan plane 42.
  • the scan plane 42 includes the peripheral scan lines 44 and 46, and an internal scan line 48 having a length r that extends outwardly from the transceiver 1OA.
  • a selected point along the peripheral scan lines 44 and 46 and the internal scan line 48 can be defined with reference to the distance r and angular coordinate values ⁇ and ⁇ .
  • the length r preferably extends to approximately 18 to 20 centimeters (cm), although any length is possible.
  • Particular embodiments include approximately seventy-seven scan lines 48 that extend outwardly from the dome 34, although any number of scan lines is possible.
  • the angular movement of the transducer can be mechanically effected and/or it can be electronically or otherwise generated.
  • the number of lines 48 and the length of the lines can vary, so that the tilt angle ⁇ sweeps through angles approximately between -60° and +60° for a total arc of approximately 120°.
  • the transceiver 1OA can be configured to generate approximately seventy-seven scan lines between the first limiting scan line 44 and a second limiting scan line 46.
  • each of the scan lines has a length of approximately 18 to 20 centimeters (cm).
  • the angular separation between adjacent scan lines 48 (FIGURE 2-2) can be uniform or non-uniform.
  • the angular separation ⁇ ⁇ and ⁇ 2 can be about 1.5°.
  • the angular separation ⁇ 1 and ⁇ % can be a sequence wherein adjacent angles can be ordered to include angles of 1.5°, 6.8°, 15.5°, 7.2°, and so on, where a 1.5° separation is between a first scan line and a second scan line, a 6.8° separation is between the second scan line and a third scan line, a 15.5° separation is between the third scan line and a fourth scan line, a 7.2° separation is between the fourth scan line and a fifth scan line, and so on.
  • the angular separation between adjacent scan lines can also be a combination of uniform and non-uniform angular spacings, for example, a sequence of angles can be ordered to include 1.5°, 1.5°, 1.5°, 7.2°, 14.3°, 20.2°, 8.0°, 8.0°, 8.0°, 4.3°, 7.8°, and so on.
  • FIGURE 2-4 is a graphical representation of a plurality of scan lines emanating from the hand-held ultrasound transceiver 1OA forming a single scan plane 42 extending through a cross-section of a section of the abdominal aorta 51.
  • the number and location of the internal scan lines emanating from the transceiver 1OA within a given scan plane 42 can thus be distributed at different positional coordinates about the axis line 11 as can be required to sufficiently visualize structures or images within the scan plane 42.
  • four portions of off-centered shadow regions 49 are exhibited as irregular regions and a cross section of the abdominal aorta 51 region of interest (ROI) is depicted as in the center.
  • Three portions can be viewable within the scan plane 42 in totality, and one can be truncated by the peripheral scan line 44.
  • FIGURE 3-1 depicts a partial schematic and partial isometric and side view of transceiver 1OB, and a scan cone array comprised of 3D-distributed scan lines in an alternate embodiment of an ultrasound system.
  • a plurality of three-dimensional (3D) distributed scan lines emanating from a transceiver that cooperatively forms a scan cone 45.
  • Each of the scan lines has a length r that projects outwardly from the transceiver 1OB of FIGURES 2-1-2-4.
  • the transceiver 1OB emits 3D-distributed scan lines within the scan cone 30 that can be one-dimensional ultrasound A-lines.
  • Other transceiver embodiments can also be configured to emit 3D-distributed scan lines.
  • these 3D-distributed A-lines define the conical shape of the scan cone 45.
  • the ultrasound scan cone 45 extends outwardly from the dome 34 of the transceiver 1OB centered about an axis line 43.
  • the 3D-distributed scan lines of the scan cone 43 include a plurality of internal and peripheral scan lines that can be distributed within a volume defined by a perimeter of the scan cone 43.
  • the peripheral scan lines 31A- 3 IF define an outer surface of the scan cone 45
  • the internal scan lines 39A-39C can be distributed between the respective peripheral scan lines 31A-31F.
  • Scan line 39B can be generally collinear with the axis 43, and the scan cone 45 can be generally and coaxially centered on the axis line 43.
  • the locations of the internal and peripheral scan lines can be further defined by an angular spacing from the center scan line 39B and between internal and peripheral scan lines.
  • the angular spacing between scan line 39B and peripheral or internal scan lines can be designated by angle ⁇ and angular spacings between internal or peripheral scan lines can be designated by angle 0.
  • the angles ⁇ j, ⁇ 2 , and ⁇ 3 respectively define the angular spacings from scan line 39B to scan lines 39A, 39C, and 3 ID.
  • angles 0 ls 0 2 , and 0 3 respectively define the angular spacings between scan line 31B and 31C, 31C and 39A, and 31D and 31E.
  • the plurality of peripheral scan lines 3 IA-E and the plurality of internal scan lines 39 A-D can be three dimensionally distributed A-lines (scan lines) that are not necessarily confined within a scan plane, but instead can sweep throughout the internal regions and along the periphery of the scan cone 45.
  • a given point within the scan cone 45 can be identified by the coordinates r , ⁇ , and 0 whose values generally vary.
  • the number and location of the internal scan lines emanating from the transceiver 1OB can thus be distributed within the scan cone 45 at different positional coordinates as required to sufficiently visualize structures or images within a region of interest (ROI) in a patient.
  • ROI region of interest
  • the angular movement of the ultrasound transducer within the transceiver 1OB can be mechanically effected, and/or it can be electronically generated.
  • the number of lines and the length of the lines can be uniform or otherwise vary, so that angle ⁇ sweeps through angles approximately between -60° between scan line 39B and 3 IA, and +60° between scan line 39B and 3 IB.
  • angle ⁇ in this example presents a total arc of approximately 120°.
  • the transceiver 1OB can be configured to generate a plurality of 3D-distributed scan lines within the scan cone 45 having a length r, in one embodiment for example, of approximately 20 to 40 centimeters (cm).
  • FIGURE 3-2 illustrates a transceiver 1OC configured with a transducer designed to provide a fan-like (e.g., having planar sides) conic scan cone 47 utilizing C- mode and B-mode ultrasound modalities.
  • the transceiver 1OC projects a series of B-mode scan planes 52A and 52B that oscillate like a pendulum between extremes within the scan cone 47.
  • the B-mode scan planes 52 A, 52B may be derived from a plurality of scan lines similar to scan lines 44, 46, and 48 of FIGURES 2-3 and 2-4.
  • the pendulum oscillating scan planes 52A and 52B can be arranged substantially at right angles to each other as depicted at axis crossing 57.
  • the oscillating scan planes 52A and/or 52B can define a series of C-scan planes 55 that vary in depth location from the transceiver dome 34.
  • the C-scan planes 55 move from the transducer vanishing point, and the B-scan planes angularly radiate from the transducer vanishing point.
  • a portion of the abdominal aorta taken as a C-mode shape is displayed on the transceiver display 16.
  • the C-scan geometry shown as scan planes 55 present a substantially square-like ultrasound area within the scan cone 47.
  • the C-Scan image information contained within scan planes 55 presents a cross-section view of the ultrasound at a particular depth probed by the transceiver 1OC.
  • the C-mode may be more representative of a portion of the abdominal aorta than the actual whole of the length of the aorta.
  • the C- Scan illustrates a cross-section view of the ultrasound at a particular depth to obtain a targeting image of the abdominal aorta.
  • the targeting image may be more of a binary image showing the lines and spaces that are inside the aorta versus those that are outside of the aorta.
  • the definition of C-mode image basically may be a plane parallel to the face of the transducer to obtain a projection image of the AAA region.
  • the C-mode acquired projection image can yield abdominal aorta information not confined to simply one a single plane parallel to the transducer surface, but multiple planes denoted as C-scans.
  • the C-mode acquired projection image may be binary, and can include a non-AAA region and an AAA region.
  • the AAA region may be presented as an interpolated shape that may be generated from one side to the opposite side, for example the left most and the right most sides of a valid segmentation, or cut, the AAA region on all planes.
  • FIGURE 3-3 illustrates a partial isometric and schematic view of an ultrasound AAA detection system 70 utilizing a transceiver probe 1OD and console 74 combination.
  • the AAA detection system 70 may be battery powered and portable and can also be referred to as the BVI9600 with AortaScreen Mode system. Other embodiments can include line power.
  • the ultrasound transceiver 1OD may be configured to send out and receive ultrasound RF signals.
  • the received RF may be transmitted to console 74.
  • the DSP in console can process the RF information to extract the information relevant to the detected feature of each line.
  • Each line is classified as being "available” or “non- available", as will be described below, based on a set threshold value.
  • the classification of information contained in each scanline can be integrated with the image processing module for accurate and optimal positioning of the transceiver for segmentation and volume measurement.
  • the transceiver 1OD can have a transceiver display 36, housing 35 and dome 34 design similar to transceivers 1OA and 1OB, and may be in signal communication to console 74 via signal cable 89.
  • the console 74 can be pivoted from console base 72.
  • the console 74 includes a display 76, detection and operation function panel 78, and select panel 80.
  • the detection and operation function provide for targeting the abdominal aorta, allow user voice annotation recording, retrieval and playback of previously recorded voice annotation files, and current and previously stored 3D and 2D scans.
  • display 76 provides a screenshot 77C with a targeting icon 79A including cross hairs centered in a cross sectional depiction of an abdominal aorta.
  • a targeting icon screenshot 77B with a plurality of directional arrows can appear and flash to guide the user to move the transceiver 1OC to center the abdominal aorta or AAA.
  • the targeting icon screenshot 77B can also appear on the display 36 of the transceiver 1OD.
  • the targeting icon screenshot 77B similarly guides the user to place the transceiver 1OD to center the abdominal aorta or AAA or other region of interest as the directional indicator panel 32 did in transceiver 1OA of FIGURE 2-1 above.
  • the transceiver 1OD can be in wireless communication via wireless signal 82 with the wireless hub 84; the output from the transceiver 1OD can be delivered to a wireless hub 84 via wireless signal port 86.
  • the wireless hub 84 can also charge batteries 88 for loading into the battery compartment (not shown) of console 74. In one embodiment, all the calculations can be performed in the imaging console 74.
  • the BVI9600-embodiment system 70 does not require a computer or network to complete the analysis and imaging processing. In other embodiments, the system 70 can utilize the wireless hub 84 as a gateway to transmit transceiver 1OD acquired imaging information in local and Internet systems similar to that described in FIGURE 4 below.
  • FIGURE 4 is a schematic illustration of a server-accessed local area network and/or an internet system in communication with a plurality of ultrasound imaging systems.
  • an ultrasound system 90 can include one or more personal computer devices 52 that can be coupled to a server 96 by a communications system 95.
  • the devices 52 can be, in turn, coupled to one or more ultrasound transceivers 1OA and/or 1OB and/or 1OC, for examples the ultrasound subsystems 94A-94D.
  • Ultrasound based images of organs or other regions of interest derived from either the signals of echoes from fundamental frequency ultrasound thereof, can be shown within scan cone 40 or 45 presented on display 54.
  • the server 96 can be operable to provide additional processing of ultrasound information, or it can be coupled to still other servers (not shown in FIGURE 4) and devices.
  • Transceivers lOA/B/C can be in wireless communication with computer 52 in sub-system 94A, in wired signal communication in sub-system 94B, in wireless communication with computer 52 via receiving cradle 50 in sub-system 94C, or in wired communication with computer 52 via receiving cradle 50 in sub-system 94D .
  • An Internet system 114 can be coupled or otherwise in communication with the ultrasound sub-systems 94A-94D.
  • a 2D real time ultrasound instrument was determined to be preferably equipped to find optimal scan locations by searching the interested area with realtime B-mode image feedback.
  • the BVI9600 model has been optimized to inform a user whether a 3D data set is valid or not, prior to collection of useful data. This was not a capability of the previous BVI6500 AAA mode, where only one B-mode image (cross section of an aorta) was provided for a user review after the 3D scan was completed.
  • the BVI6500 B-mode image did not represent the 3D data set condition.
  • the BVI9600 of the present invention investigates each scanline to determine whether the scanline contains any object information relevant to the region of interest or is just shadow blocked by air pocket.
  • the scanline may be called "Available”.
  • FIGURE 5 a plurality of measurements along a patient's midline may be taken to find the maximal measurement of the abdominal aorta.
  • gel may be applied at the midline of the body approximately 1 inch (2.5 cm) below the xiphoid (sternum); making sure that there is enough gel to cover the probe head 34.
  • FIGURE 5 shows four such exemplary scan probe locations along the midline of a patient's torso.
  • Finn pressure may be applied, noting that patients with excess tissue can require more pressure in order to obtain a good quality scan.
  • the trigger/scan button 37 of the transceiver 1OA may be then depressed.
  • An aim/scan screen appears, for example in display 36 of the transceiver (FIGURE 2.1) and/or display 76 on consol 74 and the instrument performs a 3D scan of the patient's region of interest (ROI).
  • the aiming result and B-mode image may be displayed, for example as shown in an exemplary screen shot FIGURE 10.
  • AAA ultrasound scanning a subject is required to fast overnight to have less gas in their intestines so that the images are more shadow-free due to air block. Air blocks not only degrade the B-mode image quality, but also mislead the delineation of abdominal aorta which is an advantageous factor for volume estimation accuracy. Specifically, the acoustic energy can be blocked by a layer which has high acoustic impedance difference as between two media. Examples can include a fat layer under skin or air pockets in the intestines. The presence of fat and/or air pockets prevent information relevant to the region of interest (ROI), such as an AAA wall or lumen location, to be reflected to the transducer resulting in dark or black region underneath (See FIGURE 6 as an example).
  • ROI region of interest
  • A-mode scanlines
  • FIGURE 6-1 a B-mode image of an AAA may be in the same field as a shadow region due to the air-block or fat layer, two A-mode lines 60, 65.
  • A-mode line 60 there may be an AAA between depth 85 and 120 mm, therefore A-mode line 60 has useful information such as tissue, AAA, and acoustic enhancement.
  • A-mode line 65 however there may be no information due to the acoustic block near the depth of 25mm.
  • FIGURES 6-2 and 6-3 are representative histograms of the A-mode lines 60 (FIGURE 6-2) and 65 (FIGURE 6-3) wherein, a threshold line set at 40 (67), in the exemplary histogram, is shown. The threshold line 67 may be used to determine whether the A-mode has available information or not.
  • FIGURES 7-1-7-12 and 8 illustrate how the aorta aiming guide can determine "availability" and may be used to help the user correct the position of the transceiver to scan AAA.
  • a threshold may be established, which can vary from patient to patient, but in this example the threshold is 40 (as above), and may be related to the intensity number, for example, out of 256 levels after A-law compression.
  • a searching range in this non-limiting example is between 42% of A-mode length to 100% of A- mode length, the length of which may be related to experimental testing of optimal scan depths. Given these parameters, once set, if any A mode point within the searching range is above the established threshold but below an acoustic block, the A mode line may be determined as "Available".
  • non-available scanlines can be differentiated from “available” scanlines, for example, “non-available” are shadowed.
  • A-mode lines are marked non-available by colored markings/dots underneath the B -mode image indicating the non-availability of that area.
  • colored markings/dots underneath the B-mode image represent approximately 40% unavailable scanlines, compared to >90% unavailable scanlines in FIGURE 7-6.
  • the available and non-available regions from all 12 planes of this example of the B-mode display may be then plotted in C-mode representing the pseudo- horizontal cross section of the 3D scan cone composed by the 12 sectors as shown in FIGURE 8.
  • FIGURE 8 shows the C-mode display of an "Availability Plot" 80 for
  • the "Availability Plot" 80 in C-mode depicts available regions 85 as solid regions and non-available regions 87 as white and/or hashed lines.
  • available regions 85 As solid regions
  • non-available regions 87 As white and/or hashed lines.
  • ⁇ 46% of the aggregate scanned regions are defined as "available.”
  • the percentage of available A- mode lines in the cone data, defined as "Availability" can be calculated using equation
  • the above calculation may be used to automatically determine whether the 3D full scan is appropriate for segmentation and display and without shadow within the region of interest.
  • the C-mode display 81 can be partitioned into concentric circles, in this example, an inner circle 93 and an outer circle 97, to directionally position the transceiver to an available area, .
  • Each quarter of the circle can be defined as a fan 98, from a perpendicular point defined by the cross-hairs in the middle; each overlapping fan expanding, for example, by 45° to both side in either direction .
  • Each direction can have a fan, and the fans may be overlapped for the greatest coverage.
  • FIGURE 9-2 illustrates an embodiment of a decision tree utilizing the above schema.
  • the availability may be calculated, starting within the inner circle 93. If the inner circle 93 has a greater than or equal to 75% available scanlines, the scan is a "good scan.” Otherwise, a direction may be suggested to re-position the scanhead. Repositioning may be suggested based on a criterion, generally in the direction corresponding to any inner fan and/or inner + outer fan availability determination that has the greatest calculated "availability.” In other words, the arrow guide indicates the direction of greatest availability.
  • the suggested aiming direction may be the one whose inner fan availability is maximal.
  • An arrow indicator guide may be a moderate and solid light indicator 121, 122, indicating that little movement is suggested.
  • the suggested aiming direction may be based on the direction in which both the inner and outer fan has a maximum availability.
  • the arrow indicator guide may be a flashing aggressive light indicator 124, 126, and the user may be guided in the recommended direction and orientation to achieve an optimized availability.
  • FIGURE 10 is an exemplary screen view 100 of the aiming guide of the BVI9600.
  • a B-mode image 102 and a C-mode image 103 is shown. Due to the air block in the 3D data, shown in the C-mode display 103 as non-available (by hatched lines) 87, the user is directed to move in a NW direction, i.e. towards the upper left had quadrant 110 where the greatest availability and/or "available" scanlines may be detected to collect a better data set.
  • NW direction i.e. towards the upper left had quadrant 110 where the greatest availability and/or "available" scanlines may be detected to collect a better data set.
  • the user moves the instrument around on the abdomen to measure the abdominal aorta, which may be calculated from full three-dimensional scan cone or, can be a single two dimensional plane.
  • the use of the three-dimensional scan cone removes orientation requirements, permitting the user to position the device in any orientation.
  • the user takes several three-dimensional image scans, moving along the patient's abdomen. After each scan, the volume/diameter of the section of the aorta covered by that scan may be displayed and the image may be stored if the diameter from the new scan is larger than any previous diameter.
  • the image produced in this embodiment, whether it be form two dimensional planes or three dimension scan cone, can be optimally transmitted via the internet for remote enhance processing and rendering as is illustrated in FIGURE 4.
  • an AAA detection algorithm 120 may be used to [0088]
  • the fundamental AAA detection algorithm 120 used in the BVI 9600 device begins with process block Find Initial Wall 122 using A-mode ultrasound data that incorporates data-smoothing. Find Initial Wall 122 looks for the front and back walls of the abdominal aorta illustrated and described in FIGURE 12 below. After the front and back walls are found, a line passing through the center of the abdominal aorta may be determined in the following process block Find Centroid 124.
  • Fix Initial Walls 130 refines the initial wall points, removes any outliers, and fills gaps in the detected wall location regions.
  • the Centroid 118 as shown in FIGURE 12-3, may be found by averaging the front and back walls.
  • the maximum diameter (2A) 129 of a sphere that can fit within the AAA may be calculated as the maximum distance between the front and back walls.
  • the Fix Walls algorithm 130 may be adjusted to comply with this spherical assumption to limit the segmentation of the AAA. That is, the given segmentation (from automatic detection) may be modified not to exceed the circle boundary whose diameter is the maximum of one, maximized slice.
  • FIGURE 12-1 and 12-2 B-mode images of a longitudinal section 12-1 and cross section 12-2 are shown with their exemplary segmentations for volume estimation.
  • the AAA phantom has 4 cm, and 5 cm diameter in axial and lateral directions, respectively.
  • FIGURE 12-3 is a schematic representation of the application of the limited segmentation wherein a maximum diameter may be fixed from the centroid 118.
  • Two lengths, length (2A) 129 the diameter, and length B 135, the radius, represent the maximum circle that can fit inside the AAA.
  • the mathematical steps to calculate a modified segmentation parameter are described with reference to FIGURE 11, above.
  • the volume calculation at BLOCK 134 may be the result of the aggregate calculation of volume resulting from the limited segmentation, in this example, represented by lines 136 and 138, i.e. those within the scanlines 145, 147 at points corresponding to the intersection of the front wall and the back wall of the artery being measured.
  • the diameter calculation may be derived by application mathematical relationship between volume and diameter as further discussed in relationship to FIGURE 13 below.
  • FIGURE 13 graphically represents the relationship between volume and diameter as determined using the below formulae.
  • the volume of sphere volume is:
  • V is a volume and D is a diameter
  • the mathematical functions can be calculated automatically from a lookup table of diameter values such that the diameter calculations need not be performed for each volume under consideration.
  • the above relationship can be expressed as a look-up table, where once the volume (at BLOCK 134) is determined by the limited segmentation procedure, the diameter may be automatically correlated (at a BLOCK 137) with a pre-calculated diameter according to the above relationship.
  • FIGURE 14 depicts a flowchart of the Find Initial Walls sub-algorithm of FIGURE 11.
  • Find Initial Walls 122 process may be executed on every A-mode scan line and may be subjected to averaging and low-pass filtering using, for example, a 15 or 16 sample set beginning with process block 140.
  • a local gradient at process block 142 may be computed for each sample point using a central difference formulation taken for seven samples.
  • the central difference formulation is defined by equations 4-9 (Eq. A-
  • Equation 4 The standard central difference formula is given in Equation 4:
  • the normalization factor is simply the distance between the two points. In Eq. 4 the distance separating the two means in the calculation was 1 , and in Eq. 5 the step between the two means is 2.
  • the normalization of the gradient by the step size while mathematically correct, incurs a cost in terms of operation. This operation may be neglected in the gradient calculation for the aortic wall detection algorithm with minimal effect: since the same calculation is performed for every data sample, every data sample can have the same error and thus the relative gradient values between different samples remain unchanged.
  • the gradient calculation may be expanded to three neighboring points to each side of the sample in question. This calculation is shown in Eq. 6. This calculation is simply the sum of three gradient approximations and thus the end result can be 12 times its normal value. This deviation from the true mathematical value has minimal effect since the calculation may be the same at each point and thus all the gradient values can be 12 times their usual value.
  • the full calculation is written in Eq. 7.
  • the first line shows the summation calculation to obtain the means, and the difference operations to obtain the gradient.
  • the entire sum is normalized by 15, the number of points included in each local mean.
  • the second line of the operation shows the result when the summations are simplified, and represents the maximal implementation of the calculation.
  • This calculation incurs a cost of 23 additions or subtractions, 2 floating-point multiplications, 1 floating point division, and at least 1 temporary variable. This operation cost is incurred for 91% of the data samples.
  • the cost of the calculation can be reduced by not simplifying the summations and using a running sum operation. In that manner, only one mean may need to be calculated for each point, but that mean needs to be for the i+3 point.
  • the running sum calculation uses the previous sum, and then corrects the sum by subtracting the old "left hand" end point and adding the new "right hand” end point. The operation is shown in Eq. 8. This running sum operation incurs a cost of 5 additions and subtractions.
  • Equation 9 Since the running sum was calculated for the i+3 point, all average values are available for the gradient calculation. This calculation is shown in Equation 9:
  • This equation has the same form as the one in Eq. 6 except for the normalization factor of 16.
  • This normalization factor is not a result of the gradient operation, but rather it is the normalization factor mean calculation.
  • the factor of 16 is used instead of the standard value of 15 that one would expect in a 15 -point average for this simple reason: using a factor of 16 eliminates floating-point division. If the means are normalized by 16, then the division operation can be replaced by a "righf'-shift by 4 at a significantly lower cost to the embedded system. Therefore the gradient operation has eleven additions and subtractions and one shift by 4.
  • the results from local gradient 142 are subjected to loop limit processing between blocks 144 and 158 to obtain the best front wall and back wall pair for each scan line denoted as a tissue gradient or tissue delta.
  • the best front wall and back wall pair on each line may be defined as the front wall and back wall pair for which the pixel intensity difference in the back wall gradient and front wall gradient is the maximum and the smallest local average between front wall and back wall pair is the minimum.
  • the loop limit processing begins with loop limit block 144 that receives pixel values for each sample in the detection region and subjects the pixel intensity values to determine whether the gradient is minimum at decision diamond 146. If affirmative, then the pixel values may be ascertained whether it's the best front wall-back wall (FW/BW) candidate combination at decision diamond 147. If affirmative, the FW/BW candidate pair may be saved and loop limit processing returns to limit block 144. If negative, at process block 152, the Front Wall pixel value may be saved and another back wall candidate may be sought with a subsequent return to loop limit block 152.
  • FW/BW front wall-back wall
  • sub-algorithm 122 continues to decision diamond 156 to determine whether the back wall and the gradient is maximum. If affirmative, at process block 154, a candidate BW/FW pair is established and sub-algorithm re-routes to loop limit block 144. If negative, the end of analysis for a particular FW/BW candidate occurs at loop limit block 158 either routes back to the limit loop block 144 or exits to find Centroid 124.
  • Formulations relating to Find Centroid 124 may be based on coordinate geometries described in equations 10 and 11 utilizing coordinate conversions.
  • the trigonometric functions can be calculated for a table of ⁇ values such that the cosine and sine calculations need not be performed for each of the points under consideration.
  • the closest plane can be found by finding the shortest vector from each plane to the centroid.
  • the shortest vector from a plane to a point can be the perpendicular to the projection of the centroid on the plane.
  • the projection of the centroid on the plane is defined as the dot product of the centroid vector, c, with the plane definition vector, a, divided by the length of the plane definition vector. If the plane definition vector is a unit vector, then the division operation is unnecessary. To find the perpendicular to the projection, it is sufficient to subtract the projection vector from the centroid vector as shown in Eq. 11 :
  • the length of these projections can be found by calculating the Euclidean norm for each line.
  • the Euclidean norm may be more commonly known as the length or magnitude of the vector.
  • To find the plane closest to the centroid calculate the lengths for the perpendicular to the projection of the centroid on each plane, and take the plane with the shortest of these lengths.
  • the preferred embodiment of the invention has been illustrated and described, as noted above, many changes can be made without departing from the spirit and scope of the invention. For example, it should be understood that various changes, modifications, and substitutions can be incorporated in the apparatus embodiment to achieve the ultrasonic, volumetric determination of the abdominal aorta to thus apply that to the aortic diameter reading. Accordingly, the scope of the invention is not limited by the disclosure of the preferred embodiment.

Abstract

Systems, methods, and ultrasound transceivers equipped and configured to execute analysis and extract ultrasound information related to an abdominal aortic aneurysm of a subject are described. The methods utilize algorithms to establish improved targeting of the abdominal aortic aneurysm within a region-of-interest. The targeting algorithms may be optimally applied to provide the user with real-time feedback and orientation guidance for positioning the transceiver. Additional methods utilize diameter conversion algorithms to establish the diameter of the abdominal aortic aneurysm based on conversion of the volume measurement and limited segmentation within a targeted region-of-interest of the aorta.

Description

DEVICE, SYSTEM, AND METHOD TO MEASURE ABDOMINAL AORTIC
ANEURYSM DIAMETER
COPYRIGHT NOTICE
THIS DISCLOSURE IS PROTECTED UNDER UNITED STATES AND INTERNATIONAL COPYRIGHT LAWS. © 2009 VERATHON® INCORPORATED. ALL RIGHTS RESERVED. A
PORTION OF THE DISCLOSURE OF THIS PATENT DOCUMENT CONTAINS MATERIAL WHICH IS SUBJECT TO COPYRIGHT PROTECTION. THE COPYRIGHT OWNER HAS NO OBJECTION TO THE FACSIMILE REPRODUCTION BY ANYONE OF THE PATENT DOCUMENT OR THE PATENT DISCLOSURE, AS IT APPEARS IN THE PATENT AND TRADEMARK OFFICE PATENT FILE ORRECORDS, BUT OTHERWISE RESERVES ALL COPYRIGHT RIGHTS WHATSOEVER.
CROSS REFERENCE TO RELATED APPLICATIONS
[0001] This application claims priority to, and incorporates by reference in their entireties, U.S. provisional patent application serial number 61/087,152 filed August 7, 2008 and U.S. provisional patent application serial number 61/094,003 filed September 3, 2008.
[0002] This application is also a continuation-in-part of, claims priority to, and incorporates by reference in its entirety U.S. patent application serial number 12/121,721 filed May 15, 2008, which is a continuation-in-part of U.S. patent application serial number 11/968,027 filed December 31, 2007, U.S. patent application serial number 11/926,522 filed October 27, 2007, U.S. patent application serial number 11/925,887 filed October 27, 2007, U.S. patent application serial number 11/925,896 filed October 27, 2007, U.S. patent application serial number 11/925,900 filed October 27, 2007, U.S. patent application serial number 11/925,850 filed October 27, 2007, U.S. patent application serial number 11/925,843 filed October 27, 2007, U.S. patent application serial number 11/925,654 filed October 26, 2007, and U.S. Provisional Patent Application Nos. 60/938,359 filed May 16, 2007; 60/938,371 filed May 16, 2007; and 60/938,446 filed May 16, 2007.
[0003] All of the foregoing named applications are incorporated by reference in their entireties as if fully set forth herein.
BACKGROUND OF THE INVENTION
[0004] The aorta artery in the abdomen carries blood from the heart to the abdominal region. One disorder of the abdominal aorta is known as an abdominal aortic aneurysm, which is a permanent localized dilation of the arterial wall of the abdominal aorta. When dilation of the arterial wall is greater than 1.5 times the typical, i.e. nominal, diameter, it is referred to as an aneurysm. A normal abdominal aorta is shown in FIGURE 1-1 (see U.S. Patent No. 6,905,468). FIGURE 1-2 shows a typical aortic aneurysm at. An aortic aneurysm is usually located below the renal arteries and the kidney arteries and above the aorta-iliac bifurcation. Below the aortic-iliac bifurcation are additional arteries. Abdominal aortic aneurysms are a fairly common disorder, occurring in approximately 5- 7% of the population over age 60. Since effective screening programs are not yet established, the diagnosis of AAA is frequently made at the time of rupture or impending rupture, which leads to a dramatic increase of post-operative morbidity. (Daly, et al., 2004)
[0005] Abdominal aortic aneurysms, depending upon their size, result in pressure on adjacent tissue structure and organs, causing potential embolization and/or thrombosis in those tissues/organs. Rupture of the aneurysm typically results in death, and comprises approximately 2% of all deaths in men over 60 years of age. [0006] Accurate diagnosis of an abdominal aortic aneurysm is important in preventing rupture, as well as in controlling the expansion of the aneurysm. Conventional two-dimension B-mode ultrasound scan devices are currently used to produce measurements of aortic aneurysms, both axially (longitudinally) along the aorta and laterally (radially) across the aorta. Typically, the accuracy is within three millimeters of the actual size of the aneurysm, using conventional CT or MRI processing. These conventional systems, however, are very expensive, both to purchase/lease and to maintain. Further, a trained sonographer is necessary to interpret the results of the scans. As a result, many aneurysms go undetected and/or are not consistently monitored after discovery, until rupture resulting in death of the patient.
[0007] A recent prospective study by Vidakovic, et al. (2006) sought to evaluate the diagnostic potential and accuracy in Abdominal Aortic Aneurysm (AAA) screening using an automatic bladder volume indicator (BVI) instrument. The BVI was originally designed for the estimation of post- void residual volumes. The device is inexpensive and can effectively be used after a short training. A measurement method of bladder volume is different between BVI and US, however several reports have found that BVI is as reliable as US to measure post- void residual urine. (Yucel, et ah, 2005; Byun, et al., 2003)
[0008] In the Vidakovic et al. study AAA volumes were measured in 94 patients, and compared with 2D ultrasound and CT measurement to see if these comparisons can provide a method of screening AAA within certain volume thresholds. The reported results indicated there was an 89% agreement of the diameter measurements by ultrasound (US) as compared with those made with the bladder volume indicator (BVI). Using a cut-off value for the presence of AAA of 50 ml by BVI, the BVI technique predicted AAA with a sensitivity of 94%, a specificity of 82%, a positive predictive value of 88%, and a negative predictive value of 92%. The agreement between standard US and BVI in detecting an AAA was 89%.
[0009] This study showed the potential of using the BVI volume. Compared to other portable US devices used to screen patients, the BVI is simpler for use, requires a shorter training period, and is significantly cheaper. One barrier to its adoption is that the current device does not provide automatic conversion values and/or accurate values of AAA diameter. Moreover certain impediments exist to accurate readings of the region of interest that must be overcome for accurate predictive measurements.
[0010] Hence, it would be useful to a primary care physician or emergency personnel to have a low-cost device which provides accurate information concerning aortic aneurysms by providing AAA diameter measurements, without the necessity of a trained technician to interpret the scan results. Specifically, the art fails to provide a low cost system, method, and apparatus to automatically and accurately obtain and utilize data derived from an automatic bladder volume instrument (BVI) to provide values of abdominal aortic aneurysm (AAA) diameters. SUMMARY OF THE INVENTION
[0011] An ultrasound transceiver device, system, and method to obtain, analyze, and interpret ultrasonic information from a vascular region of interest to measure the diameter of a suspected blood vessel aneurysm.
BRIEF DESCRIPTION OF THE DRAWINGS
[0012] Examples of particular embodiments of the present invention are described in detail below with reference to the following drawings: [0013] FIGURE 1-1 shows a normal abdominal aorta;
[0014] FIGURE 1-2 shows a typical abdominal aortic aneurysm;
[0015] FIGURES 2-1 thru 2-4 depict a partial schematic and a partial isometric view of a transceiver, a scan cone comprising a rotational array of scan planes, and a scan plane of the array of an ultrasound harmonic imaging system;
[0016] FIGURE 3-1 is a side elevation view of an ultrasound transceiver 1OA that includes an inertial reference unit
[0017] FIGURE 3-2 illustrates a side and partial isometric that schematically depicts an ultrasound scanner employing C-mode and B-mode ultrasound modalities;
[0018] FIGURE 3-3 illustrates a partial isometric and schematic view of an ultrasound scanner system;
[0019] FIGURE 4 is a schematic illustration of a server-accessed local area network and/or internet in communication with a plurality of ultrasound imaging systems;
[0020] FIGURE 5 shows exemplary probe locations along the midline of a patient's torso.
[0021] FIGURE 6-1 is a B-mode image of an AAA is in the same field as a shadow region ; [0022] FIGURES 6-2 and 6-3 are histograms of the A-mode lines 65 (FIGURE 6-2) and 65 (FIGURE 6-3);
[0023] FIGURES 7-1 thru 7-12 are twelve B-mode images showing shadows and non-availability;
[0024] FIGURE 8 schematically illustrates an availability plot in C-mode (view from the top of the scan cone);
[0025] FIGURE 9-1 and 9-2 illustrates an aiming direction scheme; and a decision tree for the arrow feedback from the indicator 22 of FIGURE 3-3;
[0026] FIGURE 10 illustrates a screen shot of an exemplary aiming guide;
[0027] FIGURE 11 depicts a flowchart of a aortic aneurysm detection algorithm employing fundamental ultrasound energies;
[0028] FIGURES 12-1 and 12-2 are exemplary B-mode images of longitudinal section and cross section of AAA phantom with their typical segmentations (red dotted lines) for volume estimation;
[0029] FIGURES 12-3 is a schematic diagram illustrating an exemplary embodiment of the limited segmentation;
[0030] FIGURE 13 is a graphic diagram for the conversion from volume to diameter assuming an object shape is sphere; and
[0031] FIGURE 14 depicts a flowchart of the Find Initial Walls sub-algorithm of FIGURE 11.
DETAILED DESCRIPTION OF THE PARTICULAR EMBODIMENTS
[0032] Embodiments of the invention may be described in the general context of computer-executable instructions, such as program modules, executed by one or more computer processors or other devices, computer-readable media on which such
instructions are stored, and/or the processors/devices themselves. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types. Typically the functionality of the program modules may be combined or distributed as desired in various embodiments.
[0033] Embodiments of the invention may include or otherwise utilize at least some form of computer readable media, which may be associated with one or more processors and/or memory devices. Computer readable media can be any available media that can be accessed by one or more components of such operating environment. By way of example, and not limitation, computer readable media may comprise computer storage media and communication media. Computer storage media includes volatile and nonvolatile, removable and non-removable media implemented in any method or technology for storage of information such as computer readable instructions, data structures, program modules or other data. Computer storage media includes, but is not limited to, RAM, ROM, EEPROM, flash memory or other memory technology, CD- ROM, digital versatile disks (DVD) or other optical storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other medium which can be used to store the desired information and which can be accessed by one or more components of such operating environment. Communication media typically embodies computer readable instructions, data structures, program modules or other data in a modulated data signal such as a carrier wave or other transport mechanism and includes any information delivery media. The term "modulated data signal" means a signal that has one or more of its characteristics set or changed in such a manner as to encode information in the signal. By way of example, and not limitation, communication media includes wired media such as a wired network or direct-wired connection, and wireless media such as acoustic, RF, infrared and other wireless media. Combinations of the any of the above should also be included within the scope of computer readable media. [0034] Particular embodiments are described for devices, systems, and corresponding methods encompassing ultrasound detection and the measurement of suspected abdominal aortic aneurysms. The devices, systems, and methods employ transceivers equipped to convey fundamental ultrasound frequencies, and analysis of fundamental echoes returning from a vascular region of interest (ROI). Signal processing algorithms executable by computer systems, described below, are developed to optimally extract information from fundamental ultrasound echoes delivered under A-mode, B- mode, and/or C-mode ultrasound configurations.
[0035] Disclosure below includes systems and methods to detect and measure an AAA involving transmitting ultrasound energy having at least one of a fundamental frequency to the AAA, collecting ultrasound echoes returning from the AAA and generating signals from the ultrasound echoes, and identifying within the ultrasound signals those attributable to fundamental ultrasound frequencies. Thereafter, the fundamental-frequency-derived signals undergo signal processing via computer executable program instructions to present an image of the AAA on a display, and calculate the volume and diameter of the AAA.
[0036] The signal processing applied to the transceiver echoic fundamental ultrasound signals include an algorithm having computer readable instructions for ascertaining the certainty that a given scan line traverses a AAA, a shadow region, or both a AAA and a shadow region using a grading algorithm for predicting the scan line's AAA or shadow classification.
[0037] The ultrasound transceivers or distal collection devices (DCD) are capable of collecting in vivo three-dimensional (3-D) cone-shaped ultrasound images of a patient. During the data collection process initiated by the DCD, a radio frequency pulsed ultrasound field is transmitted into the body, and the back-scattered "echoes" may be transducer-detected and presented as a one-dimensional (1-D) voltage trace, which may be also referred to as a radio frequency (RF) line. After detection of RF signal pulse envelopes, a set of 1-D data samples is interpolated to form a two-dimensional (2 -D) image. A plurality of 2-D images can be assembled to form a 3-D ultrasound image.
[0038] Particular embodiments described below include a system to detect an abdominal aortic aneurysm. The system includes an ultrasound transceiver positioned to deliver ultrasound energy and receive echoes of the ultrasound energy across a plurality of scan planes; an algorithm configured to signal process the received echoes and characterize detected signals across a plurality of scan planes, wherein a "percentage of availability" measurement may be made. Percentage of availability (POA) may be based on whether relevant information is contained within the scanplane as compared with shadows or other regions not of interest. The system includes a display for presenting a visual depiction of availability contained in a plurality of scan planes; and display for indicating positioning information of the ultrasound transceiver based on the visual depiction of availability.
[0039] The B VI9600 Transceiver and Principal of Operation:
[0040] FIGURES 2-1- 2-4 depict a partial schematic and a partial isometric view of a transceiver, a scan cone comprising a rotational array of scan planes, and a scan plane of the array of various ultrasound systems capable of collecting RF line analysis.
[0041] FIGURE 2-1 is a side elevation view of an ultrasound transceiver 1OA that includes an inertial reference unit 38, according to an embodiment of the invention. The transceiver 1OA includes a transceiver housing 35 having an outwardly extending handle 33 suitably configured to allow a user to manipulate the transceiver 1OA relative to a patient. Ultrasound transducers operating within the transceiver 1OA can be equipped to collect and ready signals for ultrasound fundamental and/or harmonic frequency analysis.
[0042] The handle 33 includes a trigger 37 that allows the user to initiate an ultrasound scan of a selected anatomical portion. The transceiver 1OA also includes a transceiver dome 34 that contacts a surface portion of the patient when the selected anatomical portion is scanned (See FIGURE 5). The dome 34 generally provides an appropriate acoustical impedance match to the anatomical portion and/or permits ultrasound energy to be properly focused as it is projected into the anatomical portion. The transceiver 1OA further includes one transducer element, or preferably an array of separately excitable ultrasound transducer elements, positioned within or otherwise adjacent to the housing 35. The transducer elements can be suitably positioned within the housing 35 or otherwise to project ultrasound energy outwardly from the dome 34, and to permit reception of acoustic reflections generated by internal structures within the anatomical portion. The one or more array of ultrasound elements can include a one- dimensional, or a two-dimensional array of piezoelectric elements that can be moved within the housing 35 by a motor. Alternatively, the array can be stationary with respect to the housing 35 so that the selected anatomical region can be scanned by selectively energizing the elements in the array.
[0043] In one embodiment of the transceiver 1OA, a directional indicator panel or aiming guide panel 32 includes a plurality of arrows that can be illuminated for initial targeting and guiding a user to access the targeting of an organ or structure within a region of interest (ROI).
[0044] In the BVI 9600 system 70 described in FIGURE 3-3 below, the directional indicator panel 32 of transceiver 1OA has a virtual equivalent in the form of a targeting icon screenshot 77B that appears on an display 76 on the console 74; both the indicator panel 32 and display 76 (displaying targeting icon 77B) can function to guide a transceiver user to place the transceiver to obtain a "good scan" of the abdominal aortic region of interest (e.g., for the detection of AAA).
[0045] In particular embodiments, if the AAA structure is centered (as indicated by reference numeral 77C), and there are no impediments to the scanlines, e.g., shadows caused by air pockets (see discussion below) from placement of the transceiver 1 OA or 1OC acoustically placed against the dermal surface at a first location of the subject, the directional arrows will be not illuminated. If the AAA is off-center or a shadow appears in the field of interest, an arrow or set of arrows can be illuminated to direct the user to reposition the transceiver 1OA, 1OC acoustically at a second or subsequent dermal location of the subject. The acoustic coupling can be achieved by liquid sonic gel applied to the skin of the patient or by sonic gel pads against which the transceiver dome 34 is placed. The directional indicator panel 32 can also be presented on the display 54 of computer 52 in imaging subsystems described in FIGURES 4 below, and/or presented on the transceiver display 36.
[0046] Turning back to FIGURE 2.1, transceiver 1OA can, in one embodiment, include an inertial reference unit that includes an accelerometer and/or gyroscope 38 positioned preferably within or adjacent to housing 35. The accelerometer 38 can be operable to sense an acceleration of the transceiver 1OA, preferably relative to a coordinate system, while the gyroscope can be operable to sense an angular velocity of the transceiver 1OA relative to the same or another coordinate system. Accordingly, the gyroscope can be of conventional configuration that employs dynamic elements, or it can be an optoelectronic device, such as the known optical ring gyroscope. In one embodiment, the accelerometer and the gyroscope can include a commonly packaged and/or solid-state device. One suitable commonly packaged device can be the MT6 miniature inertial measurement unit, available from Omni Instruments, Incorporated, although other suitable alternatives exist. In other embodiments, the accelerometer and/or the gyroscope can include commonly packaged micro-electromechanical system (MEMS) devices, which are commercially available from MEMSense, Incorporated. As described in greater detail below, the accelerometer and the gyroscope cooperatively permit the determination of positional and/or angular changes relative to a known position that is proximate to an anatomical region of interest in the patient. Other configurations related to the accelerometer and gyroscope concerning transceivers 1OA, 1OB, 1OC of FIGURES 2.-1, 3-2, and/or 3-3 equipped with inertial reference units and the operations thereof are described in co-pending U.S. Patent Application Serial No. 11/222,360 filed September 8, 2005, herein incorporated by reference.
[0047] The transceiver 1OA shown in FIGURE 2-1 includes a display 36 and/or is capable of being in signal communication with a device with a display, FIGURES 3-3 and 4) operable to allow viewing of processed results from an ultrasound scan, and/or to allow an operational interaction between the user and the transceiver 1OA. For example, the display 36 of the handheld transceiver 1OA can be configured to display alphanumeric data that indicates a proper and/or an optimal position of the transceiver 1OA relative to the selected anatomical portion. Display 36 can be used to view two- or three-dimensional images of the selected anatomical region. Accordingly, the display 36 can be a liquid crystal display (LCD), a light emitting diode (LED) display, a cathode ray tube (CRT) display, or other suitable display devices operable to present alphanumeric data and/or graphical images to a user.
[0048] In operation, to scan a selected anatomical portion of a patient, the transceiver dome 34 of the transceiver 1OA can be positioned against a surface portion of a patient that is proximal to the anatomical portion to be scanned. See for example FIGURE 5 for exemplary positioning for AAA detection. The user actuates the transceiver 1OA by depressing a trigger 37. In response, the transceiver 1OA transmits ultrasound signals into the body, and receives corresponding return echo signals that can be at least partially processed by the transceiver 1OA to generate an ultrasound image of the selected anatomical portion. In a particular embodiment, the transceiver 1 OA transmits ultrasound signals in a range that extends from approximately about two megahertz (MHz) to approximately about ten MHz. Ultrasound energies beyond 10 MHz can be utilized.
[0049] In one embodiment, the transceiver 1OA can be operably coupled to an ultrasound system that can be configured to generate ultrasound energy at a predetermined frequency and/or pulse repetition rate and to transfer the ultrasound energy to the transceiver 1OA. The system also includes a processor that can be configured to process reflected ultrasound energy that is received by the transceiver 1OA to produce an image of the scanned anatomical region. As discussed, the system generally includes a viewing device, such as a cathode ray tube (CRT), a liquid crystal display (LCD), a plasma display device, or other similar display device, that can be used to view the generated image. The system can also include one or more peripheral devices that cooperatively assist the processor to control the operation of the transceiver 1OA, such a keyboard, a pointing device, or other similar devices. In still another particular embodiment, the transceiver 1OA can be a self-contained device that includes a microprocessor positioned within the housing 35 and software associated with the microprocessor to operably control the transceiver 1OA, and to process the reflected ultrasound energy to generate the ultrasound image. Accordingly, the display 36 can be used to display the generated image and/or to view other information associated with the operation of the transceiver 1OA. For example, the information can include alphanumeric data that indicates a preferred position of the transceiver 1OA prior to performing a series of scans.
[0050] In yet another particular embodiment, the transceiver 1OA can be operably coupled to a general-purpose computer (see FIGURE 4), such as a laptop or a desktop computer that includes software that at least partially controls the operation of the transceiver 1OA, and also includes software to process information transferred from the transceiver 1OA, so that an image of the scanned anatomical region can be generated. The transceiver 1OA can also be optionally equipped with electrical contacts to make communication with receiving cradles 50 as illustrated in FIGURE 4 below. Although transceiver 1OA of FIGURE 2-1 can be used in any of the foregoing embodiments, other transceivers can also be used. For example, the transceiver can lack one or more features of the transceiver 1OA. For example, a suitable transceiver need not be a manually portable device, and/or need not have a top-mounted display, and/or can selectively lack other features or exhibit further differences.
[0051] FIGURE 2-2 is a graphical representation of a plurality of scan planes that form a three-dimensional (3D) array having a substantially conical shape. An ultrasound scan cone 40 formed by a rotational array of two-dimensional scan planes 42 projects outwardly from the dome 34 of the transceivers 1OA. Other transceiver embodiments of transceiver 1OA can also be configured to develop a scan cone 40 formed by a rotational array of two-dimensional scan planes 42. The pluralities of scan planes 40 can be oriented about an axis 43 extending through the transceivers 1 OA. One or more, or preferably each of the scan planes 42 can be positioned about the axis 43, preferably, but not necessarily at a predetermined angular position θ . The scan planes 42 can be mutually spaced apart by angles θ \ and <92. Correspondingly, the scan lines within each of the scan planes 42 can be spaced apart by angles φ j and φ 2. Although the angles θ \ and 6*2 are depicted as approximately equal, it is understood that the angles θ \ and θ 2 can have different values. Similarly, although the angles φ \ and φ 2 are shown as approximately equal, the angles φ \ and φ 2 can also have different values. Other scan cone configurations are possible. For example, a wedge-shaped scan cone, or other similar shapes can be generated by the transceiver 1OA.
[0052] FIGURE 2-3 is a graphical representation of a scan plane 42. The scan plane 42 includes the peripheral scan lines 44 and 46, and an internal scan line 48 having a length r that extends outwardly from the transceiver 1OA. Thus, a selected point along the peripheral scan lines 44 and 46 and the internal scan line 48 can be defined with reference to the distance r and angular coordinate values φ and θ . The length r preferably extends to approximately 18 to 20 centimeters (cm), although any length is possible. Particular embodiments include approximately seventy-seven scan lines 48 that extend outwardly from the dome 34, although any number of scan lines is possible.
[0053] As described above, the angular movement of the transducer can be mechanically effected and/or it can be electronically or otherwise generated. In either case, the number of lines 48 and the length of the lines can vary, so that the tilt angle φ sweeps through angles approximately between -60° and +60° for a total arc of approximately 120°. In one particular embodiment, the transceiver 1OA can be configured to generate approximately seventy-seven scan lines between the first limiting scan line 44 and a second limiting scan line 46. In another particular embodiment, each of the scan lines has a length of approximately 18 to 20 centimeters (cm). The angular separation between adjacent scan lines 48 (FIGURE 2-2) can be uniform or non-uniform. For example, and in another particular embodiment, the angular separation φ \ and φ 2 (as shown in FIGURE 2-3) can be about 1.5°. Alternately, and in another particular embodiment, the angular separation φ 1 and φ % can be a sequence wherein adjacent angles can be ordered to include angles of 1.5°, 6.8°, 15.5°, 7.2°, and so on, where a 1.5° separation is between a first scan line and a second scan line, a 6.8° separation is between the second scan line and a third scan line, a 15.5° separation is between the third scan line and a fourth scan line, a 7.2° separation is between the fourth scan line and a fifth scan line, and so on. The angular separation between adjacent scan lines can also be a combination of uniform and non-uniform angular spacings, for example, a sequence of angles can be ordered to include 1.5°, 1.5°, 1.5°, 7.2°, 14.3°, 20.2°, 8.0°, 8.0°, 8.0°, 4.3°, 7.8°, and so on.
[0054] FIGURE 2-4 is a graphical representation of a plurality of scan lines emanating from the hand-held ultrasound transceiver 1OA forming a single scan plane 42 extending through a cross-section of a section of the abdominal aorta 51. The number and location of the internal scan lines emanating from the transceiver 1OA within a given scan plane 42 can thus be distributed at different positional coordinates about the axis line 11 as can be required to sufficiently visualize structures or images within the scan plane 42. As shown, four portions of off-centered shadow regions 49 are exhibited as irregular regions and a cross section of the abdominal aorta 51 region of interest (ROI) is depicted as in the center. Three portions can be viewable within the scan plane 42 in totality, and one can be truncated by the peripheral scan line 44.
[0055] FIGURE 3-1 depicts a partial schematic and partial isometric and side view of transceiver 1OB, and a scan cone array comprised of 3D-distributed scan lines in an alternate embodiment of an ultrasound system. A plurality of three-dimensional (3D) distributed scan lines emanating from a transceiver that cooperatively forms a scan cone 45. Each of the scan lines has a length r that projects outwardly from the transceiver 1OB of FIGURES 2-1-2-4. As illustrated, the transceiver 1OB emits 3D-distributed scan lines within the scan cone 30 that can be one-dimensional ultrasound A-lines. Other transceiver embodiments can also be configured to emit 3D-distributed scan lines. Taken as an aggregate, these 3D-distributed A-lines define the conical shape of the scan cone 45. The ultrasound scan cone 45 extends outwardly from the dome 34 of the transceiver 1OB centered about an axis line 43. The 3D-distributed scan lines of the scan cone 43 include a plurality of internal and peripheral scan lines that can be distributed within a volume defined by a perimeter of the scan cone 43. Accordingly, the peripheral scan lines 31A- 3 IF define an outer surface of the scan cone 45, while the internal scan lines 39A-39C can be distributed between the respective peripheral scan lines 31A-31F. Scan line 39B can be generally collinear with the axis 43, and the scan cone 45 can be generally and coaxially centered on the axis line 43.
[0056] The locations of the internal and peripheral scan lines can be further defined by an angular spacing from the center scan line 39B and between internal and peripheral scan lines. The angular spacing between scan line 39B and peripheral or internal scan lines can be designated by angle Φ and angular spacings between internal or peripheral scan lines can be designated by angle 0. The angles Φj, Φ2, and Φ3 respectively define the angular spacings from scan line 39B to scan lines 39A, 39C, and 3 ID. Similarly, angles 0ls 02, and 03 respectively define the angular spacings between scan line 31B and 31C, 31C and 39A, and 31D and 31E.
[0057] With continued reference to FIGURE 3-1, the plurality of peripheral scan lines 3 IA-E and the plurality of internal scan lines 39 A-D can be three dimensionally distributed A-lines (scan lines) that are not necessarily confined within a scan plane, but instead can sweep throughout the internal regions and along the periphery of the scan cone 45. Thus, a given point within the scan cone 45 can be identified by the coordinates r , Φ, and 0 whose values generally vary. The number and location of the internal scan lines emanating from the transceiver 1OB can thus be distributed within the scan cone 45 at different positional coordinates as required to sufficiently visualize structures or images within a region of interest (ROI) in a patient. The angular movement of the ultrasound transducer within the transceiver 1OB can be mechanically effected, and/or it can be electronically generated. In any case, the number of lines and the length of the lines can be uniform or otherwise vary, so that angle Φ sweeps through angles approximately between -60° between scan line 39B and 3 IA, and +60° between scan line 39B and 3 IB. Thus angle Φ in this example presents a total arc of approximately 120°.
[0058] In one embodiment, the transceiver 1OB can be configured to generate a plurality of 3D-distributed scan lines within the scan cone 45 having a length r, in one embodiment for example, of approximately 20 to 40 centimeters (cm).
[0059] FIGURE 3-2 illustrates a transceiver 1OC configured with a transducer designed to provide a fan-like (e.g., having planar sides) conic scan cone 47 utilizing C- mode and B-mode ultrasound modalities. The transceiver 1OC projects a series of B-mode scan planes 52A and 52B that oscillate like a pendulum between extremes within the scan cone 47. The B-mode scan planes 52 A, 52B may be derived from a plurality of scan lines similar to scan lines 44, 46, and 48 of FIGURES 2-3 and 2-4. The pendulum oscillating scan planes 52A and 52B can be arranged substantially at right angles to each other as depicted at axis crossing 57. The oscillating scan planes 52A and/or 52B can define a series of C-scan planes 55 that vary in depth location from the transceiver dome 34. The C-scan planes 55 move from the transducer vanishing point, and the B-scan planes angularly radiate from the transducer vanishing point. For transceiver 1OC users, a portion of the abdominal aorta taken as a C-mode shape is displayed on the transceiver display 16. The C-scan geometry shown as scan planes 55 present a substantially square-like ultrasound area within the scan cone 47. The C-Scan image information contained within scan planes 55 presents a cross-section view of the ultrasound at a particular depth probed by the transceiver 1OC. The C-mode may be more representative of a portion of the abdominal aorta than the actual whole of the length of the aorta. In this depiction, the C- Scan illustrates a cross-section view of the ultrasound at a particular depth to obtain a targeting image of the abdominal aorta. The targeting image may be more of a binary image showing the lines and spaces that are inside the aorta versus those that are outside of the aorta. The definition of C-mode image basically may be a plane parallel to the face of the transducer to obtain a projection image of the AAA region. The C-mode acquired projection image can yield abdominal aorta information not confined to simply one a single plane parallel to the transducer surface, but multiple planes denoted as C-scans. In the transceives substantially similar to the BVI9600 transceiver product, the C-mode acquired projection image may be binary, and can include a non-AAA region and an AAA region. The AAA region may be presented as an interpolated shape that may be generated from one side to the opposite side, for example the left most and the right most sides of a valid segmentation, or cut, the AAA region on all planes.
[0060] FIGURE 3-3 illustrates a partial isometric and schematic view of an ultrasound AAA detection system 70 utilizing a transceiver probe 1OD and console 74 combination. The AAA detection system 70 may be battery powered and portable and can also be referred to as the BVI9600 with AortaScreen Mode system. Other embodiments can include line power. The ultrasound transceiver 1OD may be configured to send out and receive ultrasound RF signals. The received RF may be transmitted to console 74. The DSP in console can process the RF information to extract the information relevant to the detected feature of each line. Each line is classified as being "available" or "non- available", as will be described below, based on a set threshold value. The classification of information contained in each scanline can be integrated with the image processing module for accurate and optimal positioning of the transceiver for segmentation and volume measurement.
[0061] The transceiver 1OD can have a transceiver display 36, housing 35 and dome 34 design similar to transceivers 1OA and 1OB, and may be in signal communication to console 74 via signal cable 89. The console 74 can be pivoted from console base 72. The console 74 includes a display 76, detection and operation function panel 78, and select panel 80. The detection and operation function provide for targeting the abdominal aorta, allow user voice annotation recording, retrieval and playback of previously recorded voice annotation files, and current and previously stored 3D and 2D scans. [0062] As illustrated in FIGURE 3-3, display 76 provides a screenshot 77C with a targeting icon 79A including cross hairs centered in a cross sectional depiction of an abdominal aorta. Other screen shots can appear in the display 76 depending on which function key is pressed in the function panel 78. A targeting icon screenshot 77B with a plurality of directional arrows can appear and flash to guide the user to move the transceiver 1OC to center the abdominal aorta or AAA. The targeting icon screenshot 77B can also appear on the display 36 of the transceiver 1OD. The targeting icon screenshot 77B similarly guides the user to place the transceiver 1OD to center the abdominal aorta or AAA or other region of interest as the directional indicator panel 32 did in transceiver 1OA of FIGURE 2-1 above.
[0063] As illustrated in FIGURE 3-3, the transceiver 1OD can be in wireless communication via wireless signal 82 with the wireless hub 84; the output from the transceiver 1OD can be delivered to a wireless hub 84 via wireless signal port 86. The wireless hub 84 can also charge batteries 88 for loading into the battery compartment (not shown) of console 74. In one embodiment, all the calculations can be performed in the imaging console 74. The BVI9600-embodiment system 70 does not require a computer or network to complete the analysis and imaging processing. In other embodiments, the system 70 can utilize the wireless hub 84 as a gateway to transmit transceiver 1OD acquired imaging information in local and Internet systems similar to that described in FIGURE 4 below.
[0064] FIGURE 4 is a schematic illustration of a server-accessed local area network and/or an internet system in communication with a plurality of ultrasound imaging systems. In an exemplary embodiment, an ultrasound system 90 can include one or more personal computer devices 52 that can be coupled to a server 96 by a communications system 95. The devices 52 can be, in turn, coupled to one or more ultrasound transceivers 1OA and/or 1OB and/or 1OC, for examples the ultrasound subsystems 94A-94D. Ultrasound based images of organs or other regions of interest derived from either the signals of echoes from fundamental frequency ultrasound thereof, can be shown within scan cone 40 or 45 presented on display 54. The server 96 can be operable to provide additional processing of ultrasound information, or it can be coupled to still other servers (not shown in FIGURE 4) and devices. Transceivers lOA/B/C can be in wireless communication with computer 52 in sub-system 94A, in wired signal communication in sub-system 94B, in wireless communication with computer 52 via receiving cradle 50 in sub-system 94C, or in wired communication with computer 52 via receiving cradle 50 in sub-system 94D .
[0065] An Internet system 114 can be coupled or otherwise in communication with the ultrasound sub-systems 94A-94D.
[0066] Scanning and Placement of Transceiver and Determination of "Availability":
[0067] During the field evaluation of AAA scan-mode, a 2D real time ultrasound instrument was determined to be preferably equipped to find optimal scan locations by searching the interested area with realtime B-mode image feedback. The BVI9600 model has been optimized to inform a user whether a 3D data set is valid or not, prior to collection of useful data. This was not a capability of the previous BVI6500 AAA mode, where only one B-mode image (cross section of an aorta) was provided for a user review after the 3D scan was completed. The BVI6500 B-mode image did not represent the 3D data set condition.
[0068] In contrast, the BVI9600 of the present invention investigates each scanline to determine whether the scanline contains any object information relevant to the region of interest or is just shadow blocked by air pocket. In the case that the scanline has any object information, it may be called "Available". By plotting these availabilities in the aiming screen in real-time and guiding a user with aiming arrow, the user can avoid the air blocked scan.
[0069] Referring now to FIGURE 5, a plurality of measurements along a patient's midline may be taken to find the maximal measurement of the abdominal aorta. In operation, gel may be applied at the midline of the body approximately 1 inch (2.5 cm) below the xiphoid (sternum); making sure that there is enough gel to cover the probe head 34. FIGURE 5 shows four such exemplary scan probe locations along the midline of a patient's torso. Finn pressure may be applied, noting that patients with excess tissue can require more pressure in order to obtain a good quality scan. The trigger/scan button 37 of the transceiver 1OA, for example, may be then depressed. An aim/scan screen appears, for example in display 36 of the transceiver (FIGURE 2.1) and/or display 76 on consol 74 and the instrument performs a 3D scan of the patient's region of interest (ROI).
[0070] When the scan is done, the aiming result and B-mode image may be displayed, for example as shown in an exemplary screen shot FIGURE 10.
[0071] In AAA ultrasound scanning, a subject is required to fast overnight to have less gas in their intestines so that the images are more shadow-free due to air block. Air blocks not only degrade the B-mode image quality, but also mislead the delineation of abdominal aorta which is an advantageous factor for volume estimation accuracy. Specifically, the acoustic energy can be blocked by a layer which has high acoustic impedance difference as between two media. Examples can include a fat layer under skin or air pockets in the intestines. The presence of fat and/or air pockets prevent information relevant to the region of interest (ROI), such as an AAA wall or lumen location, to be reflected to the transducer resulting in dark or black region underneath (See FIGURE 6 as an example). As a result scanlines (A-mode) may be investigated to determine whether useful information is contained therein or not. As shown in FIGURE 6-1, a B-mode image of an AAA may be in the same field as a shadow region due to the air-block or fat layer, two A-mode lines 60, 65. In A-mode line 60 there may be an AAA between depth 85 and 120 mm, therefore A-mode line 60 has useful information such as tissue, AAA, and acoustic enhancement. In A-mode line 65 however there may be no information due to the acoustic block near the depth of 25mm. FIGURES 6-2 and 6-3 are representative histograms of the A-mode lines 60 (FIGURE 6-2) and 65 (FIGURE 6-3) wherein, a threshold line set at 40 (67), in the exemplary histogram, is shown. The threshold line 67 may be used to determine whether the A-mode has available information or not. [0072] FIGURES 7-1-7-12 and 8 illustrate how the aorta aiming guide can determine "availability" and may be used to help the user correct the position of the transceiver to scan AAA.
[0073] In order to guide a user to the best scanhead position and orientation, the guidelines are used. A threshold may be established, which can vary from patient to patient, but in this example the threshold is 40 (as above), and may be related to the intensity number, for example, out of 256 levels after A-law compression. A searching range, in this non-limiting example is between 42% of A-mode length to 100% of A- mode length, the length of which may be related to experimental testing of optimal scan depths. Given these parameters, once set, if any A mode point within the searching range is above the established threshold but below an acoustic block, the A mode line may be determined as "Available". Based on the established threshold, "non-available" scanlines can be differentiated from "available" scanlines, for example, "non-available" are shadowed. In FIGURES 7-1- 7-12 A-mode lines are marked non-available by colored markings/dots underneath the B -mode image indicating the non-availability of that area. In FIGURE 7-10 for example, colored markings/dots underneath the B-mode image represent approximately 40% unavailable scanlines, compared to >90% unavailable scanlines in FIGURE 7-6.
[0074] The available and non-available regions from all 12 planes of this example of the B-mode display may be then plotted in C-mode representing the pseudo- horizontal cross section of the 3D scan cone composed by the 12 sectors as shown in FIGURE 8. FIGURE 8 shows the C-mode display of an "Availability Plot" 80 for
available and non-available regions. In this example, the "Availability Plot" 80 in C-mode (i.e. the view from the top of the scan cone) depicts available regions 85 as solid regions and non-available regions 87 as white and/or hashed lines. In this example, ~ 46% of the aggregate scanned regions are defined as "available." The percentage of available A- mode lines in the cone data, defined as "Availability" can be calculated using equation
(J )
Figure imgf000024_0001
[0076] The above calculation may be used to automatically determine whether the 3D full scan is appropriate for segmentation and display and without shadow within the region of interest.
[0077] In FIGURE 9-1, the C-mode display 81 can be partitioned into concentric circles, in this example, an inner circle 93 and an outer circle 97, to directionally position the transceiver to an available area, . Each quarter of the circle can be defined as a fan 98, from a perpendicular point defined by the cross-hairs in the middle; each overlapping fan expanding, for example, by 45° to both side in either direction . Each direction can have a fan, and the fans may be overlapped for the greatest coverage.
[0078] FIGURE 9-2 illustrates an embodiment of a decision tree utilizing the above schema. For example, in order to guide a user to an optimal scanhead position and orientation, the availability may be calculated, starting within the inner circle 93. If the inner circle 93 has a greater than or equal to 75% available scanlines, the scan is a "good scan." Otherwise, a direction may be suggested to re-position the scanhead. Repositioning may be suggested based on a criterion, generally in the direction corresponding to any inner fan and/or inner + outer fan availability determination that has the greatest calculated "availability." In other words, the arrow guide indicates the direction of greatest availability. First, for example, if any inner fan has availability greater than 75%, the suggested aiming direction may be the one whose inner fan availability is maximal. An arrow indicator guide may be a moderate and solid light indicator 121, 122, indicating that little movement is suggested. Second, if there is no inner fan that meets the above exemplary criteria, the suggested aiming direction may be based on the direction in which both the inner and outer fan has a maximum availability. In this case, the arrow indicator guide may be a flashing aggressive light indicator 124, 126, and the user may be guided in the recommended direction and orientation to achieve an optimized availability.
[0079] FIGURE 10 is an exemplary screen view 100 of the aiming guide of the BVI9600. In this example, a B-mode image 102 and a C-mode image 103 is shown. Due to the air block in the 3D data, shown in the C-mode display 103 as non-available (by hatched lines) 87, the user is directed to move in a NW direction, i.e. towards the upper left had quadrant 110 where the greatest availability and/or "available" scanlines may be detected to collect a better data set.
[0080] Turning again to FIGURE 9-1 and 9-2, the program 105 of an embodiment can determine whether to show a flashing arrow or solid arrow on the indicator panel dependent on the calculation of the percentage of availability of the image. For example: [0081] Condition 1 : Current Position = > 75% available at BLOCK 110 = good scan at BLOCK 112 = no movement is indicated = start AAA detection algorithm (See FIGUREI l).
[0082] Condition 2: Any inner fan = > 75% at BLOCK 118= re-position transceiver to aim towards greatest availability = Solid arrow direction indication at BLOCK 121 (See solid arrow 122 indicated in direction) = back to Condition 1 at BLOCK I lO.
[0083] Condition 3: All Inner fans = <^25% available within inner fans = query inner + outer fans = flashing arrow direction indication at BLOCK 124 (See flashing arrow 126 indicated in direction of maximal inner + outer fan availability) = back to Condition 1 at BLOCK 110.
[0084] Using the above positioning method, the user moves the instrument around on the abdomen to measure the abdominal aorta, which may be calculated from full three-dimensional scan cone or, can be a single two dimensional plane.
[0085] The use of the three-dimensional scan cone removes orientation requirements, permitting the user to position the device in any orientation. The user takes several three-dimensional image scans, moving along the patient's abdomen. After each scan, the volume/diameter of the section of the aorta covered by that scan may be displayed and the image may be stored if the diameter from the new scan is larger than any previous diameter. The image produced in this embodiment, whether it be form two dimensional planes or three dimension scan cone, can be optimally transmitted via the internet for remote enhance processing and rendering as is illustrated in FIGURE 4.
[0086] Volume Measurement Algorithms:
[0087] A. Front and Back Walls. Once an optimal position for AAA measurement is achieved based on the "availability" calculation 105 of compared scanlines, as described above and illustrated in FIGURE 9-2, an AAA detection algorithm 120 may be used to [0088] The fundamental AAA detection algorithm 120 used in the BVI 9600 device begins with process block Find Initial Wall 122 using A-mode ultrasound data that incorporates data-smoothing. Find Initial Wall 122 looks for the front and back walls of the abdominal aorta illustrated and described in FIGURE 12 below. After the front and back walls are found, a line passing through the center of the abdominal aorta may be determined in the following process block Find Centroid 124. This center abdominal aorta line or centroid may be used as a seed from which process block Fix Initial Walls 130 utilized, as illustrated and described in FIGURE 14 below. Fix Initial Walls 130 refines the initial wall points, removes any outliers, and fills gaps in the detected wall location regions. The Centroid 118 as shown in FIGURE 12-3, may be found by averaging the front and back walls. The maximum diameter (2A) 129 of a sphere that can fit within the AAA may be calculated as the maximum distance between the front and back walls. The Fix Walls algorithm 130 may be adjusted to comply with this spherical assumption to limit the segmentation of the AAA. That is, the given segmentation (from automatic detection) may be modified not to exceed the circle boundary whose diameter is the maximum of one, maximized slice. In FIGURE 12-1 and 12-2 B-mode images of a longitudinal section 12-1 and cross section 12-2 are shown with their exemplary segmentations for volume estimation. In this example, the AAA phantom has 4 cm, and 5 cm diameter in axial and lateral directions, respectively.
[0089] FIGURE 12-3 is a schematic representation of the application of the limited segmentation wherein a maximum diameter may be fixed from the centroid 118. Two lengths, length (2A) 129 the diameter, and length B 135, the radius, represent the maximum circle that can fit inside the AAA. The mathematical steps to calculate a modified segmentation parameter are described with reference to FIGURE 11, above. The volume calculation at BLOCK 134 may be the result of the aggregate calculation of volume resulting from the limited segmentation, in this example, represented by lines 136 and 138, i.e. those within the scanlines 145, 147 at points corresponding to the intersection of the front wall and the back wall of the artery being measured. At BLOCK 137 the diameter calculation may be derived by application mathematical relationship between volume and diameter as further discussed in relationship to FIGURE 13 below.
[0090] FIGURE 13 graphically represents the relationship between volume and diameter as determined using the below formulae. The volume of sphere volume is:
Figure imgf000028_0001
[0092] where, V is a volume and D is a diameter, and
[0093] where the diameter can be derived if the volume is know by manipulation of the above relationship:
[0094] D = Ij-V (3)
V π
[0095] The mathematical functions can be calculated automatically from a lookup table of diameter values such that the diameter calculations need not be performed for each volume under consideration. For example, the above relationship can be expressed as a look-up table, where once the volume (at BLOCK 134) is determined by the limited segmentation procedure, the diameter may be automatically correlated (at a BLOCK 137) with a pre-calculated diameter according to the above relationship.
[0096] B. AAA Detection Algorithm, Finding Front and Back Abdominal Aorta Walls and Centroid:
[0097] FIGURE 14 depicts a flowchart of the Find Initial Walls sub-algorithm of FIGURE 11. Find Initial Walls 122 process may be executed on every A-mode scan line and may be subjected to averaging and low-pass filtering using, for example, a 15 or 16 sample set beginning with process block 140. Next, a local gradient at process block 142 may be computed for each sample point using a central difference formulation taken for seven samples. The central difference formulation is defined by equations 4-9 (Eq. A-
9) below:
[0098] The standard central difference formula is given in Equation 4:
Figure imgf000029_0001
[00100] This formula assumes that the function is defined at the half-index, which is usually not the case. The solution is to widen the step between the samples to 2 and arrive at the equation in 5.
[00101] Λ, = ife+1H) (5)
The normalization factor is simply the distance between the two points. In Eq. 4 the distance separating the two means in the calculation was 1 , and in Eq. 5 the step between the two means is 2. The normalization of the gradient by the step size, while mathematically correct, incurs a cost in terms of operation. This operation may be neglected in the gradient calculation for the aortic wall detection algorithm with minimal effect: since the same calculation is performed for every data sample, every data sample can have the same error and thus the relative gradient values between different samples remain unchanged.
[00102] To further amplify wall locations, the gradient calculation may be expanded to three neighboring points to each side of the sample in question. This calculation is shown in Eq. 6. This calculation is simply the sum of three gradient approximations and thus the end result can be 12 times its normal value. This deviation from the true mathematical value has minimal effect since the calculation may be the same at each point and thus all the gradient values can be 12 times their usual value. An advantage to using the three neighboring points is that more information about the edge is included in the calculation, which can amplify the strong edges of the aorta and weaken the false-edges caused by the noise process in the image. [00103] dxt = X1+3 + xi+2 + xi+l - x,_, - jc, 2 - x,_3 (6)
[00104] The full calculation is written in Eq. 7. The first line shows the summation calculation to obtain the means, and the difference operations to obtain the gradient. The entire sum is normalized by 15, the number of points included in each local mean. The second line of the operation shows the result when the summations are simplified, and represents the maximal implementation of the calculation. This calculation incurs a cost of 23 additions or subtractions, 2 floating-point multiplications, 1 floating point division, and at least 1 temporary variable. This operation cost is incurred for 91% of the data samples.
Figure imgf000030_0001
15
[00106] The cost of the calculation can be reduced by not simplifying the summations and using a running sum operation. In that manner, only one mean may need to be calculated for each point, but that mean needs to be for the i+3 point. The running sum calculation uses the previous sum, and then corrects the sum by subtracting the old "left hand" end point and adding the new "right hand" end point. The operation is shown in Eq. 8. This running sum operation incurs a cost of 5 additions and subtractions.
[00107] X1+3 =
Figure imgf000030_0002
= X1+3-, - X1+^ + xl+3+1 = x,+2 - x,_5 + JC1+10 (8)
Figure imgf000030_0003
Since the running sum was calculated for the i+3 point, all average values are available for the gradient calculation. This calculation is shown in Equation 9:
[00108] dxt = " X'-3 ~ X'-2 ~ X'-' + X'+1 + X'+2 + X'+3 (9)
[00109] This equation has the same form as the one in Eq. 6 except for the normalization factor of 16. This normalization factor is not a result of the gradient operation, but rather it is the normalization factor mean calculation. The factor of 16 is used instead of the standard value of 15 that one would expect in a 15 -point average for this simple reason: using a factor of 16 eliminates floating-point division. If the means are normalized by 16, then the division operation can be replaced by a "righf'-shift by 4 at a significantly lower cost to the embedded system. Therefore the gradient operation has eleven additions and subtractions and one shift by 4.
[00110] Adding the operational cost of the running sum calculation gives an overall cost of 16 additions and subtractions and the shift. A clear advantage in this simplification is the elimination of multiplication and division from the operation.
[00111] Returning to FIGURE 14, the results from local gradient 142 are subjected to loop limit processing between blocks 144 and 158 to obtain the best front wall and back wall pair for each scan line denoted as a tissue gradient or tissue delta. The best front wall and back wall pair on each line may be defined as the front wall and back wall pair for which the pixel intensity difference in the back wall gradient and front wall gradient is the maximum and the smallest local average between front wall and back wall pair is the minimum.
[00112] The loop limit processing begins with loop limit block 144 that receives pixel values for each sample in the detection region and subjects the pixel intensity values to determine whether the gradient is minimum at decision diamond 146. If affirmative, then the pixel values may be ascertained whether it's the best front wall-back wall (FW/BW) candidate combination at decision diamond 147. If affirmative, the FW/BW candidate pair may be saved and loop limit processing returns to limit block 144. If negative, at process block 152, the Front Wall pixel value may be saved and another back wall candidate may be sought with a subsequent return to loop limit block 152.
[00113] Returning to decision diamond 146, if the answer is negative for "Is gradient Minimum?", sub-algorithm 122 continues to decision diamond 156 to determine whether the back wall and the gradient is maximum. If affirmative, at process block 154, a candidate BW/FW pair is established and sub-algorithm re-routes to loop limit block 144. If negative, the end of analysis for a particular FW/BW candidate occurs at loop limit block 158 either routes back to the limit loop block 144 or exits to find Centroid 124.
[00114] Formulations relating to Find Centroid 124 may be based on coordinate geometries described in equations 10 and 11 utilizing coordinate conversions. The coordinate conversions are shown in Eq. 10 where 38 is the index of the broadside beam (the ultrasound ray when φ=0), φ is the index of the line, θ is the angle of the plane. The plane angle is shifted by π to ensure that the sign of the x and y coordinates match the true location in space.
Figure imgf000032_0001
[00116] The trigonometric functions can be calculated for a table of θ values such that the cosine and sine calculations need not be performed for each of the points under consideration. The closest plane can be found by finding the shortest vector from each plane to the centroid. The shortest vector from a plane to a point can be the perpendicular to the projection of the centroid on the plane. The projection of the centroid on the plane is defined as the dot product of the centroid vector, c, with the plane definition vector, a, divided by the length of the plane definition vector. If the plane definition vector is a unit vector, then the division operation is unnecessary. To find the perpendicular to the projection, it is sufficient to subtract the projection vector from the centroid vector as shown in Eq. 11 :
Figure imgf000032_0002
Figure imgf000032_0003
[00118] The length of these projections can be found by calculating the Euclidean norm for each line. The Euclidean norm may be more commonly known as the length or magnitude of the vector. To find the plane closest to the centroid, calculate the lengths for the perpendicular to the projection of the centroid on each plane, and take the plane with the shortest of these lengths. [00119] While the preferred embodiment of the invention has been illustrated and described, as noted above, many changes can be made without departing from the spirit and scope of the invention. For example, it should be understood that various changes, modifications, and substitutions can be incorporated in the apparatus embodiment to achieve the ultrasonic, volumetric determination of the abdominal aorta to thus apply that to the aortic diameter reading. Accordingly, the scope of the invention is not limited by the disclosure of the preferred embodiment.

Claims

[00120] The embodiments of the invention in which an exclusive property or privilege is claimed are defined as follows:
1. A system to detect an abdominal aortic aneurysm in a region of interest comprising: an ultrasound transceiver positioned to deliver ultrasound energy and receive echoes of the ultrasound energy across a plurality of scan planes; a processing device configured to signal process the received echoes and characterize detected signals across the plurality of scan plans, wherein a percentage of availability of scan planes scanned is calculated based on echos received from scan planes; a display configured to present a visual depiction of the percentage of availability of the region of interest contained in a plurality of scan planes; and a display configured to present positioning information of the ultrasound transceiver based on the visual depiction of the percentage of availability.
2. The system of claim 1 , wherein the processing device guides a user to position the ultrasound transceiver towards regions of interest where the calculated percentage of availability in a plurality of scanlines is 75-100% availability.
3. The system of claim 1, wherein the processing device includes at least one of a C- mode component, a B-mode component, and an A-mode component.
4. A system for abdominal aortic aneurysm evaluation and monitoring in a region of interest of an abdominal aorta, comprising: a hand-held data collection device which comprises: at least one ultrasound transceiver to deliver ultrasound energy and receive echoes of the ultrasound energy across a plurality of scan planes for obtaining three-dimensional ultrasound scan information; a display means for indicating positioning information received from the ultrasound transceiver at a first measured location of the abdominal aorta, wherein the display is a visual depiction of a calculated percentage availability of the first measured location; a guide means for guiding a user to position the ultrasound transceiver over a region of interest of the aorta based on the calculated percentage availability; a processor for converting the received scan plane information about the region of interest obtained from the transceiver; a processor for determining aorta volume information from the converted scan plan information; and a calculation circuit for calculating the diameter of the aorta at the region of interest from the aorta volume information.
5. The system of claim 4, wherein the guide means guides the user to position the ultrasound transceiver towards regions of interest where the calculated percentage of availability contained in a plurality of scans planes is 75-100% availability.
6. A method to detect and measure an abdominal aortic aneurysm comprising: transmitting ultrasound energy to the abdominal aorta; collecting ultrasound echoes returning from the abdominal aorta; generating signals from the ultrasound echoes; identifying fundamental signals from the generated signals; processing the fundamental signals using algorithms designed for fundamental signals; presenting an image of a region of interest of the abdominal aorta; calculating the volume of the abdominal aorta; and calculating the diameter of the abdominal aorta from the calculated volume of the abdominal aorta.
7. The method of claim 6, wherein a guide means guides the user to position the ultrasound transceiver to at least one position to measure the volume of the abdominal aorta.
8. Computer readable media having instructions to execute a method to detect and measure an abdominal aortic aneurysm, comprising the steps of: transmitting ultrasound energy to the abdominal aorta; collecting ultrasound echoes returning from the abdominal aorta; generating signals from the ultrasound echoes; identifying fundamental signals from the generated signals; processing the fundamental using algorithms designed for fundamental signals; presenting an image of the abdominal aorta; calculating measurements of the volume of the abdominal aorta; and calculating measurements of the diameter of the abdominal aorta from the calculated volume of the abdominal aorta.
9. The computer readable media of claim 8, wherein the computer readable media further comprise instructions to execute the steps of comparing fundamental signals returning from non available regions from available regions and calculating percentage of availability; presenting an image of non available regions and available regions to a user based on the calculated percentage of availibilty; and indicating to the user a position for an ultrasound transmitting device wherein the percentage of availability is maximized.
PCT/US2009/053201 2008-08-07 2009-08-07 Device, system, and method to measure abdominal aortic aneurysm diameter WO2010017508A1 (en)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CA2732997A CA2732997C (en) 2008-08-07 2009-08-07 Device, system, and method to measure abdominal aortic aneurysm diameter
EP09805626.0A EP2323559A4 (en) 2008-08-07 2009-08-07 Device, system, and method to measure abdominal aortic aneurysm diameter
JP2011522292A JP5658151B2 (en) 2008-08-07 2009-08-07 Apparatus, system and method for measuring the diameter of an abdominal aortic aneurysm

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
US8715208P 2008-08-07 2008-08-07
US61/087,152 2008-08-07
US9400308P 2008-09-03 2008-09-03
US61/094,003 2008-09-03

Publications (1)

Publication Number Publication Date
WO2010017508A1 true WO2010017508A1 (en) 2010-02-11

Family

ID=41663998

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2009/053201 WO2010017508A1 (en) 2008-08-07 2009-08-07 Device, system, and method to measure abdominal aortic aneurysm diameter

Country Status (5)

Country Link
US (1) US8133181B2 (en)
EP (1) EP2323559A4 (en)
JP (1) JP5658151B2 (en)
CA (2) CA2732997C (en)
WO (1) WO2010017508A1 (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015091368A1 (en) * 2013-12-20 2015-06-25 Koninklijke Philips N.V. Ultrasound imaging assembly and method for displaying ultrasound images
RU2689769C1 (en) * 2018-01-31 2019-05-28 Государственное Бюджетное Учреждение "Санкт-Петербургский Научно-Исследовательский Институт Скорой Помощи Им. И.И. Джанелидзе" Method for measuring length of proximal cervical abdominal aortic aneurysm
EP3626160A1 (en) 2018-09-24 2020-03-25 Dianova A/S A handheld imaging element with a movement sensor
US20210113194A1 (en) * 2019-10-17 2021-04-22 Verathon Inc. Systems and methods for ultrasound scanning
EP3836154A1 (en) 2019-12-13 2021-06-16 Morten Bo Søndergaard Svendsen A method and an apparatus for imaging curved surfaces
EP3839883A1 (en) 2019-12-18 2021-06-23 Morten Bo Søndergaard Svendsen A method and a system for obtaining image data of an object

Families Citing this family (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090062644A1 (en) * 2002-06-07 2009-03-05 Mcmorrow Gerald System and method for ultrasound harmonic imaging
WO2011074271A1 (en) * 2009-12-18 2011-06-23 パナソニック株式会社 Ultrasonic diagnostic device, and region-to-be-detected image display method and measurement method using same
US9585638B2 (en) 2010-03-06 2017-03-07 Russell Dewey System for classifying palpable soft tissue masses using A-mode echographs
JP6037447B2 (en) * 2012-03-12 2016-12-07 東芝メディカルシステムズ株式会社 Ultrasonic diagnostic equipment
US20140046202A1 (en) * 2012-08-13 2014-02-13 California Institute Of Technology Noninvasive systems for aortic aneurysm evaluation
US9367965B2 (en) * 2012-10-05 2016-06-14 Volcano Corporation Systems and methods for generating images of tissue
US9091628B2 (en) 2012-12-21 2015-07-28 L-3 Communications Security And Detection Systems, Inc. 3D mapping with two orthogonal imaging views
CA2974377C (en) 2015-01-23 2023-08-01 The University Of North Carolina At Chapel Hill Apparatuses, systems, and methods for preclinical ultrasound imaging of subjects
JP6792340B2 (en) * 2016-03-29 2020-11-25 ザイオソフト株式会社 Medical image processing equipment, medical image processing methods, and medical image processing programs
WO2018165448A1 (en) 2017-03-08 2018-09-13 Obma Padraic A method for identifying human joint characteristics
EP3435382A1 (en) 2017-07-27 2019-01-30 Koninklijke Philips N.V. Imaging method, controller and imaging system, for monitoring a patient post evar
CN107361790B (en) * 2017-07-31 2020-07-10 天津迈达医学科技股份有限公司 Intelligent bladder positioning method based on ultrasonic image
WO2019070812A1 (en) * 2017-10-04 2019-04-11 Verathon Inc. Multi-plane and multi-mode visualization of an area of interest during aiming of an ultrasound probe
US11278259B2 (en) 2018-02-23 2022-03-22 Verathon Inc. Thrombus detection during scanning
CN111885960B (en) * 2018-03-13 2023-10-13 韦拉索恩股份有限公司 Universal interlacing of ultrasound probes
JP6483910B1 (en) * 2018-11-28 2019-03-13 メロディ・インターナショナル株式会社 Ultrasonic inspection equipment
US11684344B2 (en) 2019-01-17 2023-06-27 Verathon Inc. Systems and methods for quantitative abdominal aortic aneurysm analysis using 3D ultrasound imaging
EP3968274A1 (en) * 2020-09-14 2022-03-16 Tata Consultancy Services Limited Method and system for asset inspection using unmanned aerial vehicles

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040054280A1 (en) * 2002-09-18 2004-03-18 Mcmorrow Gerald J. Three-dimensional system for abdominal aortic aneurysm evaluation
US20040076317A1 (en) * 1998-07-23 2004-04-22 David Roberts Method and apparatus for the non-invasive imaging of anatomic tissue structures
US20080139938A1 (en) * 2002-06-07 2008-06-12 Fuxing Yang System and method to identify and measure organ wall boundaries

Family Cites Families (218)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4005A (en) * 1845-04-22 Improvement in the manufacture of indurated rubber fabrics
US8013A (en) * 1851-04-01 Improvement in wheat-fans
US4007A (en) * 1845-04-22 Improvement in separating tannin and coloring-matter in quercitron-bark
US3613069A (en) 1969-09-22 1971-10-12 Gen Dynamics Corp Sonar system
US4431007A (en) * 1981-02-04 1984-02-14 General Electric Company Referenced real-time ultrasonic image display
FR2551611B1 (en) 1983-08-31 1986-10-24 Labo Electronique Physique NOVEL ULTRASONIC TRANSDUCER STRUCTURE AND ULTRASONIC ECHOGRAPHY MEDIA EXAMINATION APPARATUS COMPRISING SUCH A STRUCTURE
US4556066A (en) 1983-11-04 1985-12-03 The Kendall Company Ultrasound acoustical coupling pad
US4926871A (en) * 1985-05-08 1990-05-22 International Biomedics, Inc. Apparatus and method for non-invasively and automatically measuring the volume of urine in a human bladder
US4757821A (en) 1986-11-12 1988-07-19 Corazonix Corporation Omnidirectional ultrasonic probe
AU591360B2 (en) 1986-11-13 1989-11-30 National Aeronautics And Space Administration - Nasa Ultrasonic quantification of a property internal to a living subject
US4844080A (en) 1987-02-19 1989-07-04 Michael Frass Ultrasound contact medium dispenser
US4821210A (en) * 1987-04-02 1989-04-11 General Electric Co. Fast display of three-dimensional images
US5005126A (en) * 1987-04-09 1991-04-02 Prevail, Inc. System and method for remote presentation of diagnostic image information
US5058591A (en) 1987-11-10 1991-10-22 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Rapidly quantifying the relative distention of a human bladder
US5473555A (en) 1988-08-18 1995-12-05 Hewlett-Packard Company Method and apparatus for enhancing frequency domain analysis
DE68923448T2 (en) 1988-08-30 1995-12-07 Fujitsu Ltd COUPLER FOR ULTRASONIC TRANSDUCER PROBE.
US5159931A (en) 1988-11-25 1992-11-03 Riccardo Pini Apparatus for obtaining a three-dimensional reconstruction of anatomic structures through the acquisition of echographic images
US5060515A (en) 1989-03-01 1991-10-29 Kabushiki Kaisha Toshiba Image signal processing circuit for ultrasonic imaging apparatus
US5299577A (en) * 1989-04-20 1994-04-05 National Fertility Institute Apparatus and method for image processing including one-dimensional clean approximation
FR2650071B1 (en) * 1989-07-20 1991-09-27 Asulab Sa PROCESS FOR PROCESSING AN ELECTRICAL SIGNAL
US5151856A (en) 1989-08-30 1992-09-29 Technion R & D Found. Ltd. Method of displaying coronary function
EP0420758B1 (en) * 1989-09-29 1995-07-26 Terumo Kabushiki Kaisha Ultrasonic coupler and method for production thereof
US5125410A (en) * 1989-10-13 1992-06-30 Olympus Optical Co., Ltd. Integrated ultrasonic diagnosis device utilizing intra-blood-vessel probe
US5644513A (en) 1989-12-22 1997-07-01 Rudin; Leonid I. System incorporating feature-oriented signal enhancement using shock filters
US5148809A (en) 1990-02-28 1992-09-22 Asgard Medical Systems, Inc. Method and apparatus for detecting blood vessels and displaying an enhanced video image from an ultrasound scan
US6196226B1 (en) * 1990-08-10 2001-03-06 University Of Washington Methods and apparatus for optically imaging neuronal tissue and activity
US5235985A (en) 1992-04-30 1993-08-17 Mcmorrow Gerald J Automatic bladder scanning apparatus
DE9209877U1 (en) 1992-07-22 1993-11-25 Stoeger Helmut Pressure switch
US5381794A (en) * 1993-01-21 1995-01-17 Aloka Co., Ltd. Ultrasonic probe apparatus
DK0683636T3 (en) 1993-02-10 1999-10-18 Reilly Daniel Joseph O Bag for dispensing fluid material
DE69431741T2 (en) 1993-03-12 2003-09-11 Toshiba Kawasaki Kk Device for medical treatment with ultrasound
US5898793A (en) * 1993-04-13 1999-04-27 Karron; Daniel System and method for surface rendering of internal structures within the interior of a solid object
US5435310A (en) 1993-06-23 1995-07-25 University Of Washington Determining cardiac wall thickness and motion by imaging and three-dimensional modeling
US5601084A (en) * 1993-06-23 1997-02-11 University Of Washington Determining cardiac wall thickness and motion by imaging and three-dimensional modeling
JPH07136162A (en) 1993-11-17 1995-05-30 Fujitsu Ltd Ultrasonic coupler
US5465721A (en) 1994-04-22 1995-11-14 Hitachi Medical Corporation Ultrasonic diagnostic apparatus and ultrasonic diagnosis method
JP3723665B2 (en) * 1997-07-25 2005-12-07 フクダ電子株式会社 Ultrasonic diagnostic equipment
US5698549A (en) 1994-05-12 1997-12-16 Uva Patent Foundation Method of treating hyperactive voiding with calcium channel blockers
US5645077A (en) 1994-06-16 1997-07-08 Massachusetts Institute Of Technology Inertial orientation tracker apparatus having automatic drift compensation for tracking human head and other similarly sized body
US5615680A (en) * 1994-07-22 1997-04-01 Kabushiki Kaisha Toshiba Method of imaging in ultrasound diagnosis and diagnostic ultrasound system
EP0696435A3 (en) 1994-08-10 1997-03-12 Hewlett Packard Co Utrasonic probe
US5972023A (en) 1994-08-15 1999-10-26 Eva Corporation Implantation device for an aortic graft method of treating aortic aneurysm
US5526816A (en) * 1994-09-22 1996-06-18 Bracco Research S.A. Ultrasonic spectral contrast imaging
US5503152A (en) * 1994-09-28 1996-04-02 Tetrad Corporation Ultrasonic transducer assembly and method for three-dimensional imaging
US5487388A (en) * 1994-11-01 1996-01-30 Interspec. Inc. Three dimensional ultrasonic scanning devices and techniques
US6690963B2 (en) 1995-01-24 2004-02-10 Biosense, Inc. System for determining the location and orientation of an invasive medical instrument
US5575286A (en) 1995-03-31 1996-11-19 Siemens Medical Systems, Inc. Method and apparatus for generating large compound ultrasound image
US6151404A (en) 1995-06-01 2000-11-21 Medical Media Systems Anatomical visualization system
US5503153A (en) * 1995-06-30 1996-04-02 Siemens Medical Systems, Inc. Noise suppression method utilizing motion compensation for ultrasound images
US5588435A (en) 1995-11-22 1996-12-31 Siemens Medical Systems, Inc. System and method for automatic measurement of body structures
US5841889A (en) 1995-12-29 1998-11-24 General Electric Company Ultrasound image texture control using adaptive speckle control algorithm
JP3580627B2 (en) * 1996-01-29 2004-10-27 株式会社東芝 Ultrasound diagnostic equipment
US5851186A (en) 1996-02-27 1998-12-22 Atl Ultrasound, Inc. Ultrasonic diagnostic imaging system with universal access to diagnostic information and images
AU1983397A (en) * 1996-02-29 1997-09-16 Acuson Corporation Multiple ultrasound image registration system, method and transducer
US5806521A (en) 1996-03-26 1998-09-15 Sandia Corporation Composite ultrasound imaging apparatus and method
US5605155A (en) * 1996-03-29 1997-02-25 University Of Washington Ultrasound system for automatically measuring fetal head size
US5735282A (en) * 1996-05-30 1998-04-07 Acuson Corporation Flexible ultrasonic transducers and related systems
US6962566B2 (en) 2001-04-19 2005-11-08 Sonosite, Inc. Medical diagnostic ultrasound instrument with ECG module, authorization mechanism and methods of use
US6569101B2 (en) * 2001-04-19 2003-05-27 Sonosite, Inc. Medical diagnostic ultrasound instrument with ECG module, authorization mechanism and methods of use
US5846202A (en) 1996-07-30 1998-12-08 Acuson Corporation Ultrasound method and system for imaging
US6343936B1 (en) * 1996-09-16 2002-02-05 The Research Foundation Of State University Of New York System and method for performing a three-dimensional virtual examination, navigation and visualization
US7194117B2 (en) * 1999-06-29 2007-03-20 The Research Foundation Of State University Of New York System and method for performing a three-dimensional virtual examination of objects, such as internal organs
US5776063A (en) 1996-09-30 1998-07-07 Molecular Biosystems, Inc. Analysis of ultrasound images in the presence of contrast agent
US5903664A (en) * 1996-11-01 1999-05-11 General Electric Company Fast segmentation of cardiac images
US5738097A (en) * 1996-11-08 1998-04-14 Diagnostics Ultrasound Corporation Vector Doppler system for stroke screening
US6030344A (en) * 1996-12-04 2000-02-29 Acuson Corporation Methods and apparatus for ultrasound image quantification
US5782767A (en) 1996-12-31 1998-07-21 Diagnostic Ultrasound Corporation Coupling pad for use with medical ultrasound devices
US6122538A (en) 1997-01-16 2000-09-19 Acuson Corporation Motion--Monitoring method and system for medical devices
US5892843A (en) * 1997-01-21 1999-04-06 Matsushita Electric Industrial Co., Ltd. Title, caption and photo extraction from scanned document images
US6045508A (en) * 1997-02-27 2000-04-04 Acuson Corporation Ultrasonic probe, system and method for two-dimensional imaging or three-dimensional reconstruction
US5772604A (en) * 1997-03-14 1998-06-30 Emory University Method, system and apparatus for determining prognosis in atrial fibrillation
US6117080A (en) 1997-06-04 2000-09-12 Atl Ultrasound Ultrasonic imaging apparatus and method for breast cancer diagnosis with the use of volume rendering
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
US6008813A (en) 1997-08-01 1999-12-28 Mitsubishi Electric Information Technology Center America, Inc. (Ita) Real-time PC based volume rendering system
US6312379B1 (en) 1997-08-15 2001-11-06 Acuson Corporation Ultrasonic harmonic imaging system and method using waveform pre-distortion
US5945770A (en) 1997-08-20 1999-08-31 Acuson Corporation Multilayer ultrasound transducer and the method of manufacture thereof
US5928151A (en) 1997-08-22 1999-07-27 Acuson Corporation Ultrasonic system and method for harmonic imaging in three dimensions
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
US6148095A (en) 1997-09-08 2000-11-14 University Of Iowa Research Foundation Apparatus and method for determining three-dimensional representations of tortuous vessels
US6325758B1 (en) 1997-10-27 2001-12-04 Nomos Corporation Method and apparatus for target position verification
US5971923A (en) 1997-12-31 1999-10-26 Acuson Corporation Ultrasound system and method for interfacing with peripherals
US5964710A (en) 1998-03-13 1999-10-12 Srs Medical, Inc. System for estimating bladder volume
US6213949B1 (en) * 1999-05-10 2001-04-10 Srs Medical Systems, Inc. System for estimating bladder volume
JPH11267121A (en) * 1998-03-20 1999-10-05 Ge Yokogawa Medical Systems Ltd Method and device for ultrasonic photographing
US6302845B2 (en) 1998-03-20 2001-10-16 Thomas Jefferson University Method and system for pressure estimation using subharmonic signals from microbubble-based ultrasound contrast agents
US6200266B1 (en) * 1998-03-31 2001-03-13 Case Western Reserve University Method and apparatus for ultrasound imaging using acoustic impedance reconstruction
US5980459A (en) 1998-03-31 1999-11-09 General Electric Company Ultrasound imaging using coded excitation on transmit and selective filtering of fundamental and (sub)harmonic signals on receive
JP4116143B2 (en) 1998-04-10 2008-07-09 株式会社東芝 Ultrasonic diagnostic equipment
US6048312A (en) * 1998-04-23 2000-04-11 Ishrak; Syed Omar Method and apparatus for three-dimensional ultrasound imaging of biopsy needle
US6102858A (en) 1998-04-23 2000-08-15 General Electric Company Method and apparatus for three-dimensional ultrasound imaging using contrast agents and harmonic echoes
US6511325B1 (en) * 1998-05-04 2003-01-28 Advanced Research & Technology Institute Aortic stent-graft calibration and training model
JP4260920B2 (en) 1998-05-13 2009-04-30 株式会社東芝 Ultrasonic diagnostic equipment
US6511426B1 (en) * 1998-06-02 2003-01-28 Acuson Corporation Medical diagnostic ultrasound system and method for versatile processing
US6359190B1 (en) * 1998-06-29 2002-03-19 The Procter & Gamble Company Device for measuring the volume of a body cavity
US6071242A (en) * 1998-06-30 2000-06-06 Diasonics Ultrasound, Inc. Method and apparatus for cross-sectional color doppler volume flow measurement
CA2338735A1 (en) 1998-07-21 2000-02-03 Acoustic Sciences Associates Synthetic structural imaging and volume estimation of biological tissue organs
US6346124B1 (en) * 1998-08-25 2002-02-12 University Of Florida Autonomous boundary detection system for echocardiographic images
US5993390A (en) 1998-09-18 1999-11-30 Hewlett- Packard Company Segmented 3-D cardiac ultrasound imaging method and apparatus
US6126598A (en) 1998-10-01 2000-10-03 Atl Ultrasound, Inc. Ultrasonic diagnostic imaging system with adaptive spatial compounding
US6286513B1 (en) 1998-10-22 2001-09-11 Jessie L. Au Methods for treating superficial bladder carcinoma
JP2000126181A (en) 1998-10-27 2000-05-09 Mitani Sangyo Co Ltd Method for extracting and treating tumor
JP2000126182A (en) 1998-10-27 2000-05-09 Mitani Sangyo Co Ltd Tumor diagnosing method
JP2000126178A (en) 1998-10-27 2000-05-09 Mitani Sangyo Co Ltd Method of quantifying stereo surface shape and automatic identification method of malignant tumor
US6545678B1 (en) * 1998-11-05 2003-04-08 Duke University Methods, systems, and computer program products for generating tissue surfaces from volumetric data thereof using boundary traces
US6524249B2 (en) * 1998-11-11 2003-02-25 Spentech, Inc. Doppler ultrasound method and apparatus for monitoring blood flow and detecting emboli
JP2000139917A (en) * 1998-11-12 2000-05-23 Toshiba Corp Ultrasonograph
US6159150A (en) 1998-11-20 2000-12-12 Acuson Corporation Medical diagnostic ultrasonic imaging system with auxiliary processor
US6042545A (en) * 1998-11-25 2000-03-28 Acuson Corporation Medical diagnostic ultrasound system and method for transform ultrasound processing
US6272469B1 (en) 1998-11-25 2001-08-07 Ge Medical Systems Global Technology Company, Llc Imaging system protocol handling method and apparatus
WO2000039576A1 (en) 1998-12-23 2000-07-06 Image Guided Technologies, Inc. A hybrid 3-d probe tracked by multiple sensors
US6193657B1 (en) * 1998-12-31 2001-02-27 Ge Medical Systems Global Technology Company, Llc Image based probe position and orientation detection
GB9901270D0 (en) 1999-01-21 1999-03-10 Quadrant Healthcare Uk Ltd Method and apparatus for ultrasound contrast imaging
JP2000210286A (en) 1999-01-25 2000-08-02 Kasei Optonix Co Ltd Urination alarm unit
US6385332B1 (en) 1999-02-19 2002-05-07 The John P. Roberts Research Institute Automated segmentation method for 3-dimensional ultrasound
US6213951B1 (en) * 1999-02-19 2001-04-10 Acuson Corporation Medical diagnostic ultrasound method and system for contrast specific frequency imaging
US6544181B1 (en) * 1999-03-05 2003-04-08 The General Hospital Corporation Method and apparatus for measuring volume flow and area for a dynamic orifice
US6142942A (en) 1999-03-22 2000-11-07 Agilent Technologies, Inc. Ultrasound imaging system and method employing an adaptive filter
US6400848B1 (en) * 1999-03-30 2002-06-04 Eastman Kodak Company Method for modifying the perspective of a digital image
US6210327B1 (en) * 1999-04-28 2001-04-03 General Electric Company Method and apparatus for sending ultrasound image data to remotely located device
US6259945B1 (en) * 1999-04-30 2001-07-10 Uromed Corporation Method and device for locating a nerve
US6110111A (en) 1999-05-26 2000-08-29 Diagnostic Ultrasound Corporation System for quantizing bladder distension due to pressure using normalized surface area of the bladder
US6063033A (en) * 1999-05-28 2000-05-16 General Electric Company Ultrasound imaging with higher-order nonlinearities
US6575907B1 (en) 1999-07-12 2003-06-10 Biomedicom, Creative Biomedical Computing Ltd. Determination of fetal weight in utero
US6778690B1 (en) 1999-08-13 2004-08-17 Hanif M. Ladak Prostate boundary segmentation from 2D and 3D ultrasound images
US6254539B1 (en) 1999-08-26 2001-07-03 Acuson Corporation Transducer motion compensation in medical diagnostic ultrasound 3-D imaging
US6264609B1 (en) 1999-09-15 2001-07-24 Wake Forest University Ultrasound apparatus and method for tissue characterization
US7510536B2 (en) 1999-09-17 2009-03-31 University Of Washington Ultrasound guided high intensity focused ultrasound treatment of nerves
US6277073B1 (en) 1999-09-23 2001-08-21 Acuson Corporation Medical diagnostic ultrasound imaging method and system using simultaneously transmitted ultrasound beams
US6443894B1 (en) 1999-09-29 2002-09-03 Acuson Corporation Medical diagnostic ultrasound system and method for mapping surface data for three dimensional imaging
US6610013B1 (en) 1999-10-01 2003-08-26 Life Imaging Systems, Inc. 3D ultrasound-guided intraoperative prostate brachytherapy
US6440071B1 (en) 1999-10-18 2002-08-27 Guided Therapy Systems, Inc. Peripheral ultrasound imaging system
US6235038B1 (en) * 1999-10-28 2001-05-22 Medtronic Surgical Navigation Technologies System for translation of electromagnetic and optical localization systems
US6338716B1 (en) * 1999-11-24 2002-01-15 Acuson Corporation Medical diagnostic ultrasonic transducer probe and imaging system for use with a position and orientation sensor
US6466817B1 (en) * 1999-11-24 2002-10-15 Nuvasive, Inc. Nerve proximity and status detection system and method
US6350239B1 (en) * 1999-12-28 2002-02-26 Ge Medical Systems Global Technology Company, Llc Method and apparatus for distributed software architecture for medical diagnostic systems
US6491631B2 (en) 2001-01-11 2002-12-10 General Electric Company Harmonic golay-coded excitation with differential pulsing for diagnostic ultrasound imaging
JP2003521341A (en) 2000-01-31 2003-07-15 アー.ヤー. アンゲルセン、ビョルン Correction of phase plane aberration and pulse reverberation in medical ultrasound imaging
US6515657B1 (en) * 2000-02-11 2003-02-04 Claudio I. Zanelli Ultrasonic imager
US6406431B1 (en) 2000-02-17 2002-06-18 Diagnostic Ultasound Corporation System for imaging the bladder during voiding
US6494841B1 (en) 2000-02-29 2002-12-17 Acuson Corporation Medical diagnostic ultrasound system using contrast pulse sequence imaging
US6551246B1 (en) * 2000-03-06 2003-04-22 Acuson Corporation Method and apparatus for forming medical ultrasound images
US6511427B1 (en) * 2000-03-10 2003-01-28 Acuson Corporation System and method for assessing body-tissue properties using a medical ultrasound transducer probe with a body-tissue parameter measurement mechanism
US6238344B1 (en) * 2000-03-30 2001-05-29 Acuson Corporation Medical diagnostic ultrasound imaging system with a wirelessly-controlled peripheral
US6440072B1 (en) 2000-03-30 2002-08-27 Acuson Corporation Medical diagnostic ultrasound imaging system and method for transferring ultrasound examination data to a portable computing device
US6503204B1 (en) * 2000-03-31 2003-01-07 Acuson Corporation Two-dimensional ultrasonic transducer array having transducer elements in a non-rectangular or hexagonal grid for medical diagnostic ultrasonic imaging and ultrasound imaging system using same
US20020016545A1 (en) * 2000-04-13 2002-02-07 Quistgaard Jens U. Mobile ultrasound diagnostic instrument and system using wireless video transmission
US6682473B1 (en) * 2000-04-14 2004-01-27 Solace Therapeutics, Inc. Devices and methods for attenuation of pressure waves in the body
US6409665B1 (en) 2000-06-01 2002-06-25 Corey D. Scott Apparatus for applying impedence matching fluid for ultrasonic imaging
EP1162476A1 (en) * 2000-06-06 2001-12-12 Kretztechnik Aktiengesellschaft Method for examining objects with ultrasound
KR100350026B1 (en) * 2000-06-17 2002-08-24 주식회사 메디슨 Ultrasound imaging method and apparatus based on pulse compression technique using a spread spectrum signal
US6569097B1 (en) * 2000-07-21 2003-05-27 Diagnostics Ultrasound Corporation System for remote evaluation of ultrasound information obtained by a programmed application-specific data collection device
US6650927B1 (en) 2000-08-18 2003-11-18 Biosense, Inc. Rendering of diagnostic imaging data on a three-dimensional map
US6375616B1 (en) * 2000-11-10 2002-04-23 Biomedicom Ltd. Automatic fetal weight determination
US6822374B1 (en) 2000-11-15 2004-11-23 General Electric Company Multilayer piezoelectric structure with uniform electric field
IL146597A0 (en) 2001-11-20 2002-08-14 Gordon Goren Method and system for creating meaningful summaries from interrelated sets of information
US6643533B2 (en) 2000-11-28 2003-11-04 Ge Medical Systems Global Technology Company, Llc Method and apparatus for displaying images of tubular structures
US6491636B2 (en) * 2000-12-07 2002-12-10 Koninklijke Philips Electronics N.V. Automated border detection in ultrasonic diagnostic images
US6540679B2 (en) 2000-12-28 2003-04-01 Guided Therapy Systems, Inc. Visual imaging system for ultrasonic probe
US6868594B2 (en) * 2001-01-05 2005-03-22 Koninklijke Philips Electronics, N.V. Method for making a transducer
JP4614548B2 (en) 2001-01-31 2011-01-19 パナソニック株式会社 Ultrasonic diagnostic equipment
US6936009B2 (en) 2001-02-27 2005-08-30 General Electric Company Matching layer having gradient in impedance for ultrasound transducers
US6939301B2 (en) 2001-03-16 2005-09-06 Yaakov Abdelhak Automatic volume measurements: an application for 3D ultrasound
US6626836B2 (en) 2001-04-04 2003-09-30 Siemens Medical Solutions Usa, Inc. Adaptive signal processing scheme for contrast agent imaging
US6468218B1 (en) 2001-08-31 2002-10-22 Siemens Medical Systems, Inc. 3-D ultrasound imaging system and method
US7158692B2 (en) 2001-10-15 2007-01-02 Insightful Corporation System and method for mining quantitive information from medical images
US6911912B2 (en) 2001-11-08 2005-06-28 The Procter & Gamble Company Method of urinary continence training based on an objective measurement of the bladder
US7042386B2 (en) * 2001-12-11 2006-05-09 Essex Corporation Sub-aperture sidelobe and alias mitigation techniques
US6544179B1 (en) * 2001-12-14 2003-04-08 Koninklijke Philips Electronics, Nv Ultrasound imaging system and method having automatically selected transmit focal positions
KR100406098B1 (en) * 2001-12-26 2003-11-14 주식회사 메디슨 Ultrasound imaging system and method based on simultaneous multiple transmit-focusing using the weighted orthogonal chirp signals
US6674687B2 (en) 2002-01-25 2004-01-06 Navcom Technology, Inc. System and method for navigation using two-way ultrasonic positioning
US7141020B2 (en) 2002-02-20 2006-11-28 Koninklijke Philips Electronics N.V. Portable 3D ultrasound system
US6903813B2 (en) 2002-02-21 2005-06-07 Jjl Technologies Llc Miniaturized system and method for measuring optical characteristics
ATE447889T1 (en) 2002-03-15 2009-11-15 Bjorn A J Angelsen ULTRASONIC IMAGE OF OBJECTS IN MULTIPLE SCANNING PLANES
US7128711B2 (en) 2002-03-25 2006-10-31 Insightec, Ltd. Positioning systems and methods for guided ultrasound therapy systems
US6878115B2 (en) * 2002-03-28 2005-04-12 Ultrasound Detection Systems, Llc Three-dimensional ultrasound computed tomography imaging system
US6705993B2 (en) * 2002-05-10 2004-03-16 Regents Of The University Of Minnesota Ultrasound imaging system and method using non-linear post-beamforming filter
US6788620B2 (en) 2002-05-15 2004-09-07 Matsushita Electric Ind Co Ltd Acoustic matching member, ultrasound transducer, ultrasonic flowmeter and method for manufacturing the same
US6676605B2 (en) * 2002-06-07 2004-01-13 Diagnostic Ultrasound Bladder wall thickness measurement system and methods
US7520857B2 (en) * 2002-06-07 2009-04-21 Verathon Inc. 3D ultrasound-based instrument for non-invasive measurement of amniotic fluid volume
US6884217B2 (en) * 2003-06-27 2005-04-26 Diagnostic Ultrasound Corporation System for aiming ultrasonic bladder instruments
US7819806B2 (en) * 2002-06-07 2010-10-26 Verathon Inc. System and method to identify and measure organ wall boundaries
US7744534B2 (en) 2002-06-07 2010-06-29 Verathon Inc. 3D ultrasound-based instrument for non-invasive measurement of amniotic fluid volume
GB2391625A (en) * 2002-08-09 2004-02-11 Diagnostic Ultrasound Europ B Instantaneous ultrasonic echo measurement of bladder urine volume with a limited number of ultrasound beams
US20060025689A1 (en) * 2002-06-07 2006-02-02 Vikram Chalana System and method to measure cardiac ejection fraction
US8221322B2 (en) 2002-06-07 2012-07-17 Verathon Inc. Systems and methods to improve clarity in ultrasound images
US20080262356A1 (en) 2002-06-07 2008-10-23 Vikram Chalana Systems and methods for ultrasound imaging using an inertial reference unit
US8221321B2 (en) * 2002-06-07 2012-07-17 Verathon Inc. Systems and methods for quantification and classification of fluids in human cavities in ultrasound images
US20090112089A1 (en) * 2007-10-27 2009-04-30 Bill Barnard System and method for measuring bladder wall thickness and presenting a bladder virtual image
US7041059B2 (en) * 2002-08-02 2006-05-09 Diagnostic Ultrasound Corporation 3D ultrasound-based instrument for non-invasive measurement of amniotic fluid volume
US20070276247A1 (en) 2002-06-07 2007-11-29 Vikram Chalana Systems and methods for ultrasound imaging using an inertial reference unit
US7727150B2 (en) * 2002-06-07 2010-06-01 Verathon Inc. Systems and methods for determining organ wall mass by three-dimensional ultrasound
US20040127797A1 (en) 2002-06-07 2004-07-01 Bill Barnard System and method for measuring bladder wall thickness and presenting a bladder virtual image
US7450746B2 (en) 2002-06-07 2008-11-11 Verathon Inc. System and method for cardiac imaging
US20090062644A1 (en) * 2002-06-07 2009-03-05 Mcmorrow Gerald System and method for ultrasound harmonic imaging
US7611466B2 (en) 2002-06-07 2009-11-03 Verathon Inc. Ultrasound system and method for measuring bladder wall thickness and mass
US7004904B2 (en) * 2002-08-02 2006-02-28 Diagnostic Ultrasound Corporation Image enhancement and segmentation of structures in 3D ultrasound images for volume measurements
US6780152B2 (en) * 2002-06-26 2004-08-24 Acuson Corporation Method and apparatus for ultrasound imaging of the heart
AU2003282690B2 (en) 2002-10-07 2008-12-18 Best Medical International, Inc. Method and apparatus for target position verification
WO2004034694A2 (en) 2002-10-10 2004-04-22 Visualsonics Inc. High frequency high frame-rate ultrasound imaging system
US6825838B2 (en) 2002-10-11 2004-11-30 Sonocine, Inc. 3D modeling system
US6695780B1 (en) * 2002-10-17 2004-02-24 Gerard Georges Nahum Methods, systems, and computer program products for estimating fetal weight at birth and risk of macrosomia
US6628743B1 (en) 2002-11-26 2003-09-30 Ge Medical Systems Global Technology Company, Llc Method and apparatus for acquiring and analyzing cardiac data from a patient
US20040106869A1 (en) 2002-11-29 2004-06-03 Ron-Tech Medical Ltd. Ultrasound tracking device, system and method for intrabody guiding procedures
US6831394B2 (en) 2002-12-11 2004-12-14 General Electric Company Backing material for micromachined ultrasonic transducer devices
US6954406B2 (en) 2003-03-04 2005-10-11 Jones Joie Pierce Acoustical source and transducer having, and method for, optimally matched acoustical impedance
US7489299B2 (en) 2003-10-23 2009-02-10 Hillcrest Laboratories, Inc. User interface devices and methods employing accelerometers
US7090640B2 (en) * 2003-11-12 2006-08-15 Q-Vision System and method for automatic determination of a region of interest within an image
US20050135707A1 (en) 2003-12-18 2005-06-23 Turek Matthew W. Method and apparatus for registration of lung image data
EP1715897B2 (en) * 2004-01-20 2013-10-30 Sunnybrook and Women's College Health Sciences Centre High frequency ultrasound imaging using contrast agents
US20050193820A1 (en) 2004-03-04 2005-09-08 Siemens Medical Solutions Usa, Inc. Integrated sensor and motion sensing for ultrasound and other devices
US7301528B2 (en) 2004-03-23 2007-11-27 Fujitsu Limited Distinguishing tilt and translation motion components in handheld devices
US8900149B2 (en) 2004-04-02 2014-12-02 Teratech Corporation Wall motion analyzer
TWI376520B (en) 2004-04-30 2012-11-11 Hillcrest Lab Inc Free space pointing devices and methods
US7846103B2 (en) * 2004-09-17 2010-12-07 Medical Equipment Diversified Services, Inc. Probe guide for use with medical imaging systems
US7382907B2 (en) 2004-11-22 2008-06-03 Carestream Health, Inc. Segmenting occluded anatomical structures in medical images
US8167803B2 (en) 2007-05-16 2012-05-01 Verathon Inc. System and method for bladder detection using harmonic imaging
US20090105585A1 (en) * 2007-05-16 2009-04-23 Yanwei Wang System and method for ultrasonic harmonic imaging
US20090088660A1 (en) * 2007-08-29 2009-04-02 Mcmorrow Gerald System and methods for nerve response mapping

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040076317A1 (en) * 1998-07-23 2004-04-22 David Roberts Method and apparatus for the non-invasive imaging of anatomic tissue structures
US20080139938A1 (en) * 2002-06-07 2008-06-12 Fuxing Yang System and method to identify and measure organ wall boundaries
US20040054280A1 (en) * 2002-09-18 2004-03-18 Mcmorrow Gerald J. Three-dimensional system for abdominal aortic aneurysm evaluation

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See also references of EP2323559A4 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015091368A1 (en) * 2013-12-20 2015-06-25 Koninklijke Philips N.V. Ultrasound imaging assembly and method for displaying ultrasound images
US10918359B2 (en) 2013-12-20 2021-02-16 Koninklijke Philips N.V. Ultrasound imaging assembly and method for displaying ultrasound images
RU2689769C1 (en) * 2018-01-31 2019-05-28 Государственное Бюджетное Учреждение "Санкт-Петербургский Научно-Исследовательский Институт Скорой Помощи Им. И.И. Джанелидзе" Method for measuring length of proximal cervical abdominal aortic aneurysm
EP3626160A1 (en) 2018-09-24 2020-03-25 Dianova A/S A handheld imaging element with a movement sensor
WO2020064737A1 (en) 2018-09-24 2020-04-02 Dianova A/S A handheld imaging element with a movement sensor
US20210113194A1 (en) * 2019-10-17 2021-04-22 Verathon Inc. Systems and methods for ultrasound scanning
CN114554968A (en) * 2019-10-17 2022-05-27 韦拉索恩股份有限公司 System and method for ultrasound scanning
US11911220B2 (en) * 2019-10-17 2024-02-27 Verathon Inc. Systems and methods for ultrasound scanning
EP3836154A1 (en) 2019-12-13 2021-06-16 Morten Bo Søndergaard Svendsen A method and an apparatus for imaging curved surfaces
EP3839883A1 (en) 2019-12-18 2021-06-23 Morten Bo Søndergaard Svendsen A method and a system for obtaining image data of an object

Also Published As

Publication number Publication date
JP2011530352A (en) 2011-12-22
US20100036242A1 (en) 2010-02-11
CA2957778A1 (en) 2010-02-11
EP2323559A4 (en) 2016-09-21
CA2732997A1 (en) 2010-02-11
JP5658151B2 (en) 2015-01-21
CA2957778C (en) 2022-12-13
US8133181B2 (en) 2012-03-13
CA2732997C (en) 2017-03-14
EP2323559A1 (en) 2011-05-25

Similar Documents

Publication Publication Date Title
CA2732997C (en) Device, system, and method to measure abdominal aortic aneurysm diameter
CA2688778C (en) System and method for ultrasonic harmonic imaging
EP2934328B1 (en) Anatomically intelligent echocardiography for point-of-care
US7727150B2 (en) Systems and methods for determining organ wall mass by three-dimensional ultrasound
EP3080778B1 (en) Imaging view steering using model-based segmentation
US20090062644A1 (en) System and method for ultrasound harmonic imaging
CN105392428B (en) System and method for mapping the measurement of ultrasonic shear wave elastogram
EP1819279B1 (en) System and method for determining organ wall mass by three-dimensional ultrasound
US20130096430A1 (en) Ultrasonic diagnostic apparatus and ultrasonic scanning method
EP3554380B1 (en) Target probe placement for lung ultrasound
EP1679038A2 (en) Ultrasound diagnostic system and method of forming arbitrary m-mode images
US20210312652A1 (en) Automated prostate analysis system
JP6730980B2 (en) Acoustic streaming for detecting and identifying liquid pools
WO2006062867A2 (en) System and method for determining organ wall mass by three-dimensional ultrasound
Brandt Evaluation of new ultrasound techniques for clinical imaging in selected liver and vascular applications
EP2097009A2 (en) System and method for ultrasound harmonic imaging

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 09805626

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2732997

Country of ref document: CA

ENP Entry into the national phase

Ref document number: 2011522292

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 2009805626

Country of ref document: EP