Classifications

• • 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/539—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
• • 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
  • G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
  • G01S15/02—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
• • 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/52023—Details of receivers
  • G01S7/52036—Details of receivers using analysis of echo signal for target characterisation
  • G01S7/52042—Details of receivers using analysis of echo signal for target characterisation determining elastic properties of the propagation medium or of the reflective target
• • 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/52079—Constructional features
• • 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/523—Details of pulse systems

Definitions

• the present invention relates to a device for characterization of any of a ductile membrane, a surface, or sub-surface material properties.
• the invention relates to the non-contact characterization of a ductile membrane, surface or sub-surface regions by stimulating the surface or membrane using a low frequency excitation such as an audible excitation and measuring the displacement behavior with a comparatively higher frequency excitation such as ultrasound.
• the measurement of a ductile surface or membrane is traditionally done by external application of a force to develop an understanding of the reaction to the surface or membrane to the application of a force.
• an example characterization may be drawn to a spring constant over a short range of motion, or in the case of a displacement in one direction followed by a non-linear displacement in the opposite direction, for a lossy system, a hysteresis effect may also be observed.
• the membrane provides an interface and challenge surface for a surface stimulation, and the characteristics of the fluid or gel adjacent to the membrane are measured using the membrane, which otherwise does not change the characteristics being measured, as the membrane is closely coupled to the surface and underlying fluid such that displacements to the underlying fluid may be made through the surface or membrane without changing the characteristics to be measured.
• Such characterizations are widely performed using mechanical (stirring paddles) or other means in the food industry to measure the stiffness or other mechanical properties of various liquid or semi solid food items, such as those with curing or culture times. For example, in the prior art of food science, it is desired to know when yogurt has reached a particular culture level by its mechanical properties, whereas the bacterial growth interval may be determined from prior batches and a fixed interval used.
• the ripeness of a fruit or vegetable may be determined by its surface softness, elasticity, or other mechanical property, which may be performed by the application of a small force sufficient to produce measurable deflection or which causes surface or shear waves across the surface of the fruit or vegetable and measuring the minute surface or sub-surface deflections in response to that surface force.
• the eye globe may be measured by using a contact or optical method, to determine the deflection upon application of an air puff.
• the optical reflectivity characteristics of the eye may be measured as a deflection versus applied force.
• the optical interface of the eye and the various underlying structures may have resolution of motion limitations which are the consequence of the optical
• the measurement system characterizing an optically transparent structure, where the intended target (such as the cornea) and surrounding structures (such as the ocular lens and vitreous fluid) have similar optical properties and are therefore difficult to resolve from each other. It is desired to use remote sensing and remote interrogation, such as through air, to measure the elasticity of the surface of an eye or other organ.
• Ultrasound measurements are typically performed using a liquid medium as a coupling fluid, such that the ultrasound energy propagates from the transducer through the fluid as traveling acoustic waves, creating reflected ultrasound energy arising from target interfaces which have differences in acoustic index of refraction according to the transmission impedance mismatch between the coupling fluid and the target.
• the use of ultrasound has not been possible for non-contact displacement measurements resolving distances on the order of millimeter (mm) and sub-mm distance ranges because of the traditional requirement for coupling fluids between the transducer and target.
• a method for the non-contact surface characterization of a surface through remote sensing of displacement in combination with the application of a surface stimulus is desired. It is further desired to use air as the medium for stimulation and measurement, without direct contact to the surface or membrane by either the stimulus or measurement device.
• a first object of the invention is an apparatus and method for the generation of a pulsatile mechanical stimulation for application to a surface or membrane, the mechanical stimulation combined with an ultrasonic measurement of the response of a surface or membrane to the stimulation by examination of a reflected ultrasound signal, and the characterization of the reflected response to determine mechanical properties of the surface or membrane.
• a second object of the invention is the measurement of surface, sub-surface, or membrane ductility and elasticity by examination of an ultrasound signal which has reflected from the surface, sub-surface, or membrane.
• a third object of the invention is an apparatus and method for the application of a force to a food item, the force generating a surface or membrane deflection with a temporal response, including a shear wave or surface wave, the surface or membrane deflection characterized by a phase or amplitude change in an ultrasound signal which has reflected from the surface or membrane.
• a fourth object of the invention is the characterization of a substantially homogeneous fluid for elasticity or viscosity by application of a non-contact displacement force and characterization of the response to the non-contact displacement force by measuring the phase and amplitude of reflected ultrasound from transmitted ultrasound applied to the surface of the fluid.
• a fifth object of the invention is an apparatus and method for the measurement of an elasticity of an animal or human organ, whereby an air puff is applied to the organ, such as a human eye to perform tonometry, and the response of the organ is measured, where the air puff response of the organ during the response time of the relaxation from air puff results in a characteristic phase change of the reflected ultrasound energy compared to the phase of the applied ultrasound energy over a continuous wave (CW) ultrasound interrogation or a series of pulsed ultrasound interrogations.
• CW continuous wave
• the present invention may utilize an ultrasound transducer, such as a capacitive micro- machined ultrasound transducer (cMUT), which provides a diaphragm structure with a thin diaphragm which provides a good impedance match to gaseous media for propagation of ultrasonic energy through air from the displacement of the diaphragm.
• cMUT capacitive micro- machined ultrasound transducer
• the cMUT is suitable for coupling through air, and is suitable for duplication in a planar 2D hex packing array, such as on a substrate, for generation of a planar wave of desired spatial characteristic and lateral extent.
• a material which is optionally on the far surface of a membrane can be characterized by the application of a pulsatile challenge displacement force such as by an air puff, where an ultrasound transducer measures the dynamic displacement characteristics of the membrane or surface.
• the membrane transparently provides a physical barrier to the material to be characterized, and so does not significantly change the properties of the material as seen by the ultrasound transducer.
• the ultrasound transducer sends and receives ultrasound energy through a medium such as air to the surface or membrane to be characterized, the ultrasound transducer operative during an interval when a comparatively low frequency excitation source is coupled to the surface or membrane.
• the low frequency excitation source generates a subtle movement of the surface or membrane during an interval coincident with an ultrasound transmitter delivering acoustic wave ultrasound energy to the surface or membrane either in CW form or in pulsed form.
• a receiver for ultrasound reflected from the surface or membrane measures displacement of the surface or membrane as a phase change in the received signal when compared to the transmit frequency, thereby indicating a temporal displacement of the surface or membrane.
• An analysis of the temporal displacement of the surface or membrane, as measured by the phase shifts of the reflected ultrasound in response to the pneumatic excitation coupled to the surface or membrane, in combination with comparison to the temporal displacement or from templates or metrics associated with the delay in and amplitude of response between the excitation stimulus to and ultrasound response from the surface or membrane, is used to determine the mechanical characteristics of the surface or membrane.
• the mechanical characteristics measured may include ductility, elasticity, or hardness. In this manner, a non-contact measurement of the mechanical properties the fluid below the surface or membrane may be determined.
• a non-contact characterization system may comprise an excitation generator generating a non- contact displacement force directed to a remote object to be characterized.
• a non-contact characterization system may comprise an ultrasound transducer directing ultrasonic energy to a remote object; a transmit waveform generator coupled to the ultrasound transducer.
• a non- contact characterization system may comprise a receive signal processor coupled to the ultrasound transducer, and receiving ultrasound energy reflected from a remote object, the receive signal processor comparing the received ultrasound signal amplitude or phase to a transmit signal formed by the transmit waveform generator, the receive signal processor thereby forming an estimate of remote object displacement in response to the non-contact displacement force.
• the excitation generator may form a puff of air directed to a remote object.
• the puff of air may contain a non-oxidizing gas.
• the non-oxidizing gas may be at least one of: nitrogen, argon, C02, helium.
• the ultrasound transducer may be a capacitive micromachined ultrasonic transducer (cMUT).
• the transmit signal generator may be either a continuous wave (CW) or pulsed signal generator operative at a fundamental or harmonic frequency of the ultrasound transducer transmit frequency.
• the receive signal processor may perform baseband
• each phase estimate determined from a phase displacement of the received signal compared to the transmitted signal.
• the series of phase estimates may form an elasticity metric.
• Another aspect of the present disclosure provides a method for non-contact
• a non-contact force event may be provided.
• the excitation generator may provide a non-contact force to a surface or a volume of material adjacent to the surface to be characterized.
• a series of displacement measurements may be formed for each applied excitation generator non-contact force event.
• Ultrasound energy may be directed from a transducer to a surface or adjacent volume of material to be characterized.
• the transducer may receive reflected ultrasound energy from the surface or volume to be characterized.
• the phase of the transmitted ultrasound energy may be compared to the received ultrasound energy to form a displacement estimate.
• a plurality of said displacement measurements may be used to estimate an elasticity or a viscosity.
• the excitation generator non-contact force may be a puff of air released from a nozzle directed at a surface on or of a material to be characterized.
• the puff of air may be an inert gas including at least one of: Nitrogen, Argon, C02 or Helium.
• the transducer may be at least one of a cMUT or piezoelectric transducer.
• the directing ultrasound energy may be a continuous wave (CW) ultrasound energy on a transmit transducer.
• the comparing the phase of the transmitted ultrasound energy to the received ultrasound energy may form a phase displacement which may be converted to a non-contact displacement of a material to be measured.
• the estimate of elasticity or viscosity may be performed by comparing a series of measured phase displacements to a template series of measurements.
• the estimate of elasticity or viscosity may be performed by examination of a rise time or a fall time of a displacement compared to the rise time or fall time of said excitation source.
• the excitation may be a step or impulse pressure.
• the estimate of elasticity or viscosity may be performed by examination of a ringdown characteristic.
• the ringdown characteristic may be at least one of an exponential decay time or a ring cycle interval or frequency.
• the excitation may propagate to the surface or the volume to be characterized as a surface wave, a shear wave, or a combination of surface wave and shear wave.
• the surface or volume to be characterized may include a membrane over a fluid.
• the surface or volume to be characterized may include a fluid having an elasticity or viscosity.
• the surface or volume to be characterized may include an animal or human.
• the non-contact material characterization system may comprise an excitation generator.
• the excitation generator may be configured to generate a non-contact displacement force directed to an elastic surface or to a volume of material adjacent to the elastic surface.
• the non-contact material characterization system may comprise an ultrasound transducer.
• the ultrasound transducer may be configured to direct a transmit signal to the surface or to the volume.
• the ultrasound transducer may be configured to receive a reflected Doppler ultrasound signal from the surface or the volume.
• characterization system may comprise a processor coupled to the ultrasound transducer.
• the processor may be configured to determine a displacement of the surface or the volume in response to the non-contact displacement force and the reflected Doppler ultrasound signal.
• the processor may be configured to determine the displacement by comparing at least one of an amplitude or a phase of the reflected Doppler ultrasound signal to at least one of an amplitude or a phase of the transmit signal.
• the non-contact displacement force generated by the excitation generator comprises a puff of air.
• the excitation generator may be configured to generate the non-contact displacement force with a frequency of at least 20Hz.
• the non-contact displacement force generated by the excitation generator comprises an impulse excitation.
• the ultrasound transducer may comprise a capacitive micromachined ultrasonic transducer (cMUT) or a piezoelectric transducer.
• the non-contact material characterization system may comprise a waveform generator coupled to the ultrasound transducer.
• the waveform generator may comprise either a continuous wave (CW) or a pulsed signal generator.
• the waveform generator may be operative at a fundamental or a harmonic of a frequency of the transmit signal.
• the processor may be configured to determine the displacement by performing a series of comparisons of the phase of the reflected Doppler signal to the phase of the transmit signal.
• the processor is further configured to perform baseband demodulation of the reflected Doppler signal.
• the processor may be configured to determine at least one of an elasticity or a viscosity of the elastic surface or volume based on the series of phase estimates.
• the processor may be configured to determine at least one of an elasticity or a viscosity of the elastic surface or volume based on the determined displacement.
• a non-contact displacement force may be provided to the elastic surface or to a volume of material adjacent to the elastic surface (step (a)).
• Ultrasound energy may be transmitted from a transducer to the surface or a volume adjacent to the elastic surface (step (b)).
• Reflected Doppler ultrasound signal may be received at the transducer (step(c)).
• At least one of an amplitude or a phase of the transmitted ultrasound energy may be compared to at least one of an amplitude or a phase of the reflected Doppler ultrasound signal to form a displacement measurement (step (d)).
• An elasticity or a viscosity may be calculated from the displacement measurement in response to the non-contact displacement force.
• Providing the non-contact force to the elastic surface or to the volume of material adjacent to the elastic surface may comprise directing a puff of air to the elastic surface or to the volume of material adjacent to the elastic surface.
• the non-contact displacement force may comprise an impulse excitation.
• the transducer may comprise a capacitive micromachined ultrasonic transducer (cMUT) or a piezoelectric transducer. Transmitting the ultrasound energy may comprise transmitting continuous wave (CW) ultrasound energy.
• CW continuous wave
• Providing the non-contact force may comprise providing a step or impulse pressure.
• Calculating the elasticity or the viscosity may comprise analyzing a ringdown characteristic of the reflected Doppler ultrasound signal.
• the ringdown characteristic may comprise at least one of an exponential decay time, a ring cycle interval, or a ring cycle frequency.
• the non-contact displacement force may propagate to the surface or the volume as at least one of a surface wave or a shear wave.
• the surface or the volume may comprise a membrane over a fluid.
• the elastic surface may comprise tissue of an animal or a human.
• the elastic surface may comprise a solid or semi-solid food item.
• the non-contact displacement force may be provided at a frequency of at least 20Hz.
• Figure l is a block diagram of a signal processor system for estimating the
• Figure 2 is a block diagram as in figure 1 where the signal processor operates directly on received ultrasound echoes.
• Figure 3 shows waveforms for the system of figure 1.
• Figure 4A shows a plot for a sinusoidal excitation applied to deformable surface or membrane with a reflected response signal having a phase delay and amplitude level.
• Figure 4B shows a plot for a step excitation applied to a deformable surface or membrane, and a response to the step excitation having a phase delay and amplitude level.
• Figure 4C-1 shows a plot of a sinusoidal surface or membrane displacement generating more than +/-180 0 of phase shift.
• Figure 4C-2 shows the acquired data with phase wrapped from the large phase shifts of figure 4C-1.
• Figure 4C-3 shows a plot of an unwrapped phase estimate from figure 4C-2.
• Figure 5 shows a CW signal processor for continuous interrogation of a membrane or surface or membrane in response to an excitation generator.
• Figure 6 shows the waveforms for the CW system of figure 5.
• Figure 7A is a plot of a sinusoidal excitation source and associated membrane or surface or membrane displacement response.
• Figure 7B is a plot of a step excitation source and associated membrane or surface or membrane displacement response.
• Figure 1 shows a signal processor for an example embodiment of a non-contact, gas coupled surface, sub-surface, or membrane characterization system.
• a surface 150 may be a membrane or a compressible surface with elastic characteristics.
• the surface 150 is interrogated by an ultrasound beam 128 from an ultrasound transducer 160 which is optionally mounted adjacent to a source of applied non-contact force, such as a gaseous puff such as an air puff, which is conveyed through hose 122 to nozzle 124, or other means for providing a non-contact pressure stimulus to a membrane or surface such that the response of the membrane or surface to the non-contact force can be observed.
• a source of applied non-contact force such as a gaseous puff such as an air puff
• the nozzle 124 may be coupled to an interior volume which provides for coupling of dynamic pressures from surface excitation generator 120 through hose 122 to the surface 150 where the air pressures result in displacement of the surface 150.
• This displacement can also be the result of a pressure differential, such that one side of a membrane or surface is momentarily greater or lower than the other side, and the surface or membrane response characterized after the pressure differential is in steady state.
• the excitation generator 120 may generate static or dynamic pressure variations which are coupled to the surface 150.
• the excitation generator 120 may produce any suitable static or dynamic pressure modulation for displacement of the surface, including a sub-audio frequency below 20Hz, an audio frequency from 20Hz to 20KHz, or a super-audio frequency above 20KHz.
• the nature of the pressure excitation generated by the excitation generator may be an impulsive step or delta (impulse) generation, a sinusoidal pressure excitation, a square wave excitation, or any combination of these, and the excitation may be a gated burst or continuous.
• the pressure excitation may be provided with or without a static positive or negative pressure bias.
• Pressure nozzle 124 also has an adjacent ultrasound transducer 160 with electrical leads 162 and 164 coupled to transmit receive switch 118. Ultrasound transducer 160 generates ultrasound beam 128 which is directed to a central region of the surface 150.
• a controller 148 generates a variety of control signals which are distributed through the signal processor 100.
• a system reference clock 110 may be derived from a temporally stable clock source, and the reference clock 110 may also be used for demodulation of the received signal.
• System reference clock 110 is coupled to a transmit waveform generator 112 which generates a pulse train at or near the center frequency of transducer 160, transmit transducer interface 114 performs voltage level shifting and any required amplification before coupling to the transmit/receive switch 118, which couples the waveforms from transmit interface 114 to the ultrasonic transducer 160 via leads 162 and 164.
• the ultrasound transducer 160 generates and directs the ultrasonic energy in beam 128 to the surface 150.
• Reflected energy from the surface 150 is coupled from the transducer 160 back through leads 162 and 164 to the transmit/receive switch 118, where it is directed to the receive preamplifier 116, which amplifies the signal level and isolates the signal using frequency filtering to remove out-of-band frequency components, and optionally provides automatic gain control through a gain control input from controller 148.
• the output of the receive preamplifier 116 is applied to quadrature mixers 140 and 142, where a quadrature clock from clock generator 110 at the ultrasound transmitting frequency generates a quadrature output comprising an I (in- phase) baseband channel and Q (quadrature, or 90 degrees separated) baseband channel, which are coupled to identical low pass filters 136 and 138, each of which has a respective analog to digital converter 132 and 134, the output of which is stored in data buffers 144, one data buffer for each I and Q channel. It is also possible to provide multi-gate samples over time, so that if the target moves out of the range of a particular gate, the system can adaptively use samples associated with the subsequent depth.
• each transmit pulse When the received signal is mixed with the reference clock in this manner, each transmit pulse generates a single phase value per range gate sample at a particular depth, and over a series of transmit events this sequence of phase differences is used by the phase and amplitude analyzer 146 to estimate the temporal displacement of surface 150.
• the transmit clock coupled to the transducer during the transmit interval is derived from system clock 110, which is substantially at the center frequency of the transducer.
• the phase and amplitude analyzer 146 examines primarily the phase of the returned signal
• the system clock at the transmit rate, supplies synchronous clocking to quadrature mixers 140 and 142 during the receive interval to compare the receive signal phase to the system clock (at the original transmit frequency) to generate a phase difference between the transmitted pulse and the reflected pulse.
• This phase value may be compared over one or more cycles of the receive signal to establish an average phase value for that particular receive interval, and then each phase value from each receive interval is assembled to provide a continuous series of discrete estimates of surface 150 displacement, based on the wavelength of the acoustic wave and the phase value measured.
• the phase and/or amplitude analyzer 146 may operate on the amplitude of the received signal, which may be analyzed to provide information about the quality of the phase estimate made from the data (such as from signal to noise metrics), or the amplitude of the signal may be analyzed to provide a metric such as db/MHz-cm falloff, or the amplitude profile may provide an elasticity metric which indicates the elasticity or other surface characteristic based on the reflection of ultrasound energy from the surface, in response to the momentary pressure perturbation provided by the excitation generator.
• the elasticity metric is any phase or amplitude derived metric from the data presented to the amplitude and phase analyzer 146 which provides a measurement of surface elasticity, where the mobility is preferentially associated with the elasticity of the surface 150.
• Controller 148 which generates the excitation 122 also reads the output of phase and amplitude analyzer 146 over the duration of excitation generator 120 activity, and optionally the amplitude of the reflected signal, to derive a temporal response of the surface 150 to the pneumatic excitation provided through coupling tip 124.
• the pneumatic excitation may be any sub-audio, audio, or super-audio frequency or pulse as previously described.
• Figure 2 shows an alternate embodiment of the signal processor of figure 1, where the signal processor is performing direct sampling of the RF signal from the transducer, rather than using quadrature mixing to baseband of the RF signal.
• System clock 210 generates the transmit clock, which is coupled to transmit waveform generator 112.
• the operation of transmit waveform generator 112, transmit transducer interface 114, transmit receive switch 118, receive preamplifier 116, surface excitation source 120 and transducer 160 are as previously described for figure 1.
• the receive preamplifier 116 may be gain controllable, as before, with the gain determined by controller 248 to place the RF signal in optimum A/D converter 232 range.
• the output of the receive preamplifier 116 is directed to a band pass filter 236 for reduction of the noise bandwidth applied to the ADC 232, which samples at the Nyquist rate of at least 2X faster than the applied signal.
• the Nyquist sampling rate is at least 3MHz plus the skirt falloff associated with the bandwidth of the transducer 160, known in the art of signal sampling as the Nyquist sampling criteria.
• the single channel output of the ADC 232 is applied to a data buffer 244, and a signal analyzer 246 examines phase shifts in the buffered signal to determine phase changes of the RF signal to discern movement of the surface.
• the sequence of phase measurements used to form the phase measurement may be a series of measurements which are inverse-time weighted to increase the effect of recently acquired measurements, or they may be uniformly weighted over a window of phase samples.
• the use of a weighting coefficients applied to the stream of measurements over a window may provide favorable noise rejection characteristics, and weighting may be chosen to favor signals in the excitation source bandwidth to filter and reduce the effect of noise which is outside the excitation source bandwidth.
• Figure 3 shows example operation of the ultrasound processor of figure 1.
• transmit/receive events provide an estimate of the surface position as a series of phase values during a series of repeated interrogation intervals 340, each of which provides a single phase value.
• System clock waveform 302 operates continuously, and is furnished by system clock generator 110 of figure 1.
• the duration of the event interval 340 is determined by the time- of-flight from the transducer 160 to the surface 130 and back to the transducer 160 of figure 1.
• the propagation velocity of ultrasound in air is 330m/s (.33mm/ps). Accordingly, for a l.5MHz transducer, the resultant wavelength of this traveling wave in air is 0.22mm.
• the total time of flight for an ultrasound signal lOmm each direction is then 60ps, so duration 340 may be no less than 60ps in this case.
• This time of flight interval for a transmit pulse to return as a receive signal after reflection is shown as interval 343 in figure 3.
• the time of flight provides an upper limit to the pulse repetition frequency (PRF) corresponding to the sum of the transmit interval and receive interval.
• PRF pulse repetition frequency
• the transducer with a l.5MHz center frequency will have a 220u wavelength traveling in air. A displacement of the surface will result in a shortened path from the transducer to the surface, and the reflected signal from the surface back to the transducer will return with a phase shift.
• phase and amplitude analyzer observing a phase offset of 180 degrees between transmit clock and received signal compared to a datum phase offset will correspond to a 55 pm displacement of the surface.
• a transmit interval 342 for the transmission of a longer pulse train provides improved signal to noise ratio of the receive signal phase and also extends the return time of flight by the duration 342 of the transmit pulse stream, at the expense of decreased axial resolution, which may be desirable for the case of a discrete moving target such as the surface.
• transmit interval 342 is 6.6ps
• the maximum interval 340 is 66.6ps, which implies a pulse repetition frequency (PRF) of l5KHz or less.
• PRF pulse repetition frequency
• the first cycle of transmit energy reaches the surface at the same time the transducer is finishing sending the last of the transmit burst.
• the PRF will be much lower to account for the required attenuation of multi-path reflection energy which will mix with the surface reflections.
• a CW system In a CW system, separate transmit and receive transducers are used and multipath considerations can be ignored when the surface to be measured has signal energy which dominates the Doppler signal energy which is returned. It may be preferable for the system to operate in CW mode in some circumstances, and in pulsed mode in others, depending on the nature of the reflected signal energy. For pulsed mode, several parametric changes are available to improve SNR, including the use of shorter or longer bursts of transmit RF energy (307 of figure 3) which match the spatial extent of the diaphragm or fluid which is responding coherently to the excitation, such as by moving in a common direction, as a stronger reflected signal with better SNR is provided when the sample volume which is returning acoustic energy is in motion in the same direction.
• a short transmit burst 307 with spatial extent which matches the regions of diaphragm or fluid which are moving in the same direction, and to successively range gate or sample each region of the reflected signal separately to distinguish movement or displacement of one region from movement or displacement of another region, thereby treating each reflection region separately.
• This may also provide advantage where a membrane or surface moves into range of one sampled region and out of range of another sampled region (such as the spatial extent represented by waveforms 311/313 of figure 3). It is alternatively possible to provide many cycles of transmit energy to improve the phase accuracy of each measurement, particularly where a clear surface reflection boundary is present and most of the signal energy is reflected from the surface.
• the combined transmit interval and receive interval which determine the PRF may be in the repetition period range of 50ps to lms or more. As multi-path reflections may occur, it may be preferable to reduce the maximum PRF to reduce the effect of ultrasonic reflections from transmit events earlier than the current interval 340, for example.
• the maximum PRF For an ultrasound propagation velocity of 33mm/ps, when the separation from the transducer to surface is l5mm, the round trip ultrasound path requires ⁇ 90ps, and if the separation distance from transducer to surface is 20mm, the round trip path requires ⁇ l20ps.
• Transducer waveform 306 shows the transmit waveform 307 which includes bias and amplitude corrections during the transmit interval 342, and a reduced amplitude receive signal 309 from the surface.
• the received signal 309 also includes the effects of surface displacement in the form of a phase change from the system clock, which must be subtracted from any static phase value which may be present.
• Each 66ps cycle provides a single phase estimate value, which may be considered in polar coordinates using the I and Q outputs. This may be done using a range gate select a time of flight interval corresponding to the region containing a reflection from the surface to obtain each sample indicating the instantaneous phase of the surface for a particular sample from a transmit event.
• Each of the acquired values within an RX interval 344 is averaged or temporally filtered over the temporal region corresponding to the surface reflected response to reach an average phase estimate shown as 311 and 313, respectively, for I and Q waveforms 308 and 310.
• phase estimates are saved, each of these estimates spanning an extent of the Rx interval 344 and which extent corresponds to a reflection from a particular depth.
• samples of IQ are concatenated to construct a time series describing surface or sub-surface motion, since phase change over time is attributed to change in distance from the transducer.
• a succession of these sampled values are collected and compared against a surface excitation waveform which is used to form a characterization of the surface for a particular excitation waveform.
• Figure 4A shows an example sinusoidal excitation applied to a surface, such as a sinusoidal waveform 321 applied using a voice coil diaphragm displacing a volume sufficient to modulate a localized region of the membrane or surface pressure by lOOdaPa (dekapascals) p-p.
• Sub-sonic frequencies ⁇ 20Hz
• audio frequencies (20Hz-20kHz)
• super-audio frequencies >20kHz
• the sinusoidal pressure excitation 321 results in a modulation of the surface displacement, which is shown as phase plot 332, as the modulation in surface position corresponds to a change in the phase of the return signal.
• Each discrete circle of waveform 332 represents a sample point such as a polar conversion of average values for I 311 and Q 313.
• a series of audio and sub-audio tones with various cycle times 322 and measuring the phase response 330 as shown in plot 332 it is possible to estimate the characteristic such as viscosity or elasticity of the fluid behind the surface.
• an exemplar elasticity metric measurement associated with the changed density or viscosity of the fluid could be an associated change in surface or membrane response time.
• a frequency domain response of the surface may be made using a series of excitations 321 and measuring a series of surface responses 332.
• FIG. 4C-1 shows a received signal 430 with displacement-associated phase excursions which exceed l/2 (180°). Because phase excursions greater than 180° wrap to -180°, the series of samples of figure 4C-2 wrap and produce the series of samples shown, with samples of individual segments 432, 434, 436, 438, and 440. If a sufficiently high sample rate is used which exceeds the Nyquist sampling rate, it is possible to "unwrap" the samples as shown in figure 4C-3, to provide the original phase information. These techniques are well known in the prior art of Doppler signal reconstruction.
• figure 4A shows a sinusoidal excitation which may be provided in a series of such excitations to generate a phase vs. frequency response plot of the surface displacement from the series of measurements
• Figure 4B shows a time domain step response equivalent of figure 4A, where a surface step pressure excitation 362 of 50 daPa peak is applied to the surface or membrane, which generates the phase response 372 of the return signal from the surface. It is similarly possible to characterize the surface response based on a time delay 374 and amplitude response (shown as 180 degrees) for phase response plot 372, corresponding to .11/2 mm displacement.
• the phase unwrap techniques described in the series of figure 4C-1, 4C-2, 4C-3 may similarly be applied to the samples of waveform 372 of figure 4B to reconstruct phase shifts in excess of +/-180°.
• the signal processing of figure 2 operates in a similar manner as was described for figure 3, however the reflection 306 observed by the transducer is directly sampled and compared with a reference clock to determine the phase changes associated with the surface movement, for example by multiplying the reference clock with the received signal over a receive signal averaging time, and integrating this value over the duration of the receive signal to estimate a phase value for one receive interval.
• this will result in the generation of response waveform 332 from excitation source 321 interacting with the surface, as described for figure 4 A, or response waveform 372 from excitation source 362 interacting with the surface.
• the phase analyzer 146 or signal analyzer 246 is operative to detect the largest amplitude of the I or Q channel, or largest magnitude of vector generated by IQ, and associates the region of that large signal response as the surface to be characterized, subsequently sampling that region to form an estimate of phase offset from which displacement may be estimated.
• Figure 5 shows another embodiment of the invention for CW operation.
• the signal processor of figure 5 operates as in figure 1, and with the same block descriptions operative as was present in figure 1, however the transmit interface 114 is directly coupled via leads 502/504 to a transmit transducer 524 and generating transmit beam 526, which is coincident on the surface with the receive beam profile 528 of receive transducer 530, which conveys the receive signal using leads 506/508 to receive amplifier 116, where the signal processing occurs as described previously for figure 1, however, the system of figure 5 operates continuously, with the transmitter continuously transmitting, and the receiver baseband signal being continuously received. This operation is advantageous for detection of signal bandwidth which exceeds the pulsed transmit configuration described in figure 3.
• the low pass filters 570 and 572 may set to the Nyquist rate of half of the transmit clock rate (Ftx/2), or they may alternatively be set to the highest response frequency expected during measurement, not to exceed Ftx/2.
• the advantage of matching the low pass filter 570 and 572 bandwidth to the signal bandwidth is the reduction in excess noise from out-of-band noise frequency
• the systems of figures 1, 2, and 5 may require the ultrasound transceiver radiation patterns (128 of figures 1 and 2, 526 of figure 5) be visualized at the insonified surface of the membrane, so that the region of insonification which is being characterized may be easily identified visually.
• An optical illumination system (such as a collimated LED, laser, or other beam source) may be added to indicate the region being examined and characteristics by the system, such as by adding an optical source which converges or is co-linear with the radiation pattern 128 of figure 1 and 526 of figure 5.
• Figure 6 shows waveform plots for the baseband CW system of figure 5.
• the system clock 110, transmit waveform generator 112, and transmit transducer interface 114 generate a DC-biased transducer CW signal waveform 602 of figure 6, which is applied to the transmit transducer 524 of figure 5, and the receive transducer 530 of figure 5 generates receive signal 608 of figure 6.
• a DC bias may be required for a cMUT transducer, whereas a piezoelectric transducer does not require a DC bias.
• the outputs of the I and Q channel low pass filters 136 and 138, respectively, are shown as waveforms 614 and 616.
• the phase unwrapping techniques described previously may be applied to these waveforms as well, where the detected phase crosses the +/-180° boundary and wraps to the opposite boundary.
• Figures 7 A and 7B show CW output 714 for an excitation 702, and the sample points of 332 and 372 of figures 4A and 4B are not present, as the CW system of figure 5 is not subject to the baseband Nyquist sampling limitations of the pulsed dopper system of figures 2 and 3, provided that the mixer output is sampled at a sufficiently high rate to satisfy the Nyquist criteria for phase changes at the mixer output.
• the transducer types for 130 of figures 1 & 2, and 524 and 530 of figure 5A may be any of capacitive micromachined ultrasonic transducer (cMUT), or piezoelectric transducers, for example, formed with the piezoelectric material PZT.
• cMUT capacitive micromachined ultrasonic transducer
• PZT piezoelectric material
• K is the Boltzmann constant 1.38 * 10 23 ;
• T is the temperature of the system, assumed to be 300°K;
• B is the bandwidth of the sampled signal (either the bandwidth of the bandpass filter 236 of figure 2, or bandwidth of the low pass filter 136/138 of figures 1 and 5;
• R is the resistance generating the Johnson noise, typically 50 ohms.
• the noise floor for a 50 ohm system is 0.9nV /L/ ⁇ TZ. It is typically easier to perform narrowband filtering on a baseband signal such as the low pass filters 136 and 138 of figure 1 than the bandpass filter 236 of figure 2.
• a first order band pass filter 236 for a l.5MHz system might have a 3db bandwidth of 1 MHz, whereas the desired signal content is below lKHz, which is difficult to incorporate into bandpass filter 236, but simple to incorporate into low pass filter 136.
• the sample noise floor for lKHz baseband system would 28nV rms whereas the lMHz bandwidth direct sampling system would be 30x higher, or 900nV rms with the same signal energy.
• the noise factor of the system (typically governed by the first few elements in the receive chain) is managed separately, as it would scale the noise floor by the noise factor, so a 6dB noise factor would approximately double both of the above rms noise floor values.
• the invention may be practiced many different ways.
• the phase and amplitude analyzer produces an elasticity metric which is a characterization of the sequence of phase measurements from the ultrasound reflection from the surface in combination with the displacement of the surface from the surface excitation source.
• the elasticity metric which is derived from the response of the surface may provide an indication of several different phenomenon.
• the elasticity metric may indicate whether a surface adjacent to a membrane has a gaseous boundary (in which case the reflection is from the membrane itself) or fluid boundary (in which case the reflection is from both the membrane and fluid adjacent to the membrane).
• the elasticity metric may indicate, for the case of characterizing a fluid behind the membrane fluid boundary, the extents or characteristics of the fluid, or alternatively it may provide a means to measure the characteristics of an elastic fluid with or without hysteresis of response, whereby the fluid has an offset in displacement response, or "memory", such that the response behavior in one direction is similar to the response behavior in the opposite direction, but only after traveling a particular displacement distance.
• a fluid elasticity metric may be determined from the characteristic response of the surface or membrane to the surface excitation and reflected ultrasound characterization.
• the excitation generator may be a air bladder manipulated by an operator to apply a force to a membrane or surface, an air displacement generator producing alternating pressure, step pressure, or air puffs.
• the excitation generator output may be sealed to the surrounding region of the surface or unsealed and using a puff of gas such as atmospheric air or other suitable gas.
• the estimate of surface deflection may be derived from a measured estimate of velocity, an acceleration, or any other metric associated with deflection over time.
• the excitation is a step or impulse response with a rising edge, falling edge, or impulsive excitation which includes decomposed frequency content which is greater than the expected response frequency of the system, and the estimate of elasticity or viscosity is performed by examination of a ringdown characteristic.
• the ringdown characteristic comprises at least one of an exponential decay time or a ring cycle interval or frequency, such as the decomposition of a response into a ringdown characteristic, such as:
• f( ⁇ ) is the captured phase for a series of measurements
• t is the exponential decay coefficient
• f is the ring cycle frequency
• t is time.
• the excitation generator is an air puff directed through a nozzle to a membrane surface, the membrane having an underlying liquid to be characterized.
• the liquid may be a semi-solid food item such as yogurt, cheese, or other fermented item, and the elasticity measurement may provide a measurement of completion of a fermentation or other time-dependent process.
• the excitation generator is an air puff coupled to the surface of a semi-solid food item being fermented or cured, such as one using a bacterial or other time-dependent process, and the ultrasound transducer periodically receives a reflected signal from an ultrasound transmit event, where the received amplitude or phase during the relaxation time of the food item after the air puff excitation is compared to the originally transmitted signal to characterize the time-dependent displacement characteristic of the food item to derive an elasticity or other physical property based on the relaxation time change in response to the excitation.
• the excitation generator applies an impulsive pressure to the surface of a food item such as a vegetable or fruit, the ultrasound energy applied to the food item being reflected and characterized to measure the time dependent surface response of the fruit or vegetable, thereby using the reflected phase or amplitude response during a surface relaxation time to determine an elasticity or other physical property which may be correlated to the ripeness of the fruit or vegetable.
• the food item may be placed into a holder and the surface excited with a puff of gas such as air, the surface deflection response estimating ripeness or other property.
• the excitation may be a gas which is delivered at a supersonic velocity and/or at a glancing angle to the surface of the food item, or one or more food items may be placed into a chamber which has a variable pressure to measure a low frequency surface response to pressure, such as deflection vs. pressure.
• the excitation may be applied to one surface and the response measured on a different surface of the same item, such as the measurement of a propagating surface wave or a shear wave which travels through the item being characterized.
• the excitation generator applies an impulsive pressure to an animal or human organ such as the eye during an interval when ultrasound energy is applied to the eye, the phase changes of the reflected ultrasound during application of the non- contact excitation or after removal of the non-contact excitation indicating an elasticity which may be correlated to an inter-ocular pressure for measurement or diagnosis of glaucoma.
• the excitation may be an air puff containing nitrogen (N2) or other non-reactive gas such as Argon, Carbon Dioxide (C02), Helium (He), or an inert gas to prevent oxidation of the membrane of surface to be characterized.
• N2 nitrogen
• C02 Carbon Dioxide
• He Helium
• the acoustic wave velocity will have a different value compared to the velocity of air.
• the ultrasound transducer is formed by a plurality of membrane surfaces which are each electrostatically excited by a transmit voltage waveform to cause the membrane surface to form and transmit an acoustic wave, such as a capacitive micromachined ultrasound transducer (CMUT).
• the applied voltage waveform may also include a direct current (DC) component for positioning the membrane surface in a desired range of motion, while the alternating current (AC) component of the applied voltage waveform may provide movement of the membrane for generation of ultrasound energy.
• DC direct current
• AC alternating current
• the ultrasound transducer may be a ceramic piezoelectric transducer which changes mechanical shape based on the applied voltage.
• a signal processor for detection of air or fluid behind a surface, and further estimating an elasticity metric of a fluid when present, the signal processor comprising:
• an ultrasound transducer for coupling ultrasound energy to a surface
• an excitation generator producing sub-audio, audio, or super-audio excitation coupled into a surface to me measured to cause a measurable deflection in a surface
• a transmitter coupled to the ultrasound transducer during a transmit interval
• a receiver coupled to the ultrasound transducer during a receive interval which follows the transmit interval
• phase and/or amplitude analyzer comparing the phase of a transmit signal of the transmit interval to a phase and/or amplitude of a receive signal during the receive interval to estimate a surface deflection
• the signal processor deriving a metric from the phase and amplitude analyzer by comparing the surface deflection with the excitation generator output; the elasticity metric indicating whether the receive signal is a reflection from a surface which includes reflections from air or from fluid, and optionally characterizing a fluid from the surface deflection characteristic.
• a signal processor where the excitation generator generates at least one of: sinusoidal, impulse, steady state, or periodic sub-audio, audio, or super-audio excitation.
• a signal processor where the phase and amplitude analyzer is operative on received acoustic energy from the transducer at a natural center frequency of the transducer.
• a signal processor where the phase and amplitude analyzer is operative on received acoustic energy from the transducer at a baseband frequency spectrum, the baseband frequency spectrum formed by mixing the receive signal with a carrier frequency which is at substantially the center frequency of the transmitter.
• a signal processor where the transmitter generates a transmit waveform which includes an excitation voltage signal at a center frequency of the transducer during the transmit interval.
• a signal processor where the sum of the transmit interval and the receive interval is greater than 50 microseconds and less than 1 millisecond.
• a signal processor where the phase and amplitude analyzer determines a weighted or unweighted average phase with respect to a transmit clock.
• a signal processor where a metric is derived from a temporal phase change between a received signal from the transducer during the receive signal interval and a transmit clock which is operative during the receive interval.
• a signal processor where a metric is derived from a phase relationship between a mixer output baseband signal and the excitation generator output.
• a signal processor where a metric is derived from a temporal phase change in the receive signal and an excitation generator output.
• a signal processor where the ultrasound transducer generates a periodic burst of transmit signal energy.
• a signal processor where the ultrasound transducer generates continuous transmit signal energy.
• a signal processor where the phase and amplitude analyzer is operative on received signals to identify a region of first reflection from a membrane or surface, and thereafter characterizes a fluid behind an identified region as either air or liquid.
• a signal processor where, when the fluid behind an identified region is liquid, determines a surface elasticity, or fluid viscosity using a phase and amplitude response associated with a measurable deflection. [0086] A signal processor for characterizing a temporal response from a surface or membrane, the signal processor having:
• an excitation generator producing sub-audio, audio, or super-audio excitation for application to a surface or membrane to cause a displacement
• a transducer for launching acoustic waves towards a surface to be characterized and receiving reflections from the surface
• a visual indicator to allow the direction of acoustic waves from the transducer to a region of interest on the surface
• an ultrasound transmitter operative during a transmit interval and coupling a gated frequency burst to the transducer
• an ultrasound receiver operative during a receive interval and coupled to the transducer; a phase and amplitude detector comparing the phase of a transmit clock to a receive signal from the ultrasound receiver and generating a phase output;
• a response analyzer comparing the phase output to the excitation generator output, the response analyzer determining a viscosity of a fluid adjacent to surface to be characterized by comparison of the phase output and the excitation generator output.
• a signal processor where the transducer is at least one of a capacitive micro-machined ultrasound transducer (cMUT) or a piezoelectric transducer.
• cMUT capacitive micro-machined ultrasound transducer
• piezoelectric transducer a piezoelectric transducer
• a signal processor where the excitation generator is at least one of a voice coil actuator, or a moving diaphragm.
• a signal processor where the visual guide is at least one of: a laser diode, light emitting diode, or optical indicator which illuminates a region corresponding to a beam profile from the ultrasonic transducer.
• a signal processor where the ultrasound transmitter has a repetition rate of less than l5KHz.
• a signal processor where the phase and amplitude detector is a baseband mixer generating an output after a low pass filter.
• a signal processor where the phase and amplitude detector is operative at a center frequency of the transducer.
• a signal processor where the response analyzer compares the phase output and the excitation generator output over a plurality of sample points over a duration of time when the excitation generator is operative.

Landscapes

• Engineering & Computer Science (AREA)
• Physics & Mathematics (AREA)
• Radar, Positioning & Navigation (AREA)
• Remote Sensing (AREA)
• Computer Networks & Wireless Communication (AREA)
• General Physics & Mathematics (AREA)
• Acoustics & Sound (AREA)
• Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
• Ultra Sonic Daignosis Equipment (AREA)
• Length Measuring Devices Characterised By Use Of Acoustic Means (AREA)

Applications Claiming Priority (1)

Publications (2)

Family

ID=69141644

Family Applications (1)

Country Status (7)

Families Citing this family (4)

Family Cites Families (18)

• 2018
  • 2018-07-13 AU AU2018431763A patent/AU2018431763A1/en active Pending
  • 2018-07-13 KR KR1020237036251A patent/KR20230153501A/ko not_active Application Discontinuation
  • 2018-07-13 WO PCT/US2018/042138 patent/WO2020013868A1/en unknown
  • 2018-07-13 CA CA3105427A patent/CA3105427A1/en active Pending
  • 2018-07-13 CN CN201880097310.6A patent/CN112673361A/zh active Pending
  • 2018-07-13 EP EP18925875.9A patent/EP3821347A4/de active Pending
  • 2018-07-13 KR KR1020217003198A patent/KR102594614B1/ko active IP Right Grant
  • 2018-07-13 JP JP2021500897A patent/JP7249400B2/ja active Active
• 2023
  • 2023-03-17 JP JP2023042986A patent/JP2023063592A/ja active Pending

Also Published As

Similar Documents

Legal Events