CROSS REFERENCE TO RELATED APPLICATIONSThis application claims priority to, and incorporates by reference in their entireties, U.S. provisional patent application Ser. No. 61/087,152 filed Aug. 7, 2008 and U.S. provisional patent application Ser. No. 61/094,003 filed Sep. 3, 2008.
This application is also a continuation-in-part of, claims priority to, and incorporates by reference in its entirety U.S. patent application Ser. No. 12/121,721 filed May 15, 2008, which is a continuation-in-part of U.S. patent application Ser. No. 11/968,027 filed Dec. 31, 2007, U.S. patent application Ser. No. 11/926,522 filed Oct. 27, 2007, U.S. patent application Ser. No. 11/925,887 filed Oct. 27, 2007, U.S. patent application Ser. No. 11/925,896 filed Oct. 27, 2007, U.S. patent application Ser. No. 11/925,900 filed Oct. 27, 2007, U.S. patent application Ser. No. 11/925,850 filed Oct. 27, 2007, U.S. patent application Ser. No. 11/925,843 filed Oct. 27, 2007, U.S. patent application Ser. No. 11/925,654 filed Oct. 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.
All of the foregoing named applications are incorporated by reference in their entireties as if fully set forth herein.
COPYRIGHT NOTICEThis 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 or records, but otherwise reserves all copyright rights whatsoever.
BACKGROUND OF THE INVENTIONThe 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 inFIG. 1-1 (see U.S. Pat. No. 6,905,468).FIG. 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 overage 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)
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.
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.
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 al., 2005; Byun, et al., 2003)
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%.
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.
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 INVENTIONAn 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 DRAWINGSExamples of particular embodiments of the present invention are described in detail below with reference to the following drawings:
FIG. 1-1 shows a normal abdominal aorta;
FIG. 1-2 shows a typical abdominal aortic aneurysm;
FIGS. 2-1 thru2-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;
FIG. 3-1 is a side elevation view of anultrasound transceiver10A that includes an inertial reference unit
FIG. 3-2 illustrates a side and partial isometric that schematically depicts an ultrasound scanner employing C-mode and B-mode ultrasound modalities;
FIG. 3-3 illustrates a partial isometric and schematic view of an ultrasound scanner system;
FIG. 4 is a schematic illustration of a server-accessed local area network and/or internet in communication with a plurality of ultrasound imaging systems;
FIG. 5 shows exemplary probe locations along the midline of a patient's torso.
FIG. 6-1 is a B-mode image of an AAA is in the same field as a shadow region;
FIGS. 6-2 and6-3 are histograms of the A-mode lines65 (FIG. 6-2) and65 (FIG. 6-3);
FIGS. 7-1 thru7-12 are twelve B-mode images showing shadows and non-availability;
FIG. 8 schematically illustrates an availability plot in C-mode (view from the top of the scan cone);
FIG. 9-1 and9-2 illustrates an aiming direction scheme; and a decision tree for the arrow feedback from the indicator22 ofFIG. 3-3;
FIG. 10 illustrates a screen shot of an exemplary aiming guide;
FIG. 11 depicts a flowchart of a aortic aneurysm detection algorithm employing fundamental ultrasound energies;
FIGS. 12-1 and12-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;
FIGS. 12-3 is a schematic diagram illustrating an exemplary embodiment of the limited segmentation;
FIG. 13 is a graphic diagram for the conversion from volume to diameter assuming an object shape is sphere; and
FIG. 14 depicts a flowchart of the Find Initial Walls sub-algorithm ofFIG. 11.
DETAILED DESCRIPTION OF THE PARTICULAR EMBODIMENTSEmbodiments 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.
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.
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.
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 (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.
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.
The BVI9600 Transceiver and Principal of Operation:
FIGS.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.
FIG. 2-1 is a side elevation view of anultrasound transceiver10A that includes aninertial reference unit38, according to an embodiment of the invention. Thetransceiver10A includes atransceiver housing35 having an outwardly extendinghandle33 suitably configured to allow a user to manipulate thetransceiver10A relative to a patient. Ultrasound transducers operating within thetransceiver10A can be equipped to collect and ready signals for ultrasound fundamental and/or harmonic frequency analysis.
Thehandle33 includes atrigger37 that allows the user to initiate an ultrasound scan of a selected anatomical portion. Thetransceiver10A also includes atransceiver dome34 that contacts a surface portion of the patient when the selected anatomical portion is scanned (SeeFIG. 5). Thedome34 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. Thetransceiver10A further includes one transducer element, or preferably an array of separately excitable ultrasound transducer elements, positioned within or otherwise adjacent to thehousing35. The transducer elements can be suitably positioned within thehousing35 or otherwise to project ultrasound energy outwardly from thedome34, 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 thehousing35 by a motor. Alternatively, the array can be stationary with respect to thehousing35 so that the selected anatomical region can be scanned by selectively energizing the elements in the array.
In one embodiment of thetransceiver10A, a directional indicator panel or aimingguide panel32 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).
In theBVI 9600system70 described inFIG. 3-3 below, thedirectional indicator panel32 oftransceiver10A has a virtual equivalent in the form of a targetingicon screenshot77B that appears on andisplay76 on theconsole74; both theindicator panel32 and display76 (displaying targetingicon77B) 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).
In particular embodiments, if the AAA structure is centered (as indicated by reference numeral77C), and there are no impediments to the scanlines, e.g., shadows caused by air pockets (see discussion below) from placement of thetransceiver10A or10C 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 thetransceiver10A,10C 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 thetransceiver dome34 is placed. Thedirectional indicator panel32 can also be presented on thedisplay54 ofcomputer52 in imaging subsystems described inFIG. 4 below, and/or presented on thetransceiver display36.
Turning back toFIG. 2.1,transceiver10A can, in one embodiment, include an inertial reference unit that includes an accelerometer and/orgyroscope38 positioned preferably within or adjacent tohousing35. Theaccelerometer38 can be operable to sense an acceleration of thetransceiver10A, preferably relative to a coordinate system, while the gyroscope can be operable to sense an angular velocity of thetransceiver10A 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 andgyroscope concerning transceivers10A,10B,10C of FIGS.2.-1,3-2, and/or3-3 equipped with inertial reference units and the operations thereof are described in co-pending U.S. patent application Ser. No. 11/222,360 filed Sep. 8, 2005, herein incorporated by reference.
Thetransceiver10A shown inFIG. 2-1 includes adisplay36 and/or is capable of being in signal communication with a device with a display,FIGS. 3-3 and4) operable to allow viewing of processed results from an ultrasound scan, and/or to allow an operational interaction between the user and thetransceiver10A. For example, thedisplay36 of thehandheld transceiver10A can be configured to display alphanumeric data that indicates a proper and/or an optimal position of thetransceiver10A relative to the selected anatomical portion.Display36 can be used to view two- or three-dimensional images of the selected anatomical region. Accordingly, thedisplay36 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.
In operation, to scan a selected anatomical portion of a patient, thetransceiver dome34 of thetransceiver10A can be positioned against a surface portion of a patient that is proximal to the anatomical portion to be scanned. See for exampleFIG. 5 for exemplary positioning for AAA detection. The user actuates thetransceiver10A by depressing atrigger37. In response, thetransceiver10A transmits ultrasound signals into the body, and receives corresponding return echo signals that can be at least partially processed by thetransceiver10A to generate an ultrasound image of the selected anatomical portion. In a particular embodiment, thetransceiver10A 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.
In one embodiment, thetransceiver10A 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 thetransceiver10A. The system also includes a processor that can be configured to process reflected ultrasound energy that is received by thetransceiver10A 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 thetransceiver10A, such a keyboard, a pointing device, or other similar devices. In still another particular embodiment, thetransceiver10A can be a self-contained device that includes a microprocessor positioned within thehousing35 and software associated with the microprocessor to operably control thetransceiver10A, and to process the reflected ultrasound energy to generate the ultrasound image. Accordingly, thedisplay36 can be used to display the generated image and/or to view other information associated with the operation of thetransceiver10A. For example, the information can include alphanumeric data that indicates a preferred position of thetransceiver10A prior to performing a series of scans.
In yet another particular embodiment, thetransceiver10A can be operably coupled to a general-purpose computer (seeFIG. 4), such as a laptop or a desktop computer that includes software that at least partially controls the operation of thetransceiver10A, and also includes software to process information transferred from thetransceiver10A, so that an image of the scanned anatomical region can be generated. Thetransceiver10A can also be optionally equipped with electrical contacts to make communication with receivingcradles50 as illustrated inFIG. 4 below. Althoughtransceiver10A ofFIG. 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 thetransceiver10A. 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.
FIG. 2-2 is a graphical representation of a plurality of scan planes that form a three-dimensional (3D) array having a substantially conical shape. Anultrasound scan cone40 formed by a rotational array of two-dimensional scan planes42 projects outwardly from thedome34 of thetransceivers10A. Other transceiver embodiments oftransceiver10A can also be configured to develop ascan cone40 formed by a rotational array of two-dimensional scan planes42. The pluralities ofscan planes40 can be oriented about anaxis43 extending through thetransceivers10A. One or more, or preferably each of the scan planes42 can be positioned about theaxis43, preferably, but not necessarily at a predetermined angular position θ. The scan planes42 can be mutually spaced apart by angles θ1and θ2. Correspondingly, the scan lines within each of the scan planes42 can be spaced apart by angles φ1and φ2. Although the angles θ1and θ2are depicted as approximately equal, it is understood that the angles θ1and θ2can have different values. Similarly, although the angles φ1and φ2are shown as approximately equal, the angles φ1and φ2can 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 thetransceiver10A.
FIG. 2-3 is a graphical representation of ascan plane42. Thescan plane42 includes theperipheral scan lines44 and46, and aninternal scan line48 having a length r that extends outwardly from thetransceiver10A. Thus, a selected point along theperipheral scan lines44 and46 and theinternal scan line48 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-sevenscan lines48 that extend outwardly from thedome34, although any number of scan lines is possible.
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 oflines48 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, thetransceiver10A can be configured to generate approximately seventy-seven scan lines between the first limitingscan line44 and a second limitingscan line46. 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 lines48 (FIG. 2-2) can be uniform or non-uniform. For example, and in another particular embodiment, the angular separation φ1and φ2(as shown inFIG. 2-3) can be about 1.5°. Alternately, and in another particular embodiment, the angular separation φ1and φ2can 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.30, 20.2°, 8.0°, 8.0°, 8.0°, 4.3°, 7.8°, and so on.
FIG. 2-4 is a graphical representation of a plurality of scan lines emanating from the hand-heldultrasound transceiver10A forming asingle scan plane42 extending through a cross-section of a section of theabdominal aorta51. The number and location of the internal scan lines emanating from thetransceiver10A within a givenscan plane42 can thus be distributed at different positional coordinates about the axis line11 as can be required to sufficiently visualize structures or images within thescan plane42. As shown, four portions of off-centeredshadow regions49 are exhibited as irregular regions and a cross section of theabdominal aorta51 region of interest (ROI) is depicted as in the center. Three portions can be viewable within thescan plane42 in totality, and one can be truncated by theperipheral scan line44.
FIG. 3-1 depicts a partial schematic and partial isometric and side view oftransceiver10B, 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 ascan cone45. Each of the scan lines has a length r that projects outwardly from thetransceiver10B of FIGS.2-1-2-4. As illustrated, thetransceiver10B emits 3D-distributed scan lines within the scan cone30 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 thescan cone45. Theultrasound scan cone45 extends outwardly from thedome34 of thetransceiver10B centered about anaxis line43. The 3D-distributed scan lines of thescan cone43 include a plurality of internal and peripheral scan lines that can be distributed within a volume defined by a perimeter of thescan cone43. Accordingly, theperipheral scan lines31A-31F define an outer surface of thescan cone45, while theinternal scan lines39A-39C can be distributed between the respectiveperipheral scan lines31A-31F.Scan line39B can be generally collinear with theaxis43, and thescan cone45 can be generally and coaxially centered on theaxis line43.
The locations of the internal and peripheral scan lines can be further defined by an angular spacing from thecenter scan line39B and between internal and peripheral scan lines. The angular spacing betweenscan line39B 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 Ø. The angles Φ1, Φ2, and Φ3respectively define the angular spacings fromscan line39B to scanlines39A,39C, and31D. Similarly, angles Ø1, Ø2, and Ø3respectively define the angular spacings betweenscan line31B and31C,31C and39A, and31D and31E.
With continued reference toFIG. 3-1, the plurality ofperipheral scan lines31A-E and the plurality ofinternal scan lines39A-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 thescan cone45. Thus, a given point within thescan cone45 can be identified by the coordinates r, Φ, and Ø whose values generally vary. The number and location of the internal scan lines emanating from thetransceiver10B can thus be distributed within thescan cone45 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 thetransceiver10B 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° betweenscan line39B and31A, and +60° betweenscan line39B and31B. Thus angle Φ in this example presents a total arc of approximately 120°.
In one embodiment, thetransceiver10B can be configured to generate a plurality of 3D-distributed scan lines within thescan cone45 having a length r, in one embodiment for example, of approximately 20 to 40 centimeters (cm).
FIG. 3-2 illustrates atransceiver10C configured with a transducer designed to provide a fan-like (e.g., having planar sides)conic scan cone47 utilizing C-mode and B-mode ultrasound modalities. Thetransceiver10C projects a series of B-mode scan planes52A and52B that oscillate like a pendulum between extremes within thescan cone47. The B-mode scan planes52A,52B may be derived from a plurality of scan lines similar to scanlines44,46, and48 ofFIGS. 2-3 and2-4. The pendulum oscillatingscan planes52A and52B can be arranged substantially at right angles to each other as depicted at axis crossing57. Theoscillating scan planes52A and/or52B can define a series of C-scan planes55 that vary in depth location from thetransceiver dome34. The C-scan planes55 move from the transducer vanishing point, and the B-scan planes angularly radiate from the transducer vanishing point. Fortransceiver10C users, a portion of the abdominal aorta taken as a C-mode shape is displayed on the transceiver display16. The C-scan geometry shown as scan planes55 present a substantially square-like ultrasound area within thescan cone47. The C-Scan image information contained within scan planes55 presents a cross-section view of the ultrasound at a particular depth probed by thetransceiver10C. 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.
FIG. 3-3 illustrates a partial isometric and schematic view of an ultrasoundAAA detection system70 utilizing atransceiver probe10D andconsole74 combination. TheAAA detection system70 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. Theultrasound transceiver10D may be configured to send out and receive ultrasound RF signals. The received RF may be transmitted to console74. 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.
Thetransceiver10D can have atransceiver display36,housing35 anddome34 design similar totransceivers10A and10B, and may be in signal communication to console74 viasignal cable89. Theconsole74 can be pivoted fromconsole base72. Theconsole74 includes adisplay76, detection andoperation function panel78, andselect panel80. 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.
As illustrated inFIG. 3-3,display76 provides ascreenshot77C with a targetingicon79A including cross hairs centered in a cross sectional depiction of an abdominal aorta. Other screen shots can appear in thedisplay76 depending on which function key is pressed in thefunction panel78. A targetingicon screenshot77B with a plurality of directional arrows can appear and flash to guide the user to move thetransceiver10C to center the abdominal aorta or AAA. The targetingicon screenshot77B can also appear on thedisplay36 of thetransceiver10D. The targetingicon screenshot77B similarly guides the user to place thetransceiver10D to center the abdominal aorta or AAA or other region of interest as thedirectional indicator panel32 did intransceiver10A ofFIG. 2-1 above.
As illustrated inFIG. 3-3, thetransceiver10D can be in wireless communication viawireless signal82 with thewireless hub84; the output from thetransceiver10D can be delivered to awireless hub84 viawireless signal port86. Thewireless hub84 can also chargebatteries88 for loading into the battery compartment (not shown) ofconsole74. In one embodiment, all the calculations can be performed in theimaging console74. The BVI9600-embodiment system70 does not require a computer or network to complete the analysis and imaging processing. In other embodiments, thesystem70 can utilize thewireless hub84 as a gateway to transmittransceiver10D acquired imaging information in local and Internet systems similar to that described inFIG. 4 below.
FIG. 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, anultrasound system90 can include one or morepersonal computer devices52 that can be coupled to aserver96 by acommunications system95. Thedevices52 can be, in turn, coupled to one ormore ultrasound transceivers10A and/or10B and/or10C, for examples the ultrasound sub-systems94A-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 withinscan cone40 or45 presented ondisplay54. Theserver96 can be operable to provide additional processing of ultrasound information, or it can be coupled to still other servers (not shown inFIG. 4) and devices.Transceivers10A/B/C can be in wireless communication withcomputer52 in sub-system94A, in wired signal communication in sub-system94B, in wireless communication withcomputer52 via receivingcradle50 in sub-system94C, or in wired communication withcomputer52 via receivingcradle50 in sub-system94D.
AnInternet system114 can be coupled or otherwise in communication with the ultrasound sub-systems94A-94D.
Scanning and Placement of Transceiver and Determination of “Availability”:
During the field evaluation of AAA scan-mode, a2D 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.
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.
Referring now toFIG. 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 theprobe head34.FIG. 5 shows four such exemplary scan probe locations along the midline of a patient's torso. Firm 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 button37 of thetransceiver10A, for example, may be then depressed. An aim/scan screen appears, for example indisplay36 of the transceiver (FIG. 2.1) and/ordisplay76 onconsol74 and the instrument performs a 3D scan of the patient's region of interest (ROI).
When the scan is done, the aiming result and B-mode image may be displayed, for example as shown in an exemplary screen shotFIG. 10.
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 (SeeFIG. 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 inFIG. 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, twoA-mode lines60,65. InA-mode line60 there may be an AAA betweendepth 85 and 120 mm, thereforeA-mode line60 has useful information such as tissue, AAA, and acoustic enhancement. InA-mode line65 however there may be no information due to the acoustic block near the depth of 25 mm.FIGS. 6-2 and6-3 are representative histograms of the A-mode lines60 (FIG. 6-2) and65 (FIG. 6-3) wherein, a threshold line set at 40 (67), in the exemplary histogram, is shown. The threshold line67 may be used to determine whether the A-mode has available information or not.
FIGS.7-1-7-12 and8 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.
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 FIGS.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. InFIG. 7-10 for example, colored markings/dots underneath the B-mode image represent approximately 40% unavailable scanlines, compared to >90% unavailable scanlines inFIG. 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 inFIG. 8.FIG. 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) depictsavailable regions85 as solid regions andnon-available regions87 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 (1):
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.
InFIG. 9-1, the C-mode display81 can be partitioned into concentric circles, in this example, aninner circle93 and anouter circle97, to directionally position the transceiver to an available area. Each quarter of the circle can be defined as afan98, 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.
FIG. 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 theinner circle93. If theinner circle93 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. Re-positioning 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 solidlight indicator121,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 aggressivelight indicator124,126, and the user may be guided in the recommended direction and orientation to achieve an optimized availability.
FIG. 10 is anexemplary screen view100 of the aiming guide of the BVI9600. In this example, a B-mode image102 and a C-mode image103 is shown. Due to the air block in the 3D data, shown in the C-mode display103 as non-available (by hatched lines)87, the user is directed to move in a NW direction, i.e. towards the upper left hadquadrant110 where the greatest availability and/or “available” scanlines may be detected to collect a better data set.
Turning again toFIG. 9-1 and9-2, theprogram105 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:
Condition 1: Current Position=>75% available atBLOCK110=good scan atBLOCK112=no movement is indicated=start AAA detection algorithm (SeeFIG. 11).
Condition 2: Any inner fan=>75% atBLOCK118=re-position transceiver to aim towards greatest availability=Solid arrow direction indication at BLOCK121 (Seesolid arrow122 indicated in direction)=back toCondition 1 atBLOCK110.
Condition 3: All Inner fans=≦25% available within inner fans=query inner+outer fans=flashing arrow direction indication at BLOCK124 (See flashingarrow126 indicated in direction of maximal inner+outer fan availability)=back toCondition 1 atBLOCK110.
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.
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 inFIG. 4.
Volume Measurement Algorithms:
A. Front and Back Walls. Once an optimal position for AAA measurement is achieved based on the “availability”calculation105 of compared scanlines, as described above and illustrated inFIG. 9-2, anAAA detection algorithm120 may be used to
The fundamentalAAA detection algorithm120 used in theBVI 9600 device begins with process block FindInitial Wall122 using A-mode ultrasound data that incorporates data-smoothing. FindInitial Wall122 looks for the front and back walls of the abdominal aorta illustrated and described inFIG. 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 processblock Find Centroid124. This center abdominal aorta line or centroid may be used as a seed from which process block FixInitial Walls130 utilized, as illustrated and described inFIG. 14 below. FixInitial Walls130 refines the initial wall points, removes any outliers, and fills gaps in the detected wall location regions. TheCentroid118 as shown inFIG. 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. TheFix Walls algorithm130 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. InFIG. 12-1 and12-2 B-mode images of a longitudinal section12-1 and cross section12-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.
FIG. 12-3 is a schematic representation of the application of the limited segmentation wherein a maximum diameter may be fixed from thecentroid118. Two lengths, length (2A)129 the diameter, andlength B135, 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 toFIG. 11, above. The volume calculation atBLOCK134 may be the result of the aggregate calculation of volume resulting from the limited segmentation, in this example, represented bylines136 and138, i.e. those within thescanlines145,147 at points corresponding to the intersection of the front wall and the back wall of the artery being measured. At BLOCK137 the diameter calculation may be derived by application mathematical relationship between volume and diameter as further discussed in relationship toFIG. 13 below.
FIG. 13 graphically represents the relationship between volume and diameter as determined using the below formulae. The volume of sphere volume is:
where, V is a volume and D is a diameter, and
where the diameter can be derived if the volume is know by manipulation of the above relationship:
The mathematical functions can be calculated automatically from a look-up 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 BLOCK134) is determined by the limited segmentation procedure, the diameter may be automatically correlated (at a BLOCK137) with a pre-calculated diameter according to the above relationship.
B. AAA Detection Algorithm, Finding Front and Back Abdominal Aorta Walls and Centroid:
FIG. 14 depicts a flowchart of the Find Initial Walls sub-algorithm ofFIG. 11. FindInitial Walls122 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 withprocess block140. Next, a local gradient at process block142 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. 4-9) below:
The standard central difference formula is given in Equation 4:
dxi=xi+1/2−xi−1/2 (4)
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.
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.
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.
dxi=xi+3+xi+2+xi+1−xi−1−xi−2−xi−3 (6)
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.
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 “right”-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.
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.
Returning toFIG. 14, the results fromlocal gradient142 are subjected to loop limit processing betweenblocks144 and158 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 withloop limit block144 that receives pixel values for each sample in the detection region and subjects the pixel intensity values to determine whether the gradient is minimum atdecision diamond146. If affirmative, then the pixel values may be ascertained whether it's the best front wall-back wall (FW/BW) candidate combination atdecision diamond147. If affirmative, the FW/BW candidate pair may be saved and loop limit processing returns to limitblock144. If negative, atprocess block152, the Front Wall pixel value may be saved and another back wall candidate may be sought with a subsequent return toloop limit block152.
Returning todecision diamond146, if the answer is negative for “Is gradient Minimum?”, sub-algorithm122 continues todecision diamond156 to determine whether the back wall and the gradient is maximum. If affirmative, atprocess block154, a candidate BW/FW pair is established and sub-algorithm re-routes toloop limit block144. If negative, the end of analysis for a particular FW/BW candidate occurs atloop limit block158 either routes back to thelimit loop block144 or exits to findCentroid124.
Formulations relating toFind Centroid124 may be based on coordinate geometries described inequations 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.
x=(38−φ)cos(π−θ)
y=(38−φ)sin(π−θ) (10)
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.
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.