Systems and methods for reducing noise in imaging.

Abstract

To provide systems and methods for the denoising of images in the presence of broadband noise based on the detection and / or estimation of in-band noise.SOLUTION: According to various example embodiments, an estimate of broadband noise that lies within the imaging band is made by detecting or characterizing the out-of-band noise that lies outside of the imaging band. This estimated in-band noise may be employed for denoising the detected imaging waveform. According to other example embodiments, a reference receive circuit that is sensitive to noise within the imaging band, but is isolated from the imaging energy, may be employed to detect and / or characterize the noise within the imaging band. The estimated reference noise may be employed to denoise the detected in-band imaging waveform.SELECTED DRAWING: Figure 1A

Description

[Technical Field]

[0001] <Cross-references to related applications> This application claims priority to U.S. Provisional Application No. 62 / 463,431, filed February 24, 2017, entitled “SYSTEMS AND METHODS FOR NOISE REDUCTION IN IMAGING.” The entire contents of that document are incorporated herein by reference.

[0002] This disclosure relates to image processing techniques and image data processing for noise reduction. [Background technology]

[0003] Medical imaging using ultrasound and MRI relies on detecting low-amplitude signals in the radio frequency spectrum, typically ranging from 2 MHz to 200 MHz. Image quality is greatly influenced by the signal-to-noise ratio.

[0004] In intravascular ultrasound (IVUS), intracardiac echocardiography (ICE), and other forms of minimally invasive imaging, ultrasound transducers detect ultrasound signals from surrounding structures and convert their acoustic energy into electrical signals. These signals propagate along one or more conductive channels (e.g., coaxial conductors, twisted-pair conductors, flexible circuits, etc.). For a number of reasons (including cost, manufacturability, safety, biocompatibility, thermal issues, and power supply requirements), minimally invasive imaging probes that can be implanted in the body may lack amplifiers to enhance signal intensity. The electrical signals detected by minimally invasive ultrasound transducers can be very small (<10mV, generally <1mV), and much of the information about tissue structures that can be imaged by ultrasound often lies on the smaller side of the dynamic range of the detected electrical signal. The amplitude of the received ultrasound signal is limited by one or all of the following: the mechanical efficiency of the transducer, the small amplitude of the detected acoustic signal, the small size of the transducer, or attenuation along the conductors carrying the electrical signal from the transducer to the body. From this perspective, signals from minimally invasive ultrasound imaging systems tend to be very weak.

[0005] Noise can be introduced to the system from numerous sources. This includes wireless transmitters, power electronic components, transmission lines, switching transistors, and other components known in the field. Noise can be introduced through induction, direct conduction, or suboptimal isolation between electromagnetically sensitive components. Some noise may be generated by components within the imaging system itself, such as scan actuators in any electronic component of the imaging system, pulse width modulators in motor controllers, switching-mode power components, clock circuits, and transistors. Furthermore, other systems connected to the patient or in the therapeutic environment (e.g., impedance monitors, tracking systems (e.g., Carto® 3, Carto® XP, NavX® systems), temperature sensors, infusion pumps, ablation systems, ECGs, hemodynamic monitors) can introduce noise. RFID inventory control systems used in the medical field can also introduce noise.

[0006] There are approaches to reduce the amount of noise entering the ultrasonic receiving circuit of an ultrasonic imaging system. These include, for example, selecting components within the system that generate minimal RF noise, electrical isolation, shielding, proper grounding, and physically separating noise-generating components from components affected by electromagnetic noise. These approaches can be difficult to implement because noise sources may have desirable characteristics for other reasons (e.g., pulse-width modulation motor controllers are energy-efficient and have good response times), or because they are difficult to physically isolate from each other (e.g., it is desirable to place power electronic components near the imaging probe or its associated circuits).

[0007] Other approaches to reduce the impact of noise on ultrasonic signal quality (ultrasonic image quality) include filtering and image processing. Ultrasonic signals have a known bandwidth, and detected ultrasonic signals are filtered using analog or digital filtering techniques (sometimes a combination of both). Analog or digital filtering can be applied to limit a portion of the electrical signal output from the ultrasonic receiving circuit to frequency components (or their harmonics) within the operating bandwidth of the ultrasonic transducer. Selecting a filter with a narrow bandwidth and sharp cutoff can reduce the amount of noise acceptable for signals used for image generation or ultrasonic signal utilization (e.g., Doppler spectral analysis of ultrasonic signals, or evaluation of scattering flow in an ultrasonic processing field). Notch filters or comb filters are useful for removing narrowband noise within the imaging range frequency. Excessive filtering can have unwanted effects, such as reducing the amount of signal power acceptable for image generation or other uses of the ultrasonic signal. It can also negatively impact other performance aspects of the ultrasonic imaging system, such as resolution. However, if the filter's passband is too wide, the system will receive more noise.

[0008] Image processing can further reduce noise by filtering the generated image data, for example, through averaging or outlier removal. For instance, filtering can be applied within an image in the spatial domain by applying a Gaussian filter to a pixel and its adjacent pixels to smooth out blur and random noise within the image. However, this tends to degrade the spatial resolution of the image. Similarly, spatial domain filtering can be applied to imaging structures that do not move rapidly with respect to the frame repetition frequency of the imaging device. For example, pixels in an image frame are pixels at similar positions in preceding and / or succeeding frames, obtained as a result of averaging or a Gaussian filter.

[0009] A similar problem exists in MRI imaging systems, where weak signals are detected when there is noise from unwanted sources of radio frequency energy.

[0010] What is useful are methods, systems, and devices for identifying noise and dynamically removing it from one or more imaging signals.

[0011] Much noise can enter the ultrasonic receiving signal, and once it enters the system, it can be difficult to remove. This is especially true when the noise is broadband by nature, in which case some of the noise falls within the passband of the ultrasonic system. For example, in an imaging system with a transducer with a center frequency of 10 MHz and a passband of 7.5 to 12.5 MHz, the system is designed to filter out noise below 7.5 MHz and above 12.5 MHz. However, the noise within the 7.5 to 12.5 MHz band is often a considerable amount relative to the amplitude of the ultrasonic signal being detected.

[0012] Many noise sources arise as a result of fast transitions, such as when field-effect transistors or switches are turned on or off. Fast-transforming electrical signals have a very wide frequency domain, which can easily cover all or part of the passband of an ultrasonic receiving signal chain. This is especially true in power circuits and pulse-width modulation circuits where the noise has a sufficiently strong amplitude relative to the detected signal. [Overview of the project] [Means for solving the problem]

[0013] One approach to reducing broadband noise utilizes the fact that imaging energy is primarily within a selectable imaging band, and broadband noise can be detected both within and outside that imaging band. In principle, by detecting or characterizing out-of-band noise at any given time, it is possible to estimate the broadband noise that may be present within the imaging band, and then modify the detected signal to reduce the estimated in-band noise. By effectively estimating in-band noise based on out-of-band noise, it is possible to generate a signal that estimates the desired imaging energy where the estimated in-band noise is absent.

[0014] Another approach to reducing noise within the imaging band is to use a reference receiving circuit (equipped with components that can emulate resistors, capacitors, inductors, transmit lines, amplifiers, converters, inert transducers, or transducer receiving circuits) that can detect noise within the imaging band but is isolated from the imaging energy. By estimating the in-band imaging noise based on the in-band noise received by the reference receiving circuit, a signal can be generated that estimates the desired imaging energy where the estimated in-band noise is absent.

[0015] One aspect provides a method for removing noise from an imaged signal when broadband noise is present. The method is: When imaging energy is not being received, the imaging transducer receiving circuit detects an energy wave to obtain a noise characteristic identification waveform, filters the noise characteristic identification waveform to generate an in-band noise characteristic identification waveform that exists within the imaging band and an out-of-band noise characteristic identification waveform that at least a portion of which exists within the noise detection band outside the imaging band. A step of segmenting the in-band noise characteristic identification waveform and the out-of-band noise characteristic identification waveform according to one or more time windows, A step of determining the relationship between noise in the imaging band and noise in the noise detection band by processing the in-band noise characteristic identification waveform and the out-of-band noise characteristic identification waveform for at least one time window. The steps include: detecting the imaging signal with the imaging transducer receiving circuit to acquire one or more imaging waveforms; For at least one image waveform, a) A step of filtering the image waveform to generate an in-band image waveform present within the imaging band and an out-of-band noise detection image waveform present within the noise detection band. b) A step of segmenting the in-band imaging waveform and the out-of-band noise detection imaging waveform according to a time window of shape 1, c) Using the above relationship and the out-of-band noise detection imaging waveform, estimate the measurement result regarding the noise amount of the in-band imaging waveform within at least one time window. d) A step of applying noise suppression correction to a portion of the in-band imaging waveform within at least one time window in step c), It holds.

[0016] On one side, a method is provided for removing noise from an imaged signal when noise is present. The method is: When imaging energy is not being received, The steps include: detecting an energy wave using an imaging transducer receiving circuit to acquire a noise characteristic identification waveform; filtering the noise characteristic identification waveform to generate an in-band noise characteristic identification waveform present within the imaging band; A step of acquiring a reference noise characteristic identification waveform by detecting noise using a reference receiving circuit configured not to transmit imaging energy while the imaging transducer receiving circuit is detecting noise, A step of processing the in-band noise characteristic identification waveform and the reference noise characteristic identification waveform to determine the relationship between the noise in the imaging band and the noise detected by the reference receiving circuit. The steps include: detecting the imaging signal with the imaging transducer receiving circuit to acquire one or more imaging waveforms; For at least one image waveform, a) A step of filtering the imaging waveform to generate an in-band imaging waveform that exists within the imaging band, b) A step of detecting a reference noise detection waveform using the reference receiving circuit, c) A step of segmenting the in-band imaging waveform and the reference noise detection waveform according to a time window of 1 or more, d) Using the relationship and the reference noise detection waveform, estimate the measurement result regarding the noise amount of the in-band imaging waveform within at least one time window. e) A step of applying noise suppression correction to a portion of the in-band imaging waveform within at least one time window in step c), It holds.

[0017] One aspect provides a method for removing noise from an imaged signal when broadband noise is present. The method is: The steps include: detecting an energy wave with an imaging transducer receiving circuit to acquire an imaging waveform; filtering the imaging waveform to generate an in-band imaging waveform that exists within the imaging band and an out-of-band noise detection imaging waveform that exists within a noise detection band where at least a portion is outside the imaging band; A step of detecting the in-band imaging envelope of the in-band imaging waveform, Steps include detecting the out-of-band envelope of the out-of-band noise detection imaging waveform, The step of obtaining a modified outer band envelope by applying a scaling factor to the aforementioned outer band envelope, A step of obtaining a noise-corrected in-band envelope by combining the corrected out-of-band envelope and the in-band imaging envelope. It has, The scaling coefficient is selected to reduce the contribution of intraband noise to the noise-corrected intraband envelope.

[0018] One aspect provides a method for removing noise from an image signal when noise is present. The method is: The steps include: detecting an energy wave with an imaging transducer receiving circuit to acquire an imaging waveform; filtering the imaging waveform to generate an in-band imaging waveform that exists within the imaging band and an out-of-band noise detection imaging waveform that exists within a noise detection band where at least a portion is outside the imaging band; A step of obtaining a corrected waveform by applying an amplitude scaling coefficient and a frequency shift to the out-of-band noise detection imaging waveform, thereby ensuring that the corrected waveform includes frequency components present within the imaging band. A step of obtaining a noise-corrected in-band image waveform by combining the corrected waveform and the in-band image waveform, It has, The amplitude scaling coefficient is selected to reduce the contribution of in-band noise to the noise-corrected in-band imaged waveform.

[0019] One aspect provides a method for removing noise from an image signal when noise is present. The method is: The steps include: detecting an energy wave using an imaging transducer receiving circuit to acquire an imaging waveform; filtering the imaging waveform to generate an in-band imaging waveform that exists within the imaging band; A step of acquiring a reference noise detection waveform by detecting noise using a reference receiving circuit configured not to transmit imaging energy while the imaging transducer receiving circuit is detecting the noise received, A step of detecting the in-band imaging envelope of the in-band imaging waveform, A step of detecting the reference envelope of the reference noise detection waveform, The step of obtaining a modified reference envelope by applying a scaling factor to the aforementioned reference envelope, A step of obtaining a noise-corrected in-band envelope by combining the corrected reference envelope and the in-band imaging envelope. It has, The scaling coefficient is selected to reduce the contribution of intraband noise to the noise-corrected intraband envelope.

[0020] One aspect provides a method for removing noise from an image signal when noise is present. The method is: The steps include: detecting an energy wave using an imaging transducer receiving circuit to acquire an imaging waveform; filtering the imaging waveform to generate an in-band imaging waveform that exists within the imaging band; A step of acquiring a reference noise detection waveform by detecting noise using a reference receiving circuit configured not to transmit imaging energy while the imaging transducer receiving circuit is detecting the noise received, A step of adaptively filtering the reference noise detection waveform according to one or more adaptive filter parameters, A step of obtaining a noise-corrected in-band imaging waveform by combining the filtered reference noise detection waveform and the in-band imaging waveform. It has, The adaptive filter parameters are dynamically determined by processing the noise-corrected in-band imaging waveform to minimize the power of the noise-corrected in-band imaging waveform.

[0021] One aspect provides a method for removing noise from an imaged signal when broadband noise is present. The method is: The steps include: detecting an energy wave with an imaging transducer receiving circuit to acquire an imaging waveform; filtering the imaging waveform to generate an in-band imaging waveform that exists within the imaging band and an out-of-band noise detection imaging waveform that exists within a noise detection band where at least a portion is outside the imaging band; A step of detecting the in-band imaging envelope of the in-band imaging waveform, Steps include detecting the out-of-band imaging envelope of the out-of-band noise detection imaging waveform, A step of adaptively filtering the out-of-band imaging envelope according to one or more adaptive filter parameters, A step of obtaining a noise-corrected in-band imaging envelope by combining the filtered out-of-band imaging envelope and the in-band imaging envelope. It has, The adaptive filter parameters are dynamically determined by processing the noise-corrected intraband imaging envelope to minimize the power of the noise-corrected intraband imaging envelope.

[0022] One aspect provides a method for removing noise from an imaged signal when broadband noise is present. The method is: The steps include: detecting an energy wave with an imaging transducer receiving circuit to acquire an imaging waveform; filtering the imaging waveform to generate an in-band imaging waveform that exists within the imaging band and an out-of-band noise detection imaging waveform that exists within a noise detection band where at least a portion is outside the imaging band; A step of obtaining a corrected waveform by applying a frequency shift to the out-of-band noise detection imaging waveform, thereby ensuring that the corrected waveform includes frequency components present within the imaging band. A step of adaptively filtering the out-of-band noise detection imaging waveform according to one or more adaptive filter parameters, A step of obtaining a noise-corrected in-band image waveform by combining the filtered corrected waveform and the in-band image waveform. It has, The adaptive filter parameters are dynamically determined by processing the noise-corrected in-band imaging waveform to minimize the power of the noise-corrected in-band imaging waveform.

[0023] One aspect provides a method for removing noise from an imaged signal when broadband noise is present. The method is: The steps include: detecting an energy wave with an imaging transducer receiving circuit to acquire an imaging waveform; filtering the imaging waveform to generate an in-band imaging waveform that exists within the imaging band and an out-of-band noise detection imaging waveform that exists within a noise detection band where at least a portion is outside the imaging band; The steps include: selecting appropriate filter parameters for a dynamic digital filter that filters the in-band image waveform to remove in-band noise by processing the out-of-band noise detection image waveform; A step of filtering the in-band image waveform by the dynamic digital filter according to the filter parameters, It holds.

[0024] The present invention provides a method for performing noise suppression on a signal acquired by a detection system characterized by one or more noise sources. The method is as follows: When imaging energy is not being received, the imaging transducer receiving circuit detects an energy wave to obtain a noise characteristic identification waveform, filters the noise characteristic identification waveform to generate an in-band noise characteristic identification waveform that exists within the imaging band and an out-of-band noise characteristic identification waveform that exists within a noise detection band where at least a portion is outside the imaging band. A step of segmenting the in-band noise characteristic identification waveform and the out-of-band noise characteristic identification waveform according to one or more time windows, For at least one time window, process the in-band noise characteristic identification waveform and the out-of-band noise characteristic identification waveform according to a pattern recognition algorithm to identify a noise pattern that correlates with the noise in the imaging band within the noise detection band. The imaging transducer receiving circuit detects the imaging signal to acquire the imaging waveform, filters the imaging waveform to acquire the in-band imaging waveform present within the imaging band and the out-of-band noise detection imaging waveform present within the noise detection band. A step of segmenting the in-band image waveform and the out-of-band image waveform according to one or more time windows, The steps include: processing the out-of-band noise detection image waveform for at least one time window according to the pattern recognition algorithm to detect the noise pattern; When the noise pattern is detected, a step is to apply noise suppression correction to the time window of the in-band imaging waveform that is specific to the noise pattern detected in the out-of-band noise detection imaging waveform. It holds.

[0025] The present invention provides a method for performing noise suppression on a signal detected by a detection system characterized by one or more known noise sources. The method is: When imaging energy is not being received, A step of obtaining a noise characteristic identification waveform by detecting an energy wave with an imaging transducer, and filtering the noise characteristic identification waveform to generate an in-band noise characteristic identification waveform present within the imaging band. A step of acquiring a reference noise characteristic identification waveform by detecting noise using a reference receiving circuit configured not to transmit imaging energy while the imaging transducer receiving circuit is detecting the received noise, A step of segmenting the in-band noise characteristic identification waveform and the reference noise characteristic identification waveform according to one or more time windows, A step of processing the in-band noise characteristic identification waveform and the reference noise characteristic identification waveform for at least one time window to determine the relationship between the noise in the imaging band and the noise detected by the reference receiving circuit. A step of processing the in-band noise characteristic identification waveform and the reference noise characteristic identification waveform according to a pattern recognition algorithm to identify the presence of a noise pattern in the reference noise characteristic identification waveform that correlates with the noise in the in-band noise characteristic identification waveform. The imaging transducer receiving circuit detects the imaging signal and acquires the imaging waveform, and the reference receiving circuit detects a reference noise detection waveform, filters the imaging waveform, and acquires an in-band imaging waveform that exists within the imaging band. A step of segmenting the in-band imaging waveform and the reference noise detection waveform according to one or more time windows, Steps include: processing the reference noise detection waveform according to the pattern recognition algorithm for at least one time window to detect the presence of the noise pattern; and, when the noise pattern is detected, applying noise suppression correction to the time window of the in-band imaging waveform specific to the noise pattern detected within the reference noise detection waveform. It holds.

