WO2018152643A1 - Systems and methods for noise reduction in imaging - Google Patents

Systems and methods for noise reduction in imaging Download PDF

Info

Publication number
WO2018152643A1
WO2018152643A1 PCT/CA2018/050212 CA2018050212W WO2018152643A1 WO 2018152643 A1 WO2018152643 A1 WO 2018152643A1 CA 2018050212 W CA2018050212 W CA 2018050212W WO 2018152643 A1 WO2018152643 A1 WO 2018152643A1
Authority
WO
WIPO (PCT)
Prior art keywords
imaging
noise
band
waveform
detection
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Ceased
Application number
PCT/CA2018/050212
Other languages
English (en)
French (fr)
Inventor
Brian Courtney
Naimul Mefraz KHAN
Natasha ALVES-KOTZEV
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sunnybrook Research Institute
Original Assignee
Sunnybrook Research Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority to KR1020197027630A priority Critical patent/KR20190131031A/ko
Priority to JP2019546178A priority patent/JP7138113B2/ja
Priority to KR1020247026556A priority patent/KR20240125696A/ko
Priority to KR1020237024799A priority patent/KR102694702B1/ko
Priority to CN201880024819.8A priority patent/CN110537201B/zh
Priority to EP18756578.3A priority patent/EP3586303B1/en
Priority to KR1020237002915A priority patent/KR20230020570A/ko
Priority to AU2018225834A priority patent/AU2018225834B2/en
Application filed by Sunnybrook Research Institute filed Critical Sunnybrook Research Institute
Priority to CA3054383A priority patent/CA3054383C/en
Publication of WO2018152643A1 publication Critical patent/WO2018152643A1/en
Anticipated expiration legal-status Critical
Priority to JP2022140465A priority patent/JP7529732B2/ja
Priority to JP2024119350A priority patent/JP7842817B2/ja
Ceased legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/70—Denoising; Smoothing
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/5205—Means for monitoring or calibrating
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5269—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/5608—Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52077—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging with means for elimination of unwanted signals, e.g. noise or interference
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/10—Image enhancement or restoration using non-spatial domain filtering
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/10—Segmentation; Edge detection
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • A61B5/7217—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise originating from a therapeutic or surgical apparatus, e.g. from a pacemaker
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/12—Diagnosis using ultrasonic, sonic or infrasonic waves in body cavities or body tracts, e.g. by using catheters
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/46—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5269—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
    • A61B8/5276—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts due to motion
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10088—Magnetic resonance imaging [MRI]
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10132—Ultrasound image
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20048—Transform domain processing
    • G06T2207/20056—Discrete and fast Fourier transform, [DFT, FFT]

