1. Field
Embodiments of the present invention relate to medical equipment and, in particular, to ultrasound equipment.
2. Discussion of Related Art
There is a growing need for a home-based imaging system that would allow clinicians to have access to patients and be able to make diagnostic decisions remotely. Conventional imaging modalities (e.g., X-rays, computed tomography, magnetic resonance and nuclear medicine) are not portable, and they are more suitable in centralized locations, e.g., hospitals and clinics, due to their size, cost, and training required to operate them. On the other hand, ultrasound imaging is safe, non-invasive, portable, inexpensive, relatively easy to use, and capable of real-time imaging. However, current general purpose ultrasound machines are not appropriate for being used at home because they are still bulky, heavy and expensive. In addition, they need a trained specialist (e.g., sonographer) familiar with their operation and the anatomy of internal organs to scan human body and collect the images and other information for radiologists or other physicians to make diagnostic and therapeutic decisions.
An alternative to the general purpose ultrasound machine is an application-specific ultrasound machine, such as portable ultrasound machines and ultrasonic measuring devices. Many current portable ultrasound machines based on the application-specific integrated circuit (ASIC) technology tend to be relatively small, light, inexpensive, and are mainly used in small hospitals and clinics. Although the size and cost of these portable ultrasound machines have been reduced, they are still difficult to operate and expensive for a home-based imaging system. In addition, several compromises have been made, ranging from the imaging modes supported to the image quality provided. Ultrasonic measuring devices, such as bladder scanners and fetal monitors, for example, significantly reduce the cost and size compared to the portable ultrasound machines as well as general purpose ultrasound machines, but they do not provide real-time ultrasound images that are valuable for remote diagnosis, consultation and/or monitoring/screening.
In addition, since these application-specific ultrasound machines are designed based on fixed-function and hardwired design approaches such as application specific integrated circuit (or ASIC) to reduce the cost and size, they also suffer from limited flexibility, which is one of the key features required for a home-based ultrasound imaging system.
Currently, most medical ultrasound examinations are done in hospitals and clinics using general purpose or portable ultrasound machines by clinicians such as sonographers, radiologists, and physicians, for example, and are interpreted by radiologists or specially-trained physicians. Carrying out a traditional ultrasound scan by an unskilled individual at home is not currently allowed because it may lead to missing pathologies and misdiagnosis.
Another drawback of conventional ultrasound machines concerns beamforming of the received reflected signal. The role of a receive beamformer in a medical ultrasound system is to condition receive signals in order to form high-quality images. In the receive beamformer, a signal received from each individual element in an ultrasound transducer is delayed and then combined together into a single coherent signal.
Digital beamforming has been widely used in diagnostic medical ultrasound systems because it can reduce time delay errors and provide lower side lobes and better image resolution, compared with analog beamforming. In digital beamforming, receive signals are quantized using analog to digital converters (ADCs), delayed using digital circuits and then summed together.
To be able to specify a time delay with accuracy, a high speed ADC is required. From previous studies, a sampling frequency of the ADC, fs, for accurate digital beamforming is known to be
f
s≧16f0 (1)
where f0 is the center frequency of transmitted ultrasound signals. If the center frequency is 5 MHz, the ADC sampling frequency must be higher than 80 MHz, which is still very challenging to support in modern ultrasound machines even with current very large scale integration (VLSI) technology because of the number of ADCs required.
Therefore, in most medical ultrasound machines, receive signals are quantized by ADCs with the sampling frequency of 4f0 (e.g., 20 MHz) and then interpolated to simulate a sampling frequency of 16 f0. Depending on the interpolation technique used, the digital beamformers can be classified into an interpolator-based beamformer and a phase rotator-based beamformer.
In the interpolator-based beamformer, an interpolator is placed on each channel to interpolate the receive signal. While a coarse time delay is achieved by controlling the first-in first-out (FIFO) memory, a fine time delay is obtained by changing interpolator's coefficients. The interpolator-based beamforming method can achieve an accurate time delay and high contrast resolution. However, it requires a finite impulse response (FIR) filter on each channel for interpolating the receive signal. Furthermore, clock frequency as high as 16f0 might be needed during interpolation.
Alternatively, the fine time delay is obtained using a phase rotator with the assumption that the receive signal is a narrow-band signal. In the phase rotator-based beamformer, the time delay is converted into a phase value using φl=2πf0Δτl where φl and Δτl are the phase value and the time delay for the lth point in receive beamforming, respectively. In this method, the receive signals are quantized similarly as in the interpolator-based beamforming method. Then, complex baseband signals are obtained by demodulating the quantized signals. Alternatively, complex baseband signals can be derived from quantizing the demodulated receive signals. The obtained complex baseband signals are first delayed by the coarse time delay. Then, the phase of the delayed complex baseband signal is rotated by the phase value in order to compensate the phase distortion introduced by performing beamforming on the complex baseband signal.
Since the phase rotator-based beamforming performs beamforming on the complex baseband signal, it does not require interpolation filters and alleviates the high data transfer rate requirement. Thus, the phase rotator-based beamforming would be suitable for a low-cost beamforming technique compared to the interpolator-based beamforming. In addition, further reduction in the hardware complexity can be achieved by dividing the conventional phase-rotator based beamforming into two stages (i.e., coarse time delay adjustment and phase compensation) since they have different hardware requirements (i.e., high data transfer rate and computation, respectively).
Still another drawback of conventional ultrasound systems relates to the signal-to-noise ratio (SNR) and resolution of the system. For example, as described above in medical ultrasound imaging systems, electrical signals are applied to an ultrasonic transducer to generate ultrasound waves, which are then transmitted into the human body for imaging. To obtain high signal to noise ratios (SNR) and good resolution, the electrical signals typically have high peak power and short time duration. Although the time gain compensation (TGC) is applied to the receive signals, it may be difficult to obtain an appropriate SNR for an object deep inside the body due to high attenuation in soft tissues. By increasing the peak power of transmit signals, higher SNRs may be obtained, but it is not desirable because high peak power could potentially damage the ultrasonic transducer and the soft tissues underneath. Therefore, it is necessary to improve the SNRs of medical ultrasound systems without increasing the peak power of transmit signals.
Coded excitation techniques are capable of improving the SNR by increasing the average power of transmit signals instead of the peak power. In coded excitation, an elongated signal, which is encoded with high time-bandwidth (TB) product codes for increasing the average power and preserving the spatial resolution, is transmitted and then the reflected signal from the body is decoded into a short signal by pulse compression. The expected improvement in SNR from coded excitation, GSNR, is given by
GSNR=10 log10M (2)
where M represents the relative time duration of the elongated transmit signal with respect to that of the conventional short transmit signal. However, it is practically difficult to achieve the above SNR improvement due to the limited transmit power efficiency (TPE) of the encoded transmit signal, which is defined as the ratio of the transmit power available at the output and input of an ultrasonic transducer. Therefore, when selecting an encoding code, the TPE should be considered with other desirable features, such as imaging resolution and transmitter complexity.
Various TB codes, including Chirp, Golay and Barker, have been extensively examined for coded excitation. Among these codes, the Chirp codes can maximize the TPE because they can be designed to have most of their energy within the frequency bandwidth of the ultrasonic transducer. The Chirp codes are commonly weighted by a window function (e.g., Hanning and Chebyshev) to attain acceptable imaging resolution, i.e., narrow mainlobe width (MLW) and low peak sidelobe level (PSL). However, the weighted Chirp codes need a complex transmitter on each channel to amplify their arbitrary values, i.e., a linear power amplifier.
The Golay codes can provide the narrow MLW and minimal PSL with Golay sequences (i.e., +1 and −1). Although the Golay codes do not need complex power amplifiers, paired firings are needed, leading to a reduction in frame rates. In addition, for coherent summation between the complementary pairs, additional hardware is needed to store the results from pulse compression. Moreover, if there is tissue motion during paired firings, severe artifacts are introduced due to the incoherency between the complementary Golay codes.
On the other hand, the Barker codes can provide the narrow MLW and low PSL (e.g., −22 dB with the length of 13) without the need of paired firings and the complex power amplifier. For further reduction in PSL, a sidelobe suppression filter can follow matched filtering when performing pulse compression. However, the Barker codes suffer from the low TPE due to their wide frequency bandwidth that is not matched to that of the ultrasonic transducer. The low TPE results in lower sensitivity and higher temperature in the ultrasonic transducer due to high dissipated power. Therefore, it is desirable to improve the TPE of the Barker codes when using them as an elongated transmit signal in coded excitation.
As described above, the receive signal is decoded by pulse compression to improve the SNR and spatial resolution, particularly in the axial direction. Two types of pulse compression can be applied. In pre-compression, the receive signal is compressed by a pulse compressor on each channel before receive beamforming. This approach can certainly offer the effective compression of the receive signal. However, it requires multiple pulse compressors, resulting in a high complexity in pulse compression. In post-compression, the receive signals from multiple channels are coherently combined together by the receive beamformer, and then the beamformed signal is decoded by a single pulse compressor. Although the post-compression method can reduce the computational complexity in pulse compression significantly, it introduces artifacts in the images due to distortions in the elongated signals caused by dynamic receive focusing during receive beamforming. Therefore, it is desirable to develop an efficient pulse compression method for the Barker codes to achieve effective pulse compression with an acceptable computational complexity.
A home ultrasound system according to embodiments of the present invention may be capable of adapting to changing clinical needs and/or new applications, assisting non-experts to acquire clinically usable data/images, updating the examination protocols (i.e., scanning, image formation and analysis) from a remote location, supporting new clinical applications and/or adapting to changing clinical needs, supporting an efficient power management for improved portability and longer battery life, and supporting remote diagnosis, consultation and/or monitoring/screening. For some embodiments, the home ultrasound system includes a home ultrasound machine, external computing devices, local storage, central storage, and/or a remotely located diagnosis station. The home ultrasound machine may be used to scan a patient at home and acquire ultrasound data. The acquired ultrasound data may be transferred to the diagnosis station via a communication network. The home ultrasound machine may be located in a clinic, such as a local neighborhood clinic, in a physician's office, and/or in a hospital, such as in a hospital emergency room, for example. The home ultrasound machine also may be located in a vehicle, such as an aid vehicle, for example.
For some embodiments, the home ultrasound machine may include an application-specific scan head, reconfigurable hardware, a programmable processing unit, configuration memory, a power manager, a system controller, a network interface, a user interface, embedded storage, and so on.
The application-specific scan head may be changeable to support different applications. The scan head may include an identification memory that stores information associated with a type for scan head. For some embodiments, the home ultrasound machine may compare the scan head type with a scan head type specified for a particular ultrasound examination and provide an error indication if the scan head type and the scan head type specified for the particular ultrasound examination do not match.
The scan head may transmit an ultrasound signal to a patient and receive a reflected ultrasound signal. The reflected ultrasound signal may be converted to a radio frequency (RF) signal in the scan head. The scan head also may encode the transmit signal with binary phase codes, such as Barker codes, for example. The scan head may use an efficient transmit power (ETP) coding process. The scan head may encode the transmit signal with binary phase codes, such as Barker codes, for example. The ETP-coding process consists of three stages: generation of a principal transmit signal, selection of an encoding kernel, and coding the generated transmit signal with the selected encoding kernel.
For effective pulse compression with a low computational complexity, 2-stage pulse compression is applied to the receive signals. This 2-stage pulse compression consists of a pre-compressor using matched filters and a post-compressor using a single sidelobe suppression filter. In each pre-compressor, a matched filter is used to decode the receive signals coded with the Barker codes to minimize the distortion of the receive signal during receive beamforming. The decoded receive signals are combined together during receive beamforming.
A low-cost digital receive beamformer is provided by dividing a phase-rotator based beamforming into two stages (i.e., pre- and post-beamforming). In the pre-beamforming stage where high data transfer rate is needed, the appropriate complex baseband samples with coarse time delays are selected. The phase compensation requiring high computation capability is performed in the post-beamforming stage. In one embodiment, the pre-beamforming stage is implemented on the low-cost reconfigurable circuit(s) while the post-beamforming stage is implemented on the programmable digital signal processors. Thus, the cost reduction is obtained in the phase-rotator based beamforming by utilizing low-cost reconfigurable circuits and digital signal processors and taking advantage of their hardware reusability.
The reconfigurable hardware may perform pre-beamforming in which a coarse time delay adjustment may be applied to the received signal. The reconfigurable hardware also may perform the first stage of compression in which the received signal may be decoded. The reconfigurable hardware may be programmable gate array (PGA), a field programmable gate array (FPGA), a programmable logic device (PLD), or an application specific integrated circuit (ASIC). The programmable processor may be software on a digital signal processor (DSP).
The programmable processor may perform a second stage of compression in which the peak sidelobe level (PSL) of the coherent beam formed following pre- and post-beamforming is filtered. Then, the beamformed signal is filtered with the sidelobe suppression filter to reduce the PSL of the decoded receive signals. The matched filter used in the pre-compressor can be implemented by using only 2's complement adders because the Barker codes are composed of binary sequences. Thus, the matched filter can be placed in each channel without creating a large computational burden. Only a single sidelobe suppression filter, which can be implemented using complex multipliers as well as adders, is needed in the post-compressor. Therefore, the developed coded excitation technique is a cost-effective solution to improve the SNR in the medical ultrasound systems by enhancing the TPE and minimizing the artifacts from dynamic receive focusing while reducing the necessary hardware complexity.
The programmable processor also may perform back-end processing such as B-mode, spectral Doppler, color-flow and 3D processing, for example. To assist non-clinical people in acquiring clinically usable data and/or images, the programmable processing unit also may be utilized in performing assisted guidance and application-specific analysis by taking advantage of the reusability of the programmable hardware.
The reconfigurable hardware (HW) and the programmable processing unit may be capable of adapting to changing clinical needs and/or new applications by downloading the corresponding configuration information either locally (e.g., flash memory) or via the communication network. This configuration information may be stored in the configuration memory and utilized by the home ultrasound machine.
For longer battery life, the home ultrasound machine may provide an efficient power management based on transducer contract analysis. Furthermore, different levels of power saving modes are supported by adjusting the system parameters as well as by changing the display intensity.
To support remote diagnosis, consultation and monitoring/screening, the beamformed RF data, ultrasound images, and measurement values may be transferred to the diagnosis station by a communication network. The beamformed RF data may be processed by the diagnosis station utilizing signal and image processor(s) or a personal computer (PC) to generate ultrasound images and/or application-specific measurement values in the hospital/clinic.
Alternatively, the beamformed RF data may be processed by external computing devices such as a personal computer (PC), for example, using software and/or hardware at home. The acquired RF data, images and measurement values from the external computing devices may be stored locally and then directly transferred to the hospital/clinic.
In the hospital/clinic, the diagnosis station may be used to review and analyze the transferred ultrasound data and images. Alternatively, the diagnosis station may be used to optimize the parameters for the back-end processing and measurement algorithms to generate better quality images and more accurate measurements. The diagnosis station also may be utilized to generate and transfer the settings/parameters and algorithms for ultrasound data acquisition and/or signal/image processing used in the home ultrasound machine. The ultrasound images and data may be converted to standard formats, such as Digital Imaging and Communications in Medicine (DICOM) format, for example, and transferred into picture archiving and communications system (PACS) for further diagnosis and/or permanent archiving. Using the home ultrasound system according to embodiments of the present invention, ultrasound examinations could be performed at home by non-experts such as patients and their family members, for example, and then acquired ultrasound data including images can be transferred to hospitals for interpretation by radiologists and/or trained clinicians.
In the drawings, like reference numbers generally indicate identical, functionally similar, and/or structurally equivalent elements. The drawing in which an element first appears is indicated by the leftmost digit(s) in the reference number, in which:
The illustrated home ultrasound machine 102 includes an application-specific scan head 120, reconfigurable hardware 122, a programmable processing unit 124, home ultrasound system (HUS) configuration memory 126, a power manager 128, a system controller 130, a network interface 132, a user interface 134, and embedded storage 136 operatively coupled to each other. The illustrated diagnosis station 104 includes patient information (PI) manager 140, a network interface 142, a user interface 144 and a signal and image processor 146 operatively coupled to each other. The illustrated diagnosis station also is coupled to a picture archiving and communications system (PACS) 148.
The home ultrasound machine 102 may be used by non-experts (e.g., nurses, patients or their family members) to scan a patient and acquire ultrasound data. The settings and/or parameters for the ultrasound examination may be downloaded from the diagnosis station 104 located in a hospital/clinic. The home ultrasound machine 102 can be adapted to multiple applications by changing the application-specific scan head 120 and/or downloading different configuration information for the reconfigurable HW 122 and the programmable processing unit 124.
For some embodiments, the first stage of beamforming takes place in the reconfigurable HW 122 and the second stage of beamforming takes place in the programmable processing unit 124. Beamforming is used to improve the signal-to-noise-ratio (SNR) and spatial resolution by coherently summing the ultrasound signals from the scan head 120.
The acquired data from the application-specific scan head 120 may be processed in the reconfigurable HW 122 and the programmable processing unit 124, and then the generated ultrasound images and the application-specific measurement values may be transferred to a diagnosis station 104 in the hospital/clinic for remote diagnosis, consultation and/or monitoring/screening. Alternatively, the acquired data after receive beamforming can be directly transferred to the diagnosis station 104 where optimized parameters may be used for better image quality and more accurate measurements.
Alternatively, the acquired data may also be processed in the external computing devices 108, and then the generated ultrasound images and measurement values may be transferred to the diagnosis station 104. In addition, the acquired data including ultrasound images and measurements may be transferred to the central storage 112 for permanent archiving. The acquired data may be converted to standard formats, such as DICOM, for example, and then connected to the PACS 148.
The transmit information (i.e., focusing time delay and transmit power) for the home ultrasound machine 102 may be stored in the transmit pattern memory 212. This information may be initially obtained from the HUS configuration memory 126 and dynamically updated by the beamforming delay calculator (described below with reference to
The transmit information stored in the transmit pattern memory 212 may be converted to corresponding electrical signals by the pulser 206, and the electrical signals may be used to excite the transducer elements 202. The T/R switch 210 may be used to separate the transmit channel from the receive channel and to protect the receiver 208 circuitry. The multiplexer 204 may be used to multiplex the transducer elements if the total number of transducer elements is different from the number of active transducer elements used for the examination (e.g., linear array and convex array). The ultrasound acoustic signals may be generated by converting the electrical signals into the acoustic waves via transducer elements 202. The reflected acoustic signals are sensed by the transducer elements 202 and converted into radio frequency (RF) electrical signals.
In the receiver 208, the converted RF signals may be amplified in proportion to the depth or the time in order to compensate for signal attenuation (i.e., time gain compensation, TGC) after undergoing low noise amplification (LNA). After LNA and TGC, RF electrical signals may be transferred to the reconfigurable HW 122 for receive beamforming.
In the reconfigurable HW 122, the RF signals from the application-specific scan head 120 may first be digitized by the ADC 302, e.g., at the sampling frequency of 4f0 where f0 is the center frequency of the transducer elements 202. The demodulator 304 may remove the carrier frequency using a demodulation technique, such as quadrature demodulation, for example. In the quadrature demodulation, the quantized RF signal is multiplied with cos(2πf0t) and sin(2πf0t). After low pass filtering, it becomes the baseband signal of complex samples (i.e., I(t)+jQ(t)). Alternatively, the demodulator 304 may remove the carrier frequency by performing a demodulation before digitization.
For some embodiments, in the reconfigurable HW, demodulator 304 is coupled to a pre-compressor 310, which is coupled to a pre-beamformer processing unit 306.
The complex baseband signals then may undergo dynamic receive focusing in order to form high-quality images. In one embodiment, a coarse time delay adjustment may be applied to quantized in-phase and quadrature components of the complex baseband signals.
Based on the geometry of transducer elements 202 and the applied beamforming technique, the programmable processing unit 124 may compute transmit focusing time delays.
a) illustrates an estimation of the transmit pattern memory size and calculation of the transmit focusing time delay for (b) strong focusing and (c) weak focusing As shown in
where Δdtx,max is the maximum delay distance, c is the sound velocity, fs,tx is the transmit sampling frequency in the application-specific scan head 120 and N is the number of active transducer elements. In case of strong focusing (e.g., a single focal point at dztx,f), the time delays for the ith transducer element can be calculated by
where Δdtx,i is the delay distance for the ith transducer element. Similarly, the time delay for jth transducer element is given by
where Δdtx,j is the delay distance for the jth transducer element.
Strong focusing and corresponding receive beamforming techniques may be appropriate for obtaining better image quality and more accurate measurements. However, their frame rates may be limited because it is difficult to generate more than two scan lines with a single firing. On the other hand, weak focusing and following beamforming techniques can increase the frame rates significantly because they distribute acoustic energy in the transmit focal zone instead of focusing on a single point as shown in
where Δdtx,i is the delay distance for the ith transducer element and dztx,i is the distance between the center of transducer elements and the focal point for the ith transducer element within the transmit focal zone. Similarly, the time delay for jth transducer element is given by
where Δdtx,j is the delay distance for the jth transducer element and is the distance between the center of transducer elements 202 and the focal point for the jth transducer element within the transmit focal zone.
where dximg and dzimg are the location of the imaging point in the lateral and axial directions, respectively, and dxrx,ele and dzrx,ele are the location of the receive element. The receive time delay is defined by
The time delay for the post-beamformer is obtained by adding the transmit time delay (i.e., Δτtx(i,j,k,l)) and the receive time delay via
Δτtx,rx(i,j,k,l)=Δτtx(i,j,k,l)+Δτrx(i,j,k,l) (10)
The pre-beamforming delay represented as the number of samples for the kth sub-scan line, Δ{circumflex over (τ)}tx,rx(i,j,k,l), is given by
Δ{circumflex over (τ)}tx,rx(i,j,k,l)=T[Δτtx,rx(i,j,k,l)×fs,rx] (11)
where T[.] is the truncation operator to remove the fractional part and fs,rx is the ADC's sampling frequency. Similarly, the pre-beamforming delay for the k−1th sub-scan line may be obtained for the dual-beam technique.
This time delay represented in Eq. (11) may be computed in real time or may be computed beforehand and stored in memory such as the pre-beamformer delay LUT 308, for example. In one embodiment, the time delays in the pre-beamformer delay LUT 308 may be initially obtained from the HUS configuration memory 126 and dynamically updated in real time by the programmable processor(s) 124. To support various beamforming techniques, multiple pre-beamformer processing units 306 can be integrated into the reconfigurable HW 122 due to the flexibility of the reconfigurable HW 122. The configuration information for the reconfigurable HW 122 may be obtained from the HUS configuration memory 126 when the home ultrasound machine 102 is powered on. This configuration information may be downloaded from the diagnosis station 104 via the communication network 106 or may be updated locally using the external storage device 110 (e.g., flash memory) by the user. The pre-beamformed RF data may be transferred to the programmable processing unit 124 for the fine time delay adjustment and back-end processing.
For some embodiments, the post-compressor 612 is performed after post-beamforming and before signal and image processing.
The programmable processing unit 124 may perform three tasks: (1) computation of transmit and receive focusing time delays, (2) post-beamforming in which the phase compensation is applied, (3) post-compression where mismatched filtering is applied, and (4) signal and image processing for generating ultrasound images and performing assisted guidance and application-specific analysis.
For some embodiments, the signal and image processor 606 may generate the ultrasound images from the beamformed RF data.
In some embodiments, the beamforming delay calculator 608 may compute the transmit and receive focusing time delays for the application-specific scan head 120, the reconfigurable HW 122, and the post-beamformer 604 in the programmable processing unit 124. Transmit focusing time delays represented as Eqs. (4) and (5) and two types of the receive focusing time delays (i.e., the pre and post-beamforming delay) represented as Eqs. (10) and (11) may be dynamically computed in the beamforming delay calculator 608 and then transferred to the transmit pattern memory 212, the pre-beamformer delay LUT 308, and the post-beamformer delay LUT 602.
In the post-beamformer 604, the phase of the pre-beamformed complex baseband signals may be adjusted before summation in order to compensate the phase distortion introduced in phase-oration beamforming. After applying the phase compensation by the post-beamformer delay via phase rotation, the complex baseband data from all channels may be coherently combined together in a summation stage. The coherently summed data may be directly transferred to the external computing device 108, the embedded storage 136, and/or the diagnosis station 104 via the communication network 106 for further processing and display. Alternatively, the coherently summed data may be transferred to the signal and image processor 606 for further processing.
The volumetric data, i.e., 3D data 712, may be reconstructed with the scan-converted B-mode and color/power Doppler data. Alternatively, the 3D data can be obtained directly from the B-mode and color/power Doppler data. The acquired B-mode, spectral/color/power Doppler and/or 3D data may be stored in the embedded storage 136. Alternatively, the acquired B-mode, spectral/color/power and/or 3D data may be transferred to the central storage 112 and the diagnosis station 104 through the network interface 132. Additionally, the acquired B-mode, spectral/color/power Doppler and/or 3D data may be utilized in the guidance assister 708 and the application-specific image analyzer for less trained operators (e.g., nurses, patients or their family members).
The home ultrasound machine 102 will be typically used by non-experts. It may be challenging for them to acquire appropriate ultrasound images for medical purpose without any guidance.
In one embodiment, the transducer contact analysis stage 802 may detect scan lines arising from an improper transducer (i.e., transducer elements 202) contact with the underlying tissue of the patient during ultrasound examination. Those scan lines may be identified by measuring the sum of returned energy along each axial direction where the return energy can be estimated by performing inverse time gain compensation (TGC). Alternatively, the transducer contact analysis may be performed with the beamformed RF data without performing inverse TGC.
The evaluation results based on the transducer contact analysis may be indicated to the user via display and/or voice. These contact analysis could be utilized to avoid the excessive exposure of the ultrasound energy to patients. Additionally, these results could also be used to control the transmitting power of the transducer based on the status of the transducer contact. Receive beamforming and image analysis parameters could also be changed when bad contact is detected so that battery life can be extended. For example, when there is a bad transducer contact, the ultrasound machine could be switched to operate in a low-power consumption mode.
At the same time, the image quality may be quantified based on several image quality metrics in real time and the evaluation results may be indicated to the user and/or recorded as part of the image sequence. New guidance stages may be easily added to the guidance assister 708 due to its programmability and flexibility.
To support multiple clinical applications, application-specific evaluation stages may be integrated into an application-specific image analyzer 710.
For the illustrated obstetrics and gynecology application, amniotic fluid indexes can be measured, by the amniotic fluid index measurement stage 902, for example, via image segmentation based on the intensity, texture connectivity, and other information. The measurement results may be transferred to the diagnosis station 104 in the hospital/clinic via the network interface 132, and the measurement results may be used for initial diagnosis or screening of patients.
The strain measurement stage 906 may compute straining images based on deformation caused by pressure. Not only the reconstructed strain image but also analysis results, such as locations and sizes of less elastic tissues (potentially cancerous), may be transferred to the diagnosis station 104 in the hospital/clinic via the network interface 132. In addition, the bladder volume measurement stage 908 may be used to estimate the bladder volume by measuring the bladder area. The bladder region may be identified by an image segmentation algorithm and/or other information/techniques. These application-specific analyses may be conducted with ultrasound data as well as several parameters used in TGC, log compression, and/or other stages.
The system configuration parameters for various processing units in the home ultrasound machine 102 may be stored in the HUS configuration memory 126.
For some embodiments, the system configuration parameters may be sent out to the application-specific scan head 120, the reconfigurable HW 122, the programmable processing unit 124, and the power manager 128. The HUS configuration memory 126 may include information such as the initial firing sequence for the pulser 206, initial settings for the pre-beamforming LUT 308 and the post-beamforming LUT 602, configuration information for the reconfigurable HW 122, programs for the programmable processing unit 124, power management information for the power manager 128, and/or application-specific scan head identity information for the ID information memory 214.
To update examination protocols and support remote diagnosis, consultation and monitoring/screening, the home ultrasound system 102 may provide improved network interfaces.
Similar to the network interfaces 132 and 142, the user interface 134 may allow information exchanges between the user and the home ultrasound machine 102. In some embodiments, the user interface 134 may be a display for display of the ultrasound images and the current parameters and settings, a keyboard for changing certain parameters or modes in the home ultrasound machine 102, a touch screen for changing certain parameters or modes in the home ultrasound machine 102, sound and/of audio for helping the user to operate the home ultrasound machine 102, for guiding the user's to proper scanning, and for getting user's attention, a communications interface for loading/storing certain parameters in the HUS configuration memory 126 and uploading ultrasound data from the home ultrasound machine 102. Communication interfaces may include standard interfaces (e.g., USB and IEEE 1394).
In some embodiments, the home ultrasound machine provides a power management method based on transducer contact analysis.
In some embodiments the power manager 128 may provide different levels of power modes by changing several system parameters used in the application-specific scan head 120, the reconfigurable HW 122, and the programmable processing unit 124 as well as by adjusting the display intensity. The decision on the power mode to be used may be made by the user based on the trade-off between the desired image quality and the battery time left. Alternatively, this decision can be automatically made by the home ultrasound system 100 based on a predefined setting or by the diagnosis station 104 being operated by the clinician.
In a block 1204, the power manager 128 monitors the voltage across the battery source (i.e., VBAT) and the amount of current being provided by the battery source (i.e., IBAT). Using these two parameters, the power being consumed from the battery source is computed, and the energy left in the battery is estimated.
Based on the current power consumption and the battery energy left, the time left in the battery can be estimated. In a block 1206, the process 1200 determines whether the battery time left is more than the predefined threshold value. If the battery time left is not less than the predefined threshold value, the current power mode is not changed and control of the process 1200 passes to a block 1208 in which the time threshold value may be updated.
On the contrary, if the battery time left is less than the current threshold value, the user may be asked to reduce the current power consumption in the home ultrasound machine 102 by reducing the display intensity (block 1210). If the user chooses to reduce the display intensity in a block 1212, the user can reduce the display intensity by a predefined amount and then control passes to a block 1214.
In block 1214, the user may be asked to reduce the current power consumption in the home ultrasound machine 102 by degrading the image quality. If the user chooses to degrade the image quality in block 1216, the user may select a level of image quality degradation in order to prolong the battery life. Different levels of image quality degradation may be achieved by changing the system parameters used in various functional units (e.g., the application-specific scan head 120, the reconfigurable HW 122, and/or the programmable processing unit 124).
The diagnosis station 104 may be used to support remote diagnosis, consultation and/or monitoring/screening. As
For some embodiments, the PI manager 140 may be used to provide the ultrasound examination protocols. In addition, the PI manager 140 may handle the transferred patient information and ultrasound data from the home ultrasound machine 102 to confirm whether the downloaded parameters are appropriately applied during scanning and processing. After this confirmation stage, the transferred ultrasound data including images and measurement values may be transferred to the signal and image processor 146 for further processing and/or diagnosis.
The network interface 142 may handle the communications between the diagnosis station 104 in hospital/clinic and the home ultrasound machine 102.
As described above in medical ultrasound imaging systems, electrical signals are applied to an ultrasonic transducer to generate ultrasound waves, which are then transmitted into the human body for imaging. To obtain high signal to noise ratios (SNR) and good resolution, the electrical signals typically have high peak power and short time duration. Although the time gain compensation (TGC) is applied to the receive signals, it may be difficult to obtain an appropriate SNR for an object deep inside the body due to high attenuation in soft tissues. By increasing the peak power of transmit signals, higher SNRs may be obtained, but it is not desirable because high peak power could potentially damage the ultrasonic transducer and the soft tissues underneath. Therefore, it is necessary to improve the SNRs of medical ultrasound systems without increasing the peak power of transmit signals.
Coded excitation techniques are capable of improving the SNR by increasing the average power of transmit signals instead of the peak power. In coded excitation, an elongated signal, which is encoded with high time-bandwidth (TB) product codes for increasing the average power and preserving the spatial resolution, is transmitted and then the reflected signal from the body is decoded into a short signal by pulse compression. The expected improvement in SNR from coded excitation, GSNR, is given by
GSNR=10 log10M (12)
where M represents the relative time duration of the elongated transmit signal with respect to that of the conventional short transmit signal. However, it is practically difficult to achieve the above SNR improvement due to the limited transmit power efficiency (TPE) of the encoded transmit signal, which is defined as the ratio of the transmit power available at the output and input of an ultrasonic transducer. Therefore, when selecting an encoding code, the TPE should be considered with other desirable features, such as imaging resolution and transmitter complexity.
Various TB codes, including Chirp, Golay and Barker, have been extensively examined for coded excitation. Among these codes, the Chirp codes can maximize the TPE because they can be designed to have most of their energy within the frequency bandwidth of the ultrasonic transducer. The Chirp codes are commonly weighted by a window function (e.g., Hanning and Chebyshev) to attain acceptable imaging resolution, i.e., narrow mainlobe width (MLW) and low peak sidelobe level (PSL). However, the weighted Chirp codes need a complex transmitter on each channel to amplify their arbitrary values, i.e., a linear power amplifier.
The Golay codes can provide the narrow MLW and minimal PSL with Barker sequences (i.e., +1 and −1). Although the Golay codes do not need complex power amplifiers, paired firings are needed, leading to a reduction in frame rates. In addition, for coherent summation between the complementary pairs, additional hardware is needed to store the results from pulse compression. Moreover, if there is tissue motion during paired firings, severe artifacts are introduced due to the incoherency between the complementary Golay codes.
On the other hand, the Barker codes can provide the narrow MLW and low PSL (e.g., −22 dB with the length of 13) without the need of paired firings and the complex power amplifier. For further reduction in PSL, a sidelobe suppression filter can followed matched filtering when performing pulse compression. However, the Barker codes suffer from the low TPE due to their wide frequency bandwidth that is not matched to that of the ultrasonic transducer. The low TPE results in lower sensitivity and higher temperature in the ultrasonic transducer due to high dissipated power. Therefore, it is desirable to improve the TPE of the Barker codes when using them as an elongated transmit signal in coded excitation.
As described above, the receive signal is decoded by pulse compression to improve the SNR and spatial resolution, particularly in the axial direction. Two types of pulse compression can be applied. In pre-compression, the receive signal is compressed by a pulse compressor on each channel before receive beamforming. This approach can certainly offer the effective compression of the receive signal. However, it requires multiple pulse compressors, resulting in a high complexity in pulse compression. In post-compression, the receive signals from multiple channels are coherently combined together by the receive beamformer, and then the beamformed signal is decoded by a single pulse compressor. Although the post-compression method can reduce the computational complexity in pulse compression significantly, it introduces artifacts in the images due to distortions in the elongated signals caused by dynamic receive focusing during receive beamforming. Therefore, it is desirable to develop an efficient pulse compression method for the Barker codes to achieve effective pulse compression with an acceptable computational complexity.
A method and apparatus for a coded excitation technique using efficient transmit power (ETP) coding and 2-stage pulse compression according to embodiments of the present invention improve the SNR and spatial resolution in the home ultrasound system 100. To improve the transmit power efficiency (TPE), an elongated transmit signal based on the binary phase codes (e.g., Barker) may be encoded by the developed ETP coding where the frequency response of transmit signals is matched to that of the ultrasonic transducer.
In some embodiments, the ETP-coding process may include of three stages: generation of a principal transmit signal, selection of an encoding kernel, and coding the generated transmit signal with the selected encoding kernel. For effective pulse compression with a low computational complexity, 2-stage pulse compression may be applied to the receive signals. This 2-stage pulse compression includes of a pre-compressor using matched filters and a post-compressor using a single sidelobe suppression filter.
In each pre-compressor, the matched filter is used to decode the receive signals coded with the binary phase codes (e.g., Barker) to minimize the distortion of the receive signal during receive beamforming. The decoded receive signals may be combined together during receive beamforming. Then, the beamformed signal may be filtered with the sidelobe suppression filter to reduce the peak sidelobe level (PSL) of the decoded receive signals. The matched filter used in the pre-compressor may be implemented by using only 2's complement adders because Barker codes are composed of binary sequences. Thus, the matched filter may be placed in each channel without creating a large computational burden. Only a single sidelobe suppression filter, which can be implemented using complex multipliers as well as adders, may be used in the post-compressor. The developed coded excitation technique is thus a cost-effective solution to improve the SNR in the medical ultrasound systems by enhancing the TPE and minimizing the artifacts from dynamic receive focusing while reducing the necessary hardware complexity.
To improve the TPE, an elongated transmit signal based on the binary phase codes (e.g., Barker) is encoded by the developed ETP coding where the frequency response of transmit signals is matched to that of an ultrasonic transducer. The ETP-coding process consists of three stages: generation of a principal transmit signal, selection of an encoding kernel, and coding the generated transmit signal with the selected encoding kernel. For effective pulse compression with a low computational complexity, two-stage pulse compression is applied to the receive signals.
The pre-compressor 1406 may use matched filters to decode the receive signals coded with the binary phase codes (e.g., Barker) to minimize the distortion of the receive signal during receive beamforming. The decoded receive signals may be combined together in the beamformer 1408 during receive beamforming. The post-compressor 1410 may filter the beamformed signal with the sidelobe suppression filter to reduce the PSL of the decoded receive signals.
For some embodiments, the matched filter used in the pre-compressor 1406 may be implemented by using only two's complement adders because of the property of the binary phase codes. Thus, the matched filter may be placed in each channel without creating a large computational burden. Only a single sidelobe suppression filter, which can be implemented using complex multipliers as well as adders, may be used in the post-compressor 1410. Therefore, the developed coded excitation technique is a cost-effective solution to improve the SNR in the medical ultrasound system 100 by enhancing the TPE and minimizing the artifacts from dynamic receive focusing while reducing the necessary hardware complexity.
In some embodiments, the ETP-coded transmit signal stored in the ETP-coded transmit sequence memory 1402 and is transmitted through the ultrasonic transducer (i.e., the transducer elements 202). The ETP-coded transmit signal also may be utilized for decoding the receive signal in the two-stage pulse compression. The ETP-coded transmit signal may be generated based on the binary phase codes (e.g., Barker). The PSL of the Barker codes used in the present invention after matched filtering is given by
where M is the ratio of the time duration of an elongated transmit signal in coded excitation with respect to that of a conventional short transmit signal.
To improve the TPE of the elongated transmit signal based on the Barker codes, we have developed an efficient transmit power (ETP) coding technique.
as follows:
where δ(t) is the Dirac delta function.
In a block 1704, an encoding kernel may be selected. An example encoding kernel for the ETP coding is illustrated in
In a block 1706, the principal transmit signal may be encoded with the selected encoding kernel.
For some embodiments, the ETP-coded transmit signal is generated by convolving the principal transmit signal with the selected encoding kernel
e(t)=s(t)*w(t) (15)
where e(t) is the ETP-coded transmit signal, w(t) is the selected encoding kernel and * is the convolution operator.
For alternative embodiments, the external computing device 108 may generate the transmit signal, select the encoding kernel, and convolve the transmit signal with the encoding kernel to generate the transmit signal encoded with binary phase codes. The external local storage 110 may then store the transmit signal encoded with the binary phase codes, which the home ultrasound machine 102 may access via the communication network 106 and/or the network interface 132.
r(t)=a·e(t−Δτ) (16)
where a is the reflection coefficient and Δτ is the time delay corresponding to the location of the target. In the pre-compressor 1406 on each channel, a matched filter is applied to the receive signal based on the corresponding ETP-coded transmit signal as like
r(t)*e(−t)=a·e(t−Δτ)*e(−t) (17)
In case of the bi-phase rectangular window function, the ETP-coded transmit signal may be composed of the binary sequences (i.e., +1 or −1) so that the matched filter can be implemented with two's complement adders. It does not require any complex multipliers. For other window functions, the elongated transmit signal (i.e., without ETP coding) as well as the corresponding ETP-coded transmit signal can be utilized in matched filtering to remove the need for complex multipliers. In fact, the elongated transmit signal may be regarded as a subset of the ETP-coded transmit signal because it can be produced by using the single-phase rectangular window function.
After matched filtering, the pre-compressed receive signals from all the channels may be coherently combined together in the receive beamformer 1408. For reducing the sidelobes in the lateral direction, apodization can be applied as well. Thus, the beamformed receive signal, d(t), is given by
where N is the number of channels on the receiver and b(n) is the apodization coefficients.
After receive beamforming, a sidelobe suppression filter may be used in the post-compressor 1410 to remove the sidelobes in the axial direction. The sidelobe suppression filter coefficients can be obtained from the published article by Chen et al. entitled “A new algorithm to optimize Barker code sidelobe suppression filters”, IEEE Transactions on Aerospace and Electronic Systems, Vol. 26, pp. 673-677, 1990, by modeling the pre-compressed receive signal (i.e., C(f)) as the convolution between the mainlobe and sidelobe signals as follows:
C(f)=Cm(f)Cs(f) (19)
where
and Cm(f) and Cs(f) are the power spectral density of the mainlobe and sidelobe functions, respectively. By performing the inverse Fourier transform of Eq. (21), the sidelobe suppression filter coefficients are obtained. An example sidelobe suppression filter is illustrated in
Alternatively, the sidelobe suppression filter coefficients may be obtained based on the minimax and genetic optimization approaches as well. The sidelobe suppression filter has arbitrary values in its coefficients so that it may use complex multipliers as well as adders. Therefore, it is practically difficult to implement the sidelobe suppression filter with the matched filter in each channel. On the contrary, if both the matched filter and the sidelobe suppression filter are positioned after receive beamforming to reduce the hardware complexity, severe artifacts are introduced due to the distortions in the elongated transmit signals from dynamic receive focusing. However, in the two-stage pulse compression method, a matched filter is implemented in each channel to minimize the signal distortion from dynamic receive focusing while a single sidelobe suppression filter is used to reduce the sidelobes after receive beamforming.
As described above, the sidelobe suppression filter is applied to the decoded receive signal via matched filtering
where cs(t) is the sidelobe suppression filter kernel. If the matched filter has successfully reduced the distortion in the pre-compressed receive signal during receive beamforming, the convolution operation for sidelobe suppression filtering can be integrated with the summation as follows:
If the sidelobe suppression filter effectively removes the sidelobes, the convolution amongst the receive signals, the matched filter kernel and the sidelobe suppression filter kernel could be given by
e(n,t−Δτ)*e(−t)*cs(t)≈δ(n,t−Δτ) (24)
By using Eq. (24), Eq. (23) can be written by
As seen in Eq. (25), after the sidelobe suppression filter, all the information for the imaging target (i.e., the reflection coefficient and distance) can be obtained.
The pre-beamforming delay, Δ{circumflex over (τ)}(i,j,k,l) for the ith firing, the jth receive channel, the kth sub-scanline and the lth imaging point, is represented as a binary number. In the pre-beamformer processing unit 306, if the delay is ‘0’, the latch 2304 holds the data stored in it. The latch 2304 is updated by the incoming data from the demodulator 304 when the pre-beamforming delay is ‘1’. The latched content is transferred to the first-in first-out (FIFO) memory 2306 for post-beamforming. The pre-beamforming delay may also be utilized for controlling the FIFO memory 2306. In the FIFO memory 2306, only pre-beamformed complex baseband data are sequentially stacked. This allows great coarse delay memory savings and can be implemented in embedded memory of low-cost reconfigurable devices.
In some embodiments, to support multi-beam and synthetic aperture techniques, multiple pre-beamformer processing units can be integrated into the pre-beamformer.
For generating multiple scan lines simultaneously, different pre-beamforming delays are utilized in multiple pre-beamformer processing units 306. As seen in
Alternatively, multiple scan lines can be reconstructed by utilizing a single pre-beamformer processing unit with memory in place of the buffer 2402. In order to reconstruct multiple scan lines, the complex baseband data may be stored in the memory and reused in the pre-beamformer processing unit 306. In this embodiment, the multiple scan lines cannot be generated simultaneously. However, only a single pre-beamformer processing unit 306 may be used for supporting multi-beam and synthetic aperture techniques.
where d is the penetration depth. K represents the number of pre-beamformer processing units 306 or the number of scan lines for multi-beam and synthetic aperture techniques. As can be seen in the pre-beamforming delay LUT 308, only those axial points to be imaged have ‘1’, while the others have ‘0’.
For some embodiments, the size of the pre-beamformer LUT 308 may be reduced by reducing the redundancy in the control words shown in
For some embodiments, the combination of variable length coding (VLC) and run length coding (RLC) can be used to further reduce the size of the LUT 12902. When VLC/RLC encoding is used, the first lookup table (LUT 1) 2902 is coupled to a RLC/VLC decoder 2904, which may provide RLC/VLC decoding of the control words. The RLC/VLC decoder 2904 is coupled to the second lookup table (LUT 2) 2906, which may store a code book for the decoded the reduced control words.
Depending on the home ultrasound system 100 specification and complexity, the number of imaging points in the axial direction can vary. For a high-end home ultrasound system 100 or a high quality image, a large number of imaging points (e.g., 4096 points) can be used. A small number of imaging points (e.g., 512 points) can be used for a cost-effective portable home ultrasound system 100 in order to reduce the size of the FIFO memory 2306 in the pre-beamformer processing unit 306 and the computational complexity in the post-beamformer 604.
For some embodiments, the phase rotation stage 2602 receives the complex baseband signals from the reconfigurable HW 122 (or pre-beamformer) and adjusts the phase of the pre-beamformed complex baseband signal. In one embodiment, the phase rotation can be represented by
where ΔτPOS
ΔτPOS(i,j,k,l)=ΔτPOS
The phase rotation can be facilitated using complex multiplication instructions available in the programmable processing unit 124. One or more post-beamformer phase delays and/or phase compensation values may be computed in real time using a Coordinate Rotation Digital Computer (CORDIC) algorithm. Alternatively, as described earlier, the post-beamformer phase compensation values may be computed in advance, stored in the post-beamformer LUT 602, and utilized later in computing the phase rotation.
The post-beamformer LUT 602 may be organized in several ways.
The illustrated post-beamformer LUT 602 includes L rows corresponding to the number of imaging points and log2 D columns whereas D represents the precision of the computed cosine and sine values in the second post-beamforming delay LUT shown in
After the phase of the complex baseband signal is adjusted, as specified by the post-beamformer LUT 308, for example, the complex baseband data from all channels may be coherently combined together in the summation stage 2604.
Although some embodiments have been described with reference to an ultrasound machine being located in a home, embodiments are not so limited. Fore example, the ultrasound machine may be located in a clinic, such as a local neighborhood clinic, in a physician's office, and/or in a hospital, such as in a hospital emergency room, for example. The ultrasound machine also may be located in a vehicle, such as an aid vehicle, for example
In the above description, numerous specific details, such as, for example, particular processes, materials, devices, and so forth, are presented to provide a thorough understanding of embodiments of the invention. One skilled in the relevant art will recognize, however, that the embodiments of the present invention may be practiced without one or more of the specific details, or with other methods, components, etc. In other instances, structures or operations are not shown or described in detail to avoid obscuring the understanding of this description.
Reference throughout this specification to “one embodiment” or “an embodiment” means that a particular feature, structure, process, block, or characteristic described in connection with an embodiment is included in at least one embodiment of the present invention. Thus, the appearance of the phrases “in one embodiment” or “in an embodiment” in various places throughout this specification does not necessarily mean that the phrases all refer to the same embodiment. The particular features, structures, or characteristics may be combined in any suitable manner in one or more embodiments.
The terms used in the following claims should not be construed to limit embodiments of the invention to the specific embodiments disclosed in the specification and the claims. Rather, the scope of embodiments of the invention is to be determined entirely by the following claims, which are to be construed in accordance with established doctrines of claim interpretation.
This application is a divisional under 35 U.S.C. §120 of U.S. patent application Ser. No. 11/213,275, filed 26 Aug. 2005, which in turn claims priority under 35 U.S.C. §119(e) from U.S. Provisional Application No. 60/604,883, filed 27 Aug. 2004.
Number | Date | Country | |
---|---|---|---|
60604883 | Aug 2004 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 11213275 | Aug 2005 | US |
Child | 13211241 | US |