[0026] One aspect provides a method for removing noise from an image signal when noise is present. The method is: When imaging energy is not being received, the imaging transducer receiving circuit detects an energy wave to obtain a noise characteristic identification waveform, filters the noise characteristic identification waveform to generate an in-band noise characteristic identification waveform that exists within the imaging band and an out-of-band noise characteristic identification waveform that exists within a noise detection band where at least a portion is outside the imaging band. The imaging transducer receiving circuit detects the imaging signal to acquire the imaging waveform, filters the imaging waveform to generate an in-band imaging waveform present within the imaging band and an out-of-band noise detection imaging waveform present within the noise detection band. A step of performing cross-correlation between the out-of-band imaging waveform and the out-of-band noise characteristic identification waveform to determine the time delay corresponding to the maximum cross-correlation. A step of obtaining a modified in-band noise characteristic identification waveform by applying the time delay and amplitude adjustment to the in-band noise characteristic identification waveform, and subtracting the modified in-band noise characteristic identification waveform from the in-band imaging waveform, It holds.

[0027] One aspect provides a method for removing noise from an image signal when noise is present. The method is: When imaging energy is not being received, The steps include: detecting an energy wave using an imaging transducer receiving circuit to acquire a noise characteristic identification waveform; filtering the noise characteristic identification waveform to generate an in-band noise characteristic identification waveform present within the imaging band; A step of acquiring a reference noise characteristic identification waveform by detecting noise using a reference receiving circuit configured not to transmit imaging energy while the imaging transducer receiving circuit is detecting the received noise, The imaging transducer receiving circuit detects the imaging signal and acquires the imaging waveform, and the reference receiving circuit detects a reference noise detection waveform, filters the imaging waveform, and acquires an in-band imaging waveform that exists within the imaging band. A step of performing cross-correlation between the reference noise detection waveform and the reference noise characteristic identification waveform to determine the time delay corresponding to the maximum cross-correlation. A step of obtaining a modified in-band noise characteristic identification waveform by applying the time delay and amplitude adjustment to the in-band noise characteristic identification waveform, and subtracting the modified in-band noise characteristic identification waveform from the in-band imaging waveform, It holds.

[0028] One aspect provides a method for removing noise from an image signal when noise is present. The method is: A step of acquiring multiple imaging waveforms by detecting imaging signals along multiple adjacent scan lines using an imaging transducer receiving circuit, For at least two adjacent scanlines, The steps include filtering the image waveform to generate an in-band image waveform that exists within the imaging band and an out-of-band noise detection image waveform that exists within a noise detection band, at least a portion of which is outside the imaging band, A step of segmenting the in-band image waveform and the out-of-band noise detection image waveform according to a series of time windows, For at least one window, A step of processing the out-of-band noise detection image waveform to determine whether or not noise correction should be applied to the corresponding window portion of the in-band image waveform. If the in-band imaging waveform within the time window is deemed suitable for noise correction, the step is to apply noise suppression correction to the in-band imaging waveform within the time window, wherein the noise suppression correction for each sample within the window is based on one or more statistical measurement results for samples in the in-band imaging waveform of two or more adjacent windows, and each adjacent window is located within an adjacent scanline. A step of generating an image based on noise-suppressed in-band imaging waveforms corresponding to multiple scan lines, It holds.

[0029] One aspect provides a method for removing noise from an image signal when noise is present. The method is: The steps include: detecting an energy wave with an imaging transducer receiving circuit to acquire an imaging waveform; filtering the imaging waveform to generate an in-band imaging waveform that exists within the imaging band and an out-of-band noise detection imaging waveform that exists within a noise detection band where at least a portion is outside the imaging band; A step of acquiring multiple imaging waveforms by detecting imaging signals along multiple adjacent scan lines using the imaging transducer receiving circuit, A step of processing one or more out-of-band noise detection imaging waveforms to determine the period of the noise source; a step of adjusting the scan rate so that the noise does not become time-synchronized in the in-band imaging waveform corresponding to the adjacent scan line. For at least two adjacent scanlines, A step of segmenting the in-band image waveform according to a series of time windows, For at least one window, A step of applying noise suppression correction to the in-band imaging waveform within the time window, wherein the noise suppression correction for each sample in the window is based on one or more statistical measurement results for samples in the in-band imaging waveform of two or more adjacent windows, and each adjacent window is located within the corresponding adjacent scanline, A step of generating an image based on noise-suppressed in-band imaging waveforms corresponding to multiple scan lines, It holds.

[0030] One aspect provides a method for removing noise from an image signal when noise is present. The method is: For at least two adjacent scanlines, The imaging transducer receiving circuit detects the imaging signal and acquires the imaging waveform, while the reference receiving circuit detects a reference noise detection waveform, filters the imaging waveform, and acquires an in-band imaging waveform that exists within the imaging band. A step of segmenting the in-band imaging waveform and the reference noise detection waveform according to a series of time windows, For at least one window, A step of processing the reference noise detection waveform and determining whether or not noise correction should be applied to the corresponding window portion of the imaging waveform within the band. If the in-band imaging waveform within the time window is deemed suitable for noise correction, the step is to apply noise suppression correction to the in-band imaging waveform within the time window, wherein the noise suppression correction for each sample within the window is based on one or more statistical measurement results relating to samples of the in-band imaging waveform in two or more adjacent windows, and each adjacent window is located within an adjacent scanline; the step is to generate an image based on the noise-suppressed in-band imaging waveforms corresponding to each of the multiple scanlines; It holds.

[0031] Further understanding of this disclosure and its functional and useful aspects can be achieved by referring to the following detailed description and drawings.

[0032] The embodiments will be described for illustrative purposes only, with reference to the drawings. [Brief explanation of the drawing]

[0033] [Figure 1A] An example of an ultrasonic imaging system configured for noise suppression is shown. [Figure 1B] This shows an example of a conventional ultrasound receiver signal chain that processes ultrasound signals before converting them into ultrasound images. [Figure 1C] An example of an ultrasound imaging system with an internal imaging probe is shown. [Figure 1D] An example of an ultrasonic imaging system with a reference transducer for detecting intraband noise is shown. [Figure 1E] An example of an ultrasound imaging system is shown, comprising a second imaging transducer having an imaging band outside the imaging band of the first imaging transducer. The second imaging transducer is, for example, part of a circuit suitable for detecting in-band noise affecting the signal received from the first imaging transducer. [Figure 1F] An example of an ultrasonic imaging system equipped with a reference receiving circuit for detecting in-band noise is shown. The reference receiving circuit extends to a position within the imaging probe. [Figure 1G] An example of an ultrasonic imaging system equipped with a reference receiving circuit for detecting in-band noise is shown. The reference receiving circuit is located at one or more positions outside the system's imaging probe. [Figure 2A] This example demonstrates a system configuration that reduces noise on the envelope of the input waveform by suppressing estimated in-band noise. In-band noise is estimated by performing envelope detection on the out-of-band waveform and applying delay, scale, and shape adjustments before subtraction. [Figure 2B]This example demonstrates a system configuration that reduces noise on the input waveform by suppressing estimated in-band noise. In-band noise is estimated by frequency-shifting the out-of-band waveform, filtering the frequency-shifted out-of-band noise, and performing delay, scale, and shape adjustments before subtraction. [Figure 2C] An example of a system configuration for noise reduction is shown. The system includes a reference receiver circuit that detects some or all of the in-band noise detectable by the imaging transducer receiver circuit and is at least partially isolated from the imaging signal detected by the imaging transducer receiver circuit. By subtracting the noise signal detected by the reference receiver circuit from the signal received by the imaging transducer receiver circuit, the noise in the output signal can be reduced. [Figure 3A] This example shows a system configuration that reduces noise using dynamic noise cancellation. The cancellation waveform of the dynamic noise cancellation is obtained from the reference receiving circuit. [Figure 3B] This example shows a system configuration that reduces noise on the input waveform using dynamic noise cancellation. The cancellation waveform for dynamic noise cancellation is obtained by envelope detection of the out-of-band waveform. [Figure 3C] This example shows a system configuration that reduces noise on the input waveform using dynamic noise cancellation. The cancellation waveform for dynamic noise cancellation is obtained by frequency-shifting the out-of-band waveform and filtering the frequency-shifted out-of-band waveform. [Figure 4] This example demonstrates a system configuration that reduces noise on the input waveform through dynamic noise cancellation. Filtering is controlled based on feedback parameters obtained by a filter update algorithm that determines one or more filter parameters based on one or more characteristics of the noise-detected waveform. [Figure 5A]This outlines a system configuration example that reduces noise on the input waveform by applying one or more noise reduction algorithms that utilize noise parameters acquired in a first measurement stage where no imaging signal is present, and then using those noise parameters in a second measurement stage where an imaging signal is collected. [Figure 5B] An outline of another system example that detects in-band noise using a reference receiving channel is shown. [Figure 6A] This example shows a system configuration that reduces noise on the input waveform based on noise detection in the out-of-band waveform. By processing each window of the out-of-band waveform, different time windows of the in-band waveform are suppressed, and the noise window of the in-band waveform is corrected by subtracting a subtraction value that depends on the power amount within the out-of-band waveform window. [Figure 6B] An example of a scattering plot is shown. It shows the signal power of the window for the in-band waveform and the signal power of each window for the out-band waveform in the noise characteristic discrimination stage. [Figure 6C] This example shows a system configuration that reduces noise on the input waveform based on noise detection in the out-of-band waveform. By processing each window of the out-of-band waveform, different time windows of the in-band waveform are suppressed, and the noise window of the in-band waveform is corrected by subtracting a subtraction value that depends on the power amount within the out-of-band waveform window. [Figure 6D] This example shows a system configuration that reduces noise on the input waveform based on noise detection in the out-of-band waveform. By processing each window of the out-of-band waveform, different time windows of the in-band waveform are suppressed, and the noise window of the in-band waveform is corrected by multiplying it by an attenuation coefficient that depends on the power amount within the out-of-band waveform window. [Figure 6E] The chart shows a method by which different time windows of waveforms within a band are first identified as the primary signal or noise, then the noise windows around the signal windows are identified and reclassified as potentially error-prone, and finally the signal windows around the noise windows are identified and reclassified as errors. [Figure 6F] The chart shows a method by which different time windows of waveforms within a band are first identified as the primary signal or noise, then the noise windows around the signal windows are identified and reclassified as potentially error-prone, and finally the signal windows around the noise windows are identified and reclassified as errors. [Figure 6G] This example shows a system configuration that reduces noise on the input waveform based on noise detection within the filtered reference waveform measured by a reference receiving channel. Different time windows of the in-band waveform are suppressed by processing each window of the reference waveform, and the noise window of the in-band waveform is corrected by subtracting a subtraction value that depends on the power amount within the window of the filtered reference waveform. [Figure 6H] This example shows a system configuration that reduces noise on the input waveform based on noise detection within the filtered reference waveform measured by a reference receiving channel. Different time windows of the in-band waveform are suppressed by processing each window of the reference waveform, and the noise window of the in-band waveform is corrected by subtracting a subtraction value that depends on the power amount within the window of the filtered reference waveform. [Figure 6I] This example shows a system configuration that reduces noise on the input waveform based on noise detection within the filtered reference waveform measured by a reference receiving channel. Different time windows of the in-band imaging waveform are corrected by processing each window of the filtered reference waveform, and the noise window of the in-band imaging waveform is corrected by an attenuation coefficient that depends on the power amount within the window of the filtered reference waveform. [Figure 7A] This example shows a system configuration that reduces noise on the input waveform based on noise detected in one or more noise detection waveforms. At least one noise detection waveform has an out-of-band signal within the imaging band. Different time windows of the in-band waveform are denoised according to one or more patterns identified by processing one or more noise detection waveforms. [Figure 7B]This example shows a system configuration that reduces noise on the input waveform based on noise detected in one or more noise detection waveforms. At least one noise detection waveform has an out-of-band signal within the imaging band. Different time windows of the in-band waveform are denoised according to one or more patterns identified by processing one or more noise detection waveforms. [Figure 7C] This example shows a system configuration that reduces noise on the input waveform based on noise detected in a reference waveform. Different time windows of the in-band imaging waveform are denoised according to one or more patterns identified by processing one or more reference waveforms. [Figure 7D] This example shows a system configuration that reduces noise on the input waveform based on noise detected in a reference waveform. Different time windows of the in-band imaging waveform are denoised according to one or more patterns identified by processing one or more reference waveforms. [Figure 8A] This example demonstrates a system configuration that reduces noise on the input waveform based on noise detected in the out-of-band waveform. Different time windows of the in-band waveform are denoised according to the estimated in-band noise, which is time-aligned before subtraction. [Figure 8B] This example demonstrates a system configuration that reduces noise on the input waveform based on noise detected in the out-of-band waveform. Different time windows of the in-band waveform are denoised according to the estimated in-band noise, which is time-aligned before subtraction. [Figure 8C] This example shows a system configuration that reduces noise on the input waveform based on noise detected in the filtered reference waveform. Different time windows of the in-band waveform are denoised according to the estimated in-band noise that was time-aligned before subtraction. [Figure 8D] This example shows a system configuration that reduces noise on the input waveform based on noise detected in the filtered reference waveform. Different time windows of the in-band waveform are denoised according to the estimated in-band noise that was time-aligned before subtraction. [Figure 8E]This example shows a system configuration that reduces noise on the input waveform based on noise detected in the out-of-band waveform. When performing noise correction, measurement results from adjacent scanlines or duplicate scanlines are used. [Figure 8F] This example shows a system configuration that reduces noise on the input waveform based on noise detected in the filtered reference waveform. When performing noise correction, measurement results from adjacent scanlines or duplicate scanlines are used. [Figure 9] An example of a magnetic resonance imaging system configured for noise suppression is shown. [Figure 10A] Examples of images acquired using an intracardiac echocardiography system are shown. Figure 10A shows an image acquired without a noise source. [Figure 10B] The following are examples of images acquired using an intracardiac echocardiography system. Figure 10B shows an image acquired when noise generated by the electroanatomical mapping system is present. [Figure 10C] The following are examples of images acquired using an intracardiac echocardiography system. Figure 10C shows an image acquired when a noise source generated by the ablation generator is present. [Figure 11A] The images shown are those acquired when noise from the electroanatomical mapping system is present. Figure 11A shows the undenoised image. [Figure 11B] Figure 11B shows the image acquired when noise is present from the electroanatomical mapping system. [Figure 11C] The image shown is obtained after noise reduction by attenuation, when noise from the electroanatomical mapping system is present. The relaxation parameter in Figure 11C is 0.5. [Figure 11D] The image shown is obtained when noise from the electroanatomical mapping system is present, after noise reduction by attenuation. The relaxation parameter in Figure 11D is 1. [Figure 11E] The image shown is obtained when noise from the electroanatomical mapping system is present, after noise reduction by attenuation. The relaxation parameter in Figure 11E is 1.5. [Figure 12A] The image shown is obtained when noise from the ablation generator is present. Figure 12A shows the image without noise reduction. [Figure 12B] Figure 12B shows the image acquired when noise from the ablation generator is present. [Figure 13A] The image shown was acquired when noise from the magnetic tracking system was present. Figure 13A shows the image without noise reduction. [Figure 13B] The image shown is obtained when noise from the magnetic tracking system is present. Figure 13B shows the denoised image. [Figure 14] This section shows phrases used to refer to waveforms within the imaging band, noise detection band, and waveforms from the reference receiver circuit. [Modes for carrying out the invention]

[0034] Various embodiments and aspects of this disclosure are described in detail below. The following description and drawings are for illustrative purposes only and should not be construed as limitations of this disclosure. Specific details are provided to help understand the embodiments of this disclosure. However, in some cases, known or existing details are omitted for the sake of simplicity in describing the embodiments of this disclosure.

[0035] In this specification, the terms “equipped with” and “equipped with” should be interpreted as inclusive, open-ended, and non-exclusive. Specifically, in this specification and the claims, the terms “equipped with” and “equipped with” and their derivatives mean that a particular feature part, step, or component is included. These terms should not be interpreted as excluding other feature parts, steps, or components.

[0036] In this specification, the term “Example” means “serving as an example, embodiment, or explanatory example,” and should not be construed as being preferable or advantageous to any other configuration in this specification.

[0037] In this specification, the terms "approximately" and "abbreviated" mean covering variations that exist within the upper and lower limits of a value range, such as variations in characteristics, parameters, and sizes. Unless otherwise specified, the terms "approximately" and "abbreviated" mean plus or minus 25%.

[0038] Unless otherwise explicitly stated, specific scopes and groups are simplified ways of referring to each individual member of that scope or group, and the same applies to any included subranges or subgroups. Unless otherwise explicitly stated, this disclosure includes all members and combinations of subranges and subgroups.

[0039] In this specification, when the term "order of ~" is used with a quantity or parameter, it means a range from about one-tenth to about ten times that quantity or parameter.

[0040] Ultrasonic imaging relies on received echoes from a medium, and optionally transmits narrow acoustic pulses to the medium in a specific direction. In this specification, the term "scan line" refers to a line representing the spatial direction in which imaging energy is received in the medium. A 2D image is acquired by receiving echoes from multiple scan lines within the medium. The inventors have devised, developed, and tested various methods and systems to effectively reduce broadband noise from ultrasonic acquisition and / or processing systems.

[0041] Refer to Figure 1A. An example of an ultrasonic imaging system is shown. One or more ultrasonic transducers 10 are controlled to perform ultrasonic imaging across multiple scan lines 12. The transducers 10 interface with control processing hardware 100, which optionally controls transmitters 15 that generate and emit imaging energy from the transducers 10. The control processing hardware 100 is configured to receive the ultrasonic energy detected by the transducers 10. The ultrasonic energy is routed to one or more amplifiers 20 via a Tx / Rx (transmit-receive) switch 25.

[0042] The ultrasonic transducer 10 can optionally be configured to image spatial regions corresponding to multiple scan lines 12. This can be done, for example, by mechanical scanning of the transducer 10, or by electronic scanning using an image sensor array. The image sensor array can be, for example, a phase array, a ring array, a linear array, a matrix array, a curve array, etc., but is not limited to these. In the latter case, a transmitting beamformer 26 and a receiving beamformer 27 can be used to generate multiple transmitting signals and beamform multiple receiving signals.

[0043] In this specification, the term "receiving circuit" generally refers to transmitting lines (e.g., coaxial, PCB tracing, etc.), connectors, MUX / DEMUX, RX / TX switches 25, amplifiers 20, slip rings, converters, and other components known in the art.

[0044] In this specification, the term "transducer receiving circuit" includes a receiving circuit connected to one or more ultrasonic transducer elements 10 configured to receive ultrasonic signals during use.