Definitions

  • the present disclosure relates to imaging technologies and the processing of imaging data for the removal of noise.
  • Medical imaging with ultrasound and MRI imaging relies on detecting low amplitude signals in the radiofrequency spectrum, typically spanning from 2 MHz to 200 MHz. Image quality is greatly influenced by the signal-to-noise ratio.
  • IVUS intravascular ultrasound
  • ICE intracardiac echocardiography
  • the ultrasound transducer detects ultrasound signals from the surrounding structures and converts the acoustic energy into an electrical signal. This signal is then transmitted along one or more conductive channels (such as coaxial conductors, twisted pair conductors, flex circuits etc.). For many reasons, (including cost,
  • the portion of the minimally invasive imaging probe that can be inserted intracorporeally often does not contain an amplifier to boost the signal strength.
  • the electrical signals detected by minimally invasive ultrasound transducers can be very small ( ⁇ 10mV and more typically ⁇ 1 mv), and much of the information about tissue structures that can be imaged with ultrasound tends to lie in the lower portion of the dynamic range of the electrical signals that are detected.
  • the signal amplitude of a received ultrasound signal is limited by any or all of the mechanical efficiency of the transducer, the low amplitude of the acoustic signals detected, the small size of the transducer and attenuation along the conductors that carry the electrical signal from the transducer out of the body. In light of this, the signals in minimally invasive ultrasound imaging systems tend to be very weak.
  • Noise can be introduced into the system from many sources, including radio transmitters, power electronics, transmission lines, switching transistors and others known in the art. Noise can be introduced via induction or directly via conduction and suboptimal isolation of components that are sensitive to electromagnetic interference. Some of the noise may be generated by components within the imaging system itself, such as scanning actuators, pulse width modulators for motor controllers, switched mode power supplies, clocking circuits and transistors in any of the electronic components of an imaging system. Furthermore, other systems coupled to a patient or in the procedural environment, such as impedance monitors, tracking systems (like those found in Carto ® 3, Carto ® XP or NavXTM systems), temperature sensors, infusion pumps, ablation systems, ECG and hemodynamic monitors can introduce noise.
  • impedance monitors like those found in Carto ® 3, Carto ® XP or NavXTM systems
  • tracking systems like those found in Carto ® 3, Carto ® XP or NavXTM systems
  • temperature sensors like those found in Carto ® 3, Carto
  • RFID inventory control systems used in some clinical areas can also introduce noise.
  • Several approaches are directed at reducing the amount of noise that enters into the ultrasound receive circuitry of ultrasound imaging systems, including selection of components within the system that generate minimal RF noise, electrical isolation, shielding, proper grounding, and physically separating noise-generating components from components that are susceptible to electromagnetic noise. These approaches are often difficult to implement, as the sources of the noise often have preferred characteristics for other reasons (i.e. pulse width modulated motor controllers are energy efficient and have good response times) or are difficult to physically isolate from one another (i.e. it may be desirable to have power electronics in close proximity to the imaging probe or its associated circuitry).
  • Ultrasound signals typically have a known bandwidth and the detected ultrasound signal may be filtered using either analog or digital filtering techniques (often a combination of the two). Analog or digital filtering can be applied to limit the portions of the electrical signal output from the ultrasound receive circuitry to those portions whose frequencies lie within the operational bandwidth (or harmonics thereof) of the ultrasound transducer. Selecting filters with narrow bandwidths and sharp cutoffs can reduce the amount of noise that is allowed into the signals used to generate images or otherwise make use of the ultrasound signals (such as for Doppler
  • Notch or comb filters are helpful in removing narrowband noise within the imaging range of frequencies. Overly aggressive filtering can have the unwanted effect of reducing the amount of signal power that gets accepted for generating images or for other use of the ultrasound signals. It may also negatively impact other performance aspects of an ultrasound imaging system, such as resolution. However, if the passband of the filters is too large, then more noise is accepted into the system.
  • Image processing can further reduce the noise by filtering the image data generated, such as by averaging or removing outlier values.
  • filtering can be applied within the image in the spatial domain by applying a Gaussian filter to a pixel and its neighboring pixels in order to blur or smoothen out any random noise in the image.
  • this tends to reduce the spatial resolution of the image.
  • spatial domain filtering can be applied in the structures being imaged that do not move rapidly with respect to the frame repetition frequency of the imaging modality.
  • a pixel in an image frame can be the average or Gaussian- filtered result of the pixels at similar positions in one or more preceding and / or trailing frames.
  • noise enter into the ultrasound receive signal chain and can become difficult to remove once they enter the system, especially if they are broadband in nature, wherein a portion of the noise lies within the passband of the ultrasound system.
  • the system may be designed to heavily filter out any portions of the noise that are less than 7.5MHz and any portions of the noise that are more than 12.5MHz.
  • the amplitude of the noise within the 7.5- 12.5MHz bandpass may frequently be appreciable relative to the amplitude of the ultrasound signal that is being detected.
  • An electrical signal with rapid transients in it has a very broad frequency domain representation that can easily span all or a portion of the passband of the ultrasound receive signal chain. This is particularly true of power supplies or pulse width modulation circuits where the noise can have a strong enough amplitude to compete with the signal being detected.
  • broadband noise exploits the fact that imaging energy predominantly lies within a selective imaging band but that broadband noise can be detected both within the imaging band and outside of the imaging band.
  • detecting or characterizing noise outside of the imaging band at any point in time one can estimate the broadband noise that might lie within the imaging band and alter the detected signal to reduce the estimated in-band noise.
  • By effectively creating an estimate of the in-band noise based on out-of-band noise one can generate a signal that estimates the desired imaging energy in the absence of the estimated in-band noise.
  • reference receive circuits comprising resistors, capacitors, inductors, transmission lines, amplifiers, transformers, inactivated transducers or components that can emulate a transducer receive circuit
  • a method of denoising imaging signals detected in the presence of broadband noise comprising:
  • characterization waveform and the out-of-band noise characterization waveform to determine a relationship between noise in the imaging band and noise in the noise-detection band; detecting imaging signals with the imaging transducer receive circuit thereby obtaining one or more imaging waveforms;
  • a method of denoising imaging signals detected in the presence of noise comprising:
  • detecting noise with a reference receive circuit configured to avoid transduction of imaging energy while detecting noise received by the imaging transducer receive circuit, thereby obtaining a reference noise characterization waveform
  • a method of denoising imaging signals detected in the presence of broadband noise comprising: detecting energy waves with an imaging transducer receive circuit, thereby obtaining an imaging waveform, and filtering the imaging waveform to generate an in-band imaging waveform residing within an imaging band and an out-of-band noise-detection imaging waveform residing within a noise- detection band that lies, at least in part, beyond the imaging band;
  • scaling factor is selected to reduce a contribution of in- band noise in the noise-corrected in-band envelope.
  • a method of denoising imaging signals detected in the presence of noise comprising:
  • the amplitude scaling factor is selected to reduce a contribution of in-band noise in the noise-corrected in-band imaging waveform.
  • a method of denoising imaging signals detected in the presence of noise comprising:
  • detecting energy waves with an imaging transducer receive circuit thereby obtaining an imaging waveform, and filtering the imaging waveform to generate an in-band imaging waveform residing within an imaging band; and detecting noise with a reference receive circuit configured to avoid transduction of imaging energy while detecting noise received by the imaging transducer receive circuit, thereby obtaining a reference noise-detection waveform;
  • scaling factor is selected to reduce a contribution of in- band noise in the noise-corrected in-band envelope.
  • a method of denoising imaging signals detected in the presence of noise comprising: detecting energy waves with an imaging transducer receive circuit, thereby obtaining an imaging waveform, and filtering the imaging waveform to generate an in-band imaging waveform residing within an imaging band; and detecting noise with a reference receive circuit configured to avoid transduction of imaging energy while detecting noise received by the imaging transducer receive circuit, thereby obtaining a reference noise-detection waveform;
  • adaptive filter parameters are actively determined by processing the noise-corrected in-band imaging waveform to minimize the power of the noise-corrected in-band imaging waveform.
  • a method of denoising imaging signals detected in the presence of broadband noise comprising: detecting energy waves with an imaging transducer receive circuit, thereby obtaining an imaging waveform, and filtering the imaging waveform to generate an in-band imaging waveform residing within an imaging band and an out-of-band noise-detection imaging waveform residing within a noise- detection band that lies, at least in part, beyond the imaging band;
  • the adaptive filter parameters are actively determined by processing the noise-corrected in-band imaging envelope to minimize the power of the noise-corrected in-band imaging envelope.
  • a method of denoising imaging signals detected in the presence of broadband noise comprising: detecting energy waves with an imaging transducer receive circuit, thereby obtaining an imaging waveform, and filtering the imaging waveform to generate an in-band imaging waveform residing within an imaging band and an out-of-band noise-detection imaging waveform residing within a noise- detection band that lies, at least in part, beyond the imaging band;
  • a method of denoising imaging signals detected in the presence of broadband noise comprising: detecting energy waves with an imaging transducer receive circuit, thereby obtaining an imaging waveform, and filtering the imaging waveform to generate an in-band imaging waveform residing within an imaging band and an out-of-band noise-detection imaging waveform residing within a noise- detection band that lies, at least in part, beyond the imaging band;
  • a detection system characterized by one or more noise sources comprising:
  • a detection system characterized by one or more known noise sources comprising:
  • detecting noise with a reference receive circuit configured to avoid transduction of imaging energy while detecting noise received by the imaging transducer receive circuit, thereby obtaining a reference noise characterization waveform
  • processing the in-band noise characterization waveform and the reference noise characterization waveform to determine a relationship between noise in the imaging band and noise detected by the reference receive circuit
  • segmenting the in-band imaging waveform and the reference noise- detection waveform according to one or more time windows; for at least one time window, processing the reference noise-detection waveform according to the pattern recognition algorithm to detect the presence of the noise pattern;
  • a method of denoising imaging signals detected in the presence of noise comprising:
  • a method of denoising imaging signals detected in the presence of noise comprising:
  • detecting noise with a reference receive circuit configured to avoid transduction of imaging energy while detecting noise received by the imaging transducer receive circuit, thereby obtaining a reference noise characterization waveform
  • a method of denoising imaging signals detected in the presence of noise comprising:
  • the in-band imaging waveform within the time window is deemed to be suitable for noise correction, applying a denoising correction to the in-band imaging waveform within the time window, wherein the denoising correction for each sample in the window is based on one or more statistical measures associated with samples in the in-band imaging waveforms from two or more adjacent windows, each adjacent window residing in a respective adjacent scan line; and
  • a method of denoising imaging signals detected in the presence of noise comprising:
  • the denoising correction for each sample in the window is based on one or more statistical measures associated with samples in the in-band imaging waveforms from two or more adjacent windows, each adjacent window residing in a respective adjacent scan line; and generating an image based on denoised in-band imaging waveforms respectively associated with the plurality of scan lines.
  • a method of denoising imaging signals detected in the presence of noise comprising:
  • the in-band imaging waveform within the time window is deemed to be suitable for noise correction, applying a denoising correction to the in-band imaging waveform within the time window, wherein the denoising correction for each sample in the window is based on one or more statistical measures associated with samples in the in-band imaging waveforms from two or more adjacent windows, each adjacent window residing in a respective adjacent scan line; and
  • FIG. 1A shows an example of an ultrasound imaging system
  • FIG. 1 B shows an example of a conventional ultrasound receive signal chain for the processing of an ultrasound signal before conversion to an ultrasound image.
  • FIG. 1C shows an example of an ultrasound imaging system including an intracorporeal imaging probe.
  • FIG. 1 D shows an example of an ultrasound imaging system including a reference transducer for detecting in-band noise.
  • FIG. 1 E shows an example of an ultrasound imaging system including a second imaging transducer having an imaging band that lies outside the imaging band of a first imaging transducer.
  • the second imaging transducer may be part of a circuit suitable for detecting in-band noise that affects the signals received from first imaging transducer.
  • FIG. 1 F shows an example of an ultrasound imaging system including a reference receive circuit for detecting in-band noise, where the reference receive circuit extends to a location within the imaging probe.
  • FIG. 1G shows an example of an ultrasound imaging system including a reference receive circuit for detecting in-band noise, where the reference receive circuit is located in one or more portions of the system that are external to the imaging probe.
  • FIG. 2A illustrates an example system configuration for noise reduction on the envelope of an input waveform via the estimation and suppression of the estimated in-band noise, where the in-band noise is estimated by performing envelope detection of an out-of-band waveform, followed by delay, scale and shape adjustment prior to subtraction.
  • FIG. 2B illustrates an example system configuration for noise reduction on an input waveform via the estimation and subtraction of in-band noise, where the in-band noise is estimated by frequency shifting an out-of-band waveform, filtering the frequency-shifted out-of-band waveform, followed by delay, scale and shape adjustment prior to subtraction.
  • FIG 2C illustrates an example of a system configuration for noise reduction with a reference receive circuit that is sensitive to some or all of the in-band noise that the imaging transducer receive circuit is sensitive to, but is at least partially isolated from the imaging signals detected by the imaging transducer receive circuit.
  • the subtraction of the noise signals detected by the reference receive circuit from the signals received by the imaging transducer receive circuit reduce the noise in the output signal.
  • FIG. 3A illustrates an example of a system configuration for noise reduction via active noise cancellation, where a cancellation waveform for active noise cancellation is obtained from a reference receive circuit.
  • FIG. 3B illustrates an example of a system configuration for noise reduction on an input waveform via active noise cancellation, where a cancellation waveform for active noise cancellation is obtained by envelope detection of an out-of-band waveform.
  • FIG. 3C illustrates an example of a system configuration for noise reduction on an input waveform via active noise cancellation, where a cancellation waveform for active noise cancellation is obtained by frequency shifting an out-of-band waveform and filtering the frequency-shifted out-of- band waveform.
  • FIG. 4 illustrates an example of a system configuration for noise reduction on an input waveform via filtering of an in-band waveform, where the filtering is controlled based on a feedback parameter obtained by a filter update algorithm that determines one or more parameters of the filter based on one or more characteristics of a noise-detection waveform.
  • FIG. 5A schematically illustrates an example system configuration for noise reduction on an input waveform based on noise parameters obtained during a first measurement stage in the absence of an imaging signal, and by the application of one or more noise reduction algorithms that utilize the aforementioned noise parameters during a second measurement stage when imaging signals are collected.
  • FIG. 5B schematically illustrates an alternative example system in which a reference receive channel is employed to detect in-band noise.
  • FIGS. 6A and 6C illustrate example system configurations for noise reduction on an input waveform based on detection of noise in an out-of-band waveform, in which different time windows of an in-band waveform are suppressed, based on the processing of a respective window of an out-of- band waveform, and where noise windows of the in-band waveform are corrected by subtracting with a subtrahend value dependent on the amount of power within the window of the out-of-band waveform.
  • FIG. 6B illustrates an example scatter plot showing signal power of windows of an in-band waveform versus signal power of respective windows of an out-of-band waveform in a noise characterization stage.
  • FIG. 6D illustrates an example system configuration for noise reduction on an input waveform based on detection of noise in an out-of-band waveform, in which different time windows of an in-band waveform are suppressed, based on the processing of a respective window of an out-of- band waveform, and where noise windows of the in-band waveform are corrected by multiplication with an attenuation factor dependent on the amount of power within the window of the out-of-band waveform.
  • FIG. 6E and 6F show charts pertaining to a method in which different time windows of an in-band waveform are initially identified as being predominantly signal or noise, after which noise windows surrounded by signal windows are identified as being likely erroneous and are reclassified, and conversely after which signal windows surrounded by noise windows are identified as being erroneous and are reclassified as noise.
  • FIGS. 6G and 6H illustrate an example system configuration for noise reduction on an input waveform based on detection of noise in a filtered reference waveform measured with a reference receive channel, in which different time windows of an in-band waveform are suppressed, based on the processing of a respective window of a reference waveform, and where noise windows of the in-band waveform are corrected by subtracting with a subtrahend value dependent on the amount of power within the window of the filtered reference waveform.
  • FIG. 61 illustrates an example system configuration for noise reduction on an input waveform based on noise measured in a filtered reference waveform measured with a reference receive channel, in which different time windows of an in-band imaging waveform are corrected, based on the processing of a respective window of the filtered reference waveform, and where noise windows of the in-band imaging waveform are corrected by an attenuation factor dependent on the amount of power within the window of the filtered reference waveform.
  • FIG. 7A and 7B illustrate an example system configuration for noise reduction on an input waveform based on noise detected in one or more noise-detection waveforms, of which at least one noise-detection waveform comprises signal that is out-of-band from the imaging band.
  • Different time windows of an in-band waveform undergo noise reduction according to one or more patterns identified by processing one or more noise-detection waveforms.
  • FIG. 7C and 7D illustrate an example system configuration for reduction on an input waveform based on noise detected in a reference waveform, in which different time windows of an in-band imaging waveform undergo noise reduction according to one or more patterns identified by processing one or more reference waveforms.
  • FIGS. 8A and 8B illustrate an example system configuration for noise reduction on an input waveform based on noise detected in an out-of-band waveform, in which different time windows of an in-band waveform undergo noise reduction according to the estimated in-band noise that is temporally aligned prior to reduction.
  • FIGS. 8C and 8D illustrate an example system configuration for noise reduction on an input waveform based noise detected in a filtered reference waveform, in which different time windows of an in-band waveform undergo noise reduction according to the estimated in-band noise that is temporally aligned prior to reduction.
  • FIG. 8E shows an example system configuration for noise reduction on an input waveform based on noise detected in an out-of-band waveform, in which measures from adjacent or replicate scan lines are employed when performing noise correction.
  • FIG. 8F shows an example system configuration for noise reduction on an input waveform based on noise detected in a filtered reference waveform, in which measures from adjacent or replicate scan lines are employed when performing noise correction.
  • FIG. 9 shows an example of a magnetic resonance imaging system configured for noise suppression.
  • FIGS. 10A-C show example images obtained using an intra-cardiac echo system showing (A) an image obtained in the absence of a noise source; (B) an image obtained in the presence of noise generated via an electroanatomic mapping system; and (C) an image obtained in the presence of noise generated from an ablation generator.
  • FIGS. 11A-B show images obtained in the presence of noise from an electroanatomic mapping system, without (A) and with (B) noise reduction.
  • FIGS. 11 C-E show images obtained in the presence of noise from an electroanatomic mapping system after noise reduction by attenuation, where the relaxation parameter was set as 0.5 (C), 1 (D) and 1.5(E).
  • FIGS. 12A-B show images obtained in the presence of noise from an ablation generator, without (A) and with (B) the application of a noise reduction method.
  • FIGS. 13A-B show the images in obtained in the presence of noise from a magnetic tracking system, without (A) and with (B) the application of a noise reduction method.
  • FIG. 14 shows the phrases used to refer to waveforms in the imaging band, noise detection band and waveforms from a reference receive circuit
  • the terms “about” and “approximately” are meant to cover variations that may exist in the upper and lower limits of the ranges of values, such as variations in properties, parameters, and dimensions. Unless otherwise specified, the terms “about” and “approximately” mean plus or minus 25 percent or less.
  • any specified range or group is a shorthand way of referring to each and every member of a range or group individually, as well as each and every possible sub-range or sub -group encompassed therein and similarly with respect to any sub-ranges or sub-groups therein.
  • the present disclosure relates to and explicitly incorporates each and every specific member and combination of sub-ranges or sub-groups.
  • the term "on the order of”, when used in conjunction with a quantity or parameter, refers to a range spanning approximately one tenth to ten times the stated quantity or parameter.
  • Ultrasound imaging relies on receiving echoes from a medium, optionally after sending a narrow acoustic pulse out in the medium in a particular direction.
  • scan line refers to a line representing a spatial direction in the medium from which imaging energy is to be received.
  • a 2D image is obtained by receiving echoes from a plurality of scan lines within the medium.
  • FIG. 1A an example ultrasound imaging system is shown, in which one or more ultrasound transducers 10 are controlled to perform ultrasound imaging across a plurality of scan lines 12.
  • the transducer 10 interfaces with control and processing hardware 100, which optionally controls a transmitter 15 for the generation and emission of imaging energy by the transducer 10.
  • the control and processing hardware 100 is configured to receive ultrasound energy signals detected by the transducer 10, which are routed, typically via a Tx/Rx (transm it-receive) switch 25, to one or more amplifiers 20.
  • the ultrasound transducer(s) 10 may optionally be configured to image a spatial region associated with a plurality of scan lines 12, for example, via mechanical scanning of the transducer 10, or, for example, via electronic scanning via the use of an array of imaging elements, such as, but not limited to, a phased array, ring array, linear array, matrix array or curvilinear array.
  • a transmit beamformer 26 and receive beamformer 27 may be employed to generate a plurality of transmit signals and to beamform a plurality of received signals.
  • receiver circuit generally refers to components such as a transmission line (e.g. coax, PCB tracings, others), connectors, mux/demux, RX/ TX switches 25, amplifiers 20, slip rings, transformers and other components known in the art.
  • transducer receive circuit may include a receive circuit connected to one or more ultrasound transducer elements 10 configured to receive ultrasound signals at the time of use.
  • ultrasound receive signal chain includes a receive circuit, but can include additional components such as analog-to- digital converters (ADCs) and further digital processing components and/or processing logic, including, but not limited to, noise removal processing module 150, before the signal enters the process of being converted into an image (such as via scan conversion) and subsequent image processing.
  • ADCs analog-to- digital converters
  • noise removal processing module 150 includes, but not limited to, noise removal processing module 150, before the signal enters the process of being converted into an image (such as via scan conversion) and subsequent image processing.
  • channel may refer to conductive electrical circuits, wireless channels, optical channels, or other signal paths.
  • an imaging receive channel which denotes the path traversed by detected imaging signals, is shown at 13 in FIG. 1A.
  • the system may employ a single receive channel per transducer, or several receive channels per transducer (such as may be the case for an array transducer where there may a channel for each piezoelectric transducer element or groups of piezoelectric transducer elements in the array).
  • ASICs and other devices may be used along the signal receive chain to multiplex signals along a channel from more than one piezoelectric transducer elements.
  • the system may be configured to control the transducer 10 to optionally transmit energy to a medium, and to detect imaging energy within an imaging frequency band (henceforth referred to as an "imaging band").
  • the imaging band may constitute a single continuous frequency band, or two or more frequency intervals (such as in harmonic imaging), within which imaging energy is detected. Imaging energy or noise within the imaging band are henceforth referred to as being "in-band”.
  • Figure 14 is a representation of the grouping and terminology used to describe the various waveforms.
  • the system may also be configured to detect, via one or more channels connected to the transducer, energy in one or more additional frequency bands, where at least one frequency band lies, at least in part, beyond the imaging band.
  • These one or more additional frequency bands are henceforth referred to as "detection bands”.
  • a waveform that lies, at least in part, beyond an imaging band is henceforth referred to as being "out-of-band".
  • a detection band may reside within the imaging band.
  • Noise- detection bands may be either out-of-band or within-band.
  • At least one detection band may be selected such that the signal-to-noise ratio within the detection band is substantially less than the signal-to-noise ratio in the imaging band when the transducer is used in an imaging mode (i.e. when the transducer detects imaging energy).
  • the detection band may lie outside the full-width, half-maximum bandwidth of the imaging band or another bandwidth corresponding to a threshold below the maximum strength of the signal employed.
  • imaging waveform refers to a waveform (analog or digitally sampled) that is obtained from an imaging transducer receive circuit when the imaging transducer is receiving or is expected to be receiving imaging energy.
  • in-band imaging waveform refers to an imaging waveform (analog or digitally sampled) that lies in the imaging band.
  • An in-band imaging waveform is expected to include imaging energy and may also include unwanted noise energy.
  • an in-band imaging waveform is processed to remove noise energy for the generation of a denoised image.
  • detection-band imaging waveform refers to a waveform obtained from an imaging transducer receive channel and residing within one or more noise detection bands.
  • a detection-band imaging waveform may be out-of-band or within-band.
  • a within-band noise-detection imaging waveform may be employed, for example, in order to confirm the presence of noise within the imaging band. More specifically, a "within-band noise-detection imaging waveform” may be employed to confirm that a noise source having a noise component outside of the imaging band also has a noise component within the imaging band.
  • a detection-band imaging waveform that lies, at least in part, outside the imaging band is referred to as an "out-of-band noise-detection waveform".
  • the system may be configured to be in a "noise- characterization mode", during which the transducer 10 does not transmit energy to a medium, and does not detect imaging energy from the medium.
  • noise-characterization waveform refers to a waveform obtained when the imaging transducer is not receiving imaging energy.
  • in-band noise-characterization waveform refers to a waveform that resides in the imaging band, obtained from an imaging transducer receive channel when the imaging transducer is not receiving imaging energy.
  • detection-band noise-characterization waveform refers to a waveform that resides in a noise detection band, obtained from an imaging transducer receive channel when the imaging transducer is not receiving imaging energy.
  • a detection-band noise characterization waveform that lies, at least in part, outside the imaging band is referred to as an "out-of-band noise characterization waveform”.
  • a detection-band noise characterization waveform that lies entirely within the imaging band is referred to as a "within-band noise characterization waveform”.
  • baseline noise-characterization waveform refers to a waveform obtained when the imaging transducer is not receiving imaging energy, and when a selected noise source is expected to be off (i.e. absent of producing noise), such that the baseline noise-characterization waveform provides a baseline for the selected noise source.
  • a baseline noise- characterization waveform that lies in an imaging band is referred to as an "in- band baseline noise-characterization waveform”.
  • a baseline noise- characterization waveform that lies in a noise-detection band is referred to as a "detection-band baseline noise-characterization waveform”.
  • an optional reference receive circuit 1 1 may be provided that includes a receive circuit configured not to receive reflected ultrasound signals during imaging, while being capable of detecting noise energy similar to the noise that gets coupled into one or more transducer receive circuits during imaging.
  • a reference receive circuit may employ one or more components of a transducer receive circuit (for example, a reference receive circuit and a transducer receive circuit may utilize different channels of an amplifier or an ADC).
  • the system may be configured to detect noise within the imaging band via one or more reference receive channels, optionally connected to a reference ultrasound transducer (not shown) that is acoustically isolated or inactivated such that it does not transduce reflected ultrasound waves but is sensitive to the noise received by the imaging transducer receive circuit.
  • the one or more imaging transducers 10 and the one or more reference transducers need not be oriented in a common spatial direction.
  • reference waveforms The signals received by the one or more reference transducer receive circuits or reference electrical receive circuits (on a reference receive channel) are henceforth referred to as reference waveforms.
  • a reference waveform is predominantly noise and not imaging energy.
  • reference waveform refers to a waveform obtained from one or more reference receive channels.
  • a reference waveform may be filtered to reside within the imaging band and/or outside the imaging band.
  • reference noise-detection waveform refers to a reference waveform obtained from a reference receive channel when the imaging transducer is receiving or is expected to be receiving imaging energy.
  • reference noise-characterization waveform refers to a reference waveform obtained from a reference receive channel when the imaging transducer is not receiving imaging energy.
  • the system may optionally be configured to suppress noise using a combination of detection-band waveforms and reference waveforms.
  • the phrase "noise-detection waveform” refers to either a reference waveform or a detection-band waveform.
  • a noise-detection waveform is referred to as a "noise- detection imaging waveform”.
  • the noise-detection waveform is referred to as a "noise-detection characterization waveform”.
  • FIG. 1 A shows a single transducer element
  • transducers with multiple piezoelectric elements may be employed.
  • a plurality of transducer elements may be controlled as a phased array or linear array or 2D array.
  • the transducer may not be limited to one that transmits imaging energy for the purpose of producing multi-dimensional 2D cross sectional images or 3D volumes (including 4D imaging datasets comprising 3D images over time), but may include transducers used for Doppler assessment of flow, transducers used as ultrasound beacons (e.g. as described in US Patent Publication No.
  • the transducer may not be limited to one that both transmits and receives imaging energy as shown in FIG. 1A, but may include transducers that receive ultrasonic energy from a medium that has been excited by other means, such as an optical energy (photoacoustic imaging), or by a separate ultrasound transducer. Further, although FIG.
  • the scan lines may be unidirectional, such as in M-mode imaging or during certain Doppler modalities, such as pulse-wave or continuous wave Doppler to assess flow.
  • a single transducer receive channel is configured to receive imaging energy within the imaging band, and to coincidentally also receive additional energy within one or more noise- detection bands, of which at least one comprises out-of-band noise.
  • one or more imaging transducer receive channels may be employed to receive imaging energy within the imaging band, and one or more transducer receive channels may be employed to receive additional energy within one or more noise-detection bands of which at least one comprises out-of-band noise.
  • one or more reference receive channels may be employed to receive noise energy (i.e. reference noise-detection waveforms) while being isolated from the imaging energy in the imaging band.
  • the reference receive channel may be filtered in a manner similar to the imaging transducer receive channel by using an imaging band pass filter.
  • the reference receive channel may not be filtered at all, or may have different filters other than an imaging band pass filter to better facilitate estimation of noise within the imaging band.
  • the control and processing hardware 100 may include, for example, one or more processors 1 10, memory 1 15, a system bus 105, one or more input/output devices 120, and a plurality of optional additional devices such as communications interface 135, data acquisition interface 140, display 125, and external storage 130.
  • control and processing hardware 100 may be provided as an external component that is interfaced to a processing device.
  • an optional transmit beamformer 26 and an optional receive beamformer 27 may be included as a component of control and processing hardware 100 (as shown within the dashed line), or may be provided as one or more external devices.
  • Some aspects of the present disclosure can be embodied, at least in part, in software, which, when executed on a computing system, configures the computing system as a specialty-purpose computing system that is capable of performing the signal processing and noise reduction methods disclosed herein, or variations thereof. That is, the techniques can be carried out in a computer system or other data processing system in response to its processor, such as a microprocessor, CPU or GPU, executing sequences of instructions contained in a memory, such as ROM, volatile RAM, non-volatile memory, cache, magnetic and optical disks, cloud processors, or other remote storage devices. Further, the instructions can be downloaded into a computing device over a data network, such as in a form of a compiled and linked version.
  • the logic to perform the processes as discussed above could be implemented in additional computer and/or machine readable media, such as discrete hardware components as large-scale integrated circuits (LSI's), application-specific integrated circuits (ASIC's), or firmware such as electrically erasable programmable read-only memory (EEPROM's) and field-programmable gate arrays (FPGAs).
  • LSI's large-scale integrated circuits
  • ASIC's application-specific integrated circuits
  • firmware such as electrically erasable programmable read-only memory (EEPROM's) and field-programmable gate arrays (FPGAs).
  • a computer readable medium can be used to store software and data which when executed by a data processing system causes the system to perform various methods.
  • the executable software and data can be stored in various places including for example ROM, volatile RAM, non-volatile memory and/or cache. Portions of this software and/or data can be stored in any one of these storage devices.
  • a machine-readable medium includes any mechanism that provides (i.e., stores and/or transmits) information in a form accessible by a machine (e.g., a computer, network device, personal digital assistant, manufacturing tool, any device with a set of one or more processors, etc.).
  • Examples of computer-readable media include but are not limited to recordable and non-recordable type media such as volatile and non-volatile memory devices, read only memory (ROM), random access memory (RAM), flash memory devices, floppy and other removable disks, magnetic disk storage media, optical storage media (e.g., compact discs (CDs), digital versatile disks (DVDs), etc.), network attached storage, cloud storage, among others.
  • the instructions can be embodied in digital and analog communication links for electrical, optical, acoustical or other forms of propagated signals, such as carrier waves, infrared signals, digital signals, and the like.
  • the phrases “computer readable material” and “computer readable storage medium” refer to all computer-readable media, except for a transitory propagating signal per se.
  • noise reduction filter based on noise sensed in the environment.
  • one or more of waveforms, data, filter parameters and other pertinent information described in the example embodiments below may be transmittable to a network and assessed remotely for further analysis and / or optimization of the noise reduction filter implementation. Once optimized, the noise reduction filter algorithms and / or parameters can then be transmitted to the system to enable improved noise reduction.
  • the example control and processing hardware 100 includes an imaging processing module 145 and a noise suppression module 150.
  • the image processing engine 145 may be configured or programmed to execute known image processing methods, such as scan conversion.
  • the noise reduction could occur in a post-processing fashion.
  • the data on a transducer receive channel or a reference receive channel could be digitized and stored before or after any filtering, envelope detection, shifting, shape / phase or delay adjustments, signal characterization, attenuation, subtraction or other steps in the described embodiments of the present invention.
  • FIG. 1 B illustrates an example of steps that may be employed by the control and processing hardware 100 and receive channel to process a detected an in-band imaging waveform from an imaging transducer receive channel prior to image generation.
  • the detected waveform from an imaging transducer receive circuit may be amplified 201 and filtered 202 prior to analog-to-digital conversion 203.
  • a band-pass filter 200 (which may employ multiple pass bands and stop bands) is employed to filter the detected waveform and retain the signal in the imaging band.
  • the envelope of the filtered waveform is then generated through an envelope detector 210.
  • the resulting envelope-detected waveform is then optionally decimated or expanded 220 and provided to the image processing module 230 for the generation of an image.
  • the example control and processing hardware 100 includes one or more noise suppression modules 150, which includes instructions for processing detected data (e.g. raw RF data, envelope data, or image data) to reduce a contribution of noise, according to noise reduction algorithms described in detail below.
  • the noise suppression algorithms disclosed herein may be employed to remove or reduce noise at several potential steps during the processing flow shown in FIG. 1 B, based on processing one or more noise-detection waveforms. In the case of systems using array transducers, noise suppression may occur either before or after beamforming (or both).
  • noise reduction of imaging data is performed based on measurements of signal energy (power, amplitude, intensity, or other measures of signal strength) or waveform patterns of the noise-detection waveform (such as an out-of-band noise-detection imaging waveform detected via an imaging receive channel or a reference noise-detection waveform detected on a reference receive channel).
  • a noise-detection waveform which in many embodiments is substantially absent of imaging energy, may contain noise that is correlated with or coincidental with the noise within the imaging band (or imaging bands).
  • One or more relationships between the imaging band noise and the noise detected by the noise-detection waveform may be employed for the correction of imaging band signals for the removal or reduction (e.g. suppression) of imaging band noise.
  • FIG. 1 C an example system is shown for ultrasound imaging of a region with an intracorporeal imaging probe 350 which connects via patient interface module (PI M) 300 to a control and processing hardware 100.
  • the intracorporeal ultrasound imaging device may be configured to receive acoustic imaging energy from a one-dimensional, two-dimensional or three-dimensional region, optionally via mechanical or electronic scanning.
  • the imaging probe 350 comprises an imaging assembly 353 remote from its proximal end with an electrical and/or optical channel 354 that passes through an optional conduit 354 along at least a portion of its length, and a connector 351 at its proximal end.
  • an imaging assembly 353 generally refers to a component or collection of components of the imaging probe 350 with which imaging energy (e.g.
  • the imaging assembly may optionally include one or more emitters of imaging energy, and includes at least one receiver of imaging energy.
  • the imaging assembly may contain an ultrasonic imaging transducer 10 that is both an emitter and receiver of acoustic energy.
  • the ultrasonic imaging transducer may be mounted on an imaging assembly that is optionally attached or otherwise connected to a rotatable conduit or shaft (e.g. a torque cable) 352 housed within a hollow sheath of an intracorporeal ultrasound imaging probe to facilitate mechanical scanning.
  • Optional PIM 300 facilitates transmission of signals within any wires or conduits to the appropriate image processing unit 100 via a PI M cable 320, such as when the imaging probe 350 does not connect directly to the control and processing hardware 100.
  • the PI M may incorporate one or more amplifiers 20 to amplify the signals from one or more transducer receive channels.
  • the PI M may optionally incorporate a motor drive unit 301 that imparts rotational motion to a rotatable conduit 354.
  • Motor drive unit 301 may include slip rings, rotary transformers or other components that couple the signals of probe 350 to control and processing hardware 100, thus allowing the imaging conduit to rotate while the PI M cable 320 does not.
  • the PIM 300 may also optionally incorporate a pullback mechanism 302 or a reciprocating push-pull mechanism to facilitate longitudinal translation of the imaging assembly 353.
  • Such longitudinal translation of the imaging assembly may occur in conjunction with the longitudinal translation of an external shaft that surrounds the imaging conduit, or may occur within a relatively stationary external shaft.
  • One or more reference receive circuits that detect noise correlated with the noise detected by an imaging receive circuit may be useful in suppressing in-band noise.
  • the following are example implementations of reference receive circuits for noise reduction of an imaging signal in an ultrasound imaging system.
  • FIG. 1 D illustrates an example embodiment in which the imaging transducer is replicated by one or more non-imaging reference transducers 361 that resides within the imaging probe 350.
  • the reference transducer has its own electrical channel 360 that passes through the optional conduit 352 and connector 351.
  • the reference transducer 361 may be coated with epoxy or with some other acoustic damping material 362 so that it is acoustically isolated from receiving imaging energy.
  • the piezoelectric may be de-poled to render the piezoelectric inactive, or may otherwise be substituted with a substrate that is not acoustically sensitive.
  • This implementation may be extended to array transducers, where there may be a plurality of ultrasound transducer elements that are configured to receive acoustic imaging energy.
  • FIG. 1 E illustrates an example embodiment in which the imaging probe contains two or more imaging transducers, each with a separate electrical channel, 370 and 354.
  • the two or more imaging transducers may be sensitive to receive acoustic imaging energy at substantially non-overlapping spectral bandwidths.
  • the first transducer may be configured to receive acoustic energy for frequencies around 10MHz
  • a second imaging transducer may be configured to receive imaging energy for frequencies around 40MHz.
  • the 40MHz band of the 10MHz-transducer may act as a reference noise channel for the 40MHz-transducer
  • the 10MHz band of the 40MHz-transducer may act as a reference noise channel for the 10MHz-transducer.
  • FIG. 1 F illustrates an example embodiment in which the imaging transducer receive channel is replicated by a reference receive circuit that resides within the imaging probe.
  • the reference receive circuit may optionally have some combination of resistors, inductors, capacitors and / or other components, configured such that electrical impedance of the reference electric circuit 381 is matched to the impedance of the imaging transducer receive circuit or such that the sensitivity of the reference electric circuit to noise is rendered more similar to the sensitivity of the transducer receive channel to noise.
  • the advantage of such embodiments is that they may be less expensive, easier to manufacture and easier to miniaturize some of the components of reference electric circuit 381 by not requiring an actual ultrasound transducer.
  • a portion of the reference receive circuit may serve additional purposes, such as transmitting energy to drive an actuator (such as, but not limited to, a magnetic actuator), or carrying a signal (including, but not limited to temperature, pressure or current generated from an electromagnetic field for position sensing).
  • an actuator such as, but not limited to, a magnetic actuator
  • a signal including, but not limited to temperature, pressure or current generated from an electromagnetic field for position sensing.
  • FIG. 1 G illustrates an example embodiment in which the reference receive circuit 391 terminates within the PI M where, similar to the imaging receive channel, it is subject to noise received by the motor drive unit and the PI M cable.
  • Reference noise detection may optionally be used in combination with detection-band (out-of-band or within-band) noise detection to further reduce noise of an in-band imaging waveform from an imaging transducer receive channel.
  • a detection-band waveform is employed as the noise-detection waveform
  • the noise-detection waveform may be less costly to produce than those that employ the use of a reference circuit or reference transducer to generate a reference waveform, as the former does not require the physical implementation of a reference channel, such as in cases where an imaging probe, or part thereof is not used repeatedly across different patients.
  • a noise-detection waveform from a reference channel may be more effective in some imaging systems at reducing noise, as it can provide information about noise that resides within the imaging band, whereas an out- of-band noise-detection waveform does not provide a direct estimate of the in- band noise and instead relies on noise whose in-band properties can be predicted based, at least in part, on its out-of-band properties. It is also noted that the noise estimates obtained using either a detection-band waveform or a reference noise-detection waveform may be used to reduce noise in more than one imaging channel.
  • a phased array transducer where there are a plurality of piezoelectric elements, it is possible to use a single reference receive channel or a single out-of-band noise-detection waveform to estimate noise that might be collected by all or a subset of the piezoelectric elements, and thus apply the same noise estimation scheme to the signals collected from all or a subset of the piezoelectric elements.
  • noise measurements and/or noise characteristics may be determined from the measurements in the noise-detection waveform in order to increase the signal-to-noise ratio within the imaging band.
  • noise measurements include any one or more of:
  • the estimation of noise characteristics may be performed when an imaging transducer is not receiving imaging energy (e.g. after ultrasound energy from the most recent emission of an ultrasound pulse is expected to have been extinguished from the
  • the estimation of noise characteristics may be performed during imaging, when imaging energy is expected to be detected (e.g. when the transducer is in image acquisition mode).
  • noise characteristics in which noise characteristics are measured in the absence of imaging energy, such noise characteristics may be updated intermittently in order to adapt to and compensate for time-dependent changes in the noise characteristics.
  • noise reduction embodiments involve the measurement and use of energy from a noise-detection waveform while the imaging transducers are receiving imaging energy
  • alternative embodiments may employ measurement from a noise-detection waveform that are obtained during time periods in which there is an absence of imaging energy, or a combination of both.
  • FIG. 5A illustrates an example of such an embodiment, in which noise characterization is performed while the transducer receive channel is not receiving imaging energy, and where the resulting noise characterization can be employed to enable noise reduction of in-band imaging waveforms acquired while the transducer receive channel is receiving imaging energy.
  • the noise characterization step would occur prior to the acquisition and processing of imaging data that is noise reduced, but with appropriate recording of the imaging data, the information gathered during noise characterization could be employed post-hoc on the recorded imaging data.
  • energy is detected in a transducer receive circuit during a first time window when it is either anticipated or known that at least one transducer receive circuit is not receiving imaging energy, such that the waveform detected by an imaging transducer receive channel 13 is deemed to be noise 405.
  • the detected waveform is filtered at 200 and 410 to generate an in-band noise
  • the in-band noise characterization waveform and the detection-band noise characterization waveform are processed to characterize their noise properties, as shown at 420.
  • the noise characterization 420 may be employed, for example, to generate characteristic parameters 430 that characterize the noise. Examples of suitable noise characterization parameters are provided in the forthcoming example embodiments.
  • energy may be detected in a transducer receive circuit during an additional baseline noise characterization stage when it is either anticipated or known that at least one transducer receive circuit is not receiving imaging energy, and a certain noise source is selectively known to be off, such that the waveform detected by an imaging transducer receive channel 13 is deemed to be baseline noise for a selected noise source.
  • the detected baseline noise characterization waveform is filtered at 200 and 410 to generate an in-band baseline noise characterization waveform 407 and a detection-band baseline noise characterization waveform.
  • noise parameters 430 may include parameters obtained during a baseline noise characterization stage.
  • Characteristic noise parameters may be calculated before or during an imaging session, or may be retrieved from a pre-stored database located on a local or remote storage drive (network drive, cloud, etc.). After having characterized the noise in the absence of imaging energy, the characteristic parameters 430 may be employed to perform noise reduction of the in-band imaging waveform 437 detected while the
  • transducers are receiving imaging energy.
  • the waveform detected from an imaging transducer receive channel during imaging 435, containing imaging energy and noise, is filtered at 200 and 410 to generate an in-band imaging waveform 437 and a detection-band imaging waveform 438. Therefore, the in- band imaging waveform 437 contains detected imaging energy and noise, and the detection-band imaging waveform 438 contains information associated with the likely presence of noise in the in-band imaging waveform 437.
  • noise characterization stage may then be employed for the detection and/or estimation of noise 440 within the in-band imaging waveform, and to perform noise suppression 500 of the in-band imaging waveform.
  • suitable noise characterization parameters are how they are used to reduce noise are provided in the forthcoming example embodiments.
  • FIG. 5B illustrates an alternative embodiment in which a reference receive channel, such as the reference receive channel described with reference to FIG. 1 A, is employed to detect noise and generate a reference noise characterization waveform 406.
  • This reference noise characterization waveform is filtered at 202 to generate a filtered reference noise
  • the reference channel filter could be an imaging band pass filter.
  • the reference channel filter may be different from an imaging band pass filter.
  • the reference noise characterization waveform 409, and optionally the in-band noise characterization waveform 407, are processed at step 420 to provide noise characterization parameters 430.
  • a reference waveform 436 is detected and optionally filtered to generate a filtered reference noise-detection waveform 439.
  • the characteristic parameters 430 obtained during the noise characterization stage may then be employed for the detection and/or estimation of noise 440 within the in-band imaging waveform 437, and to perform noise suppression 500 of the in-band imaging waveform 437.
  • both the detection-band imaging waveform 438 (shown in FIG. 5A) and the filtered reference noise-detection waveform 439 are processed to provide information about the likely presence of noise in the in-band imaging waveform 437.
  • noise suppression may be achieved by processing the in-band imaging waveform using one or more of the following methods: subtracting the estimated noise from the signal in the imaging band; attenuation of the estimated noise energy by multiplying the signal in the imaging band with an attenuation factor; and filtering the signal in the imaging band.
  • a subtrahend value may be proportional to the amount of power detected within the out-of-band noise-detection imaging waveform.
  • an attenuation factor may be inversely proportional to a measure associated with the amount of noise in the in-band imaging waveform, thereby attenuating portions of the in-band imaging waveform associated with noise.
  • Noise characterization may be performed once, or may alternatively be performed multiple times, or continuously.
  • noise characterization may be performed intermittently over time (e.g. at periodic or aperiodic intervals) in order to adapt to, and to compensate for, time-dependent changes in the noise characteristics.
  • the noise characterization waveforms collected for noise characterization may digitized and collected in multiple discrete arrays, such as arrays that are each long enough to store imaging data along a single scan line in ultrasound imaging, or they may be collected in a more continuous fashion as one or more data streams that get stored into a large array, a circular buffer or other data structure.
  • noise characterization may be user initiated (e.g. by pressing a button), for example, at the start of an imaging session or when the user observes or suspects a degradation in image quality.
  • noise characterization may be triggered, such as either automatically or after prompting the user, when the absence of imaging energy is detected.
  • periods of absence of imaging energy may be detected when the relative energy between an in- band imaging waveform and a noise-detection imaging waveform is unchanged for a prescribed period of time and lies within a predefined range.
  • Periods of absence of imaging energy may also be determined, for example, when the energy in the in-band imaging waveform after noise correction at 501 is below a predefined threshold, indicating the absence of imaging energy.
  • the noise characterization step 420 may also be useful in alerting the user or the system that the noise profile has changed in a manner that may cause the noise reduction algorithm to adversely affect the system (if a new isolated source of noise is detected, the noise suppression module 500 may erroneously cause suppression of the in-band signal or be less effective at suppressing in-band noise).
  • a noise-detection waveform e.g. the out-of-band noise-detection waveform 438 of FIG. 5A or a reference noise-detection waveform 439 of FIG. 5B
  • an optional noise monitoring module could be employed that monitors the characteristics (such as peak energy, power, frequency content slope, skew, kurtosis, histogram or some other metrics) of the noise-detection waveform. If the characteristics of a noise-detection waveform change (e.g. if the peak energy exceeds a threshold value), the noise monitoring module can communicate with other parts of the system (such as via a message, interrupt, alarm or other) to alert that the noise content has changed.
  • characteristics such as peak energy, power, frequency content slope, skew, kurtosis, histogram or some other metrics
  • an error value may be evaluated in a noise characterization stage, where noise suppression is performed on the in-band imaging waveform in the absence of imaging energy and the error value is the energy of the in-band imaging waveform after noise correction. If the error value exceeds a pre-defined threshold, an alert is generated. An alert may prompt a re-characterization of the noise, or the system may choose to ignore one or more out-of-band noise-detection bands or one or more reference receive channels in its noise removal algorithm.
  • Noise characterization 420 may also be useful in determining the noise sources in the environment. Noise sources may be determined, for example, by a pattern recognizer, such as one described in step 570 of FIG. 7A
  • Information on the noise sources can be used, for example, in order to access a database (local or networked) to select parameters for noise suppression or in order to determine the sequence of noise reduction methods to be used. For example, it may be preferably to remove periodic noise first (as described in detail in Embodiment 6Fig 8A-D), followed by less periodic noise.
  • noise characterization may be useful at detecting the type of electroanatomic mapping system being used during an ablation procedure, or detecting the activation and deactivation of an ablation catheter, such as a catheter that uses radiofrequency energy to perform ablation to treat arrhythmias. This may be achieved, for example, by a pattern
  • the one or more noise-detection waveforms may further be capable of detecting the duration or the relative intensity or frequency of the ablation energy being applied. Such information could be useful to an intracardiac imaging system, as it may facilitate annotation of an imaging dataset with information about when a noise source, such as an ablation catheter, was activated during a procedure.
  • Noise suppression may be performed once, multiple times, intermittently or continuously. For example, noise suppression may be user initiated. Alternatively, noise suppression may be performed intermittently over time (e.g. at periodic or aperiodic intervals) in order to compensate for time-dependent noise source. In another example, noise suppression may be controlled by the external device that emanates the noise. For example, noise suppression may be enabled or disabled by the controls of an RF ablation generator, such that noise suppression is performed when RF energy is being delivered.
  • Embodiment 1 Noise Reduction Based on Suppressing Envelope Detected Out-of-Band Noise with Optional Amplitude, Shape and Delay Correction
  • a suppression operator may be a subtractor that subtracts the estimated noise from the in-band imaging waveform 437.
  • the suppression operator may be an attenuator that attenuates the in-band imaging waveform with an attenuation factor that is derived from the estimated noise.
  • One or more transducer receive channels are employed to detect imaging energy, where the energy detected includes both the imaging band and the noise-detection band.
  • the waveforms can be digitally sampled, split (or copied) and filtered thereby obtaining a sampled in-band imaging waveform 437 and a sampled out-of-band noise-detection imaging waveform 438.
  • the sampled waveforms may be detected as a set of samples that are received in a time window (listening window). For example, in the case of ultrasound imaging, the listening window may occur immediately or shortly after pulsing an ultrasound transducer such that it emits energy into the adjacent environment. Pulsing could correspond to sending out one or more pulses.
  • the imaging waveform 435 is filtered (digitally or analog) using an imaging band pass filter 200 and a noise-detection-band filter 410 that spans frequencies outside the imaging band. Envelope detection may then be performed on the in-band imaging waveform 437 and out-of-band noise-detection imaging waveform 438, as shown at 210 and 41 1 , respectively.
  • the out-of-band noise- detection imaging waveform 438 is employed to reduce noise of the in-band imaging waveform 437.
  • the amplitude of the envelope- detected out-of-band noise-detection imaging waveform 438 is optionally scaled via an amplitude adjustment factor, as shown at 510, in order to compensate for differences in the noise power within the imaging band and the noise-detection band.
  • the amplitude adjustment factor may be determined based on the power spectrum of the noise, as determined in the absence of imaging energy, i.e. in a noise characterization stage.
  • the amplitude adjustment factor may be selected, or modified, by an operator, in order to provide a desired level of noise reduction, or determined after cross correlation at 580 (described below).
  • the envelope-detected out-of-band imaging waveform may be temporally dilated, compressed or shape-adjusted using some other linear or non-linear temporal scaling function, as shown at 510, in order to compensate for differences in the shape of the noise waveforms between the imaging band and the noise-detection band.
  • the two bandpass filters 200 and 410 may not transform the input waveforms similarly. Either the band pass filters or the properties of the noise itself may result in an offset of the noise as it propagates through the band pass filters. In the absence of a delay correction, the noise may be erroneously shifted prior to suppression, which can negatively impact the noise reduced signal.
  • delay adjustment may be achieved by calculating a cross-correlation between the in-band imaging waveform and the out-of-band imaging waveform, and aligning the waveforms at the point where the cross-correlation is maximum.
  • cross-correlation can be employed to determine a time delay correction value for correcting the relative temporal misalignment of the envelopes of the in-band imaging waveform and the out-of-band imaging waveform.
  • the time delay correction value and the amplitude correction value may be calculated using a plurality of sampled in-band imaging waveforms and (co-incidental) sampled out-of-band noise-detection imaging waveforms (which may be referred to as "arrays"), or in one or more time windows in a sampled in-band imaging waveform and (co-incidental) a sampled out-of-band noise-detection imaging waveform.
  • out-of-band noise may be
  • the dependence of the power on frequency among the plurality of noise-detection bands may be employed to select a suitable amplitude adjustment for estimating the noise power that is present in the imaging band for suppression at step 525.
  • the average noise power within multiple noise-detection bands may be fitted to a functional dependence on frequency, such as a linear fit, in order to estimate the noise power within the imaging band. This functional dependence on frequency may be determined in the absence of imaging energy i.e. in a noise characterization stage 420.
  • FIG. 2B illustrates an alternative example embodiment of a subtractive or signal-attenuating noise correction method in which the out-of-band noise- detection imaging waveform is frequency-shifted prior to delay and amplitude adjustment.
  • a frequency shift operation (for example, by multiplying by a complex exponential) is performed on an out-of-band noise- detection imaging waveform, shifting the spectrum of the out-of-band noise- detection imaging waveform so that it lies within, or overlaps with the imaging band.
  • the frequency shift operation may be performed such that the center frequency of the frequency-shifted out-of-band noise-detection imaging waveform coincides, or is approximately equal to a frequency within the imaging band where it is anticipated or known that a portion of the in-band noise resides.
  • Frequency shifting may be advantageous over the envelope detection embodiment illustrated in FIG. 2A, because frequency shifting may result in better correlation between the noise in the imaging band and the noise-detection band, which may result in better noise suppression.
  • envelope detection 210 may be applied to the output signal 501 after it has undergone noise reduction to obtain a noise-reduced signal envelope 520.
  • FIG. 2C an example method is illustrated in which in- band noise, detected via a reference receive channel (using a reference receive circuit) is employed to reduce noise in the in-band imaging waveform via a suppression operator (i.e. subtractor or attenuator).
  • a suppression operator i.e. subtractor or attenuator.
  • One or more imaging transducer receive channels are employed to receive imaging energy, and one or more reference receive channels are employed to receive noise energy (i.e. reference waveforms) that is anticipated to correlate with the noise energy received by the imaging transducer receive channel.
  • the waveforms may be digitally sampled, thereby obtaining sampled in-band imaging waveforms and sampled filtered reference noise-detection waveforms.
  • the noise suppression could be performed with analog electronics, such as by using an analog signal adder with the input into the adder from the reference receive circuit being inverted in the delay, scale and shape adjustment block 510, thus resulting in subtraction of the estimated noise.
  • the suppression can be embodied as an amplifier with a time-varying gain, wherein the gain is modulated by the noise detected in the reference receive channel.
  • the input waveforms from an imaging transducer receive channel and the reference receive channel are filtered (digitally or analog) using an imaging band pass filter 200, and optional reference channel filter 202 thereby providing the in-band imaging waveform and a filtered reference noise-detection waveform, respectively.
  • reference channel filter 202 may be similar to the imaging band pass filter 200.
  • Envelope detection is then optionally performed on the filtered signals, as shown at 210 and 41 1 .
  • the reference noise-detection waveform, measured by the reference receive channel is employed to reduce noise in the in-band imaging waveform.
  • the amplitude of the filtered reference noise- detection waveform (or its envelope) is optionally scaled via an amplitude adjustment factor, as shown at 510, in order to compensate for differences in the noise power within the filtered reference noise-detection waveform and the in-band imaging waveform.
  • the amplitude adjustment factor may be determined based on the power spectrum of the noise, as determined in the absence of imaging energy i.e. in a noise characterization stage.
  • the amplitude adjustment factor may be selected by an operator in order to provide a desired level of noise reduction or determined after cross correlation at 580.
  • delay adjustment may be achieved by calculating a 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.
  • the time delay correction value and the amplitude correction value may be calculated using a plurality of in-band imaging waveforms and (co-incidental) reference noise-detection waveforms or one or more windows of an in-band imaging waveform and a (co-incidental) reference noise-detection waveform.
  • Embodiment 2 Noise Reduction using reference noise-detection waveform from a reference receive channel as Input to Adaptive Filter
  • FIG. 3A illustrates an example embodiment of a noise correction method in which an adaptive filter is employed, in an active noise control (ANC) scheme, by applying a noise reducing correction to the in-band imaging waveform based on a reference noise-detection waveform, where the reference noise-detection waveform is correlated with the noise detected by an imaging transducer receive circuit.
  • the reference noise-detection waveform is filtered at 202 and the waveform from the imaging transducer receive channel is filtered at 200.
  • the in-band imaging waveform and the reference noise-detection waveform are filtered within the same band (e.g. 7-13 MHz for an exemplary intracardiac
  • An adaptive filter is a linear filter that has a transfer function controlled by variable parameters and a means to adjust those parameters according to an optimization algorithm.
  • Adaptive filters are typically digital finite-impulse- response (FIR) or infinite-impulse-response (M R) filters.
  • An active noise control (ANC) scheme is provided for the primary input which receives a signal (S) from the signal source that is corrupted by the presence of a noise (N) that is uncorrelated with the signal.
  • the reference input receives noise (N r ) that is uncorrelated with the signal but is correlated in some way with the primary input noise (N).
  • the reference noise passes through an adaptive filter to produce an output noise (N est imat e ) that is an estimate of the primary input noise (N).
  • the noise estimate is subtracted from the corrupted signal to produce an estimate of the noise reduced signal (S est imat e ).
  • the adaptive filter actively adjusts its coefficients to minimize the output power E[s est imate 2 ]. Since the signal S is uncorrelated with N and N r , while noise N is correlated with noise N r , minimizing the total output power maximizes the signal-to noise ratio.
  • Minimization algorithms such as a stochastic Least Mean Squares (LMS) algorithm or the deterministic Recursive Least Squares (RLS) algorithm may be used to find filter coefficients that minimize the output noise power.
  • the reference noise is measured via a reference receive channel using a reference receive circuit that is isolated from the imaging energy.
  • the primary input is obtained by applying an imaging band pass filter 200 to an input waveform from an imaging transducer receive channel 13.
  • a reference input N r is obtained by applying a reference channel filter 202 to the reference noise-detection waveform obtained via a reference receive channel.
  • FIGS. 3B, and 3C illustrate example embodiments of noise correction methods in which an adaptive filter 540 is employed to apply a noise reducing correction to the in-band imaging waveform using an out-of-band noise-detection imaging waveform.
  • the primary input is obtained by applying an imaging band pass filter
  • an out-of-band noise-detection imaging waveform is obtained by applying a detection band filter 410 to an input waveform from an imaging transducer receive channel.
  • the out-of-band noise-detection imaging waveform is demodulated via envelope detection 41 1 (in a manner similar to the embodiment shown in FIG. 2A) to obtain a reference input for ANC.
  • the out-of-band imaging waveform is frequency-shifted (at 530) to the imaging band (e.g. 7-13 MHz) in a manner similar to the embodiment shown in FIG. 2B and filtered using an imaging band pass filter 203 to obtain a reference input for ANC.
  • the imaging band e.g. 7-13 MHz
  • Embodiment 3 Noise Reduction Based on Frequency Shift, using Detection-Band Waveform as Input to Variable Filter of In-Band
  • FIG. 4 illustrates an example embodiment of a noise correction method in which a dynamic filter 550 is employed to filter the in-band imaging waveform, where the dynamic filter is controlled by a filter update algorithm 560 that updates filter coefficients after processing an out-of-band noise- detection imaging waveform that includes out-of-band noise, and optionally, a within-band noise detection imaging waveform that includes noise within the all or part of the imaging band.
  • a filter update algorithm 560 that updates filter coefficients after processing an out-of-band noise- detection imaging waveform that includes out-of-band noise, and optionally, a within-band noise detection imaging waveform that includes noise within the all or part of the imaging band.
  • the input waveform is separately filtered with an imaging band pass filter 200 and a noise-detection band-pass filter 410, thereby generating an in-band imaging waveform 437 and at least one out-of- band noise-detection imaging waveform that includes out-of-band noise 438.
  • the one or more out-of-band noise-detection imaging waveforms are processed by a filter-update algorithm at 560.
  • the filter-update algorithm analyzes the out-of-band noise and may evaluate signal characteristics, such as by performing a Fourier transform on a waveform array and identifying spectral maxima and the frequencies at which they occur.
  • the filter update algorithm may use this information to control the coefficients of a dynamic digital filter that filters the imaging waveform at 550.
  • the present method may be employed to reduce noise in a signal containing harmonic noise.
  • harmonic noise may be generated from a switching rectangular pulse source where, in frequency domain, the spacing of spectral lines is dependent on the pulse repetition frequency. If the pulse-width-modulation source generates noise that spans the 3-40MHz band, where the imaging band lies within 7-13 MHz, and the dynamic filter 550 is a comb or multiple-notch filter, the filter update algorithm 560 may process signals from the noise-detection band (e.g.
  • filter update algorithm 560 and / or dynamic filter 550 may optionally also probe the a within-band noise-detection imaging waveform to confirm the presence of in-band noise at one or more selected sub-bands within the imaging band prior to removing or reducing noise. For example, if the imaging band lies within 7-13 MHz and the filter update algorithm recognizes that there is harmonic noise at 15 MHz, 18 MHz, 21 MHz, 24 MHz and 24 MHz (integer multiples in the 15-25 MHz range), then the filter update algorithm may set the dynamic filter to filter out signals at 9 MHz and 12 MHz (integer multiples of 3MHz within the imaging band).
  • such a filter may optionally only be applied if the presence of signal at 9 and 12 MHz is greater than expected relative to other signals within the imaging band are confirmed.
  • a noise characterization step may be performed, in the absence of imaging energy, to determine whether or not harmonic noise is present within the imaging band.
  • Embodiment 4 Noise Reduction Based on Pattern Recognition
  • FIG. 7A and FIG. 7B illustrate an example embodiment of a noise correction method in which pattern recognition is employed to detect noise and to perform noise reduction on in-band imaging waveforms.
  • matched sets of patterns in detection-band noise- characterization waveforms (at least one of which is lies, at least in part, outside the imaging band) and related in-band noise-characterization waveforms are initially identified during a noise characterization period in the absence of receiving imaging energy, as shown in FIG 7A. Having correlated noise-detection patterns with in-band noise patterns via the noise
  • these correlations may be employed during imaging to perform noise reduction of in-band imaging waveforms, based on the identification of patterns in the one or more detection-band imaging waveforms.
  • a noise characterization stage energy is detected within both an imaging band (in- band) and a noise-detection band in the absence of imaging energy (such as during a non-imaging noise characterization stage), thereby obtaining correlated measurements of in-band and detection-band noise- characterization waveforms.
  • Samples from the in-band noise-characterization waveform and detection-band noise-characterization waveform are recorded as pairs of arrays, where an array pair refers to a sampled in-band waveform and a second corresponding sampled detection-band noise-detection waveform recorded at the same time.
  • the detection-band noise-characterization array and the in-band noise characterization array may be windowed at 566 and 565.
  • the windows may be sliding windows, with optional overlap.
  • the windows may be centered around a peak noise amplitude, or be time-locked to a noise amplitude threshold.
  • the windows may also be conditioned to reduce artifacts induced by windowing, such as by applying a window function, such as a Hamming window, Blackman window or other window functions well known in the art of signal processing.
  • the array data (or windows thereof) may be processed to identify the presence of one or more noise patterns at 570.
  • one or more detection-band noise- characterization arrays are processed to identify waveform patterns associated with patterns in the in-band noise-characterization waveforms.
  • the pattern recognizer in step 570 extracts features from the detection-band characterization arrays and uses a predictive model to classify the features into noise 'classes'.
  • the extracted features could be statistical features (including, but not limited to, variance, standard deviation, power, skewness, and kurtosis) in time domain, frequency domain (e.g. peak frequency), time- frequency domain (e.g. wavelet coefficients).
  • the choice of features to be extracted may be made beforehand using feature selection algorithms such as forward selection or backward elimination methods.
  • Extracted features are fed to the predictive model in step 570, which may be trained to identify a pattern in a detection-band waveform and assign the pattern to a noise class.
  • machine learning methods may be used to train the predictive model to recognize patterns in the detection-band noise-characterization array using the extracted features.
  • the predictive model may include an unsupervised learning model (such as k-means clustering), or a supervised learning model (such as a linear classifier, an artificial neural network or a nearest-neighbor classifier). Supervised learning may be used if prior information about noise sources are known, for example, the sources and sequence of the noise patterns may be known beforehand and noise class labels may be assigned to a waveform pattern in the detection-band noise-characterization waveform.
  • the predictive model in step 570 may also accept as input class weights or a priori probabilities. The higher the a priori probability or weight of a class, the more likely it is to be recognized.
  • a database shown in step 575 may store in-band noise- characterization waveform patterns that are known to be co-incidental with detection-band noise-characterization waveform patterns.
  • the database may store exemplary or average temporal in-band noise- characterization waveform patterns, paired with features of the co-incidental detection-band noise characterization patterns and noise class labels.
  • the detection-band noise-characterization waveform patterns and their coincidental in-band noise-characterization waveform patterns may be determined on a per-window basis or otherwise. Further noise
  • characterization may be performed in the temporal domain, or in the spatial domain after image generation at steps 230 and 231 , in which case spatial features may also be extracted in step 570.
  • the database may be of any suitable format used in computing, such a lookup tables.
  • these correlations may be employed during imaging to perform noise reduction of in- band imaging waveforms, based on the identification of patterns in the one or more detection-band imaging waveforms.
  • an in-band imaging waveform is obtained by applying an imaging-band bandpass filter 200 to a waveform detected by an imaging transducer receive channel while the transducer is receiving imaging energy, and optionally performing envelope detection 210.
  • a detection-band imaging waveform is obtained by applying a noise-detection bandpass filter 410 to a waveform detected by an imaging transducer receive channel, and optionally performing envelope detection 41 1 .
  • the detection-band imaging waveforms and the in-band imaging waveforms may be windowed at 566 and 565, similar to the windowing step in the noise characterization stage.
  • the waveforms may be sampled and represented as arrays.
  • Features may be extracted from the detection-band imaging array (or windows thereof), similar to the feature extraction step in the noise- characterization stage. Extracted features, and optionally class weights, are employed by the pattern recognizer 570, trained in the noise-characterization stage (described above), to identify the presence of one or more patterns in the detection-band imaging waveform.
  • the period (i.e. repetition frequency) of a noise pattern may be used to adjust the a priori probability of that pattern class (e.g. for a Bayes classifier) or to adjust the weight of that class (e.g. for a Support Vector Machine) while applying a pattern classification algorithm in 570.
  • the repetition interval of a pattern is known, the pattern is expected to be present at given times with a higher probability.
  • the a priori probability or weight of that class could be adjusted to be higher at those times, increasing the likelihood that the pattern classifier will recognize that noise pattern.
  • This repetition interval may be determined in the noise characterization stage and stored in a database 578, or may be loaded from a pre-stored database (local, networked, cloud storage).
  • Patterns in the detection-band imaging waveform that are identified in step 570 as being associated with one or more noise classes are then employed to generate noise corrections to the in-band imaging waveform (e.g. the in-band imaging array). These corrections may be generated based on finding a correlated in-band pattern in step 575, where matched sets of features of detection-band noise-characterization waveform patterns, in-band noise-characterization waveform patterns and noise class labels are stored in a searchable database or other classification scheme.
  • an in-band noise correction is generated on a per-window basis, and subtracted from the in- band imaging waveform at 525, on a per-window basis, optionally after a delay and/or amplitude adjustment and/or shape adjustment 510 that temporally aligns the in-band noise pattern retrieved from the database 575 or other classification scheme with the in-band waveform.
  • the detection-band noise- characterization array is first processed to extract one or more features which are then stored.
  • the temporal intervals at which a given noise pattern is detected may also be determined at this stage and stored in 578.
  • the corresponding correlated temporal pattern in the in-band noise- characterization array is also stored in 575 (for example, in a look-up-table).
  • the detection-band imaging arrays for one or more detection-band imaging waveforms predominantly contain noise, and are processed via the same feature extraction process.
  • a weight vector which assigns weights (or a priori probabilities) for each class of noise patterns, may optionally be obtained.
  • the repetition frequency of each pattern may be loaded from the database 578 created in a noise characterization stage or from a pre-stored database.
  • the weight for each class may be adjusted dynamically so that it is dependent on that pattern's repetition frequency, the time instance when that pattern was previously detected, and the certainty with which that pattern was previously detected.
  • features extracted from the detection-band imaging waveform, and optionally class weights are again fed to the trained predictive model (trained in the noise- characterization stage), which may identify a noise pattern in the detection- band imaging waveform and assign it a class.
  • a corresponding and correlated in-band noise pattern for the noise class is then obtained from the database in 575, where, for example, in-band noise waveform patterns and features of detection-band waveform patterns for each noise class may be stored in during the noise characterization stage (e.g. in-band temporal waveforms stored in the look up table).
  • the in-band noise waveform pattern extracted from the class comparison could be, for example, an average of all coincidental in-band noise patterns for the current noise class, or the in-band noise pattern whose co-incidental detection-band pattern features are closest to the features of the current detection-band imaging array, for example, determined through a nearest-neighbour calculation.
  • This co-incidental pattern is then subtracted from the input after amplitude and delay adjustment to obtain a noise reduced in-band imaging waveform.
  • FIG. 7A and FIG. 7B show an example implementation in which a single detection-band waveform (including energy residing, at least in part, beyond the imaging band) is generated by a single detection band pass filter 410 from which input to pattern recognizer 570 is derived.
  • multiple detection-band waveforms may be generated by multiple detection band pass filters, of which at least one detection-band waveform is out-of- band.
  • one or more detection-band waveforms may carry noise within all or a portion of the imaging band (i.e. within-band noise-detection waveforms).
  • within-band data may be useful for the pattern recognizer 570 to confirm that the noise predicted by the out-of-band noise-detection waveforms in fact exists in the imaging band (either during the noise characterization stage, or during imaging).
  • within-band noise-detection imaging waveforms may also be employed by the pattern recognizer to identify a noise source.
  • a noise source associated with the imaging sub-band may be identified.
  • a peak filter centered at 8MHz may be used to obtain a within-band noise detection waveform within an imaging band ranging from 7 to 13 MHz, and a detection band-pass filter with a pass band of 15 to 25MHz may be used to obtain an out-of-band noise-detection waveform.
  • the pattern recognizer may be able to adjust its weights to preferentially detect a specific noise source (i.e. noise class).
  • the system may better select a correlated in- band noise pattern from the database at 573 to remove that 8 MHz peak than if it solely relied on information that was out-of-band to the imaging band.
  • image data e.g. B-mode image data
  • image data may be processed to determine spatial noise patterns instead of processing time domain (e.g. RF or envelope- detected) signals.
  • processing time domain e.g. RF or envelope- detected
  • FIG. 7B decimation (220 and 221 ) and B-mode image line generation (230 and 231 ) is performed prior to 570 and 575.
  • 2D imaging windows may be used to detect spatial patterns, such as for B-mode image data.
  • texture features may be extracted in the spatial domain (e.g. gray level co-occurrence matrices), frequency domain (e.g. Fourier spectrum measurements), or spatial frequency domain (e.g. energy of 2D wavelet coefficients).
  • FIG. 7D an alternative example embodiment is shown in which a reference noise-detection waveform is employed, instead of the detection-band imaging waveform of FIG. 7B, when performing pattern recognition during imaging.
  • a reference noise-characterization waveform may be employed during the initial pattern recognition stage that is performed in the absence of imaging signal.
  • the algorithms or schema described above, with reference to FIG. 7A and FIG. 7B, may thus be adapted to the present example embodiment by replacing the detection-band noise-characterization waveform (and associated array measures) with the reference noise-characterization waveform (shown in FIG. 7C), and replacing the detection-band imaging waveform (and associated array measures) with the reference noise-detection waveform (shown in FIG. 7D).
  • Embodiment 5 Noise Reduction Based on Relative Energy Measures
  • noise reduction is performed by selectively attenuating a windowed portion of an in-band imaging waveform, based on criteria that are assessed according to measurements from one or more detection-band imaging waveforms, at least one of which is an out-of- band imaging waveform.
  • Attenuating may refer, for example, to subtracting a derived subtrahend value from the envelope of the windowed in-band imaging waveform, and/or multiplying the windowed in-band imaging waveform or its envelope with an attenuation factor, where the subtrahend value and/or the attenuating factor are determined from measurements on the noise-detection imaging waveforms or reference noise-detection waveforms.
  • a noise characterization stage energy is detected within both an imaging band (in- band) and a noise-detection band in the absence of imaging energy (such as during a non-imaging noise characterization stage), thereby obtaining correlated measurements of in-band noise and detection-band noise.
  • At least one noise-detection band is out-of-band.
  • Samples from the in-band noise- characterization waveform and the detection-band noise-characterization waveform are recorded as pairs of arrays, where an array pair refers to a sampled in-band waveform and a second corresponding sampled detection- band waveform recorded at the same time.
  • the noise characterization stage may include an additional stage, referred to as a baseline noise characterization stage, when it is either known or anticipated that the imaging transducer receive circuit is not receiving imaging energy and not receiving noise energy.
  • a baseline noise characterization stage may be obtained by applying an imaging band pass filter 200 to an input waveform detected from an imaging transducer receive channel 13 in the absence of receiving imaging energy and in the absence of receiving noise energy, and optionally detecting an envelope of the filtered waveform at 210.
  • a detection-band baseline noise- characterization array may be obtained by applying a detection-band filter 410 to an input waveform from an imaging transducer receive channel 13, and optionally detecting an envelope of the filtered waveform at 41 1.
  • the in-band baseline noise characterization array and the detection-band baseline noise- characterization array, measured in the absence of imaging energy and absence of noise energy, are denoted as Gi and Gn, respectively.
  • a given array pair may optionally be segmented according to a plurality of time windows, as illustrated in FIG. 6A at 565 and 566, to obtain windowed array pairs, denoted as Gi w and Gn w .
  • the windows may be sliding windows, with optional overlap between adjacent windows.
  • the maximum power within a windowed in-band baseline noise-characterization array may be denoted as Ti.
  • the maximum power within a windowed out-of-band baseline noise characterization array may be denoted as Tn.
  • the noise characterization stage includes a stage when the imaging transducer receive circuit is not receiving imaging energy but is anticipated to receive noise energy.
  • an in-band noise characterization array 407 may be obtained by applying an imaging band pass filter 200 to an input waveform detected from an imaging transducer receive channel 13 in the absence of receiving imaging energy, and optionally detecting an envelope of the filtered waveform at 210.
  • a detection-band noise-characterization array may be obtained by applying a noise-detection band filter 410 to an input waveform from an imaging transducer receive channel 13, and optionally detecting an envelope of the filtered waveform at 41 1 .
  • the in-band noise characterization array 407, and the -detection-band noise characterization array 408, measured in the absence of imaging energy, are denoted as Ci and Cn, respectively (as shown in FIG. 6A).
  • a given array pair may optionally be segmented according to a plurality of time windows, as illustrated in FIG. 6A at 565 and 566, to obtain windowed array pairs, denoted as Ci w and Cn w .
  • the windows may be sliding windows, with optional overlap between adjacent windows.
  • the windows may be centered around a peak amplitude of the noise waveforms in one or more detection bands and / or in the imaging band.
  • the time windows may be time-locked to noise onsets determined when the amplitude of noise in one or more detection bands and / or the imaging band exceeds predefined thresholds.
  • the threshold may be proportional to parameters Tn and/or Ti obtained during a baseline noise characterization stage.
  • the in-band and detection-band noise-characterization array pairs Ci w and Cn w are processed to obtain one or more measures associated with the energy in the imaging band and a noise-detection band for each time window, in order to characterize the relative intensity of the noise within the two bands. For example, as shown in FIG. 6A, for each pair of windows Ciw and Cn w the power in the imaging band and the power in the noise-detection-band may be calculated at 570 and 572.
  • Noise-characterizing measurements may optionally be calculated only from select windows, where selection criteria may be assessed according to in-band energy measurements, and optionally detection-band energy measurements. For example, only windows for which the in-band power exceeds a predefined threshold may be selected for obtaining noise- characterizing measurements, as shown in FIG. 6B.
  • the threshold may be proportional to Ti obtained during a baseline noise characterization stage. In another example, only windows for which the out-of-band power exceeds a predefined threshold may be selected for obtaining noise-characterizing measurements.
  • the threshold may be derived from Tn obtained during a baseline noise characterization stage.
  • maximum and minimum values may refer to either upper and lower percentiles, or true maximum and minimum values.
  • 98 th and 2 nd percentiles values may be used instead of the maximum and minimum values.
  • Other statistical thresholds such as the 95 th and 5 th percentiles, 90 th and 10 th percentiles and 80 th and 20 th percentiles or others may be used to represent the maximum and minimum values for characterization purposes.
  • One or more noise measurements may be calculated from the per- window energy measurements in a noise characterization stage, and may be used to define a relationship between power in the imaging band and power in the detection band in the presence of a noise source.
  • one or more pairs of in-band and detection-band power values may be selected as inflection points for generating a piece-wise linear function to define the relationship between power in the imaging band and power in the detection band in the presence of a noise source, as shown in Fig 6B.
  • a piece-wise linear function defining the relationship between the in-band power and the detection-band power in the presence of a noise source may be generated based on maximum and minimum power values obtained in a noise characterization stage, as shown in FIG. 6B.
  • Minimum and maximum in-band power values from a noise characterization stage may be evaluated (e.g. as absolute maximum/minimum values or using statistical measures) and denoted as Pimin and Pimax, respectively.
  • a set of windows whose in-band power falls within a preselected range relative to Pimin e.g. within a percentile range
  • the minimum detection-band power may be selected as Pn m in.
  • a set of windows whose in-band power falls within a preselected range relative to Pimax may be identified, and, from among the identified set of windows, the minimum detection-band power may be selected as Pn ma x.
  • Power pairs (Pn m in, Pimin) and (Pnmax, Pimax) may be used for the fitting of a function defining an estimated relationship between in-band and detection-band power. It is to be understood that the example implementation is just one non-limiting example of selecting values of in-band and detection-band powers to provide suitable fitting points and/or a functional relationship between in-band and detection- band power, and other methods may be alternatively employed.
  • the ratio of the power in the in-band noise-characterization waveform to the power in the detection-band noise-characterization waveform may be calculated on a per-window basis, and the maximum ratio R off across a plurality of windows (an example use of this quantity is described below when determining whether or not to apply a noise correction during imaging) may be obtained.
  • f(Pn) may be a piece-wise linear function whose slope, intercept and/or inflection points are defined by the points (Pn m in , Pimin) and (Pn ma x , Pimax) calculated in a noise characterization stage.
  • f(Pn) may be a non-linear polynomial, or a combination of one or more linear or non-linear polynomials.
  • f() may be a set of values defined for one or more ranges of Pn values.
  • f(Pn) may be assigned a value Pi a for ai ⁇ Pn ⁇ a 2
  • f(Pn) may be assigned a value Pib for bi ⁇ Pn ⁇ b 2 and so forth, where [ai , a 2 ] and [bi , b 2 ], and so forth, are non- overlapping intervals of Pn.
  • a set of windows of the detection-band noise characterization array whose power lies between ai and a 2 may be identified, and the associated windows of the in-band noise characterization array may be employed to determine a value for Pi a for the interval [ai , a 2 ].
  • Pia may be a representative power value (such as maximum, mean, median or some other measure) calculated from in-band power measurements of all windows whose detection-band power, Pn, lies within the range [ai , a 2 ]..
  • Noise measurements obtained in a noise characterization stage may be employed to perform noise reduction of the in-band imaging waveform obtained during imaging, as shown in FIGS. 6C and FIG 6D.
  • the input waveforms filtered as shown at 200 and 410 to provide the in-band imaging waveform 407 and the detection-band imaging waveform 408, optionally after envelope detection at 210 and 41 1 .
  • the waveforms may be sampled to obtain in-band imaging arrays and detection-band imaging arrays.
  • the term "in-band imaging array” is employed to refer to a sampled in-band imaging waveform.
  • the term “detection-band imaging array” is employed to refer to a sampled detection- band imaging waveform.
  • a set of arrays may be recorded, where each array may be respectively associated with a given scan line.
  • a first in- band imaging array may be associated with a first scan line
  • a second in-band imaging array may be associated with a second scan line, and so forth.
  • the in-band and detection-band imaging arrays are denoted as Qie, and Qne, respectively, where ⁇ is an index identifying a given period of acquisition, such as one corresponding to a scan line.
  • the arrays may be windowed, as shown at 565 and 566 in FIG. 6C, such as using the same window properties as those employed to window the noise characterization arrays in the characterization stage.
  • the in-band imaging arrays and detection-band imaging arrays, temporally segmented according to the windows, are denoted as Qie.w and Qne.w, respectively, where the subscript w is an integer denoting the window number.
  • Qii,io refers to the 10th window portion of the in-band imaging array corresponding to the 1 st scan line.
  • power or another suitable energy measurement may be calculated, as shown at 570 and 572.
  • In-band and detection-band power values are calculated as P(Qie,w) and P(Qne, w ), respectively. These energy measurements may then be used for suppressing noise based on measurements obtained in a noise characterization stage.
  • noise may be suppressed by subtracting a power value from the envelope of the in-band imaging array Qie.w at 525.
  • Noise-detection band power P(Qne, w ) may be used to estimate a noise energy, Pi N e, w , within the in-band imaging array window Qie.w based on the function f() obtained in a noise characterization stage. For example, the noise energy Pi N e,w within Qie, w may be estimated at
  • the estimated in-band noise Pi N e,w may optionally be scaled by multiplication with a scaling factor ⁇ , where 0 ⁇ 1 .
  • the scaled estimated noise Pi N ,w may optionally be clamped below an upper limit, such as, but not limited to 0.8xP(Qie,w), 0.9xP(Qie,w) or 1xP(Qie,w) to obtain a subtrahend value.
  • the subtrahend value may be subtracted from the envelope of Qie.w at 525 to obtain a noise reduced waveform envelope.
  • noise may be suppressed by multiplying the elements of the in-band imaging array Qie.w by an attenuating factor at 526.
  • the in-band power in an imaging window Qie.w may be calculated as P(Qie, w ).
  • Noise-detection band power P(Qne,w) may be used to estimate a noise energy, Pi N e, w , within the in-band imaging window Qie.w based on the function f() obtained in a noise
  • Pi N e,w may optionally be scaled by multiplication with a scaling factor ⁇ , where 0 ⁇ 1.
  • An attenuating factor may be selected to be proportional to
  • scaling factor ⁇ may be selected to lie between zero and unity.
  • the determination of ⁇ may be dependent on the depth of the window within the waveform (thus corresponding to a depth within the imaged tissue).
  • the parameter ⁇ may optionally be user-controlled.
  • the attenuating factor may optionally be clamped below an upper limit, such as, but not limited to, unity, 0.95, 0.9, or 0.8.
  • the attenuating factor may additionally or alternatively be clamped above a lower limit, such as, but not limited to, 0, 0.01 , 0.05, or 0.1 .
  • the attenuating factor may be multiplied with the array Qie.w at 526 or its envelope to obtain a noise reduced array.
  • the system may optionally be configured to group array pairs Cn w and Ciw into one or more categories, referred to as classes.
  • one or more detection-band noise-characterization waveforms may be processed by a pattern recognizer in order to identify one or more classes of noise patterns and to assign a class to temporal windows associated with an identified noise pattern.
  • a set of windows of the detection-band noise-characterization waveform belonging to a class k may be selected and denoted as fw k).
  • In-band and out-of-band power measures for windows within set ⁇ w_k ⁇ may be used to derive a functional relationship fk() specific to the class k, using methods similar to those described in the preceding paragraphs.
  • the system may be configured to employ a pattern recognizer to identify noise patterns when the imaging transducer circuit is receiving imaging energy.
  • a pattern recognizer may be processed by a pattern recognizer in order to identify one or more classes of noise patterns and to assign a class to temporal windows associated with an identified noise pattern.
  • the estimated in-band noise power may be derived from a functional relationship f() specific to the identified class. For example, for windows belonging to class 1 , function fi() may be used to derive an estimated in-band noise (i.e .
  • noise in a given window Qie.w of the in-band imaging waveform may be reduced by subtracting a subtrahend value from the envelope of Qie.w, or by multiplying Qie.w with an attenuation factor, where the subtrahend value or the attenuation factor are derived from Pi N e,w.
  • noise classified by the pattern recognizer as class 1 may originate from an electroanatomic mapping system and noise classified as class 2 may originate from an ablation energy generator.
  • fi() could be used to estimate the in-band noise generated by the electroanatomic mapping system based on the power in the noise-detection band when the noise is recognized by the system to come from the mapping system
  • f2() could be used to estimate the in-band noise generated by the ablation generator based on the power in the noise- detection band when the noise is recognized by the system to come from the ablation generator.
  • systems that emanate noise may be monitored to determine which functional relationship to use for noise estimation and suppression.
  • the controls of an ablation generator may be monitored so that a binary gating signal is enabled when the ablation generator is actively generating energy (and associated noise). This gating signal may be used to determine the time periods when function f2() is to be used for in-band noise estimation.
  • Energy measurements may optionally be employed to estimate whether or not the in-band imaging array is likely to exhibit a low signal-to- noise ratio.
  • energy measurements may be employed to classify windows as to whether or not to apply a noise reduction correction (via subtraction or multiplication with an attenuating factor).
  • the decision about whether or not to apply a noise reduction correction for a given window, w can be made based on the power detected in the noise-detection band.
  • Noise correction is applied for Qie.w if P(Qne, w ) exceeds a predefined threshold.
  • the threshold may be obtained in a noise characterization stage.
  • the decision about whether or not to apply a noise reduction correction for a given window can be made based on ratio of the power in the imaging band relative to the power in the noise-detection band. For example, in a noise-characterization stage, for each window of the pairs of in-band noise-characterization arrays and detection-band noise- characterization arrays, the ratio of the power in the imaging band to the power in the detection band may be obtained.
  • the representative maximum ratio across all windows, denoted R 0ff may be used as a threshold to decide whether or not to apply noise reduction correction when the imaging transducer receive circuit is receiving imaging energy, as described below.
  • the decision about whether or not to apply a noise reduction correction for a given window, w, of an in-band imaging waveform can be made based on ratio of the power in the imaging band relative to the power in the noise- detection band, denoted as R 0n w .
  • R° n w is compared to yR 0 ", where ⁇ is a relaxation parameter and R 0ff is the
  • R° n w is determined to be greater than yR 0 ", then it is estimated that the signal- to-noise ratio of the in-band signal is sufficiently high and a noise reducing correction is not applied. Conversely, if R° n w is determined to be less than or equal to YR° ff , then it is estimated that the signal-to-noise is sufficiently low to warrant the application of a noise reduction correction.
  • the value of ⁇ may be employed to adjust the sensitivity to signal-to- noise, and may be used, in some cases, as an adjustment factor for cases in which weak portions of the signal would otherwise be suppressed. Lowering the value of ⁇ will lower the threshold for applying noise reduction, thereby reducing the number of windows that undergo noise reduction and allowing more imaging energy (and noise) to persist in the final output.
  • ⁇ may be selected to lie between zero and unity.
  • the determination of ⁇ may be dependent on the depth of the window within the waveform (thus
  • the parameter ⁇ may optionally be user-controlled.
  • the user can reduce the value of this parameter in order to lessen the effect of noise reduction.
  • any suitable noise reduction or suppression method may be employed to reduce noise of the in- band imaging array within the window. It will be understood that a wide variety of noise reduction corrections may be applied, such as, but not limited to, corrections involving subtraction as shown in FIG. 6C and/or multiplication with an attenuation factor as shown in FIG. 6D.
  • the noise suppression can cause erroneous noise reduction on a per-window basis.
  • some windows that contain a small amount of imaging energy may inadvertently undergo noise reduction based on an erroneous determination of a low signal-to-noise ratio in a window (i.e. false window classification). This can result in some small image "holes" in a surrounding homogenous signal region of an image, or residual noise pixels in a surrounding low-noise region of an image.
  • the status of adjacent windows i.e. windows in the spatial neighborhood
  • adjacent windows in one or more adjacent arrays may be employed to assess whether or not noise reduction by amplitude attenuation of the given window should be performed.
  • a given window is identified as not being suitable for noise reduction, yet one or more adjacent windows are identified as being suitable for noise reduction, the window may be flagged as being likely misclassified.
  • the status of the given window may be overridden and the given window may instead by identified as being suitable for noise reduction, such that noise reduction is applied to the given window.
  • a given window is identified as being suitable for noise reduction, yet adjacent windows within adjacent arrays are identified as not being suitable for noise reduction, then the initial determination of the status of the given window may be overridden such that the given window is instead identified as not being suitable for noise reduction and is flagged as a window that is likely misclassified.
  • Samples of misclassified windows may be replaced by samples of one or more non-noisy neighboring windows (i.e. replaced by or interpolated from samples of adjacent windows identified as not being suitable for noise reduction as per the aforementioned methods) optionally after performing delay and amplitude adjustments.
  • FIG. 6E An example implementation of this method is illustrated in FIG. 6E, where adjacent windows in adjacent arrays are interrogated to determine whether or not the classification of the current window is consistent with its surrounding windows in adjacent arrays. Since two windows (array 2, windows 4 and 5) classified as containing primarily noise (marked “N” ) are surrounded by windows (marked "S") classified as having a sufficiently high signal power to avoid the need for noise reduction, the two "N" windows may be replaced with samples from neighboring "S" windows (for example, by copying or by interpolating with amplitude adjustment and/or shape
  • window 4 of array 5 initially marked “S”
  • window 4 of array 5 may be reclassified as "N", shown in FIG. 6F, since preceding and proceeding windows of the array are marked "N".
  • the present example embodiment may be adapted to process image data as opposed to the processing of time-domain signals.
  • a plurality of in-band and detection-band imaging arrays representing a plurality of adjacent scan lines, may be acquired and post-processed to obtain in-band image and detection-band image frames.
  • the in-band and detection-band image pixels are denoted as Bie.d and Bne.d, respectively, where ⁇ is denotes a scan line and d is the depth.
  • the detection- band image may be used to evaluate attenuation (i.e. via subtraction or multiplication) values for each pixel in the in-band image.
  • Attenuation values may be obtained from a corresponding pixel in a detection-band image on a per-pixel basis (i.e. the value at Bne.d may be used to attenuate the pixel intensity at Bie.d) or by processing a region-of-interest in the detection-band image that would correspond to a local spatial neighborhood of an in-band image pixel (e.g. the values of a 3x3 neighborhood (e.g. in polar or Cartesian co-ordinates) around Bne.d may be used to attenuate the pixel intensity at Bie.d).
  • FIGS. 6G, 6H and 6J alternative example
  • FIG. 6G a reference receive channel is used when performing noise characterization (FIG. 6G) and for the determination of a suitable subtrahend value (FIG. 6H) or attenuation factor (FIG. 6J) for noise suppression during imaging.
  • a reference noise characterization waveform 409 is employed, instead of the detection-band noise
  • the reference channel filter could be an imaging band pass filter. Alternatively, if the noise estimation benefits from input of out-of-band noise, then the reference channel filter may be different from an imaging band pass filter.
  • a reference noise-detection waveform 439 is employed, instead of the detection-band imaging waveform of FIG. 6B, to determine subtrahend values and apply a noise reduction by subtraction.
  • a reference noise-detection waveform 439 is employed, instead of the detection-band imaging waveform of FIG. 6B, to determine and apply an attenuation factor for noise reduction.
  • 6A to 6F may thus be adapted to the present example embodiment by replacing the detection-band noise-characterization waveform 408 (and associated power measures) with the reference noise- characterization waveform 409 shown in FIG. 6G, and replacing the detection- band imaging waveform 438 (and associated array and power measures) with the reference noise-detection waveform 439 shown in FIG. 6H and FIG. 6I.
  • Embodiment 6 Noise Reduction of Pseudo-Periodic Noise Sources
  • noise reduction is performed by estimating and subtracting in-band noise, where the noise is expected to originate from a pseudo-periodic noise source or a pseudo-periodic sequence of noise sources.
  • the in-band noise is estimated based on measurements made during a noise characterization stage when an imaging transducer receive circuit is not receiving imaging energy.
  • waveforms are detected within both an imaging band and a noise-detection band in the absence of imaging energy (e.g. when an ultrasound transducer is not receiving imaging energy) and sampled, thereby obtaining a pair of co- incidental in-band and detection-band (out-of-band) noise characterization arrays.
  • an in-band noise-characterization waveform 590 may be obtained by applying an imaging band pass filter 200 to an input waveform detected from an imaging transducer receive channel 13, and optionally detecting an envelope 210 of the filtered data.
  • a detection-band noise-characterization waveform 595 may be obtained by applying a noise- detection bandpass filter 410 to an input waveform from an imaging transducer receive channel (where at least one noise-detection band comprises signal from outside of the imaging band), and optionally detecting an envelope 41 1 of the filtered data.
  • the in-band noise-characterization waveforms and detection-band noise characterization waveforms may be sampled to obtain in-band noise- characterization arrays and detection-band noise characterization arrays, denoted as Ci and Cn, respectively.
  • Ci (and Cn) should capture one or more periods of a periodic noise source.
  • the correlation between the detection-band characterization array and a detection-band imaging array may be used to estimate adjustment parameters for subtracting the in-band characterization array from the in-band imaging array. .
  • the input waveforms from an imaging transducer receive circuit 13 are filtered as shown at 200 and 410 to provide in-band imaging waveforms and detection-band imaging waveforms . These waveforms may be sampled (before or after performing envelope detection) to obtain an in-band imaging array 537 and a detection- band imaging array 538.
  • the in-band and detection-band imaging arrays are denoted as Qi and Qn, respectively.
  • each in-band imaging array is processed for noise reduction using the in-band noise-characterization array Ci to subtract noise from the in-band imaging array Qi.
  • the in-band noise characterization array Ci should be temporally aligned with the in-band imaging array Qi such that noise is co-incidental. Such alignment is possible in the case of a periodic noise source that generates noise in the imaging band that is correlated with noise in a noise-detection-band.
  • the temporal alignment may be achieved, for example, by segmenting the in-band and -detection-band imaging arrays Qi and Qn into a plurality of time windows at 565 and 566 (as described in the preceding example embodiments).
  • the windows should preferably be long enough to capture one or more periods of the periodic noise source.
  • the imaging arrays, temporally segmented according to the windows, are denoted as Qi w and Qn w , where the subscript w is an integer denoting the window number.
  • temporal alignment may be achieved on a per-window basis.
  • the temporal alignment may be achieved by calculating, within each window, the cross- correlation between the detection-band noise characterization array Cn and the detection-band imaging array Qn w , and selecting the relative time delay ⁇ corresponding to the maximum cross-correlation, as shown at 580. Due to the co-incidental relation between the noise in the imaging band and noise in the noise-detection band, this time delay ⁇ can also be applied to align the in- band noise characterization array Ci relative to the in-band imaging array Qiw, on a per-window basis as shown at 510. A scaling factor may also be applied to the aligned in-band noise characterization array.
  • a windowed portion of the aligned in-band noise characterization array which is denoted by Ci w is then subtracted from the in-band imaging array Qiw, resulting in a noise-reduced in-band imaging array, Qi w optionally after having taken the absolute value post-subtraction or applying a floor function in order to eliminate negative values. This process may then be repeated for each additional window for which noise reduction is desired.
  • a scaling factor may be applied when subtracting the aligned windowed segment of the in-band noise characterization array Ci w from the in-band imaging array Qi w .
  • FIGS. 8C and 8D an alternative example
  • FIG. 8C a reference noise characterization waveform 596 is employed, instead of the detection-band noise-characterization waveform 595 of FIG. 8A, when performing noise characterization.
  • FIG. 8D a filtered reference noise-detection waveform 439 is employed, instead of the detection-band imaging waveform of FIG. 8B, to determine and apply the amplitude adjustment.
  • the methods described above, with reference to FIGS. 8A and 8B, may thus be adapted to the present example embodiment by replacing the detection-band noise characterization waveform 408 with the reference noise-characterization waveform 409 shown in FIG. 8C, and replacing the detection-band imaging waveform 438 (and associated arrays ) with the reference noise-detection waveform 439 shown in FIG. 8D.
  • Embodiment 7 Noise Reduction using a Plurality of Scans by Changing Scan Rate
  • the system may be configured such that a set of two or more in-band imaging arrays associated with imaging energy from the same scan line or with scan lines with substantial spatial overlap are obtained.
  • the imaging energy within the set of in-band imaging arrays will have redundant temporal / depth dependency.
  • Averaging (or performing some other statistical processing, such as evaluating the minimum value) the set of redundant in- band imaging arrays may suppress noise if the noise itself is not time-locked to the trigger that prompts the imaging transducer receive circuit to start receiving imaging energy for each scan line. For example, if the pulse repetition frequency of the voltage pulse that excites an imaging ultrasound transducer is 200us, a periodic noise that repeats every 2us and will always have a component at Ous, 2us, 4us and so forth for each in-band imaging array.
  • a first imaging array will have noise components at Ous, 2us, 4us and so forth, and a second imaging array will have noise components at 1 us, 3us, 5us and so forth. Noise may be suppressed by averaging two successive redundant in- band imaging arrays, optionally after performing envelope detection.
  • noise in a detection-band waveform may be used to determine the period of one or more in-band noise sources.
  • the system may be prompted to adjust its scan rate so that imaging scan period is not an integer multiple of the period of a noise source.
  • an auto-correlation function may be used to detect periodicity in a detection-band imaging waveform while an imaging transducer receive circuit is receiving imaging energy.
  • the period of one or more noise sources may be determined from a detection-band noise-characterization waveform or reference noise-characterization waveform in a noise
  • characterization stage when the imaging transducer is not receiving imaging energy, or may be loaded from a pre-stored database. The system may then be prompted to adjust its scan rate so that imaging scan period is not an integer multiple of the period of a noise source.
  • in-band imaging waveforms obtained from a plurality of scan lines are sampled after performing envelope detection to obtain a set of in-band imaging arrays.
  • An in-band imaging array is denoted as Qie, where ⁇ is a scan line.
  • the system may be configured to suppress noise by averaging (or performing another statistical measurement, such as taking the minimum) across a plurality of arrays associated with adjacent scan lines or scan lines with significant spatial overlap.
  • a sample Qie[k] may be replaced by the average of [Qie-i[k], Qie[k], Qie+i[k]].
  • a sample Qie[k] may be selectively retained after performing some other numerical analysis on the set [Qie-i [k], Qie[k], Qie+i[k]] and determining whether the sample warrants noise reduction. For example, if for [Qie-i[k], Qie[k], Qie+i[k]] the minimum value is greater than half the maximum value, the range of sample values may not be large enough to warrant noise reduction by averaging and that the sample is therefore left unchanged.
  • Such a scheme may also be useful in MRI imaging in the presence of a periodic noise source, by ensuring that RF excitation pulses are initiated at times that do not correlate with the timing or periodicity of noise sources in the local environment. Additional feature: Selectively Choosing Samples to Undergo Noise Reduction when using a Plurality of Scans
  • the scan rate is adjusted as described above, and noise reduction is performed by selectively replacing portions of an in-band imaging array based on statistical measurements (such as average or minimum) from arrays from a plurality of adjacent scan lines, where there is sufficient overlap in the scan region associated with the adjacent scan lines. Only segments of the in-band imaging array that are assessed as noisy using detection-band measurements are replaced.
  • 410 to provide in-band waveforms and detection-band-waveforms may be sampled after performing envelope detection to obtain in-band imaging arrays and detection-band imaging arrays.
  • Data corresponding to a plurality of scan lines is recorded as sets of arrays for each scan line.
  • a set of arrays are obtained for each scan line, where one in-band imaging array is obtained in an imaging band for each scan line, and at least one out-of-band imaging array is obtained for each scan line.
  • the in-band and detection-band imaging arrays in a pair are denoted as Qie and Qne, respectively, where ⁇ is a scan line.
  • the detection-band imaging arrays are segmented into a plurality of windows, each containing J samples of Qne.
  • An array segmented according to windows is denoted as Qne.w, where the subscript w is an integer denoting the window number, ⁇ is an index denoting the scan line, and Qne.w contains samples [ Qn 0 [k w ], Qne[k w +1 ], ... Qne[k w +J-1 ], where k w is an index of the first sample in the window.
  • Each window Qne.w in the detection band is assessed for the presence of or absence of noise, as shown at 600. If noise (determined, for example, by measurements of waveform energy such as peak, RMS, etc.) exceeds a threshold, the window is deemed noisy.
  • the threshold may be selected in a noise characterization stage.
  • all the co-incidental in- band samples i.e. Qie[k w ], Qie[k w +i], ... Qie[k w +j-i] may be identified as samples that are suitable for noise reduction.
  • a given sample Qie[k] may be associated with more than one window. There may be instances when the sample is associated with both noisy and noise-free windows. In these cases, the system may be configured to pool noise assessments from multiple out-of- band windows before determining if the sample is suitable for noise reduction.
  • Samples deemed suitable for noise reduction may be replaced at 620 using statistical measures (such as duplicate values, average, minimum, and the like) from samples from arrays associated with adjacent scan lines, computed as shown at 610. For example, if sample Qie[k] is deemed suitable for noise reduction, and the system is configured to group arrays from 3 scan lines, sample Qie[k] may be replaced by the minimum of [Qie-i[k], Qie[k], Referring now to FIG. 8F, an alternative example embodiment is shown in which a reference waveform employed, instead of the detection-band imaging waveform of FIG. 8E, when performing noise reduction during imaging. The methods described above, with reference to FIG. 8D, may thus be adapted to the present example embodiment by replacing the detection- band imaging waveform (and associated arrays) with the filtered reference noise-detection waveform.
  • embodiments disclosed herein may employ frequency-domain or time- frequency domain processing during one or more steps.
  • steps 570 and 572 instead of using maximum power and ratio of maximum power in the imaging band and noise-detection band, short- term Fourier Transform or wavelet transforms can be performed on the in- band waveform and the out-of-band waveform on a per window or per array basis. Analysis of the transform coefficients (e.g. average, mean square, and the like) can then be used to characterize noise or detect windows when noise is present or absent.
  • the transform coefficients for the current window can be attenuated.
  • an inverse transform can be performed on the attenuated frequency-domain signal to obtain a noise reduced time-domain signal.
  • FIG. 2A step 510
  • FIGS. 8B and 8D step 580
  • Fourier transform algorithms may be used for efficient computation of cross-correlation.
  • machine learning algorithms may be used to classify noise patterns in step 570.
  • These patterns may be defined by frequency domain and/or time- frequency domain features, which will require frequency domain or time- frequency domain processing of time-series waveforms.
  • FIG. 9 an alternative example system is illustrated in which the signals that undergo noise reduction are obtained from a magnetic resonance (MR) system.
  • the example system includes a magnetic resonance scanner 50 that employs a main magnet 52 to produce a main magnetic field BO, which generates a polarization in a patient 60 or the examined subject.
  • the example system includes gradient coils 54 for generating magnetic field gradients.
  • a reception coil 58 detects the MR signals from patient 60.
  • the reception coil 58 can also be used as a transmission coil.
  • body coil 56 may be employed to radiate and/or detect radio frequency (RF) pulses.
  • the RF pulses are generated by an RF unit 65, and the magnetic field gradients are generated by a gradient unit 70.
  • the manner by which MR signals are detected using the sequence of RF pulses and magnetic field gradients, and how MR images are
  • a reference receive circuit may comprise a coil that is in the same room as the scanner, but is not located immediately adjacent to the imaged sample (such as a patient) from which the sought-after MRI signals are emitting. Electromagnetic noise that is traveling near the MRI machine will be detected by both the imaging reception coil 58 and the reference receive circuit.
  • the coil of a reference receive circuit might be oriented and positioned such that it is likely to receive some of the same noise as the imaging reception coil, but be distant enough from the imaged sample such that there is negligible imaging energy detected by the reference receive circuit.
  • the reference receive circuit coil may be tuned to have the same bandwidth as the imaging reception coil 58, or it may have a different bandwidth that is still able to collect noise signals in the environment that are correlated with noise that might be coupled into the imaging reception coil.
  • the reference receive circuit may further comprise a collection of receive coils, such as 3 coils whose alignments are orthogonal to each other. This would allow the collection of electromagnetic noise in a manner that a weighted sum of the noise collected in each of the 3 coils might more closely match the noise collected in the imaging reception coil, thus taking into account the directionality of dominant sources of electromagnetic noise in the MRI environment.
  • the MR system can have additional units or components that are not shown for clarity, such as, but not limited to, additional control or input devices, and additional sensing devices, such as devices for cardiac and/or respiratory gating.
  • additional control or input devices such as devices for cardiac and/or respiratory gating.
  • additional sensing devices such as devices for cardiac and/or respiratory gating.
  • the various units can be realized other than in the depicted separation of the individual units. It is possible that the different components are assembled into units or that different units are combined with one another.
  • Various units (depicted as functional units) can be designed as hardware, software or a combination of hardware and software.
  • control and processing hardware 100 obtains magnetic resonance images of patient 60 according to a suitable pulse sequence.
  • Control and processing hardware 100 is interfaced with magnetic resonance imaging scanner 50 for receiving acquired images and for controlling the acquisition of images.
  • Control and processing hardware 100 receives image data from RF unit 65 and processes the imaging data according to the methods described below.
  • Control and processing hardware 100 may be programmed with a set of instructions which when executed in the processor causes the system to perform one or more methods described in the disclosure in order to reduce noise in signals obtained from the magnetic resonance imaging system.
  • control and processing hardware 100 may be programmed with instructions in the form of a set of executable image processing modules, such as, but 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, image acquisition and image processing modules may be implemented using algorithms known to those skilled in the art for pulse sequence generation, image acquisition, and image reconstruction, respectively.
  • RF data is received from RF coils 56 and/or 58, and optionally one or more reference receive circuits.
  • Data may be sampled and filtered to obtain an in-band waveform.
  • a reference waveform via a reference receive circuit or a noise-detection-band waveform measured via filtering of RF from coils 56 and / or 58 are collected.
  • One or more noise suppression methods described in FIGS. 2-8 may be employed for noise suppression at 100.
  • the pulse generation module establishes the sequence of RF pulses and magnetic field gradients depending on the desired imaging sequence, and the image acquisition module stores the MR signals detected by the coils 56 and/or 58 in raw data space.
  • the image processing module 145 processes the acquired optionally noise-suppressed RF data to perform image reconstruction of an MR image.
  • the ability to estimate in-band noise and improve the imaging signal by removing the estimated noise from the imaging signal is provided.
  • Example 1 Noise Reduction of Unknown Noise Source via Attenuation Factor (example of Embodiment 5)
  • the present example involved the collection of ultrasound data using an intra-cardiac echo (ICE) system, in the presence of two noise sources.
  • the transducer was configured to detect ultrasound energy at frequency of 9 MHz.
  • Two band pass 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 is beyond the frequency range of the emitted ultrasound.
  • RF radio-frequency
  • the first noise source was an electroanatomic mapping system (Carto ® 3).
  • the system has an electromagnetic tracking module and an impedance- based tracking module, for which patches to measure impedance and estimate device position are attached to a patient. These patches can couple a significant amount of noise into the imaging band of the ICE images.
  • a heart phantom was used in a saline bath. The electrodes from the impedance patches were submerged in the bath.
  • the second noise source was from powering on an ablation generator that was connected to the Carto ® 3 console. The noise generated from this second noise source was determined to likely be noise propagating from the ablation generator, through the Carto ® 3 console and through the patch electrodes into the saline bath when the generator was powered on.
  • FIG. 10A shows an ultrasound image collected in the absence of either noise source
  • FIGS. 10B and 10C show the effect of the first and second noise sources, respectively, on the imaging quality.
  • Noise reduction of ultrasound waveforms detected by the ultrasound transducer of the ICE console was performed according to an implementation of the method illustrated in FIGS. 6A and 6D.
  • in-band baseline noise characterization array denoted as Gi.
  • the in-band baseline noise characterization array was obtained from sampling an in-band waveform for 125us at 200MS/s.
  • the noise-detection band was configured as the frequency band spanning approximately 15 to 25 MHz.
  • 512 waveforms (each 125us in duration sampled at 200MS/S) were obtained. Accordingly, 512 in-band and detection-band noise-characterization waveform pairs were employed to estimate noise characteristics.
  • Ciw and Cn w the power in the imaging band and the power in the noise-detection-band were calculated. Only windows whose in-band power was greater than threshold Ti were selected for further noise characterization.
  • the window with the near- minimum noise-detection band power (20 th percentile within the set) was chosen as w_min, and power values P(Cn w _min) and P(Ci w _min) were calculated and denoted as Pn m in and Pimin, respectively.
  • Pn w is the power in the detection band. It will be understood that the present implementation is provided to illustrate an example algorithm for selecting a functional relationship between the relative powers in the imaging and detection bands, and that other functional forms may alternatively be employed.
  • the function f() was employed to perform noise reduction of in-band imaging waveforms obtained during imaging (i.e. when the pulser was periodically emitting imaging energy, and ultrasound energy was being received by the transducer receive circuit).
  • Waveforms in the imaging band and the noise-detection band were sampled to obtain pairs of in-band imaging arrays and detection-band imaging arrays, denoted as Qie and Qne.
  • ⁇ is an index identifying a period of acquisition corresponding to a scan line.
  • the imaging arrays were temporally segmented according to the windows, denoted as Qie.w and Qne.w, where the subscript w is an integer denoting the window number.
  • powers P(Qie,w) and P(Qn 0 ,w) were calculated, and denoted as Pie.w and Pne.w respectively.
  • out-of-band imaging power Pne.w was used to estimate a noise energy, Pi N e, w , within the in-band imaging window Qie.w .
  • Pi N e,w was calculated as f(Pne,w), where ⁇ was set to unity and parameters m, c, Pnmax and Pn m in , were obtained in the noise characterization stage.
  • Pi N o,w was clamped below Pie.w.
  • An attenuating factor was calculated as [Pi e , w _ ⁇ ⁇ , ⁇ ]/ Pie,w ⁇ The attenuating factor was clamped above 0.02. The attenuating factor was multiplied with the array Qie.w to obtain a noise reduced array Qie.w. Qiefor all 512 scan lines were enveloped and processed by an image generator (step 230).
  • FIGS. 1 1 A and 1 1 B show the images obtained in the presence of the first noise source, without (A) and with (B) the application of the present noise reduction algorithm.
  • FIGS. 12A and 12B show the images obtained in the presence of the second noise source, without (A) and with (B) the application of the present noise reduction algorithm.
  • a clear improvement in the signal-to-noise ratio was observed in both cases.
  • FIGS. 1 1 A and 1 1 B show the images obtained in the presence of the first noise source, without (A) and with (B) the application of the present noise reduction algorithm.
  • 1 1 C, 1 1 D, and 1 1 E show the images obtained in the presence of the first noise source with the application of the present noise reduction method when the relaxation parameter ⁇ of function f() was set as 0.5 (C), unity (D) and 1 .5 (E).
  • noise reduction of ultrasound waveforms detected by the ultrasound transducer of the ICE console was performed according to an implementation of the method illustrated in FIGS. 8A and 8B.
  • the data for this example was collected with an ICE console with the presence of an electromagnetic tracker.
  • one waveform (125us in length sampled at 200MS/S) was obtained when collecting data in raw/RF mode while the imaging transducer is not receiving imaging energy.
  • Data was saved as a pair of arrays (i.e. sampled waveforms), the first array consisting of an in-band noise-charactehzation waveform (7-13 MHz) sampled after having performed envelope detection, and a corresponding array consisting of an out-of-band noise-charactehzation waveform in the 15- 25MHzband sampled after having performed envelope detection.
  • These in- band and detection band noise-charactehzation arrays are denoted as Ci and, Cn, respectively.
  • the correlated in-band noise and out-of-band noise was then used to perform noise reduction of the in-band imaging waveform during the acquisition of imaging data i.e. while the ultrasound transducer is receiving ultrasound energy.
  • While imaging a pair of arrays were obtained for several periods of acquisition corresponding, each period corresponding to a scan line, ⁇ . For each array pair, one array was obtained in the imaging band (7-13MHZ), and one array was obtained in the noise-detection-band (15-25MHz) . Each array was 125us in length sampled at 200MS/S.
  • the arrays referred to as the in- band imaging and out-of-band imaging arrays, and denoted as Qie and Qne, respectively, were measured during imaging, i.e. while the transducer was receiving imaging energy.
  • Each in-band imaging array was processed for noise reduction using the in-band noise characterization array to reduce noise from the in-band imaging array.
  • the in-band and out-of-band imaging arrays Qie and Qne were segmented into a plurality of time windows.
  • the window size was 800 samples (in this case 4 ⁇ at a rate of 200MS/S), with a 75% overlap between adjacent windows.
  • the imaging arrays, temporally segmented according to the windows, are denoted as Qie.w and Qne.w, where the subscript w is an integer denoting the window number.
  • Temporal alignment between Qne.w and Cn was achieved on a per- window basis.
  • the temporal alignment was achieved by calculating , the cross-correlation between the out-of-band noise characterization array Cn and the out-of-band imaging array Qne.w, and selecting the relative time delay ⁇ corresponding to the maximum cross-correlation. Due to the correlation between the in-band noise and the out-of-band noise, this time delay ⁇ was also applied to align the in-band noise characterization array relative to the windowed in-band imaging array. The aligned in-band noise characterization array is then windowed (denoted by Ci w ).
  • a scaling factor was applied to the in-band noise characterization array Ciw before it was subtracted from the in-band imaging array Qie.w.
  • FIGS. 13A and 13B show the images in obtained in the presence of the noise source, without (A) and with (B) the application of the present noise reduction method. A clear improvement in the signal-to-noise ratio

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Signal Processing (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biophysics (AREA)
  • Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Public Health (AREA)
  • Animal Behavior & Ethology (AREA)
  • Molecular Biology (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Artificial Intelligence (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Psychiatry (AREA)
  • Physiology (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
PCT/CA2018/050212 2017-02-24 2018-02-23 Systems and methods for noise reduction in imaging Ceased WO2018152643A1 (en)

Priority Applications (11)

Application Number Priority Date Filing Date Title
KR1020237002915A KR20230020570A (ko) 2017-02-24 2018-02-23 이미징에서의 노이즈 감소를 위한 시스템 및 방법
KR1020247026556A KR20240125696A (ko) 2017-02-24 2018-02-23 이미징에서의 노이즈 감소를 위한 시스템 및 방법
KR1020237024799A KR102694702B1 (ko) 2017-02-24 2018-02-23 이미징에서의 노이즈 감소를 위한 시스템 및 방법
CN201880024819.8A CN110537201B (zh) 2017-02-24 2018-02-23 用于成像降噪的系统和方法
EP18756578.3A EP3586303B1 (en) 2017-02-24 2018-02-23 Systems and methods for noise reduction in imaging
KR1020197027630A KR20190131031A (ko) 2017-02-24 2018-02-23 이미징에서의 노이즈 감소를 위한 시스템 및 방법
JP2019546178A JP7138113B2 (ja) 2017-02-24 2018-02-23 イメージングにおけるノイズを削減するシステムおよび方法
AU2018225834A AU2018225834B2 (en) 2017-02-24 2018-02-23 Systems and methods for noise reduction in imaging
CA3054383A CA3054383C (en) 2017-02-24 2018-02-23 Systems and methods for noise reduction in imaging
JP2022140465A JP7529732B2 (ja) 2017-02-24 2022-09-05 イメージングにおけるノイズを削減するシステムおよび方法
JP2024119350A JP7842817B2 (ja) 2017-02-24 2024-07-25 イメージングにおけるノイズを削減するシステムおよび方法

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US201762463431P 2017-02-24 2017-02-24
US62/463,431 2017-02-24

Publications (1)

Publication Number Publication Date
WO2018152643A1 true WO2018152643A1 (en) 2018-08-30

Family

ID=63252434

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CA2018/050212 Ceased WO2018152643A1 (en) 2017-02-24 2018-02-23 Systems and methods for noise reduction in imaging

Country Status (8)

Country Link
US (6) US10482582B2 (https=)
EP (1) EP3586303B1 (https=)
JP (3) JP7138113B2 (https=)
KR (4) KR20240125696A (https=)
CN (1) CN110537201B (https=)
AU (1) AU2018225834B2 (https=)
CA (1) CA3054383C (https=)
WO (1) WO2018152643A1 (https=)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111507931A (zh) * 2019-01-14 2020-08-07 阿里巴巴集团控股有限公司 一种数据处理方法和装置
CN113567959A (zh) * 2020-07-27 2021-10-29 北京一径科技有限公司 重频信号的检测方法及装置、处理设备及存储介质
CN115587291A (zh) * 2022-09-26 2023-01-10 华中科技大学 一种基于裂纹超声散射矩阵的去噪表征方法及系统
CN117052380A (zh) * 2023-10-10 2023-11-14 四川宏大安全技术服务有限公司 一种无线压力测量装置及方法
CN117991172A (zh) * 2024-04-03 2024-05-07 山东德源电力科技股份有限公司 一种具有故障识别功能的电压互感器
CN120454643A (zh) * 2025-04-24 2025-08-08 苏州阿珐特机电工程有限公司 一种振荡器相位噪声抑制的方法及装置

Families Citing this family (44)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10482582B2 (en) * 2017-02-24 2019-11-19 Sunnybrook Research Institute Systems and methods for noise reduction in imaging
US10705166B2 (en) * 2017-05-18 2020-07-07 Synaptive Medical (Barbados) Inc. Transmit coil frequency response correction for magnetic resonance imaging
CN107440740B (zh) * 2017-07-21 2021-06-25 无锡海斯凯尔医学技术有限公司 介质粘弹性定量方法及装置
US11255588B2 (en) * 2018-08-03 2022-02-22 Hoshizaki America, Inc. Ultrasonic bin control in an ice machine
EP3633534A1 (en) * 2018-10-04 2020-04-08 Thales Dis France SA A connected device adapted to measure at least a physical quantity
CN109360173B (zh) * 2018-11-23 2021-09-03 珠海医凯电子科技有限公司 基于改进的方差的彩色多普勒血流图像降噪方法
CN109639931B (zh) * 2018-12-25 2022-08-19 努比亚技术有限公司 一种拍照降噪的方法、移动终端以及计算机可读存储介质
JP7336760B2 (ja) * 2019-02-04 2023-09-01 国立大学法人富山大学 超音波断層像生成方法、超音波断層像生成装置、およびプログラム
WO2020219915A1 (en) * 2019-04-24 2020-10-29 University Of Virginia Patent Foundation Denoising magnetic resonance images using unsupervised deep convolutional neural networks
CN110430340B (zh) * 2019-06-28 2020-08-25 北京大学 一种脉冲阵列信号的降噪方法与系统
TWI734252B (zh) * 2019-11-08 2021-07-21 立積電子股份有限公司 雷達及雷達回波訊號的背景成分更新方法
WO2021181242A1 (en) * 2020-03-09 2021-09-16 Marvell Asia Pte Ltd Noise mitigation in an automotive ethernet network
CN112328949B (zh) * 2020-10-26 2024-02-27 中科上声（苏州）电子有限公司 一种汽车发动机主动降噪系统的参考信号生成方法及装置
CN114494028B (zh) * 2020-11-12 2022-12-09 生物岛实验室 粒子束成像降噪方法及装置
CN112698256B (zh) * 2020-12-07 2023-08-29 深圳航天科技创新研究院 降低磁共振成像设备电磁噪声的主动降噪系统的降噪方法
EP4264254A4 (en) 2020-12-16 2024-10-30 Evident Canada, Inc. 3D IMAGE ACCENTUATION FOR DEFECT DETECTION
CN112826494B (zh) * 2020-12-30 2023-05-23 上海联影医疗科技股份有限公司 Mr设备振动和声学噪声消减方法、系统、装置及存储介质
GB202101599D0 (en) * 2021-02-05 2021-03-24 Imperial College Innovations Ltd Processing an arterial Doppler ultrasound waveform
US12178134B2 (en) * 2021-02-10 2024-12-24 Cirrus Logic Inc. Driver circuitry
EP4302258A1 (en) * 2021-04-01 2024-01-10 Huawei Technologies Co., Ltd. Noise reconstruction for image denoising
CN113537147B (zh) * 2021-08-09 2022-04-12 桂林电子科技大学 一种基于光照补偿的夜间车道线检测方法
JP2023041375A (ja) * 2021-09-13 2023-03-24 キヤノン株式会社 情報処理装置、情報処理方法及びプログラム
EP4152651B1 (en) * 2021-09-17 2025-09-03 Nxp B.V. Signal processor for a radio receiver
DE102021210499B3 (de) * 2021-09-21 2023-02-23 Siemens Healthcare Gmbh Verfahren und Vorrichtung zur Unterdrückung von bei Magnetresonanzaufnahmen emittierten elektrischen und/oder magnetischen Feldern
US12169526B2 (en) * 2021-09-24 2024-12-17 Microsoft Technology Licensing, Llc Generating and presenting a text-based graph object
JP7577045B2 (ja) * 2021-10-14 2024-11-01 Hoya株式会社 電子内視鏡システム
JP7690387B2 (ja) * 2021-12-08 2025-06-10 富士フイルム株式会社 超音波診断装置および超音波診断装置の制御方法
EP4460716A1 (en) * 2022-01-08 2024-11-13 Koninklijke Philips N.V. Physiology sensing intraluminal device with index for spectral flow assessment, and associated devices, systems, and methods
EP4209179A1 (en) 2022-01-08 2023-07-12 Koninklijke Philips N.V. Physiology sensing intraluminal device with index for spectral flow assessment and associated devices, systems and methods
WO2023146976A1 (en) * 2022-01-26 2023-08-03 Boston Scientific Scimed, Inc. Reducing catheter rotation motor pwm interference with intravascular ultrasound imaging
US12138126B2 (en) * 2022-05-05 2024-11-12 Fujifilm Sonosite, Inc. Detecting electromagnetic emissions on ultrasound systems
CN114748098B (zh) * 2022-05-13 2025-05-23 青岛海信医疗设备股份有限公司 一种对血流进行滤波的方法、装置及存储介质
JP7775141B2 (ja) * 2022-05-16 2025-11-25 富士フイルム株式会社 超音波時系列データ処理装置及び超音波時系列データ処理プログラム
CN115163052B (zh) * 2022-06-18 2023-07-28 杭州丰禾石油科技有限公司 超声井径的参数测量方法和超声井径随钻测井装置
KR102584250B1 (ko) * 2022-06-22 2023-10-05 주식회사 에어스메디컬 데이터 처리 방법, 컴퓨터 프로그램 및 장치
US12598101B2 (en) * 2022-06-30 2026-04-07 Intel Corporation Crest factor reduction cancellation pulse having a reduced length
JP2024054939A (ja) * 2022-10-06 2024-04-18 コニカミノルタ株式会社 超音波診断装置、超音波画像生成方法及びプログラム
CN116743269B (zh) * 2023-04-10 2025-08-26 北京邮电大学 基于自适应级联经验模态分解的降噪方法、系统及设备
US20250099080A1 (en) * 2023-09-25 2025-03-27 Fujifilm Sonosite, Inc. Suppressing interference artifacts in ultrasound
KR102785610B1 (ko) * 2024-01-19 2025-03-25 (주)윈스테크 초음파 유량 측정에서 수신 신호를 교정하는 장치 및 방법
CN117930370B (zh) * 2024-03-20 2024-06-07 厦门伟卓智慧海洋科技有限公司 海洋磁力场测量及成像方法、装置、设备及介质
CN118010849B (zh) * 2024-04-10 2024-06-14 丰宁满族自治县启源建筑有限公司 一种高速公路桥涵损伤检测方法及系统
KR102774655B1 (ko) * 2024-07-15 2025-03-04 (주) 앤에스알시 Nv 양자센서를 이용한 나노버블 검출장치 및 방법
CN119268998B (zh) * 2024-12-10 2025-03-21 济宁三联钢结构有限公司 一种钢构架的性能测试装置、测试方法及系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070189635A1 (en) * 2006-02-08 2007-08-16 Anja Borsdorf Method for noise reduction in imaging methods
JP2008111062A (ja) 2006-10-31 2008-05-15 Toray Ind Inc 熱可塑性樹脂組成物およびその製造方法
JP2009091756A (ja) 2007-10-05 2009-04-30 Yamaha Livingtec Corp 連動引戸
US20090136104A1 (en) * 2007-11-27 2009-05-28 Hajian Arsen R Noise Reduction Apparatus, Systems, and Methods
US20120197113A1 (en) 2011-01-31 2012-08-02 Sunnybrook Health Sciences Centre Ultasonic probe with ultrasonic transducers addressable on common electrical channel
US20160045184A1 (en) 2013-03-15 2016-02-18 Colibri Technologies Inc. Active localization and visualization of minimally invasive devices using ultrasound

Family Cites Families (56)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS61274495A (ja) * 1985-05-29 1986-12-04 Nippon Hoso Kyokai <Nhk> Ｃｎ比測定方法および装置
JPS62270138A (ja) * 1986-05-20 1987-11-24 富士通株式会社 超音波診断装置
JPH0712358B2 (ja) * 1988-11-16 1995-02-15 株式会社東芝 超音波診断装置
JP2770391B2 (ja) * 1989-03-31 1998-07-02 株式会社日立メディコ 超音波撮像装置
JP2797616B2 (ja) * 1990-03-16 1998-09-17 松下電器産業株式会社 雑音抑圧装置
US5099848A (en) * 1990-11-02 1992-03-31 University Of Rochester Method and apparatus for breast imaging and tumor detection using modal vibration analysis
JPH04303435A (ja) * 1991-03-29 1992-10-27 Shimadzu Corp 超音波診断装置
JPH05161645A (ja) * 1991-12-11 1993-06-29 Fuji Electric Co Ltd 超音波診断装置
US5525906A (en) * 1994-07-11 1996-06-11 General Electric Detection and elimination of wide bandwidth noise in MRI signals
US5935074A (en) * 1997-10-06 1999-08-10 General Electric Company Method and apparatus for automatic tracing of Doppler time-velocity waveform envelope
WO2001024167A1 (en) * 1999-09-30 2001-04-05 Fujitsu Limited Noise suppressor
JP4542258B2 (ja) * 2000-12-21 2010-09-08 アロカ株式会社 超音波診断装置
US6678209B1 (en) 2001-11-21 2004-01-13 Luc Peng Apparatus and method for detecting sonar signals in a noisy environment
JP4219611B2 (ja) * 2002-05-01 2009-02-04 独立行政法人理化学研究所 雑音除去システムおよび雑音除去方法
EP1573993B1 (en) * 2002-10-31 2010-01-13 ZTE Corporation A method and system for broadband predistortion linearizaion
FR2848742B1 (fr) * 2002-12-17 2005-02-04 Thales Sa Procede de correction et boucle de correction d'un signal numerique complexe
US8326621B2 (en) 2003-02-21 2012-12-04 Qnx Software Systems Limited Repetitive transient noise removal
US6814701B1 (en) 2003-04-16 2004-11-09 Aloka Co., Ltd. Method and apparatus for ultrasound diagnostic imaging
JP4575014B2 (ja) * 2004-03-30 2010-11-04 オリンパス株式会社 超音波診断装置
US7513873B2 (en) 2004-10-08 2009-04-07 Supertex, Inc. Low-noise ultrasound method and beamformer system for doppler processing
CN101107608B (zh) 2005-02-16 2011-01-12 埃克森美孚上游研究公司 通过采样其它频率处的噪声估算一个频率处的噪声的方法和装置
CN100525713C (zh) 2005-06-10 2009-08-12 深圳迈瑞生物医疗电子股份有限公司 声谱图包络曲线提取方法
US7946990B2 (en) 2005-09-30 2011-05-24 Siemens Medical Solutions Usa, Inc. Ultrasound color flow imaging at high frame rates
JP4277848B2 (ja) 2005-11-21 2009-06-10 トヨタ自動車株式会社 車載受信装置
CN101138507B (zh) 2006-09-05 2010-05-12 深圳迈瑞生物医疗电子股份有限公司 多普勒血流声音信号的处理方法及装置
JP5162923B2 (ja) * 2007-02-27 2013-03-13 株式会社日立製作所 超音波撮像装置
US8711249B2 (en) * 2007-03-29 2014-04-29 Sony Corporation Method of and apparatus for image denoising
US7532079B2 (en) * 2007-06-18 2009-05-12 Nanoamp Solutions, Inc. (Cayman) Digital tuning of crystal oscillators
JP5459975B2 (ja) * 2008-04-22 2014-04-02 日立アロカメディカル株式会社 超音波診断装置
JP2010240131A (ja) * 2009-04-06 2010-10-28 Toshiba Corp 超音波診断装置及び超音波診断装置制御方法
US8666092B2 (en) 2010-03-30 2014-03-04 Cambridge Silicon Radio Limited Noise estimation
US8970217B1 (en) * 2010-04-14 2015-03-03 Hypres, Inc. System and method for noise reduction in magnetic resonance imaging
US8660226B2 (en) 2011-01-19 2014-02-25 General Electric Company Systems and methods for multichannel noise reduction
JP2013055983A (ja) * 2011-09-07 2013-03-28 Konica Minolta Medical & Graphic Inc 圧電センサー、超音波探触子及び超音波画像診断装置
CN103747736B (zh) * 2011-09-07 2016-02-17 株式会社岛津制作所 图像处理装置以及具备该图像处理装置的放射线摄影装置
WO2013065515A1 (ja) * 2011-11-01 2013-05-10 富士フイルム株式会社 放射線撮影装置、放射線撮影システム及び放射線撮影方法
US9332239B2 (en) * 2012-05-31 2016-05-03 Apple Inc. Systems and methods for RGB image processing
US9105078B2 (en) * 2012-05-31 2015-08-11 Apple Inc. Systems and methods for local tone mapping
JP6049314B2 (ja) * 2012-06-11 2016-12-21 キヤノン株式会社 放射線撮像装置及び画像処理方法
JP6091181B2 (ja) * 2012-11-21 2017-03-08 キヤノン株式会社 放射線画像撮影装置及びその制御方法、並びに、プログラム
BR112015023723A2 (pt) * 2013-03-15 2017-07-18 Colibri Tech Inc métodos para formar uma exibição visual compósita, para calcular mudanças para uma posição e/ou orientação de um dispositivo de formação de imagem e para identificar e realçar uma região em um volume
WO2014144492A1 (en) * 2013-03-15 2014-09-18 Drs Rsta, Inc. Method of shutterless non-uniformity correction for infrared imagers
US9693754B2 (en) * 2013-05-15 2017-07-04 Acist Medical Systems, Inc. Imaging processing systems and methods
US9070405B2 (en) * 2013-09-30 2015-06-30 Seagate Technology Llc Signal noise rejection
CN105594196B (zh) * 2013-10-02 2017-10-31 奥林巴斯株式会社 摄像装置和摄像方法
CN105593698B (zh) * 2013-10-07 2019-05-31 阿西斯特医疗系统有限公司 用于血管内成像的信号处理
US10317502B2 (en) * 2014-03-31 2019-06-11 Koninklijke Philips N.V. Magnetic resonance imaging with RF noise detection coils
CN107408196B (zh) * 2014-12-19 2021-03-19 顶级公司 从图像中提取特征的方法
JP6671109B2 (ja) * 2015-05-15 2020-03-25 キヤノン株式会社 放射線撮影装置、放射線撮影方法、及びプログラム
JP6609119B2 (ja) * 2015-06-03 2019-11-20 キヤノン株式会社 放射線撮影装置、放射線撮影方法、放射線撮影システム、及びプログラム
JP6491545B2 (ja) * 2015-06-12 2019-03-27 キヤノン株式会社 画像処理装置、放射線撮影装置、画像処理方法、プログラム、および記憶媒体
JP6570330B2 (ja) * 2015-06-12 2019-09-04 キヤノン株式会社 画像処理装置、放射線撮影装置、画像処理方法、プログラム、および記憶媒体
JP6538287B2 (ja) * 2016-02-08 2019-07-03 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． 皮膚検出に対する装置、システム及び方法
KR20170107252A (ko) * 2016-03-15 2017-09-25 삼성전자주식회사 이미지 데이터의 노이즈 제거 방법
JP7005206B2 (ja) * 2016-07-26 2022-01-21 キヤノンメディカルシステムズ株式会社 超音波診断装置及び超音波イメージングプログラム
US10482582B2 (en) * 2017-02-24 2019-11-19 Sunnybrook Research Institute Systems and methods for noise reduction in imaging

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070189635A1 (en) * 2006-02-08 2007-08-16 Anja Borsdorf Method for noise reduction in imaging methods
JP2008111062A (ja) 2006-10-31 2008-05-15 Toray Ind Inc 熱可塑性樹脂組成物およびその製造方法
JP2009091756A (ja) 2007-10-05 2009-04-30 Yamaha Livingtec Corp 連動引戸
US20090136104A1 (en) * 2007-11-27 2009-05-28 Hajian Arsen R Noise Reduction Apparatus, Systems, and Methods
US20120197113A1 (en) 2011-01-31 2012-08-02 Sunnybrook Health Sciences Centre Ultasonic probe with ultrasonic transducers addressable on common electrical channel
US20160045184A1 (en) 2013-03-15 2016-02-18 Colibri Technologies Inc. Active localization and visualization of minimally invasive devices using ultrasound

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111507931A (zh) * 2019-01-14 2020-08-07 阿里巴巴集团控股有限公司 一种数据处理方法和装置
CN111507931B (zh) * 2019-01-14 2023-04-18 阿里巴巴集团控股有限公司 一种数据处理方法和装置
CN113567959A (zh) * 2020-07-27 2021-10-29 北京一径科技有限公司 重频信号的检测方法及装置、处理设备及存储介质
CN115587291A (zh) * 2022-09-26 2023-01-10 华中科技大学 一种基于裂纹超声散射矩阵的去噪表征方法及系统
CN115587291B (zh) * 2022-09-26 2023-05-30 华中科技大学 一种基于裂纹超声散射矩阵的去噪表征方法及系统
CN117052380A (zh) * 2023-10-10 2023-11-14 四川宏大安全技术服务有限公司 一种无线压力测量装置及方法
CN117052380B (zh) * 2023-10-10 2024-01-02 四川宏大安全技术服务有限公司 一种无线压力测量装置及方法
CN117991172A (zh) * 2024-04-03 2024-05-07 山东德源电力科技股份有限公司 一种具有故障识别功能的电压互感器
CN117991172B (zh) * 2024-04-03 2024-06-11 山东德源电力科技股份有限公司 一种具有故障识别功能的电压互感器
CN120454643A (zh) * 2025-04-24 2025-08-08 苏州阿珐特机电工程有限公司 一种振荡器相位噪声抑制的方法及装置

Also Published As

Publication number Publication date
EP3586303A4 (en) 2020-12-30
JP2022172285A (ja) 2022-11-15
US10482582B2 (en) 2019-11-19
KR20240125696A (ko) 2024-08-19
US11769230B2 (en) 2023-09-26
JP2024156743A (ja) 2024-11-06
KR20230116072A (ko) 2023-08-03
EP3586303A1 (en) 2020-01-01
CN110537201A (zh) 2019-12-03
US20240404022A1 (en) 2024-12-05
US20230401676A1 (en) 2023-12-14
US11538137B2 (en) 2022-12-27
KR20190131031A (ko) 2019-11-25
KR20230020570A (ko) 2023-02-10
EP3586303B1 (en) 2025-08-20
US20200082509A1 (en) 2020-03-12
US12056855B2 (en) 2024-08-06
CA3054383A1 (en) 2018-08-30
US20180253830A1 (en) 2018-09-06
JP7138113B2 (ja) 2022-09-15
JP7529732B2 (ja) 2024-08-06
AU2018225834A1 (en) 2019-10-03
US10902564B2 (en) 2021-01-26
US20210217141A1 (en) 2021-07-15
AU2018225834B2 (en) 2022-12-01
KR102694702B1 (ko) 2024-08-12
JP2020508168A (ja) 2020-03-19
CA3054383C (en) 2023-07-11
JP7842817B2 (ja) 2026-04-08
CN110537201B (zh) 2024-04-09
US20230117435A1 (en) 2023-04-20

Similar Documents

Publication Publication Date Title
US12056855B2 (en) Systems and methods for noise reduction in imaging
US8597196B2 (en) Cardiac signal processing apparatus and cardiac signal processing method
JP6055862B2 (ja) コレジスタ・イメージングのための血管内超音波システム
EP2640274B1 (en) Interference reduction and signal to noise ratio improvement for ultrasound cardiac ablation monitoring
EP3542346B1 (en) Adaptive ringdown subtraction for coronary and peripheral intravascular ultrasound (ivus)
US12558056B2 (en) Ultrasound based blood flow velocity measurements
RU2588285C2 (ru) Уменьшение помех и улучшение отношения сигнал-шум для ультразвукового мониторинга кардиальной абляции
JP2004337461A (ja) 超音波診断装置
Thakur et al. Comparison of ECG Signal Denoising Algorithms in Fir and Wavelet Domains

Legal Events

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

Ref document number: 18756578

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 3054383

Country of ref document: CA

ENP Entry into the national phase

Ref document number: 2019546178

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

ENP Entry into the national phase

Ref document number: 20197027630

Country of ref document: KR

Kind code of ref document: A

ENP Entry into the national phase

Ref document number: 2018756578

Country of ref document: EP

Effective date: 20190924

ENP Entry into the national phase

Ref document number: 2018225834

Country of ref document: AU

Date of ref document: 20180223

Kind code of ref document: A

WWG Wipo information: grant in national office

Ref document number: 2018756578

Country of ref document: EP