[0045] In this specification, the term “ultrasonic receiving signal chain” includes a receiving circuit, but may also include other components, such as an analog-to-digital converter (ADC) and digital processing components and / or processing logic circuits. This includes, but is not limited to, a noise reduction processing module 150 before the signal is converted into an image (e.g., via scan conversion) and subjected to subsequent image processing.

[0046] In this specification, the term "channel" refers to a conductive electrical circuit, a wireless channel, an optical channel, or other signal path. For example, Figure 1A shows an imaging reception channel 13 that indicates the path through which the detection imaging signal intersects. The system may use one or more reception channels per transducer (for example, in the case of an array transducer, where there is a channel for each electrostatic transducer element or group thereof). Signals can also be multiplexed along channels from one or more electrostatic transducer elements using an ASIC or other device along the signal reception chain.

[0047] In “imaging mode,” the system can be configured to control the transducer 10 to optionally transmit energy to the medium and detect imaging energy within an imaging frequency band (referred to as the “imaging band”). The imaging band may consist of a single continuous frequency band or two or more frequency intervals for detecting imaging energy. Imaging energy or noise within the imaging band is considered to be “in the band.”

[0048] Figure 14 shows the groupings and terminology used to represent various waveforms.

[0049] The system can also be configured to detect energy in one or more other frequency bands, where at least one frequency band at least partially exceeds the imaging band, via one or more channels connected to the transducer. These one or more other frequency bands are called “detection bands.” Waveforms that at least partially exceed the imaging band are called “out-of-band.” Detection bands may also be within the imaging band. Waveforms that are entirely within the imaging band and have frequency components within the entire imaging band or within a subband of the imaging band are called “in-band.” Noise detection bands can be out-of-band or in-band. When using a transducer in imaging mode (i.e., when the transducer detects imaging energy), at least one detection band can be selected so that the signal-to-noise ratio within the detection band is substantially less than the signal-to-noise ratio of the imaging band. For example, the detection band may be outside the full width at half maximum of the imaging band, or outside other bands corresponding to a threshold below the maximum intensity of the signal used.

[0050] In this specification, "imaging waveform" refers to the waveform (analog or digitally sampled) obtained from the imaging transducer receiving circuit when the imaging transducer is receiving or is assumed to be receiving imaging energy.

[0051] In this specification, the phrase "in-band image waveform" refers to an image waveform (analog or digitally sampled) that falls within the imaging band. The in-band image waveform is assumed to contain imaging energy and unwanted noise energy. In the embodiments of this disclosure, the in-band image waveform is processed to remove noise energy in order to generate a noise-free image.

[0052] In this specification, the phrase "detection band imaging waveform" refers to a waveform acquired from the imaging transducer receiving circuit and located within one or more detection bands. Detection band imaging waveforms can be either out-of-band or in-band. For example, an in-band noise detection imaging waveform can be used to confirm the presence of noise within the imaging band. More specifically, an "in-band noise detection imaging waveform" can be used to confirm that a noise source having noise components outside the imaging band also has noise components within the imaging band. A detection band imaging waveform in which at least a portion is outside the imaging band is called an "out-of-band noise detection waveform."

[0053] Refer to Figure 1A. When transducer 10 is not transmitting energy to the medium and is not detecting imaging energy from the medium, the system can be configured in "noise characteristic identification mode".

[0054] In this specification, the phrase "noise characteristic identification waveform" refers to the waveform acquired when the imaging transducer is not receiving imaging energy.

[0055] In this specification, the phrase "intraband noise characteristic identification waveform" refers to a waveform that exists within the imaging band and is acquired from the imaging transducer receiving channel when the imaging transducer is not receiving imaging energy.

[0056] In this specification, the phrase "detection band noise characteristic identification waveform" refers to a waveform obtained from the imaging transducer receiving channel when the imaging transducer is located within the noise detection band and is not receiving imaging energy. A detection band noise characteristic identification waveform in which at least a portion is outside the imaging band is called an "out-of-band noise characteristic identification waveform." A detection band noise characteristic identification waveform in which the entirety is within the imaging band is called an "in-band noise characteristic identification waveform."

[0057] In this specification, the phrase "baseline noise characteristic identification waveform" refers to the waveform acquired when the imaging transducer is not receiving imaging energy and the selected noise source is assumed to be OFF (i.e., not generating noise), thereby providing a baseline for the selected noise source. The baseline noise characteristic identification waveform within the imaging band is referred to as the "in-band baseline noise characteristic identification waveform." The baseline noise characteristic identification waveform within the noise detection band is referred to as the "detection band baseline noise characteristic identification waveform."

[0058] Refer to Figure 1A. An optional reference receiver circuit 11 can be provided. This circuit is configured not to receive reflected ultrasonic signals during imaging and includes a receiver capable of detecting noise-like noise energy coupled with one or more transducer receiver circuits during imaging. The reference receiver circuit can use components from one or more transducer receiver circuits (for example, the reference receiver circuit and the transducer receiver circuits can utilize different channels of the amplifier or ADC).

[0059] In one embodiment, the system can be configured to detect noise within the imaging band via one or more reference receiving channels. The reference receiving channels are optionally connected to acoustically isolated or deactivated reference ultrasonic transducers (not shown) that detect noise received by the imaging transducer receiving circuit without converting reflected ultrasonic waves. The one or more imaging transducers 10 and the one or more reference transducers do not need to be oriented in a common spatial direction.

[0060] The signal received by one or more reference transducer receiving circuits or reference electrical receiving circuits (on a reference receiving channel) is called the reference waveform. The reference waveform is mainly noise, not imaging energy.

[0061] In this specification, the phrase "reference waveform" refers to a waveform obtained from one or more reference receiving channels. The reference waveform may be filtered to be within and / or outside the imaging band.

[0062] In this specification, the phrase "reference noise detection waveform" refers to a reference waveform obtained from a reference receiving channel when the imaging transducer is not receiving or is assumed not to be receiving imaging energy.

[0063] In this specification, the phrase "reference noise characteristic identification waveform" refers to the reference waveform obtained from the reference receiving channel when the imaging transducer is not receiving imaging energy.

[0064] The system can optionally be configured to suppress noise using a combination of a detection band waveform and a reference waveform. In this specification, the phrase "noise detection waveform" refers to either the reference waveform or the detection band waveform. When the system is in imaging mode and the imaging transducer receiving circuit is receiving or is assumed to be receiving imaging energy, the noise detection waveform is referred to as the "noise detection imaging waveform." When the system is in noise characteristic identification mode and the imaging transducer receiving circuit is not receiving imaging energy, the noise detection waveform is referred to as the "noise detection characteristic identification waveform."

[0065] Although Figure 1A shows a single transducer element, it should be understood that the embodiment shown in Figure 1A is merely a non-limiting configuration example, and transducers having multiple electrostatic elements can be used. For example, in one embodiment, multiple transducer elements may be controlled as a phase array, a linear array, or a 2D array. Furthermore, the transducers are not limited to those that transmit imaging energy to generate multidimensional 2D cross-sectional images or 3D images (including 4D imaging datasets with 3D images over time), but also include transducers used for Doppler evaluation of flow, transducers used as ultrasonic beacons (e.g., as described in U.S. Patent Publication 2016 / 0045184, titled “Active localization and visualization of minimally invasive devices using ultrasound,” which is incorporated in its entirety by reference), or ultrasonic transducers used to detect the position of moving elements (e.g., as described in U.S. Patent Publication 2012 / 0197113, titled “Ultrasonic probe with ultrasonic transducers addressable on common electrical channel,” which is incorporated in its entirety by reference).

[0066] The transducer is not limited to one that transmits and receives imaging energy, as shown in Figure 1A, but also includes transducers that receive ultrasonic energy from a medium excited by other means. For example, optical energy (photoacoustic imaging) or another ultrasonic transducer. Furthermore, although Figure 1A shows a configuration for imaging spatial regions corresponding to multiple scan lines in different directions, the scan lines may be unidirectional. For example, this is the case when evaluating flow with M-mode imaging, Doppler imaging equipment, pulsed waves, or continuous wave Doppler.

[0067] In one embodiment, one transducer receiving channel is configured to receive imaging energy within the imaging band and simultaneously receive other energy within one or more noise detection bands. This other energy includes at least out-of-band noise. In another embodiment, one or more imaging transducer receiving channels can be used to receive imaging energy within the imaging band, and one or more transducer receiving channels can be used to receive other energy within one or more noise detection bands, at least part of which includes out-of-band noise. In another embodiment, one or more reference receiving channels can be used to receive noise energy (i.e., a reference noise detection waveform) while isolating it from imaging energy within the imaging band. The reference receiving channels can be filtered in the same way as the imaging transducer receiving channels by using an imaging bandpass filter. Alternatively, in the embodiment, the reference receiving channels may not be filtered at all, or a filter other than an imaging bandpass filter may be used to improve noise estimation within the imaging band.

[0068] The control processing hardware 100 includes, for example, one or more processors 110, memory 115, system bus 105, one or more input / output devices 120, and multiple optional devices (e.g., communication interface 135, data acquisition interface 140, display 125, external storage device 130).

[0069] The system example shown in Figure 1A is a non-limiting embodiment and should be understood as not being intended to limit the system to the components shown. For example, the system may include one or more other processors and memory devices. Furthermore, one or more components of the control processing hardware 100 may be provided as external components connected to the processing device. For example, as shown, the optional transmitting beamformer 26 and optional receiving beamformer 27 may be components of the control processing hardware 100 (shown by dotted lines) or may be provided as one or more external devices.

[0070] Some aspects of this disclosure can be implemented, at least in part, by software. This software is configured such that, when executed by a computer system, the computer system becomes a purpose-specific computer capable of implementing the signal processing and noise suppression methods (or variations thereof) described herein. That is, the technology can be implemented in a computer system or other data processing system as a response to a processor. The processor is, for example, a microprocessor, CPU, or GPU that executes a sequence of instructions in memory. The memory is, for example, ROM, volatile RAM, non-volatile memory, cache, magnetic disk, optical disk, cloud processor, or other remote storage device. Furthermore, these instructions can be downloaded to a computer device via a data network, for example, in the form of a compiled linked version. Alternatively, the logic for implementing the above process can be implemented on another computer and / or on a machine-readable medium. For example, discrete hardware components such as large-scale integrated circuits (LSIs) and purpose-specific integrated circuits (ASICs), or firmware such as electrically erasable programmable read-only memory (EEPROM) and field-programmable gate arrays (FPGAs).

[0071] Computer-readable media can store software and data that, when executed by a data processing system, perform various actions. Executable software and data can be stored in various locations, such as ROM, volatile RAM, non-volatile memory, and / or cache. This software and data portion can be stored in any of these storage devices. Generally, machine-readable media have mechanisms (i.e., storage and / or transmission) that provide information in a form accessible to machines (e.g., computers, network devices, personal digital assistants, manufacturing tools, any device with one or more processor sets, etc.).

[0072] Examples of computer-readable media include, but are not limited to, read-only memory (ROM), random access memory (RAM), flash memory devices, floppy disks and other removable disks, magnetic disk storage media, optical storage media (e.g., compact discs (CDs), digital general-purpose discs (DVDs), etc.), network storage, cloud storage, and other recordable and non-recordable media. Instructions can be implemented over digital and analog communication links, such as electrical, optical, acoustic and other forms of propagating signals, carrier waves, infrared signals, and digital signals. In this specification, the phrases “computer-readable material” and “computer-readable storage media” refer to all computer-readable media except for transient propagating signals.

[0073] Many of the embodiments described herein utilize the adjustment of a noise reduction filter based on environmentally measured noise. In the embodiments, one or more of the waveforms, data, filter parameters, and other relevant information described in the following embodiments can be transmitted to a network for remote evaluation and analysis and / or optimization of the noise reduction filter. After optimization, the noise reduction filter algorithm and / or parameters can be transmitted to the system to improve the noise reduction effect.

[0074] As shown in Figure 1A, the control processing hardware 100 includes an imaging processing module 145 and a noise suppression module 150. The image processing module 145 can be configured or programmed to perform known image processing methods (e.g., scan conversion).

[0075] While some of these embodiments are described in a manner that enables real-time noise reduction, it should be understood that noise reduction may also be performed as a post-processing step. For example, data on a transducer receiving channel or a reference receiving channel can be digitized and stored before or after the steps described in the embodiments of the present invention, such as filtering, envelope detection, shifting, shape / phase or delay adjustment, signal characteristic identification, attenuation, subtraction, and other steps.

[0076] Figure 1B shows an example of the steps by which the control processing hardware 100 and the receiving channel process the in-band imaging waveform detected from the imaging transducer receiving channel before image generation. The detected waveform from the imaging transducer receiving circuit is amplified (201) and filtered (202) before analog-to-digital conversion (203). After digitization, the detected waveform is filtered using a bandpass filter 200 (which may use multiple passbands and stopbands) to retain the signal within the imaging band. The envelope of the filtered waveform is generated by an envelope detector (210). The resulting envelope detection waveform is optionally decimated or expanded (220) and provided to the image processing module for image generation (230).

[0077] Refer to Figure 1A. The control processing hardware 100 includes one or more noise reduction modules 150. The noise reduction module 150 has commands to process the detected data (e.g., raw RF data, envelope data, image data) to remove the effects of noise according to the noise reduction algorithm described below. As described below, noise is removed or suppressed in the steps of the processing flow shown in Figure 1B, based on processing one or more noise detection waveforms using the noise suppression algorithm of this specification (shown as the noise suppression module 150 in Figure 1A). In the case of a system using an array transducer, noise suppression may be performed before or after beamforming (or both). In the embodiments described below, noise reduction of image data (including, but not limited to, raw waveforms, sampled waveforms, envelope waveforms, Fourier transformed signals, and processed image data) is performed based on signal energy measurement (power, amplitude, intensity, and other signal intensity measurement results) or waveform pattern measurement of the noise detection waveform (e.g., an out-of-band noise detection image waveform detected on a reference receiving channel). In many embodiments, the noise detection waveform has substantially no imaging energy and contains noise that is correlated with or coexists with the noise in the imaging band. One or more relationships between the imaging band noise and the noise detected by the noise detection waveform can be used to correct the imaging band signal and remove or reduce (e.g., suppress) the imaging band noise.

[0078] <Implementation of a reference receiving circuit for in-vivo imaging> Refer to Figure 1C, which shows an example of an ultrasound imaging system. It comprises an in-vivo imaging probe 350 connected to control processing hardware 100 via a patient interface module (PMI) 300. The in-vivo ultrasound imaging device can be configured to receive acoustic imaging energy from one-dimensional, two-dimensional, or three-dimensional regions. This is optionally done by mechanical or electronic scanning.

[0079] The imaging probe 350 comprises an imaging component 353 located away from the proximal end, an electrical and / or optical channel 354 passing through an optical conduit 354 along at least a portion of its length, and a connector 351 at the proximal end. For the purposes of this disclosure, the imaging component 353 is generally a component or set of components of the imaging probe 350 for detecting imaging energy (e.g., an acoustic or optical signal) to image a region adjacent to the imaging component. The imaging component optionally comprises one or more emitters for imaging energy and at least one receiver for imaging energy. For example, the imaging component comprises an ultrasonic imaging transducer 10 which is both an emitter and receiver for acoustic energy. The ultrasonic imaging transducer can be mounted on the imaging component. The imaging component can optionally be mounted or connected to a rotatable conduit or shaft (e.g., a torque cable) 352 housed within a hollow sheath of the internal ultrasound imaging probe to perform mechanical scanning.

[0080] The optional PIM300 enables signal transmission via a PIM cable 320 to a suitable imaging processing unit 100 within any wire or conduit, for example, when the imaging probe 350 is not directly connected to the control processing hardware 100. The PIM includes one or more amplifiers 20, which amplify signals from one or more transducer receiving channels. The PIM optionally includes a motor drive unit 301 that provides rotational motion to a rotatable conduit 354. The motor drive unit 301 includes slip rings, a rotary transducer, and other components that connect the probe 350's signal to the control processing hardware 100, allowing the imaging conduit to rotate while the PIM cable 320 does not. The PIM300 optionally includes a pullback mechanism 302 or a push-pull mechanism that allows the imaging component 353 to move longitudinally. The longitudinal movement of the imaging component occurs with the longitudinal movement of the external shaft surrounding the imaging conduit, or occurs within a relatively stable external shaft.

[0081] Many electrical components within an imaging system can pick up unwanted energy from environmental noise sources. Examples of such components include the imaging component 353, the imaging conduit 352, the motor drive unit 301, and the PIM cable 320. One or more reference receiving circuits that detect noise correlated with the noise detected by the imaging receiving circuit are useful for suppressing in-band noise. The following is an example of a reference receiving circuit for noise reduction of imaging signals in an ultrasonic imaging system.

[0082] Figure 1D shows an embodiment in which an imaging transducer is replicated to one or more non-imaging reference transducers 361 within the imaging probe 350. The reference transducer has its own electrical channel 360 passing through an optical conduit 352 and a connector 351. The reference transducer 361 can be coated with epoxy or other acoustically attenuating material 362, thereby acoustically isolating it from the received imaging energy. Alternatively, the piezoelectric element can be depolarized and inactivated, or replaced with an acoustically responsive substrate. This implementation can be extended to an array transducer with multiple ultrasonic transducer elements configured to receive acoustic imaging energy. One or more elements in the array can be made acoustically inresponsive, thereby not reflecting the acoustic imaging energy and functioning as a reference and transducer receiving circuit.

[0083] Figure 1E shows an embodiment in which the imaging probe comprises two or more imaging transducers, each having separate electrical channels 370 and 354. The two or more imaging transducers have sensitivity to receive acoustic imaging energy in substantially non-overlapping spectral bandwidths. For example, the first transducer is configured to receive acoustic energy around a frequency of 10 MHz, and the second transducer is configured to receive acoustic energy around a frequency of 40 MHz. The 40 MHz band of the 10 MHz transducer acts as a reference noise channel for the 40 MHz transducer, and similarly, the 10 MHz band of the 40 MHz transducer acts as a reference noise channel for the 10 MHz transducer.

[0084] Figure 1F shows an embodiment in which the imaging transducer receiving channel is replicated by a reference receiving circuit within the imaging probe. The reference receiving circuit may optionally comprise a combination of resistors, inductors, capacitors, and / or other components. These are configured such that the electrical impedance of the reference electrical circuit 381 matches the impedance of the imaging transducer receiving circuit, or so that the noise sensitivity of the reference electrical circuit is similar to that of the transducer receiving channel. The advantages of this embodiment are that it is low-cost, easy to manufacture, and allows for miniaturization of the components of the reference electrical circuit 381, as it does not require an actual ultrasonic transducer. Furthermore, parts of the reference receiving circuit can act for other purposes, such as transmitting energy to drive an actuator (e.g., a magnetic actuator), or carrying signals (including, but not limited to, temperature, pressure, or current generated from a magnetic field for position detection). Thus, by using the reference receiving circuit for purposes other than collecting a reference noise detection waveform to reduce noise in the imaging signal, miniaturization is facilitated, costs are reduced, and / or functionality is improved.

[0085] Figure 1G shows an embodiment in which the reference receiver circuit 391 is terminated within the PIM. Within the PIM, as with the imaging receiver channel, noise is received by the motor drive unit and the PIM cable. Only a portion of the transducer receiver circuit is replicated. Reference noise detection may be optionally used in conjunction with detection band (out-of-band or in-band) noise detection to further reduce noise in the in-band imaging waveform from the imaging transducer receiver channel.

[0086] Embodiments in which the detection band waveform is used as the noise detection waveform have lower manufacturing costs than embodiments in which the reference waveform is generated using a reference circuit or reference transducer. This is because the former does not require the physical implementation of a reference channel, for example, when the imaging probe or part thereof is not repeatedly used between subjects. Noise detection waveforms from a reference channel may be more effective in noise reduction imaging systems. This is because the waveform can provide information about noise within the imaging band, while out-of-band noise detection waveforms do not provide a direct estimate of in-band noise, but instead rely on noise that allows for the estimation of in-band characteristics based at least partially on out-of-band characteristics.

[0087] Noise estimation results obtained using a detection band waveform or a reference noise detection waveform can be used to reduce noise within one or more imaging channels. For example, in a phase array transducer with multiple electrostatic elements, the noise that can be collected from all or a subset of the electrostatic elements can be estimated using one reference receiving channel or one out-of-band noise detection waveform. Therefore, the same noise estimation method can be applied to signals collected from all or a subset of the electrostatic elements.

[0088] <Noise Measurement> The following sections of this disclosure describe several embodiments of performing denoising of in-band image data based on out-of-band noise detection, reference channel noise detection, or both out-of-band noise detection (which may also be supplemented by in-band noise detection) and reference channel noise detection.

[0089] As described in the following examples, the noise measurement results and / or noise characteristics can be determined from the measurement results of the noise detection waveform to increase the signal-to-noise ratio within the imaging band. Non-limiting examples of noise assumptions are as follows: • Energy measurement of noise detection waveforms and in-band imaging waveforms (amplitude, root mean square of amplitude, mean power); • Energy measurement in two or more noise detection bands; • Temporal, spectral, and / or time-frequency characteristics of the noise detection waveform; including spatial or spatial-temporal patterns, features, or parameters that represent or characterize patterns in images generated using the noise detection waveform; • The time characteristics, spectral characteristics, and / or time-frequency characteristics of the noise detection imaging waveform, and patterns simultaneously present within the in-band imaging waveform, including features representing these patterns; spatial or space-time patterns in the image generated using the noise detection waveform, and spatial patterns simultaneously present in the image generated from the in-band imaging waveform from the imaging transducer channel, including features or parameters representing or characterizing these patterns; Filter parameters determined or controlled using parameters acquired based on waveform characteristics (e.g., energy), and the spectral interval of harmonic peaks detected in the noise detection waveform.

[0090] In the embodiment, noise characteristics can be estimated when the imaging transducer is not receiving imaging energy (e.g., after it is assumed that the ultrasonic energy from the most recent emission of the ultrasonic pulse has dissipated from the image). Alternatively, noise characteristics can be estimated during imaging when it is assumed that imaging energy is being detected (e.g., when the transducer is in image acquisition mode). In embodiments where noise characteristics are measured when there is no imaging energy, the noise characteristics can be intermittently updated to adapt to and compensate for time-dependent changes in the noise characteristics.

[0091] <Noise reduction using measurements where imaging energy is absent> One embodiment of noise reduction involves measuring and using energy from a noise detection waveform while the imaging transducer is receiving imaging energy, while another embodiment uses measurement results from a noise detection waveform acquired during periods when imaging energy is absent, or uses a combination of these.

[0092] Figure 5A shows an embodiment of this system. Noise characteristic identification is performed while the transducer receiving channel is not receiving imaging energy, and the resulting noise characteristics are used to denoise the in-band imaging waveform acquired while the transducer receiving channel is not receiving imaging energy. The noise characteristic identification step is usually performed before acquiring and processing the denoised image data, but by appropriately recording the image data, the information collected in noise characteristic identification can be used retrospectively on the recorded image data.

[0093] According to the method of this embodiment, in a first time window in which it is assumed or known that at least one transducer receiving circuit is not receiving imaging energy, energy is detected in the transducer receiving circuit, and the waveform detected by the imaging transducer receiving channel 13 is considered to be noise 405. The detected waveform is filtered at 200 and 410, thereby generating an in-band noise characteristic identification waveform 407 and a detected band noise characteristic identification waveform 408.

[0094] The in-band noise characteristic identification waveform and the detected band noise characteristic identification waveform are processed to identify the noise characteristics (shown in 420). For example, using the noise characteristic identification 420, characteristic parameters 430 for characteristically identifying noise can be generated. Suitable examples of noise characteristic identification parameters are described in the following embodiments.

[0095] Optionally, if it is predicted or known that at least one transducer receiver circuit is not receiving imaging energy, and a specific noise source is known to be selectively OFF, and the waveform detected by the imaging transducer receiver channel 13 is considered the baseline of the selected noise source, then in another baseline noise characteristic identification stage, the transducer receiver circuit can detect the energy. By filtering the detected baseline noise characteristic identification waveform (200 and 410), an in-band baseline noise characteristic identification waveform 407 and a detected band baseline noise characteristic identification waveform are generated. It should be understood that the noise parameters 430 may include parameters obtained in the baseline noise characteristic identification stage.

[0096] Characteristic identification noise parameters can be calculated before or during an imaging session, or retrieved from a pre-stored database located on a local or remote storage device (such as a network drive or cloud).

[0097] After characterizing the noise when there is no imaging energy, the characteristic parameter 430 can be used to denoise the in-band imaging waveform 437 detected by the transducer while receiving imaging energy. The waveform detected from the imaging transducer receiving channel during imaging (435) contains imaging energy and noise, and by filtering this (200 and 410), the in-band imaging waveform 437 and the detected band imaging waveform 438 are generated. Therefore, the in-band imaging waveform 437 contains the detected imaging energy and noise, and the detected band imaging waveform 438 contains information about noise that may be present in the in-band imaging waveform 437. Using the characteristic parameter 430 acquired in the noise characteristic identification stage, noise in the in-band imaging waveform is detected and / or estimated (440), and noise suppression (500) of the in-band imaging waveform is performed. A suitable example of noise suppression using the noise characteristic identification parameter is described in the following embodiment.

[0098] Figure 5B shows another embodiment in which noise is detected using a reference receiving channel (for example, the reference receiving channel described in Figure 1A) to generate a reference noise characteristic identification waveform 406. The reference noise characteristic identification waveform is filtered (202) to generate a filtered reference noise characteristic identification waveform 409. In one embodiment, the reference channel filter is an imaging bandpass filter. Alternatively, if the noise estimation is based on an out-of-band noise input, the reference channel filter may be something other than an imaging bandpass filter. The reference noise characteristic identification waveform 409, and optionally an in-band noise characteristic identification waveform 407, are processed in step 420 to provide noise characteristic identification parameters 430.

[0099] During imaging, a reference waveform 436 is detected and optionally filtered to generate a filtered reference noise detection waveform 439. Using characteristic parameters 430 acquired in the noise characteristic identification stage, noise in the in-band imaging waveform 437 is detected and / or estimated (440), and noise suppression (500) of the in-band imaging waveform 437 is performed. In another embodiment, both the detected band imaging waveform 438 (shown in Figure 5A) and the filtered reference noise detection waveform 439 are processed to provide information about potential noise present in the in-band imaging waveform 437.

[0100] In the embodiment, noise suppression can be achieved by processing the in-band image waveform using one or more of the following methods: subtracting estimated noise from the signal within the imaging band; attenuating the estimated noise energy by multiplying the signal within the imaging band by an attenuation coefficient; or filtering the signal within the imaging band. For example, the subtraction amount is proportional to the amount of power detected in the out-of-band noise detection image waveform. Alternatively, the attenuation coefficient is inversely proportional to the measurement result corresponding to the amount of noise in the in-band image waveform, thereby attenuating the noise-related portion of the in-band image waveform.

[0101] Noise characteristic identification (e.g., as shown at 420 in Figures 5A-5B) may be performed once, or it may be performed multiple times consecutively. For example, noise characteristic identification may be performed intermittently (e.g., periodically or irregularly) to adapt to and compensate for time-dependent changes in noise characteristics.

[0102] When noise characteristic identification is performed, the noise characteristic identification waveforms collected for noise characteristic identification may be digitized and collected in multiple individual arrays (for example, arrays that are sufficiently long and can store image data along one scan line of ultrasound imaging), or they may be collected in a more continuous manner as one or more data streams to be stored in a large array, a circular buffer, or other data structure.

[0103] In the embodiment, noise characteristic identification may be initiated by the user, for example, at the start of the imaging session or when the user detects or suspects a deterioration in image quality (e.g., by pressing a button).

[0104] In another embodiment, noise characteristic identification may be automatically triggered or prompted by the user when the absence of imaging energy is detected. For example, the period of absence of imaging energy can be detected when the relative energy between the in-band imaging waveform and the noise-detected imaging waveform remains unchanged and within a specified range over a certain period of time. The period of absence of imaging energy can be detected, for example, when the energy of the in-band imaging waveform after noise correction (501) is below a specified threshold (which suggests the absence of imaging energy).

[0105] The noise characteristic identification step 420 can also be used to alert the user or the system if the noise profile has changed in a way that could adversely affect the noise suppression algorithm (for example, if a new noise source is detected, the noise suppression module 500 may incorrectly suppress in-band signals or reduce the effectiveness of in-band noise suppression). For example, the noise detection waveform (e.g., the out-of-band noise detection waveform 438 in Figure 5A or the reference noise detection waveform 439 in Figure 5B) can be processed to determine if the noise characteristics have changed. For example, when noise suppression is being performed, a noise monitoring module can optionally be used to monitor the characteristics of the noise detection waveform (e.g., peak energy, power, frequency content loop, skew, kurtosis, histogram, or other metrics). If the characteristics of the noise detection waveform change (e.g., if the peak energy exceeds a threshold), the noise monitoring module can communicate with other parts of the system (e.g., by message, interrupt, alarm, etc.) to alert that the noise content has changed.

[0106] In another embodiment, the error value in the noise characteristic identification stage can be evaluated. In this case, when imaging energy is absent, noise suppression is performed on the in-band imaging waveform, and the error value is the energy of the in-band imaging waveform after noise correction. If the error value exceeds a specified threshold, an alert is generated. The alert prompts the system to re-characterize the noise, or the system may ignore one or more out-of-band noise detection bands or one or more reference receiving channels in the noise reduction algorithm.

[0107] The noise characteristic identifier 420 can also be used to determine noise sources in the environment. Noise sources can be determined, for example, by a pattern recognizer as shown in step 570 of Figure 7A (detailed in Embodiment 4). Information about noise sources can be used, for example, to access a database (local or on a network) to select parameters for noise suppression or to determine the sequence of noise suppression methods to be used. For example, it is desirable to first remove periodic noise (detailed in Embodiment 6, Figures 8A to 8D), and then remove noise with low periodicity.

[0108] Alternatively, noise characteristic recognition can be used to detect the type of electroanatomical mapping system used in the ablation process, or to detect the activation and deactivation of an ablation catheter that uses radio frequency energy to perform ablation and treat arrhythmias. This can be achieved, for example, by a pattern recognition system as shown in step 570 of Figure 7A. One or more noise detection waveforms can further detect the duration, relative intensity, and frequency of the ablation energy being used. This information is useful for in-vivo imaging systems because it allows for tagging of imaging datasets with information about when the noise source (e.g., the ablation catheter) was activated.

[0109] Noise suppression (e.g., 500 in Figures 5A-5B) may be performed once, multiple times, intermittently, or continuously. For example, noise suppression may be initiated by the user. Alternatively, time-dependent noise sources may be compensated for by intermittently performing noise suppression (e.g., periodically or aperiodicly). As an alternative, noise suppression may be controlled by an external device that emits noise. For example, noise suppression may be initiated or stopped by the control of the RF ablation generator so that noise suppression is performed when RF energy propagates.

[0110] The following embodiments describe one imaging waveform and one noise detection waveform. It should be understood that these embodiments can be extended to multiple imaging waveforms and / or multiple noise detection waveforms.

[0111] <Embodiment 1: Noise suppression based on envelope detection out-of-band noise suppression (optional: amplitude, shape, and delay correction)> Refer to Figure 2A. An example of a method for suppressing noise in the in-band imaging waveform 437 via a suppression operator 525 using an out-of-band noise detection waveform 438 is shown. In one embodiment, the suppression operator is a subtractor that subtracts estimated noise from the in-band imaging waveform 437. In another embodiment, the suppression operator is an attenuator that attenuates the in-band imaging waveform by an attenuation coefficient derived from the estimated noise. Imaging energy is detected using one or more transducer receiving channels, and the detected energy includes both the imaging band and the noise detection band. By digitally sampling, separating (or copying), and filtering the waveform, the sampled in-band imaging waveform 437 and the sampled out-of-band noise detection imaging waveform 438 can be obtained. The sampled waveform can be detected as a sample set received within a time window. For example, in the case of ultrasound imaging, the listening window is immediately after or slightly after the ultrasound transducer pulses out, thereby allowing energy to be irradiated to the adjacent environment. The pulse output corresponds to one or more pulse transmissions.

[0112] In the embodiment shown in Figure 2A, the image waveform 435 is filtered (digitally or analogically) using an imaging bandpass filter 200 and a noise detection band filter 410 having frequencies extending outside the imaging band. Next, envelope detection is performed on the in-band image waveform 437 and the out-of-band noise detection image waveform 438. This is shown in 210 and 411.

[0113] In the embodiment shown in Figure 2A, the noise in the in-band imaging waveform 437 is suppressed using the out-of-band noise detection imaging waveform 438. Before the suppression process, the amplitude of the envelope-detected out-of-band noise detection imaging waveform 438 is optionally scaled (as shown in 510), thereby compensating for the difference between the noise power inside the imaging band and the noise detection band. In one embodiment, the amplitude adjustment coefficient is determined based on the power spectrum of the noise determined when there is no imaging energy (i.e., the noise characteristic identification stage). In another embodiment, the amplitude adjustment coefficient may be selected or changed by an operator to perform the desired level of noise suppression, or it may be determined after the cross-correlation 580 (described later).

[0114] Prior to the suppression process, the envelope-detected out-of-band imaging waveform may be temporally expanded, compressed, or reshaped using other linear or nonlinear time scaling functions (as shown in 510) to compensate for the difference in noise waveform shape between the imaging band and the noise detection band.

[0115] As shown in 510, it is also useful to apply delay correction to the out-of-band noise detection imaging waveform detected by envelope detection before noise suppression. For example, the two bandpass filters 200 and 410 do not transform the input waveform in the same way. The characteristics of the bandpass filter or the noise itself may result in noise offset because these propagate through the bandpass filter. Without delay correction, the noise shifts to produce an error before suppression, which adversely affects the noise suppression signal.

[0116] In one embodiment, delay adjustment can be performed by calculating the cross-correlation between the in-band and out-of-band imaged waveforms and aligning the waveforms at the point where the cross-correlation is maximum. In other words, the cross-correlation can be used to determine a time delay correction value to correct the relative time difference between the envelopes of the in-band and out-of-band imaged waveforms.

[0117] Time delay correction values ​​and amplitude correction values ​​can be calculated using multiple post-sampling in-band imaging waveforms and (simultaneously existing) post-sampling out-band noise detection imaging waveforms (referred to as an "array"), or within one or more time windows of the post-sampling in-band imaging waveforms and (simultaneously existing) post-sampling out-band noise detection imaging waveforms.

[0118] In one embodiment, out-of-band noise can be queryed as multiple noise detection bands, and the power dependence with respect to frequency between the multiple noise detection bands can be used to select an appropriate amplitude adjustment to estimate the noise power present in the imaging band for noise suppression in step 525. For example, the noise power in the imaging band can be estimated by fitting the average noise power across the multiple noise detection bands to a function dependence with respect to frequency (e.g., linear fitting). This function dependence with respect to frequency can be determined when there is no imaging energy (i.e., noise characteristic identification stage 420).

[0119] <Another embodiment: Noise suppression based on the use of frequency-shifted out-of-band noise detection waveforms (optional: amplitude, shape, and delay correction)> Figure 2B shows another embodiment of a subtractive or attenuated noise correction method in which the out-of-band noise detection imaging waveform is frequency-shifted before delay and amplitude adjustment. As shown in 530, a frequency shift operation (e.g., by multiplying by a complex exponent) is performed on the out-of-band noise detection imaging waveform to shift the spectrum of the out-of-band noise detection imaging waveform so that it is within the imaging band or overlaps with the imaging band.

[0120] The frequency shift calculation 530 can be performed such that the center frequency of the frequency-shifted noise detection waveform matches or is approximately equal to the center frequency of the imaging band. For example, if the imaging band is 7-13 MHz, the center frequency fc1 is 10 MHz. If the noise detection band is 15-25 MHz, the center frequency fc2 of the noise detection band is 20 MHz. Therefore, the frequency shift calculation can be performed such that the out-of-band noise detection imaging waveform is shifted by fc2-fc1=10 MHz. Alternatively, the frequency shift can be performed such that the center frequency of the frequency-shifted out-of-band noise detection imaging waveform matches or is approximately equal to the frequency within the imaging band (where some in-band noise is assumed to exist or is known).

[0121] After the frequency shift, a bandpass filter 203 in a separate stage is performed to filter out the effect of the sum frequency (in the above example, the sum frequency is fc2 + 10MHz = 30MHz). The frequency shift is useful in the envelope detection embodiment shown in Figure 2A. This is because the frequency shift results in a better correlation of noise between the imaging band and the noise detection band, thereby improving noise suppression. In the embodiment shown in Figure 2B, after noise suppression, the envelope detection 210 is applied to the output signal 501 to obtain a noise-suppressed signal envelope 520.

[0122] <Another embodiment based on the use of a reference noise signal> Refer to Figure 2C. An example of a method for suppressing noise in the in-band image waveform via a suppression operator (i.e., a subtractor or attenuator) using in-band noise detected via a reference receiving channel (using a reference receiving circuit) is shown. Imaging energy is received using one or more imaging transducer receiving channels, and noise energy (i.e., a reference waveform) is received using one or more reference receiving channels. This noise energy is expected to correlate with the noise energy received by the imaging transducer receiving channels.

[0123] The waveform is digitally sampled to obtain the sampled in-band image waveform and the sampled filtered reference noise detection waveform. Alternatively, noise suppression may be performed using analog electronic components. For example, the input from the reference receiving circuit to the analog signal adder is inverted by the delay, scale, and shape adjuster block 510, thereby subtracting the estimated noise. In another embodiment of analog signal suppression, the suppression process can be implemented as an amplifier with a time-varying gain. This gain is modulated by the noise received by the reference receiving channel.

[0124] In the embodiment shown in Figure 2C, the input waveforms from the imaging transducer receiving channel and the reference receiving channel are filtered (digitally or analogously) using the imaging bandpass filter 200 and an optional reference channel filter 202, thereby providing the filtered reference noise detection waveform along with the in-band imaging waveform. As described above, the reference channel filter 202 is the same as the imaging bandpass filter 200. Envelope detection is optionally performed on the filtered signal as shown in 210 and 411. In the embodiment of Figure 2C, the reference noise detection waveform (measured by the reference receiving channel) is used to suppress noise in the in-band imaging waveform. Before subtracting from the in-band imaging waveform (or its envelope) or attenuating the in-band imaging waveform (or its envelope), the amplitude of the filtered reference noise detection waveform (or its envelope) is optionally scaled by an amplitude adjustment coefficient (shown in 510), thereby compensating for the difference in noise power between the filtered reference noise detection waveform and the in-band imaging waveform. In one embodiment, the amplitude adjustment coefficient can be determined based on the power spectrum of the noise determined when imaging energy is absent (i.e., at the noise characteristic identification stage). In another embodiment, the amplitude adjustment coefficient can be selected by an operator to provide a desired level of noise suppression, or determined after cross-correlation 580.

[0125] As shown in 510, it is useful to apply delay correction to the envelope-detected and filtered reference noise detection waveform before subtraction. In one embodiment, delay adjustment can be achieved by calculating the cross-correlation between the in-band imaging waveform and the filtered reference noise detection waveform and aligning the waveforms at the point where the cross-correlation is maximum. Similar to the previous embodiment, the time delay correction value and amplitude correction value can be calculated using multiple in-band imaging waveforms and (simultaneously existing) reference noise detection waveforms, or within one or more windows of the in-band imaging waveforms and (simultaneously existing) reference noise detection waveforms.

[0126] <Embodiment 2: Noise suppression using a reference noise detection waveform from a reference receiving channel as input to an adaptive filter> Figure 3A shows an example of a noise correction method using an adaptive filter in a dynamic noise control (ANC) scheme, by applying noise suppression correction to the in-band imaging waveform based on a reference noise detection waveform. The reference noise detection waveform is correlated with the noise detected by the imaging transducer receiving circuit. The reference noise detection waveform is filtered (202), and the waveform from the imaging transducer receiving channel is filtered (200). In a preferred embodiment, the in-band imaging waveform and the reference noise detection waveform are filtered within the same band (e.g., 7-13 MHz in an example of an in-vivo echocardiography system).

[0127] An adaptive filter is a linear filter with a transfer function controlled by variable parameters, and is a means of adjusting its parameters according to an optimization algorithm. Adaptive filters are typically digital finite impulse response (FIR) filters or intermittent impulse response (IIR) filters. Dynamic noise control (ANC) is applied to the main input. The main input receives a signal (S) from a signal source that is corrupted by the presence of noise (N) that is not correlated with the signal. The reference input receives noise (Nr) that is not correlated with the signal but is somehow correlated with the main input noise (N). The reference noise passes through the adaptive filter to obtain output noise (Nr), which is an estimate of the main input noise (Nr). estimate) is generated. The noise estimation result is subtracted from the corrupted signal to produce a noise suppression signal (S estimate The adaptive filter dynamically adjusts its coefficients to generate an estimated result of the output power E[S]. estimate 2 Minimize the signal-to-noise ratio. Since the signal S is uncorrelated with N and Nr, and the noise N is correlated with the noise Nr, the signal-to-noise ratio can be maximized by minimizing the total output power. A minimization algorithm (e.g., stochastic least squares (LMS) or deterministic iterative least squares (RLS) algorithm) can be used to find the filter coefficients that minimize the output noise power.

[0128] In the example of the ANC method shown in Figure 3A, a reference receiving circuit isolated from the imaging energy is used to measure the reference noise via the reference receiving channel. The main input is obtained by applying the imaging bandpass filter 200 to the input waveform from the imaging transducer receiving channel 13. The reference input Nr is obtained by applying the reference channel filter 202 to the reference noise detection waveform obtained via the reference receiving channel.

[0129] <Another embodiment: Dynamic noise cancellation using out-of-band noise> Unlike the dynamic noise control method described above, which detects and processes a reference noise detection waveform and an in-band imaging waveform within overlapping frequency bands (and common frequency bands), Figures 3B and 3C show an embodiment of a noise correction method in which an adaptive filter 540 applies noise suppression correction to the in-band imaging waveform using an out-of-band noise detection imaging waveform.

[0130] The main input is obtained by applying an imaging bandpass filter 200 to the input waveform from the imaging transducer receiving channel and acquiring the envelope 210. In the embodiment shown in Figures 3B and 3C, the out-of-band noise detection imaging waveform is obtained by applying a detection band filter 410 to the input waveform from the imaging transducer receiving channel.

[0131] In Figure 3B, the out-of-band noise detection image waveform is demodulated via envelope detection 411 (using the same method as in the embodiment of Figure 2A), thereby obtaining the reference input for ANC.

[0132] In Figure 3C, the out-of-band imaging waveform is frequency-shifted (530) to the imaging band (e.g., 7-13 MHz) using the same method as in the embodiment of Figure 2B, and filtered using the imaging bandpass filter 203, thereby obtaining the reference input for ANC.

[0133] <Embodiment 3: Noise suppression based on frequency shift using the detected band waveform as the input to a variable filter of the in-band waveform> Figure 4 shows an embodiment of a noise correction method that filters an in-band image waveform using a dynamic filter 550. The dynamic filter is controlled by a filter update algorithm 560 that updates the filter coefficients after processing an out-of-band noise detection image waveform that includes out-of-band noise (optionally, an in-band noise detection image waveform that includes noise inside all or part of the imaging band). Similar to the embodiments described above (Figures 2B and 3B), as shown in Figure 4, the input waveform is filtered separately by an imaging bandpass filter 200 and a noise detection bandpass filter 410, thereby generating an in-band image waveform 437 and at least one out-of-band noise detection image waveform (including out-of-band noise 438). One or more out-of-band noise detection image waveforms are processed by the filter update algorithm 560.

[0134] The filter update algorithm analyzes out-of-band noise and evaluates signal characteristics. This is done, for example, by performing a Fourier transform on the waveform array to identify the maximum spectrum and the frequency at which it occurs. The filter update algorithm can use this information to control the coefficients of the dynamic digital filter (550) that filters the imaged waveform.

[0135] In one embodiment, the present method can be used to suppress noise in signals containing high-frequency noise. For example, the high-frequency noise originates from a switching rectangular pulse source where the spacing of spectral lines in the frequency domain depends on the pulse repetition frequency. A pulse width modulation source generates noise in the 3-40 MHz band, the imaging band is 7-13 MHz, and the dynamic filter 550 is a comb filter or a multi-notch filter. The filter update algorithm 560 processes the signal from the noise detection band (e.g., evaluates the spacing and position of spectral lines observed in 15-25 MHz) and uses this information to control the stopband of the dynamic in-band filter 550 to remove or suppress harmonic noise.

[0136] In addition to updating the dynamic filter using out-of-band noise, the filter update algorithm 560 and / or the dynamic filter 550 can optionally probe the in-band noise detection imaging waveform to confirm the presence of in-band noise in one or more selected subbands within the imaging band before denoising or suppressing it. For example, if the imaging band is 7–13 MHz and the filter update algorithm detects the presence of harmonic noise at 15 MHz, 18 MHz, 21 MHz, and 24 MHz (integer multiples of 15–25 MHz), the filter update algorithm sets the dynamic filter to filter the signals at 9 MHz and 12 MHz (integer multiples of 3 MHz within the imaging band). In one embodiment, such filters may optionally be applied only if the 9 MHz and 12 MHz signals are stronger than expected relative to the signals within the imaging band. In another embodiment, a noise characteristic identification step may be performed when there is no imaging energy to determine whether or not harmonic noise is present within the imaging band.

[0137] <Embodiment 4: Noise suppression based on pattern recognition> Figures 7A and 7B show an embodiment of a noise correction method that uses pattern recognition to detect noise and perform noise suppression on in-band image waveforms. In this embodiment, as shown in Figure 7A, a matched pattern set of detected band noise characteristic identification waveforms (at least one of which is at least partially outside the imaging band) and the associated in-band noise characteristic identification waveforms are first identified during a noise characteristic identification period in which no received imaging energy exists. After correlating the noise detection pattern with the in-band noise pattern via the noise characteristic identification stage, noise suppression of the in-band image waveform is performed on the in-band image waveform at imaging time using these correlations based on the pattern identification of one or more detected band image waveforms.

[0138] According to the first stage of this method, in the noise characteristic identification stage, when imaging energy is absent (for example, in the non-imaging noise characteristic identification stage), energy is detected within both the imaging band (in-band) and the noise detection band, thereby obtaining correlated measurement results for the in-band and detection band noise characteristic identification waveforms. Samples from the in-band noise characteristic identification waveform and the detection band noise characteristic identification waveform are recorded as an array pair. An array pair is a simultaneously recorded post-sampling in-band waveform and a corresponding post-sampling detection band noise detection waveform.

[0139] The detected band noise characteristic identification array and the intraband noise characteristic identification array may be windowed (566 and 565). The window may be, for example, a sliding window and may optionally overlap. Optionally, the window may be centered around the peak noise amplitude or time-fixed with respect to the noise amplitude threshold. The window may be tuned to suppress any adverse effects caused by windowing. This can be done, for example, by applying a Hamming window, a Blackman window, or other window functions well known in the field of signal processing. The array data (or its window) can be processed to identify the presence of one or more noise patterns (570).

[0140] Refer to Figure 7A. One or more detection band noise characteristic recognition arrays are processed to identify waveform patterns corresponding to the patterns of in-band noise characteristic recognition waveforms. The pattern recognizer 570 extracts features from the detection band characteristic recognition arrays and classifies these features into noise "classifications" using a predictive model. The extracted features are, for example, statistical features (including, but not limited to, variance, standard deviation, power, skewness, and kurtosis) in the time domain, frequency domain (e.g., peak frequency), and time-frequency domain (e.g., wavelet coefficients). The selection of features to be extracted can be done in advance using feature selection algorithms such as variable augmentation or backward selection.

[0141] The extracted features are fed into the predictive model in step 570. The predictive model can be trained to identify patterns in the detected band waveform and assign patterns to noise classification. For example, a machine learning method can be used to train the predictive model and use the extracted features to recognize patterns in the detected band noise characteristic discrimination array. The predictive model can be an unsupervised learning model (e.g., k-means clustering) or a supervised learning model (e.g., a linear classifier, an artificial neural network, or a nearest neighbor classifier). If prior information about the noise source is known, a learning model can be used. For example, the noise source and the sequence of noise patterns may be known in advance, and noise classification labels can be assigned to the waveform patterns of the detected band noise characteristic discrimination waveform. The predictive model in step 570 can accept input classification weights or prior probabilities. Higher prior probabilities or classification weights result in higher recognizability.

[0142] The database in step 575 can store intraband noise characteristic waveform patterns that are known to coexist with the detected band noise characteristic waveform pattern. For example, the database can store samples or averages of intraband noise characteristic waveform patterns paired with the coexisting detected band noise characteristic waveform patterns and noise classification labels. The detected band noise characteristic waveform pattern and its coexisting intraband noise characteristic waveform patterns can be determined, for example, for each window. Further noise characteristic recognition may be performed in the time domain, or in the spatial domain after image generation in steps 230 and 231. In this case, spatial features are also extracted in step 570. The database can be in any suitable format used by the computer, for example, a lookup table.

[0143] By correlating the detected band noise pattern with the in-band waveform noise pattern via the noise characteristic identification stage described above, noise suppression of the in-band imaged waveform can be performed during imaging based on pattern identification of one or more detected band imaged waveforms using these correlations.

[0144] Refer to Figure 7B. An in-band image waveform is obtained by applying an imaging bandpass filter 200 to the waveform detected by the imaging transducer receiving channel while the transducer is receiving imaging energy. Envelope detection 210 is performed as an option. A detected band image waveform is obtained by applying a noise detection bandpass filter 410 to the waveform detected by the imaging transducer receiving channel. Envelope detection 411 is performed as an option. The detected band image waveform and the in-band image waveform can be windowed (566 and 565) in the same way as the windowing step of the noise characteristic identification stage. The waveforms can be sampled and presented as an array.

[0145] Similar to the feature extraction step in the noise characteristic recognition stage, features can be extracted from the detection band imaging array (or its window). The pattern recognizer 570 (described above), trained in the noise characteristic recognition stage, uses the extracted features and, optionally, classification weights to identify the presence of one or more patterns in the detection band imaging waveform.

[0146] The pattern classification algorithm 570 can be applied, and the prior probability of the pattern classification (e.g., a Bayesian classifier) ​​or the weight of the classification (e.g., that of a support vector machine) can be adjusted using the period (i.e., repetition frequency) of the noise pattern. A higher prior probability or classification weight increases the likelihood of recognition. If the repetition interval of the pattern is known, the pattern is expected to exist with a high probability at a given time. The prior probability or classification weight can be adjusted to be high at that time to increase the likelihood that the pattern classifier will recognize the noise pattern. This repetition interval can be determined in the noise characteristic identification stage and stored in the database 578, or read from a pre-stored database (local, network, or cloud storage).

[0147] Noise corrections are generated for intraband imaging waveforms (e.g., intraband imaging arrays) using the detection band imaging waveform patterns identified in step 570 as corresponding to one or more noise classifications. These corrections can be generated based on the discovery of correlated intraband patterns in step 575, and the detection band noise characteristic identification waveform patterns, the matching set of features between the intraband noise characteristic identification waveform patterns and noise classification labels, and the noise classification labels are stored in a searchable database or other classification mechanism.

[0148] In the example method shown in Figure 7B, intraband noise correction is generated for each window and subtracted from the intraband image waveform 525 for each window. Optionally, this is done after delay and / or amplitude and / or shape adjustment 510, which temporally align the intraband noise patterns obtained along with the intraband waveforms from the database 575 or other classification mechanisms.

[0149] In one embodiment, during the noise characteristic identification stage (where no received imaging energy exists), the detected band noise characteristic identification array is first processed to extract and store one or more features. The time interval for detecting a given noise pattern can also be determined and stored in this stage (578). The corresponding correlation time patterns of the in-band noise characteristic identification array are also stored (575) (for example, stored in a lookup table).

[0150] In this example, in the imaging stage where denoising is performed, the detection band imaging array of one or more detection band waveforms is largely noisy and is processed through the same feature extraction process. A weight vector (assigning a weight or prior probability to each noise classification pattern) may be optionally obtained. The repetition frequency of each pattern can be read from database 578 generated in the noise characteristic identification stage or from a pre-stored database. The weight of each classification can be dynamically adjusted to depend on the repetition frequency of the pattern, the time when the pattern was previously detected, and the accuracy of the previous detection. Features (and optionally classification weights) extracted from the detection band waveforms are supplied to a trained prediction model (trained in the noise characteristic identification stage). The prediction model identifies noise patterns in the detection band waveforms and assigns classifications to them. The corresponding in-band noise patterns for the noise classifications are obtained from database 575. Database 575 stores, for example, features of in-band noise waveform patterns and detection band waveform patterns for each noise classification in the noise characteristic identification stage (e.g., in-band time waveforms stored in a lookup table). The intraband noise waveform pattern extracted from the classification comparison is, for example, the average of all concurrently existing intraband noise patterns in the current noise classification, or the intraband noise pattern (determined, for example, by nearest neighbor calculation) whose concurrently existing detection band pattern features are closest to the features of the current detection band imaging array. This concurrently existing pattern is subtracted from the input after adjusting for amplitude and delay, thereby obtaining a noise-suppressed intraband imaging waveform.

[0151] Figures 7A and 7B show an embodiment in which a single detection bandpass filter 410 leads the input to the pattern recognizer 570 to generate a single detection band waveform (in which at least some of the energy is outside the imaging band). Alternatively, multiple detection bandpass filters (in which at least one detection band waveform is outside the band) may be used to generate multiple detection band waveforms.

[0152] In addition to at least one detected band waveform that conveys out-of-band noise, one or more detected band waveforms may also convey noise in all or part of the imaging band (i.e., in-band noise detection waveforms). Such in-band data is useful when the pattern recognizer 570 confirms that the noise predicted by the out-of-band noise detection waveform actually exists within the imaging band (at the noise characteristic identification stage or during imaging).

[0153] For example, the pattern recognizer can identify the noise source using the in-band noise detection imaging waveform. For example, if a sub-band of the imaging band is substantially different from one or more other sub-bands of the imaging band, or if it is substantially different from the net energy within the imaging band, the noise source corresponding to the imaging sub-band can be identified. For example, using a center 8 MHz peak filter, an in-band noise detection waveform within an imaging band of 7 - 13 MHz can be obtained, and using a detection band pass filter with a pass band of 15 - 25 MHz, an out-of-band noise detection waveform can be obtained. If the 8 MHz in-band noise detection waveform detects an increased energy relative to the energy of the 15 - 25 MHz out-of-band noise detection waveform, the pattern recognizer can adjust its weight to detect a specific noise source (i.e., noise classification). If relying only on information outside the imaging band, the system can select a correlated in-band noise pattern from the database 573 and remove its 8 MHz peak.

[0154] The above embodiments have been described in the context of detecting the time pattern of the raw signal or the envelope-processed signal, but it should be understood that the above algorithms are also applicable to embodiments using image data. For example, instead of processing image data (e.g., B-mode image data) to process a time-domain (e.g., RF or envelope detection) signal, a spatial noise pattern can be determined. These alternative embodiments are shown by the dashed lines in Figure 7B, and decimation (220 and 221) and B-mode image line generation (230 and 231) are performed before 570 and 575.

[0155] Instead, when processing an image in the spatial domain, a spatial pattern can be detected using a 2D image window. For example, it is B-mode image data. In B-mode data, texture features are extracted in the spatial domain (e.g., gray level co-occurrence matrix), frequency domain (e.g., Fourier spectrum measurement), or spatial frequency domain (e.g., energy of 2D wavelet coefficients).

[0156] Referring to FIG. 7D, another embodiment is shown in which a reference noise detection waveform is used instead of the detection band imaging waveform of FIG. 7B when performing pattern recognition during imaging. Similarly, in FIG. 7C, a reference noise characteristic identification waveform is used in the initial pattern recognition stage performed when no imaging signal is present. Therefore, the above algorithms or methods shown in FIGS. 7A and 7B are adjusted as in this embodiment by replacing the detection band noise characteristic identification waveform (and corresponding array measurement) with the reference noise characteristic identification waveform (shown in FIG. 7C) and replacing the detection band imaging waveform (and corresponding array measurement) with the reference noise detection waveform (shown in FIG. 7D).

[0157] <Embodiment 5: Noise Suppression Based on Relative Energy Measurement> In this embodiment, noise suppression is performed by selectively attenuating the window portion of the in-band imaging waveform according to conditions evaluated based on measurement results from one or more detection band imaging waveforms (at least one of which is an out-of-band imaging waveform). The attenuation is, for example, subtracting the derived subtraction value from the envelope of the windowed in-band imaging waveform and / or multiplying the windowed in-band imaging waveform or its envelope by an attenuation coefficient. The subtraction value and / or attenuation coefficient are determined from the measurement results of the noise detection imaging waveform or the reference noise detection waveform.

[0158] In the first stage of this example, in the noise characteristic recognition stage, when imaging energy is absent (for example, in the non-imaging noise characteristic recognition stage), energy is detected in both the imaging band (in-band) and the noise detection band, thereby obtaining the correlation measurement result between in-band noise and detected band noise. At least one noise detection band is out-of-band. Samples of the in-band noise characteristic recognition waveform and the detected band noise characteristic recognition waveform are recorded as array pairs. An array pair is a simultaneously recorded post-sampling in-band waveform and its corresponding post-sampling detected band waveform.

[0159] Optionally, the noise characteristic identification stage may include another stage. This is called the baseline noise characteristic identification stage and is performed when it is known or predictable that the imaging transducer receiving channel is not receiving imaging energy and noise energy. As shown in Figure 6A, when imaging energy and noise energy are not being received, the in-band baseline noise characteristic identification array can be obtained by applying the imaging bandpass filter 200 to the input waveform detected from the imaging transducer receiving channel 13 (optionally by detecting the envelope of the filtered waveform (210)). The detected band baseline noise characteristic identification array can be obtained by applying the detected band filter 410 to the input waveform from the imaging transducer receiving channel 13 (optionally by detecting the envelope of the filtered waveform (411)). The in-band baseline noise characteristic identification array and the detected band baseline noise characteristic identification array are measured when imaging energy and noise energy are absent and are represented as Gi and Gn, respectively. The given array pair can optionally be segmented according to multiple time windows, as shown in 565 and 566 of Figure 6A, to obtain a window array pair (Gi w and Gn w) can be obtained. The window is, for example, a sliding window and optionally overlaps with adjacent windows. In the baseline noise characteristic identification stage, one or more noise measurement results can be calculated from the energy measurement results for each window. For example, the maximum power in the windowed band baseline noise characteristic identification array is shown as Ti. Similarly, the maximum power in the windowed band baseline noise characteristic identification array is shown as Tn.

[0160] The noise characteristic identification stage includes a stage in which the imaging transducer receiving circuit is not receiving imaging energy but is assumed to be receiving noise energy. See Figure 6A. When imaging energy is not received, the in-band noise characteristic identification array 407 can be obtained by applying the imaging bandpass filter 200 to the input waveform from the imaging transducer receiving channel 13 (optionally by detecting the envelope of the filtered waveform (210)). The detected band noise characteristic identification array can be obtained by applying the noise detection band filter 410 to the input waveform from the imaging transducer receiving channel 13 (optionally by detecting the envelope of the filtered waveform (411)). The in-band noise characteristic identification array 407 and the detected band noise characteristic identification array 408 are measured when imaging energy is absent and are shown as Ci and Cn, respectively (as shown in Figure 6A).

[0161] Optionally, a given array pair can be segmented according to multiple time windows, as shown in 565 and 566 of Figure 6A, to create a windowed array pair (Ci w and Cn wIt can be obtained (shown as). The window is, for example, a sliding window, and optionally overlaps with adjacent windows. In another example, the window has a center frequency around the peak amplitude of the noise waveform of one or more detection bands and / or imaging bands. In another example, the time window is fixed in time to the occurrence of noise determined when the noise amplitude of one or more detection bands and / or imaging bands exceeds a specified threshold. For example, the threshold is proportional to the parameters Tn and / or Ti obtained in the baseline noise characteristic identification stage.

[0162] Within the band and the detection band noise characteristic identification array pair Ci w and Cn w are processed to obtain one or more measurement results corresponding to the energies of the imaging band and the noise detection band of each time window, thereby characterizing the relative intensity of the noise in the two bands. For example, as shown in FIG. 6A, for each pair of windows Ci w and Cn w , the power of the imaging band and the power of the noise detection band are calculated (570 and 572).

[0163] The noise characteristic identification measurement results may optionally be calculated only from the selected windows. The selection conditions can be evaluated according to the in-band energy measurement results and optionally the detection band energy measurement results. For example, only the windows with in-band power exceeding a specified threshold are selected to obtain the noise characteristic identification measurement results (shown in FIG. 6B). The threshold is proportional to Ti obtained in the baseline noise characteristic identification stage. In another example, only the windows with out-of-band power exceeding a specified threshold are selected to obtain the noise characteristic identification measurement results. The threshold may be derived from Tn obtained in the baseline noise characteristic identification stage.

[0164] In the above example, it should be understood that the maximum and minimum values ​​refer to the upper and lower percentiles, or the true maximum and minimum values. For example, the 98th and 2nd percentiles can be used instead of the maximum and minimum values. Other statistical thresholds (e.g., the 95th and 5th percentiles, the 90th and 10th percentiles, the 80th and 20th percentiles, etc.) can also be used to represent maximum and minimum values ​​for characteristic identification.

[0165] From the energy measurement results for each window in the noise characteristic identification stage, one or more noise measurement results can be calculated and used to define the relationship between the power of the imaging band and the power of the detection band when a noise source is present. For example, one or more pairs of in-band and detection band power values ​​can be selected as inflection points for generating a piecewise linear function to define the relationship between the power of the imaging band and the power of the detection band when a noise source is present (shown in Figure 6B).

[0166] In one example, a piecewise linear function defining the relationship between in-band power and detected band power in the presence of a noise source can be generated based on the maximum and minimum power at the noise characteristic identification stage (shown in Figure 6B). The maximum and minimum in-band power values ​​from the noise characteristic identification stage are evaluated (e.g., using absolute maximum / minimum values ​​or statistical measurement results), and Pi min and Pi max This is shown as follows. In one embodiment, the in-band power is Pi min For this, identify window sets that are within a pre-selected range (e.g., percentile range), and from the identified window sets, select the minimum detected band power Pn min Select as such. Similarly, if the in-band power is Pi max For this, identify the window set that is within the pre-selected range, and from the identified window set, select the minimum detected band power Pn max Select as: Power pair (Pn min Pi min ) and (Pn max Pi maxUsing [[ID=]], a function that defines the estimated relationship between the in-band power and the detected band power can be fitted. The examples are only non-limiting examples for selecting the in-band power and the detected band power to provide appropriate fitting points and / or for providing the functional relationship between the in-band power and the detected band power, and it should be understood that other methods can be used.

[0167] Optionally, the ratio of the power of the in-band noise characteristic identification waveform to the power of the detected band noise characteristic identification waveform can be calculated for each window, and the maximum ratio R between multiple windows off (An example of the use of this quantity is when determining whether to apply noise correction during imaging, which will be described below) can be obtained.

[0168] One or more relationships f(pn) between the power of the in-band noise characteristic identification waveform (Pi) and the power of the detected band noise characteristic identification waveform (Pn) can be obtained. For example, as shown in FIG. 6B, f(Pn) is a piecewise linear function whose slope, intercept, and / or inflection point are calculated at the noise characteristic identification stage (Pn min , Pi min ) and (Pn max , Pi max ). In another example, f(Pn) is a non-linear polynomial or a combination of a linear polynomial and one or more non-linear polynomials. In another example, f() is a set of values defined for a range of one or more Pn values. For example, f(Pn) is Pia when a1 ≤ Pn < a2, Pib when b1 ≤ Pn < b2, and so on. [a1, a2] and [b1, b2] are non-overlapping intervals of Pn. A set of windows of the detected band noise characteristic identification array whose power is between a1 and a2 can be identified, and Pia in the interval [a1, a2] can be determined using the corresponding windows of the in-band noise characteristic identification array. For example, Pia is a representative power value (e.g., maximum, minimum, average, median, or other statistical value) calculated from the in-band power measurement results of all windows where the detected band power Pn is within the interval [a, a2].

[0169] Noise measurement results obtained in the noise characteristic identification stage can be used to suppress noise in the in-band imaging waveform acquired during imaging (shown in Figures 6C and 6D). The input waveform is filtered (200 and 410), and optionally envelope detection is performed (210 and 411) to provide the in-band imaging waveform 407 and the detected band imaging waveform 408. The waveform is sampled to obtain the in-band imaging array and the detected band imaging array.

[0170] In this embodiment, the term "in-band imaging array" refers to the sampled in-band imaging waveform. The term "detection band imaging array" refers to the sampled detection band imaging waveform. The array set is recorded, and each array corresponds to a given scan line. For example, the first in-band imaging array corresponds to the first scan line, the second in-band imaging array corresponds to the second scan line, and so on. The in-band and detection band imaging arrays are each Qi θ and Qn θ This is shown as follows: θ is an index that identifies a given acquisition period (e.g., corresponding to one scanline).

[0171] As shown in Figures 6C, 565 and 566, the array can be windowed using the same window characteristics used when windowing the noise characteristic recognition array in the characteristic recognition stage. The in-band imaging array and the detection band imaging array are segmented temporally according to the window, and each Qi θ,w and Qn θ,w This is shown as follows. The sign w is an integer indicating the window number. For example, Qi 1,10 This refers to the tenth window portion of the in-band imaging array corresponding to the first scan line. For each window portion of the in-band imaging array and the detection band imaging array, calculate the power or other appropriate energy measurement results (570 and 572). Calculate the in-band and detection band power values ​​for each window, and P(Qi) for each. θ,w ) and P(Qn θ,wThese are shown as follows. Using these energy measurement results, noise suppression is performed based on the measurement results obtained in the noise characteristic identification stage.

[0172] In one embodiment, as shown in Figure 6C, the in-band imaging array Qi θ,w Noise is suppressed by subtracting the power value (525) from the envelope. Noise detection band power P(Qn θ,w Using ) based on the function f() obtained in the noise characteristic recognition stage, the in-band imaging array window Qi θ,w Noise energy P^i Nθ,w It is possible to estimate this. For example, Qi θ,w Internal noise energy P^i Nθ,w is, P^i Nθ,w =f(P(Qn θ,w It can be estimated as (575). Estimated intraband noise P^i Nθ,w Optionally, it can be scaled by multiplying by a scaling factor β (0 ≤ β ≤ 1). Scaled estimated noise βP^i Nθ,w As an option, clamp it to below the upper limit (for example, 0.8 × P(Qi θ,w ), 0.9 × P(Qi θ,w ), 1.0 × P(Qi θ,w (However, this is not limited to) the subtraction value can be obtained. Qi θ,w The noise-suppressed waveform envelope is obtained by subtracting a subtraction value (525) from the envelope.

[0173] In another example shown in Figure 6D, the in-band imaging array Qi θ,w Noise can be suppressed by multiplying the element by an attenuation coefficient (526). During imaging, the imaging window Qi θ,w The power within the band is P(Qi θ,w It can be calculated as: Noise detection band power P(Qn θ,w Using ) based on the function f() obtained in the noise characteristic recognition stage, the in-band imaging window Qi θ,w Internal noise energy P^i Nθ,w It is possible to estimate this. For example, Qi θ,w Internal noise energy P^iNθ,w is, P^i Nθ,w =f(P(Qn θ,w It can be estimated as (575). Estimated intraband noise P^i Nθ,w It can be optionally scaled by multiplying by a scaling factor β (0 ≤ β ≤ 1). The damping factor is [P(Qi θ,w )-βP^i Nθ,w ] / P(Qi θ,w It can be selected in proportion to ).

[0174] In the embodiment, the scaling coefficient β can be selected within the range of 0 to 1. The attenuation of ultrasonic energy decreases the imaging energy over time, and the determination of β depends on the window depth in the waveform (and therefore corresponds to the depth in the imaged tissue). The parameter β may optionally be user-controlled. The attenuation coefficient can optionally be clamped to an upper limit (e.g., 0.95, 0.9, 0.8, etc., but not limited to these). In addition to or instead of this, the attenuation coefficient can be clamped to a lower limit (e.g., 0, 0.01, 0.05, 0.1, etc., but not limited to these). The attenuation coefficient is set to Array Qi θ,w Alternatively, by multiplying it by its envelope, a noise-suppressed array can be obtained.

[0175] <Pattern recognition for suppressing classification-specific noise> In the noise characteristic identification stage (for example, Figure 6A), the system can optionally be configured to group array pairs Cnw and Ciw into one or more categories (referred to as classifications). Refer again to the pattern recognizer in step 570 of Figure 7A. One or more detected band noise characteristic identification waveforms can be processed by the pattern recognizer to identify one or more noise pattern classifications and assign classifications to time windows corresponding to the identified noise patterns. A window set of detected band noise characteristic identification waveforms belonging to classification k is selected and shown as {w_k}. Using the in-band and out-of-band power measurement results of the windows in set {w_k}, a functional relation f specific to classification k can be determined in the same manner as above. k () can be derived.

[0176] The system can be configured to identify a noise pattern using a pattern recognizer when the imaging transducer circuit is receiving imaging energy. Using a method similar to that described in step 570 of FIG. 7B, one or more detected band imaging waveforms can be processed by a pattern recognizer to identify one or more noise pattern classifications and assign a classification to a time window corresponding to the identified noise pattern. The in-band noise power can be derived from a function relationship f() specific to the identified classification. For example, for a window belonging to classification 1, the in-band noise (i.e., P^i Nθ,w =f1(P(Qn θ,w ))) can be derived. Similarly, for a window belonging to classification 2, the in-band noise can be derived using the function f2(). f1(), f2(), f3(), etc. are obtained in the noise characteristic identification stage. As described above, the noise within the window Qi θ,w of the in-band imaging waveform can be suppressed by subtracting a subtraction value from the envelope of Qi θ,w or by multiplying by an attenuation factor and Qi θ,w . The subtraction value or attenuation factor is derived from P^i Nθ,w . As a specific example, the noise classified as classification 1 by the pattern recognizer is from an electroanatomic mapping system, and the noise classified as classification 2 is from an ablation energy generator. Therefore, when the system recognizes that the noise has arrived from the mapping system, based on the power of the noise detection band, the in-band noise generated by the electroanatomic mapping system can be estimated using f1(). Also, when the system recognizes that the noise has arrived from the ablation generator, based on the power of the noise detection band, the in-band noise generated by the ablation generator can be estimated using f2().

[0177] In an alternative example, a noise-emitting system may be monitored to determine the functional relationships used for noise estimation and suppression. For instance, the control of an ablation generator may be monitored to enable a binary gate signal when the ablation generator is generating energy (and noise). This gate signal can then be used to determine the time period over which the function f2() should be used for in-band noise estimation.

[0178] <Selective implementation of noise suppression> The energy measurement results can be optionally used to estimate whether the in-band imaging array has a low signal-to-noise ratio. In other words, the energy measurement results can be used to classify whether or not the window should be noise-suppressed (via subtraction or multiplication of the attenuation coefficient).

[0179] In one embodiment, the decision of whether or not to apply noise suppression correction to window w can be made based on the power detected in the noise detection band. P(Qn θ,w If the noise exceeds the specified threshold, Qi noise correction will be disabled. θ,w This is applied to [the specified value]. The threshold can be obtained in the noise characteristic identification stage.

[0180] In another example, the decision of whether or not to apply noise suppression correction to a given window can be made based on the ratio of the imaging band power to the noise detection band power. For example, in the noise characteristic identification stage, the ratio of the imaging band power to the detection band power is obtained for each window of the pair of in-band noise characteristic identification arrays and detection band noise characteristic identification arrays. The representative maximum ratio (R) across all windows is obtained. Off This value can be used as a threshold to determine whether or not to apply noise suppression correction when the imaging transducer receiving circuit is receiving imaging energy, as described above.

[0181] The decision of whether or not to apply noise suppression correction to the window w of the in-band imaging waveform is based on the ratio (R) of the imaging band power to the noise detection band power.On w This can be implemented based on ). In one example, R On w γR Off Compare with γ, where R is the relaxation parameter. Off This is the representative maximum ratio calculated in the noise characteristic discrimination stage. On w γR Off If R is greater than this, the signal-to-noise ratio of the in-band signal is estimated to be sufficiently large, and noise suppression correction is not applied. Conversely, R On w γR Off The signal-to-noise ratio is estimated to be sufficiently small if the following conditions are met, and noise suppression correction is applied.

[0182] The value of γ can be used to adjust the sensitivity of the signal-to-noise ratio, and in some cases, the value of γ can be used as an adjustment coefficient when weak parts of the signal are suppressed. Lowering the value of γ lowers the threshold for applying noise suppression, thereby reducing the number of windows for noise suppression and allowing more imaging energy (and noise) to be retained in the final output.

[0183] In the embodiment, γ can be selected between zero and a unit value. In ultrasound, the imaging energy decreases over time due to the attenuation of ultrasonic energy, and the determination of γ depends on the window depth in the waveform (i.e., the depth in the tissue being imaged). The parameter γ may be optionally user-controlled. For example, if tissue or other structures in the image are unwanted or excessively attenuated, the user can reduce this parameter value to weaken the noise suppression effect.

[0184] In a window identified for noise suppression, noise in the in-band imaging array within the window can be suppressed using any appropriate denoising or noise suppression method. It should be understood that various noise suppression corrections can be applied. Examples include, but are not limited to, subtraction correction shown in Figure 6C and / or multiplication by the attenuation coefficient shown in Figure 6D.

[0185] Noise suppression can lead to incorrect noise reduction for each window. For example, windows with low imaging energy may be negatively affected by noise suppression based on incorrectly determining a low signal-to-noise ratio within the window (i.e., incorrect window classification). This can result in small image "holes" in the peripheral uniform areas of the image, or residual noise pixels in the peripheral low-noise areas of the image. In one embodiment, it is possible to determine whether a window determined to be suitable for noise suppression should actually have the process applied to it using the state of adjacent windows (i.e., spatial adjacency). If a given window is determined to be suitable for noise suppression by the above method, it is possible to evaluate whether noise suppression should be performed using the attenuation coefficient of the given window, using adjacent windows in one or more adjacent arrays (i.e., arrays corresponding to adjacent scan lines).

[0186] For example, if a given scanline is determined to be unsuitable for noise suppression, and one or more adjacent windows are determined to be suitable for noise suppression, the window can be flagged as potentially misclassified. The state of the given window is overwriteable, and the given window is identified as suitable for noise suppression, and noise suppression is applied to that window. Conversely, if a given window is determined to be suitable for noise suppression, and adjacent windows in the adjacent array are determined to be unsuitable for noise suppression, the initial determination of the state of the given window is overwritten, and the given window is flagged as unsuitable for noise suppression and potentially misclassified. Samples from misclassified windows can be replaced by samples from one or more noise-free adjacent windows (i.e., replaced or interpolated by samples from adjacent windows determined to be unsuitable for noise suppression according to the method described above). This is optionally done after delay and amplitude adjustments have been performed.

[0187] Figure 6E shows an example implementation of this method. The adjacent windows within the adjacent array are examined to determine whether the classification of the current window is consistent with the surrounding windows within the adjacent array. Two windows (array 2, windows 4 and 5) classified as mainly containing noise (marked "N") are surrounded by a window (marked "S") classified as having a sufficiently high signal-to-noise ratio. As shown in Figure 6F, the two "N" windows are replaced with samples from the adjacent "S" windows (e.g., by copying or interpolation by amplitude adjustment and / or shape adjustment), and noise suppression by noise estimation is not performed on these windows. In this embodiment, two adjacent windows on either side of the array are used to check the state of a given window, but any number of windows may be used in other embodiments.

[0188] When evaluating whether or not to change the classification of a given window, one or more temporally adjacent windows before and after the given window can also be used. For example, in Figure 6E, window 4 of array 5 is initially marked as "S," and in Figure 6F, it is reclassified as "N." This is because the preceding and succeeding windows of the array are marked as "N."

[0189] Although this embodiment performs signal processing before image processing, this embodiment can be applied to image data processing rather than time-domain signal processing. For example, multiple in-band image arrays and detected band image arrays representing multiple adjacent scan lines can be acquired, and frames of the in-band image and detected band image can be obtained by post-processing. The in-band and detected band image pixels are respectively Bi θ,d and Bn θ,d This is expressed as follows: θ represents the scanline and d represents the depth. Using the detected band image, the attenuation value (i.e., subtraction or multiplication) of each pixel in the in-band image can be evaluated. These attenuation values ​​can be obtained pixel by pixel from the corresponding pixels in the detected band image (i.e., Bn θ,d Using the value in Bi θ,dpixel intensity in can be attenuated), or can be obtained by processing a region of interest in the detection band image corresponding to the local spatial adjacent positions of the in-band image pixels (e.g., Bn θ,d Using the values of the surrounding 3×3 adjacent positions (e.g., polar coordinates or Cartesian coordinates), Bi θ,d pixel intensity in can be attenuated).

[0190] Refer to FIGS. 6G, 6H, and 6J. Separate embodiments are shown: another embodiment (FIG. 6G) using a reference reception channel during noise characteristic identification, another embodiment (FIG. 6H) for determining an appropriate subtraction value for noise suppression during imaging, and another embodiment (FIG. 6J) for determining an attenuation coefficient. In FIG. 6G, during noise characteristic identification, instead of the detection band noise characteristic identification waveform of FIG. 6A, a reference noise characteristic identification waveform 409 is used. The reference channel filter is, for example, an imaging band pass filter. Instead of this, when performing noise estimation from the input of out-of-band noise, the reference channel filter is not an imaging band pass filter. In FIG. 6H, instead of the detection band imaging waveform of FIG. 6B, a reference noise detection waveform 439 is used to determine a subtraction value, and noise suppression is applied by subtraction. Similarly, in FIG. 6I, instead of the detection band imaging waveform of FIG. 6B, a reference noise detection waveform 439 is used to determine and apply an attenuation coefficient for noise suppression. Therefore, the methods of FIGS. 6A to 6F can be applied to this embodiment by replacing the detection band noise characteristic identification waveform 408 (and the corresponding power measurement results) with the reference noise characteristic identification waveform 409 of FIG. 6G, and by replacing the detection band imaging waveform 438 (and the corresponding array and power measurement results) with the reference noise detection waveform 439 shown in FIGS. 6H and 6I.

[0191] <Embodiment 6: Noise Suppression of Quasi-Periodic Noise Sources> In this embodiment, noise suppression is performed by estimating and subtracting in-band noise. The noise is assumed to originate from a quasi-periodic noise source or a quasi-periodic sequence of noise sources. The in-band noise is estimated based on the measurement results in the noise characteristic identification stage when the imaging transducer receiving circuit is not receiving imaging energy.

[0192] In the first stage of this embodiment, waveforms within the imaging band and noise detection band are detected and sampled when imaging energy is absent (e.g., when the ultrasonic transducer is not receiving imaging energy), thereby obtaining pairs of simultaneously existing in-band and detection band (out-of-band) noise characteristic identification arrays.

[0193] As shown in Figure 8A, an intraband noise characteristic identification waveform 590 can be obtained by applying an imaging bandpass filter 200 to the input waveform detected from the imaging transducer receiving channel 13, and optionally by detecting the envelope of the filtered data (210). A detected band noise characteristic identification waveform 595 can be obtained by applying a noise detection bandpass filter 410 to the input waveform detected from the imaging transducer receiving channel 13 (at least one noise detection band includes signals from outside the imaging band), and optionally by detecting the envelope of the filtered data (411). The intraband noise characteristic identification waveform and the detected band noise characteristic identification waveform are measured when no imaging energy is present, and these can be sampled to obtain an intraband noise characteristic identification array and a detected band noise characteristic identification array. These are represented as Ci and Cn, respectively. It is assumed that Ci (and Cn) capture one or more periods of a periodic noise source.

[0194] When in-band and detected-band noise characteristic recognition arrays are acquired in the noise characteristic recognition stage, adjustment parameters for subtracting the in-band characteristic recognition array from the in-band imaging array can be estimated using the correlation between the detected-band characteristic recognition array and the detected-band imaging array.

[0195] As shown in Figure 8B, during imaging, the input waveform from the imaging transducer receiver circuit 13 is filtered (200 and 410) to provide an in-band imaging waveform and a detected band imaging waveform. These waveforms are sampled (before or after envelope detection) to obtain the in-band imaging array 537 and the detected band imaging array 538. The in-band and detected band imaging arrays are represented as Qi and Qn, respectively.

[0196] According to this method, noise is subtracted from the in-band imaging array Qi using the in-band noise characteristic identification array Ci, and each in-band imaging array is processed for noise suppression. However, in order to perform noise suppression by subtraction, it is necessary to time-align the in-band noise characteristic identification array Ci with the in-band imaging array Qi so that the noise coexists. This alignment is possible in the case of a periodic noise source that generates noise in the imaging band correlated with the noise in the noise detection band.

[0197] Time alignment can be achieved, for example, by segmenting the in-band and detection-band imaging arrays Qi and Qn into multiple time windows 565 and 566 (as described in the above embodiment). The windows should preferably be long enough to capture one or more periods of the periodic noise source. The imaging arrays time-segmented according to the windows are Qi w and Qn w It is represented as follows. The sign w is an integer indicating the window number.

[0198] In one embodiment, time alignment can be achieved for each window. In this embodiment, time alignment is achieved within each window by the detection band noise characteristic identification array Cn and the detection band imaging array Qn w This can be achieved by calculating the cross-correlation between and and by selecting the relative time delay τ corresponding to the maximum cross-correlation (580). Due to the coexistence relationship between noise in the imaging band and noise in the noise detection band, this time delay τ is applied to the in-band noise characteristic identification array Cn and the in-band imaging array Qi wIt can be aligned for each window (510). A scaling factor can also be applied to the aligned in-band noise characteristic identification array.

[0199] The window part of the in-band noise characteristic identification array after alignment ( T w shown as) is subtracted from the in-band imaging array Qi w to obtain the noise-suppressed in-band imaging array Qi w This is performed as an option after taking the absolute value after subtraction or applying a floor function to remove negative values. This process is repeated for each window for which noise suppression is attempted.

[0200] In one embodiment, when adjacent windows overlap, a scaling factor may be applied when subtracting the aligned window segment of the in-band noise characteristic identification array from the in-band imaging array Qi w For example, the subtraction (and optionally modulo) can be calculated by T w =Qi Qi w =Qi w -α T w where α is the scaling factor for windowing, α = 1 - β, and β is the overlap factor. For example, when β = 0.75, the scaling factor is α = 0.25. This embodiment is provided to illustrate an example of a method for scaling the subtraction component, and it should be understood that other functional forms may be adopted.

[0201] Refer to Figures 8C and 8D. Another embodiment is shown in which a reference receiving channel is used to determine appropriate amplitude adjustment when performing noise characteristic identification and when suppressing noise during imaging. In Figure 8C, when performing noise characteristic identification, the reference noise characteristic identification waveform 596 is used instead of the detection band noise characteristic identification waveform 595 in Figure 8A. Similarly, in Figure 8D, the filtered reference noise detection waveform 439 is used instead of the detection band imaging waveform 595 in Figure 8B to determine and apply amplitude adjustment. Therefore, the methods described in Figures 8A and 8B can be applied to this embodiment by replacing the detection band noise characteristic identification waveform 408 with the reference noise characteristic identification waveform 409 in Figure 8C, and by replacing the detection band imaging waveform 438 (and the corresponding array) with the reference noise detection waveform 439 in Figure 8D.

[0202] <Embodiment 7: Noise suppression using multiple scans by changing the scan rate> The system can be configured to acquire two or more sets of in-band imaging arrays corresponding to imaging energies from the same scanline, or to scanlines that substantially overlap in space. The imaging energies within the set of in-band imaging arrays have redundant time / depth dependencies. By averaging the redundant in-band imaging arrays (or performing other statistical processing such as evaluating the minimum value), noise can be suppressed if the noise itself is not time-fixed relative to the trigger that causes the imaging transducer receiving circuit to begin receiving imaging energy on each scanline. For example, if the pulse repetition frequency of the voltage pulse that excites the imaging ultrasonic transducer is 200 μs, the periodic noise that repeats every 2 μs will always have components such as 0 μs, 2 μs, 4 μs, etc., in each in-band imaging array. However, if the pulse repetition frequency is adjusted to 199 μs, the first imaging array will have noise components such as 0 μs, 2 μs, 4 μs, etc., and the second imaging array will have noise components such as 1 μs, 3 μs, 5 μs, etc. Noise can be suppressed by averaging two consecutive redundant residualband imaging arrays, and optionally after performing envelope detection.

[0203] In this embodiment, the period of one or more in-band noise sources can be determined using the noise in the detected band waveform (out-of-band). The system is prompted to adjust the scan rate so that the imaging scan period is not an integer multiple of the noise source period. For example, an automatic correlation function can be used to detect the periodicity of the detected band imaging waveform while the imaging transducer receiving circuit is receiving imaging energy. Alternatively, the period of one or more noise sources can be determined from the detected band noise characteristic identification waveform or a reference noise characteristic identification waveform during the noise characteristic identification stage when the imaging transducer is not receiving imaging energy, or it can be read from a pre-stored database. The system is prompted to adjust the scan rate so that the imaging scan period is not an integer multiple of the noise source period.

[0204] After determining the optimal scan rate and adjusting it, the in-band imaging waveforms acquired from multiple scan lines are sampled after envelope detection to obtain a set of in-band imaging arrays. The in-band imaging arrays are then Qi θ This is expressed as follows: θ is the scan line. Samples from the in-band imaging array are taken using Qi θ Denoted as [k], k=1··K is the sample index, where K is the number of samples in the array. The system can be configured to suppress noise by averaging (or performing other statistical measurements such as taking the minimum value) multiple arrays corresponding to adjacent scan lines or scan lines with strong spatial overlap. For example, if the system is configured to group three in-band imaging arrays, sample Qi θ [k] is [Qi θ-1 [k], Qi θ [k], Qi θ+1 [k]] is replaced with the mean. Optionally, sample Qi θ [k] is [Qi θ-1 [k], Qi θ [k], Qi θ+1[k]] is selectively retained after performing other numerical analyses on it and determining whether or not to suppress noise in the sample. For example, [Qi θ-1 [k], Qi θ [k], Qi θ+1 In [k], if the minimum value is greater than half the maximum value, the sample values ​​may not be large enough to suppress noise through averaging, and are therefore maintained without modification.

[0205] This method is useful in MRI imaging where periodic noise sources are present. This is because the RF excitation pulse is reliably output at a time that does not correlate with the timing or periodicity of the noise source in the local environment.

[0206] <Another component: Selecting samples to suppress noise when using multiple scans> In this embodiment, noise suppression is performed by adjusting the scan rate as described above and selectively replacing portions of the in-band imaging array based on statistical measurements (e.g., average or minimum) from an array of multiple adjacent scan lines. The overlap with adjacent scan lines in the scan region is assumed to be sufficiently large. Only segments of the in-band imaging array that are evaluated as having noise using the detection band measurement results are replaced.

[0207] As shown in Figure 8E, the input waveform is filtered (200 and 410) to provide in-band and detected band waveforms, and these are sampled after envelope detection to obtain the in-band imaging array and detected band imaging array.

[0208] Data corresponding to multiple scanlines is recorded as an array set for each scanline. An array set is acquired for each scanline. For each scanline, one in-band imaging array is acquired within the imaging band, and at least one out-of-band imaging array is acquired for each scanline. The in-band and detection band imaging arrays are each Qi θ and Qn θThis is expressed as follows: θ is the scan line. Samples from the in-band imaging array and the detection band imaging array are Qi, respectively. θ [k] and Qn θ This is represented as [k]. k=1··K is the sample index, and K is the number of samples in the array.

[0209] The detection band imaging array is segmented into multiple windows, each corresponding to Qn θ It contains J samples. The array segmented according to the window is Qn θ,w It is expressed as follows: the sign w is an integer representing the window number, θ is an integer representing the scanline, and Qn θ,w This is a sample [Qn θ [kw], Qn θ [kw+1],···Qn θ Includes [kw+J-1]]. kw is the index of the first sample within the window.

[0210] Each window Qn within the detection band θ,w The presence or absence of noise is evaluated (600). If the noise (determined by the measurement results of waveform energy such as peaks and RMS) exceeds the threshold, the window is considered to have noise. The threshold can be selected in the noise characteristic identification stage.

[0211] According to this embodiment, each detection band array window Qn is classified as having noise. θ,w Regarding the concurrently existing intraband samples (i.e., Qn θ [kw], Qn θ [kw+1],···Qn θ [kw+J-1]) is identified as a sample suitable for noise suppression.

[0212] When windows overlap, given sample Qi θ[k] may be associated with two or more windows. There are cases where a sample is associated with both a noisy window and a noise-free window. In this case, the system can be configured to pool noise evaluations from multiple out-of-band windows before determining whether the sample is suitable for noise suppression.

[0213] Using statistical values ​​(e.g., duplicate value, mean, minimum, etc.) from array samples corresponding to adjacent scan lines (calculated in 610), samples deemed suitable for noise suppression are replaced (620). For example, sample Qi θ If [k] is considered suitable for noise suppression and the system is configured to group arrays from three scan lines, then sample Qi θ [k] is [Qi θ-1 [k], Qi θ [k], Qi θ+1 This can be replaced with the minimum value of [k].

[0214] Refer to Figure 8F. Another embodiment is shown in which a reference waveform is used instead of the detection band imaging waveform in Figure 8E when noise suppression is performed during imaging. Therefore, by replacing the detection band imaging waveform (and the corresponding array) with the filtered reference noise detection waveform, the method described in Figure 8D can be applied to this embodiment.

[0215] <Time-domain vs. frequency-domain processing> The above embodiments have been described in the context of time-domain processing. However, many embodiments herein can utilize frequency-domain or time-frequency-domain processing in one or more steps. For example, in steps 570 and 572 of Figures 6A and 6C, instead of using the minimum power ratio of the imaging band to the noise detection band along with the maximum power, short-time Fourier transforms or wavelet transforms can be performed on the in-band and out-of-band waveforms motion-by-motion or array-by-array. Using analysis of the transformation coefficients (e.g., mean, mean square, etc.), noise can be characterized or windowed, whether or not noise is present. When noise is detected, instead of attenuating the time-domain signal (526), ​​the transformation coefficients of the current window can be attenuated. Then, an inverse transform can be performed on the attenuated frequency-domain signal to obtain a noise-suppressed time-domain signal.

[0216] In embodiments where cross-correlation is required between two time-continuous waveforms, such as Figure 2A (step 510), Figure 8B, and Figure 8D, and other embodiments, a Fourier transform algorithm can be used to efficiently calculate the cross-correlation.

[0217] Furthermore, as explained in Figures 7A-7D, noise patterns can be classified in step 570 using machine learning algorithms. These patterns can be defined by frequency domain and / or time-frequency domain characteristics. These require frequency domain or time-frequency domain processing of time-continuous waveforms.

[0218] <Generalization to applications other than ultrasound> Although the above examples of systems and methods for suppressing image noise have been described in the context of ultrasound imaging, it should be understood that the embodiments described herein are applicable to a variety of imaging devices, systems, and methods.

[0219] Another example of an imaging system to which noise suppression can be applied according to the above embodiment is a magnetic resonance imaging system. See Figure 9. Another system example is shown in which a signal to be noise-suppressed is obtained from a magnetic resonance (MR) system. This system has a magnetic resonance scanner 50 that generates a main magnetic field B0 using a main magnet 52. The main magnetic field B0 causes polarization in the patient 60 or subject. The system has a gradient coil 54 that generates a magnetic field gradient. A receiving coil 58 detects the MR signal from the patient 60. The receiving coil 58 can also be used as a transmitting coil. Alternatively, a body coil 56 can be used to emit and / or detect radio frequency (RF) pulses. The RF pulses are generated by an RF unit 65, and the magnetic field gradient is generated by a gradient unit 70. The manner in which the MR signal is detected using the RF pulses and magnetic field gradient, and the method of reconstructing the MR image, are known to those skilled in the art.

[0220] The reference receiving circuit may include a coil located in the same room as the scanner, but it is not positioned in close proximity to the imaging sample (e.g., patient) emitting the MRI signal. Electromagnetic noise propagating near the MRI machine is detected by both the imaging receiving coil 58 and the reference receiving circuit. The coil of the reference receiving circuit is positioned and oriented to receive the same noise as the imaging receiving coil, but is far enough away from the imaging sample that the imaging energy detected by the reference receiving circuit can be ignored.

[0221] The reference receiving circuit coil is tuned to have the same bandwidth as the imaging receiving coil 58, or to have a different bandwidth capable of collecting ambient noise signals correlated with the noise received by the imaging receiving coil.

[0222] The reference receiving circuit further includes a group of receiving coils, for example, three coils aligned orthogonally to each other. This allows for the collection of electromagnetic noise in such a manner that the weighted sum of the noise collected by the three coils more closely matches the noise collected by the imaging receiving coil. This allows for consideration of the directionality of the dominant electromagnetic noise source in the MRI environment.

[0223] It should be understood that the MR system may have other units or components not shown for clarity. These may include, but are not limited to, control devices, input devices, and detection devices (e.g., cardiac and / or respiratory synchronization devices). Furthermore, each unit may be implemented in a manner other than the individual units shown. Different components can be incorporated into a unit, or different units can be combined. Each unit (shown as a functional unit) can be designed as hardware, software, or a combination thereof.

[0224] In the system shown in Figure 9, the control processing hardware 100 acquires magnetic resonance images of the patient 60 according to an appropriate pulse sequence. The control processing hardware 100 has an interface with a magnetic resonance imaging scanner 50 that receives the acquired images and controls image acquisition. The control processing hardware 100 receives image data from the RF unit 65 and processes the image data according to the method described above.

[0225] The control processing hardware 100 can be programmed by a set of instructions that, when executed by the processor, cause the system to perform one or more of the above methods. This suppresses noise in the signals acquired from the magnetic resonance imaging system. For example, as shown in Figure 9, the control processing hardware 100 can be programmed by instructions in the form of a set of executable image processing modules. Examples include, but are not limited to, a pulse sequence generation module (not shown), an image acquisition module (not shown), an image processing module 145, and a noise suppression module 150. The pulse sequence generation module, image acquisition module, and image processing module can be implemented using pulse sequence generation, image acquisition, and image reconstruction algorithms, respectively, that are known to those skilled in the art. RF data is received from RF coils 56 and / or 58, and optionally from one or more reference receiving coils. The data is sampled and filtered to obtain an in-band waveform. A reference waveform via the reference receiving circuit, or a noise detection band waveform obtained by filtering the RF from coils 56 and / or 58, is collected. One or more noise suppression methods described in Figures 2 to 8 can be used for noise suppression (100). The pulse generation module generates a sequence of RF pulses and a magnetic field gradient based on the desired imaging sequence. The image acquisition module stores the MR signals detected by coils 56 and / or 58 in the raw data space. The image processing module 145 processes the acquired noise-suppressed RF data to perform image reconstruction of the MRI image.

[0226] This system provides a function to estimate in-band noise by detecting noise correlated with noise in the bandwidth of the imaging signal (via a reference receiving channel or noise detection band), and a function to improve the imaging signal by removing the estimated noise from the imaging signal.

[0227] This improves the signal-to-noise ratio (SNR) in cages that shield MRI systems from ambient noise, or improves the operation of MRI systems in open / unshielded environments that are more susceptible to noise.

[0228] <Examples> The following examples are provided to enable those skilled in the art to understand and implement embodiments of the disclosure. They are for illustrative purposes only and should not be construed as limiting the scope of the disclosure.

[0229] <Example 1: Noise suppression of an unknown noise source using an attenuation coefficient (Example of Embodiment 5)> This example concerns the acquisition of ultrasound data using an intracardiac echocardiography (ICE) system in the presence of two noise sources. The transducer was configured to detect ultrasound energy at a frequency of 9 MHz. Two bandpass filters were used in parallel to separate the radio frequency (RF) signal into an imaging band of 7–13 MHz and a noise detection band of 15–25 MHz, which exceeds the frequency range of the emitted ultrasound.

[0230] The first noise source is the electroanatomical mapping system (Carto® 3). The system has an electromagnetic tracking module and an impedance-based tracking module, and patches are attached to the patient to estimate impedance and device position. These patches can link a large amount of noise to the imaging bands of the ICE image. In this experiment, a cardiac dynamics phantom was used in a saline bath. Electrodes from the impedance patches were immersed in the water bath. The second noise source is the ablation generator connected to the Carto 3 console. The noise generated from this second noise source was assumed to be noise propagating from the ablation generator to the saline solution via the Carto 3 console and patch electrodes when the generator is ON.

[0231] Figure 10A shows an ultrasound image collected when no noise source is present. Figures 10B and 10C show the effects of the first and second noise sources on image quality, respectively.

[0232] Noise suppression of the ultrasonic waveform detected by the ultrasonic transducer of the ICE console was performed according to the methods shown in Figures 6A and 6D.

[0233] In the first baseline noise characteristic recognition stage of this experiment, when there was no received imaging energy or received noise energy, the energy within the imaging band was detected. This allowed us to obtain an in-band baseline noise characteristic recognition array (shown as Gi). The in-band baseline noise characteristic recognition array was obtained by sampling an in-band waveform for 125 μs at 200 MS / s. Using slides and superimposed windows (window size = 64 samples, 20% overlap), the in-band power measurement results were calculated for each window. The representative maximum power (90th percentile of the entire window) was calculated, and this value was assigned to the threshold Ti.

[0234] In the second noise characteristic discrimination stage of this experiment, when no received imaging energy was present (i.e., when the transducer was not supplied with a voltage pulse and was not receiving ultrasonic energy), energy was detected in both the imaging band and the noise detection band. This allowed us to obtain correlation measurements between the in-band noise and the noise detection band noise (Ci and Cn). The noise detection band was configured as a frequency band of approximately 15-25 MHz.

[0235] In this embodiment of the ICE system, 512 waveforms were acquired (each with a period of 125 μs and a sampling rate of 200 MS / s). Therefore, noise characteristics were estimated using 512 pairs of in-band and detected-band noise characteristic identification waveforms. For both the imaging band and the noise detection band, sliding and superimposed windows (window size = 64, 20% overlap) were used to determine the in-band and detected-band noise characteristic identification arrays (Ci w and Cn w The window pair was obtained. w is an integer representing the window number. w and Cn w For each pair, the power of the imaging band and the power of the noise detection band were calculated. Only windows where the intraband power was greater than the threshold Ti were selected for noise characteristic identification.

[0236] Statistical noise power measurement results were calculated in noise characteristic identification. A window set of the in-band noise characteristic identification array was selected where the power is in the 96th to 99th percentile of the in-band noise characteristic identification array power value. Within this set, the window with the noise detection band power near the minimum value (20th percentile in the set) was selected as w_max, and the power value P(Cn w_max ) and P(Ci w_max ) was calculated (Pn each) max and Pi max (As shown). Similarly, a window set of the in-band noise characteristic recognition array was selected where the in-band power is in the 1st to 5th percentile of the in-band noise characteristic recognition array power value. Within this set, the window with the noise detection band power near the minimum value (20th percentile in the set) was selected as w_min, and the power value P(Cn w_min ) and P(Ci w_min ) was calculated (Pn each) min and Pi min (To be shown as)

[0237] As shown in Figure 6B, point (Pn max Pi max ) and point (Pn min Pi min The straight line passing through ) has a gradient m and a y-intercept c. The relaxation parameter β is set to a unit value. The function f is given as follows:

number

[0238] The function f() was used to suppress noise in the in-band image waveform acquired during imaging (i.e., when the pulsar periodically emits imaging energy and the transducer receiving circuit receives ultrasonic energy).

[0239] Sample the waveforms of the imaging band and the noise detection band to create a pair of in-band imaging arrays and detection band imaging arrays (Qi θ and Qn θ ) was obtained. θ is an index that identifies the acquisition period corresponding to the scan line. The imaging array was segmented temporally according to the window, and each segment was Qi θ,w and Qn θ,w This is shown as follows: w is an integer indicating the window number. For each window portion of the in-band imaging array and the out-of-band imaging array, the power P(Qi θ,w ) and P(Qn θ,w ) was calculated (each using Pi θ,w and Pn θ,w (To be shown as)

[0240] For each window, out-of-band imaging power Pn θ,w Using the in-band imaging window Qi θ,w Internal noise energy P^i Nθ,w We estimated P^i Nθ,w is f(Pn θ,w ) was calculated as follows. β was set to a unit value, and the parameters m, c, and Pn were used in the noise characteristic identification stage. max、 Pn min Obtained P^i Nθ,w to Pi θ,w The following was clamped. The damping coefficient was set to [Pi θ,w -P^i Nθ,w ] / Pi θ,w The calculation was performed as follows. The damping coefficient was clamped to 0.02 or higher. The damping coefficient was set to array Qi θ,w Multiply by the noise attenuation array Q i θ,w The image was obtained by the image generator (230) for all 512 scan lines. Q i θ It was processed by including it.

[0241] When the first and second noise sources were present, noise in the in-band imaged waveform was suppressed using the above algorithm. Figures 11A and 11B show images acquired when the first noise source was present. Figure 11A shows the result with this noise suppression algorithm applied, while Figure 11B shows the result without it. Figures 12A and 12B show images acquired when the second noise source was present. Figure 12A shows the result with this noise suppression algorithm applied, while Figure 12B shows the result without it. In both cases, the signal-to-noise ratio is clearly improved (approximately 6 dB improvement). As an alternative example, Figures 11C, 11D, and 11E show images acquired by applying this noise suppression method when the first noise source was present. The relaxation parameter β of the function f() is 0.5 (Figure 11C), unit value (Figure 11D), and 1.5 (Figure 11E).

[0242] <Example 2: Noise suppression of periodic noise source using delay correction (Embodiment 6)> In this example, noise suppression of the ultrasonic waveform detected by the ultrasonic transducer of the ICE console was performed according to the method shown in Figures 8A and 8B. Data for this example was collected by the ICE console when an electromagnetic tracker was present. A pseudo-periodic noise pattern was generated in the ultrasonic image data by observing the control unit of the electromagnetic tracker. This is shown in Figure 13A.

[0243] In this embodiment of the ICE system, when the imaging transducer was not receiving imaging energy and data acquisition was performed in raw / RF mode, one waveform (length 125 μsm, sampling 200 MS / s) was acquired.

[0244] The data was stored as array pairs (i.e., sampled waveforms). The first array has a corresponding array containing in-band noise characteristic identification waveforms (7-13 MHz) sampled after envelope detection, and out-of-band noise characteristic identification waveforms in the 15-25 MHz band sampled after envelope detection. These in-band and detected band noise characteristic identification arrays are shown as Ci and Cn.

[0245] By using the correlation between in-band and out-band noise, noise suppression of the in-band image waveform was performed during image data acquisition (i.e., when the ultrasonic transducer is receiving ultrasonic energy).

[0246] During imaging, array pairs were acquired at multiple acquisition periods corresponding to the scanline θ. For each array pair, one array was acquired in the imaging band (7-13 MHz) and one array in the noise detection band (15-25 MHz). Each array had a length of 125 μs and a sampling rate of 200 MS / s. The in-band and out-of-band imaging arrays were each Qi-based. θ and Qn θ This was expressed as such, and measured during imaging (i.e., when the transducer was receiving imaging energy).

[0247] Each in-band imaging array was processed for noise suppression using an in-band noise characteristic identification array, thereby suppressing noise from the in-band imaging arrays.

[0248] Qi Intraband and Out-of-Band Imaging Array θ and Qn θ The image was segmented into multiple time windows. In this embodiment, the window size was 800 samples (4 μs at a rate of 200 MS / s in this example), with 75% overlap between adjacent windows. The imaging array, segmented temporally according to the windows, was processed using Qi θ,w and Qn θ,w This is expressed as follows: w is an integer representing the window number.

[0249] Qn θ,w Time alignment between Cn and Qn was performed for each window. The time alignment was performed between the out-of-band noise characteristic identification array Cn and the out-of-band imaging array Qn. θ,wThis was achieved by calculating the cross-correlation between and and selecting the relative time delay τ corresponding to the maximum cross-correlation. Due to the correlation between in-band noise and out-of-band noise, a time delay τ was applied to align the in-band noise characteristic recognition array with the windowed in-band imaging array. The aligned in-band noise characteristic recognition array was then windowed ( T w (This is how it is expressed.)

[0250] Intraband imaging array Qi θ,w Before subtracting from, the in-band noise characteristic identification array T w A scaling factor was applied to it. Qi θ,w =Qi θ,w -α T w Accordingly, noise suppression in-band imaging array Qi θ,w The following was calculated. α = 0.25 is the scaling factor to handle windowing. Negative values ​​were replaced with 0.

[0251] Using the method described above, noise in data acquired during imaging in the presence of the Aurora® electromagnetic tracking system (Northern Digital Inc.), which acts as a noise source, was suppressed. Figures 13A and 13B show images acquired when the noise source is present. Figure 13A shows the result with this noise suppression method applied, while Figure 13B shows the result without it. When the noise suppression method was implemented, a clear improvement in the signal-to-noise ratio (approximately 5 dB) was observed.

[0252] The specific embodiments described above are illustrative, and it should be understood that these embodiments are capable of various modifications and alternative forms. Furthermore, it should be understood that the claims are not intended to limit the scope to any particular form, but rather to cover all modifications, equivalents, and alternative forms included in the gist and scope of this disclosure.

Claims

1. A method for removing noise from an image signal when broadband noise is present, If imaging energy is not received, the imaging transducer receiving circuit detects the energy wave to obtain a noise characterization waveform, and filters the noise characterization waveform using an imaging bandpass filter and a noise detection bandpass filter, respectively, to generate an in-band noise characterization waveform that exists within the imaging band and an out-of-band noise characterization waveform that is at least partially within the noise detection band outside the imaging band; A step of processing the in-band noise characterization waveform and the out-of-band noise characterization waveform to determine the functional relationship between the noise in the imaging band and the noise in the noise detection band; A step of obtaining one or more imaging waveforms by detecting the imaging signal with the imaging transducer receiving circuit; For at least one of the aforementioned imaging waveforms: A step of filtering the imaging waveform using the imaging bandpass filter and the noise detection bandpass filter, respectively, to generate an in-band imaging waveform that exists within the imaging band and an out-of-band noise detection imaging waveform that exists within the noise detection band; A step of estimating noise present within the imaging band using the above-mentioned functional relationship and the out-of-band noise detection imaging waveform; A step of reducing noise in the imaging waveform within the band by applying noise reduction correction to at least a portion of the imaging waveform within the band using the estimated noise within the imaging band; A method of having.

2. The method according to claim 1, wherein the noise reduction correction is determined and applied for each time window of at least two time windows of the in-band imaging waveform.

3. The step of detecting the imaging signal includes the step of detecting the imaging signal using the imaging transducer receiving circuit along a plurality of scan lines, and the method further includes A step of obtaining a plurality of imaging waveforms, each associated with the plurality of scan lines; A step of generating an image based on denoised in-band imaging waveforms associated with each of the plurality of scan lines; Having, The method according to claim 1.

4. The method according to claim 1, wherein the noise reduction correction is applied within a given time window when a measured value related to the amount of noise in the in-band imaging waveform within a given time window is greater than a threshold.

5. The method according to claim 4, wherein the measured value relating to the amount of noise in the in-band imaging waveform is based on the functional relationship and the amount of power of the in-band noise characterizing waveform.

6. The method according to claim 4, wherein, before applying the noise reduction correction, each time window having the measured value exceeding the threshold is re-evaluated with respect to a spatially adjacent time window corresponding to a spatially adjacent scan line, thereby replacing a sample initially estimated to be associated with the presence of noise with a sample from a spatially adjacent time window where the adjacent time window is estimated to be associated with the absence of noise.

7. The method according to claim 4, wherein, before applying the noise reduction correction, each time window having the measured value exceeding the threshold is re-evaluated with respect to a temporally adjacent time window, thereby re-evaluating a time window initially estimated to be associated with the absence of noise as being associated with the presence of noise if the adjacent time window is estimated to be associated with the presence of noise.

8. The method according to claim 4, wherein the threshold is determined based on baseline measurements performed when imaging energy and noise are absent.

9. The method according to claim 1, wherein the functional relationship between the noise in the imaging band and the noise in the noise detection band is associated with the measured relative energy within the in-band noise characterization waveform and the out-of-band noise characterization waveform.

10. The noise reduction correction method according to claim 1, wherein the noise reduction correction includes a reduction value.

11. The method according to claim 10, wherein the noise reduction correction is determined and applied for each time window, and the reduction value associated with a given time window increases with the amount of energy detected in the out-of-band noise detection imaging waveform in the given time window.

12. The noise reduction correction method according to claim 1, wherein the noise reduction correction includes multiplication by an attenuation coefficient.

13. The method according to claim 12, wherein the noise reduction correction is determined and applied for each time window, and the attenuation coefficient related to a given time window decreases with the amount of noise in the in-band image waveform for the given time window, thereby attenuating the portion of the in-band image waveform related to the noise.

14. The method according to claim 1, wherein the noise reduction correction is determined and applied for each time window, and the noise reduction correction applied to a predetermined time window depends on the depth of the predetermined time window in the in-band imaging waveform.

15. The method according to claim 1, further comprising the step of repeating noise characterization to re-establish the functional relationship between the noise in the imaging band and the noise in the noise detection band.

16. The method according to claim 15, wherein the noise characterization is repeated according to user input.

17. The method according to claim 15, wherein the noise characterization is automatically repeated when the absence of imaging energy is detected.

18. The method according to claim 1, further comprising the step of detecting a change in the noise within the noise detection band by monitoring one or more parameters related to the out-of-band noise detection imaging waveform.

19. The method according to claim 18, further comprising the step of generating a warning or message when a change in at least one parameter is detected.

20. The above method further, A step of performing noise suppression on the in-band imaging waveform while no imaging energy is present; A step of generating an error value based on the energy of the in-band imaging waveform after noise correction; If the error value exceeds a pre-selected threshold, the step of generating an alert to repeat the noise characterization; Having, The method according to claim 1.

21. The method according to claim 1, wherein the imaging transducer receiving circuit includes an ultrasonic transducer.

22. The method according to claim 1, wherein the imaging transducer receiving circuit includes a coil for detecting a magnetic field.

Images