CA2021509A1 - Method of processing an electrical signal - Google Patents
Method of processing an electrical signalInfo
- Publication number
- CA2021509A1 CA2021509A1 CA002021509A CA2021509A CA2021509A1 CA 2021509 A1 CA2021509 A1 CA 2021509A1 CA 002021509 A CA002021509 A CA 002021509A CA 2021509 A CA2021509 A CA 2021509A CA 2021509 A1 CA2021509 A1 CA 2021509A1
- Authority
- CA
- Canada
- Prior art keywords
- signal
- set forth
- frequency
- transformation
- envelope
- 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.)
- Abandoned
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0858—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving measuring tissue layers, e.g. skin, interfaces
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B17/00—Measuring arrangements characterised by the use of infrasonic, sonic or ultrasonic vibrations
- G01B17/02—Measuring arrangements characterised by the use of infrasonic, sonic or ultrasonic vibrations for measuring thickness
Abstract
METHOD OF PROCESSING AN ELECTRICAL SIGNAL
ABSTRACT OF THE DISCLOSURE
This method for measuring the distance e separating two faces (3, 4) of an object (1) employs the ultrasonic wave diffused back by said faces and received by a sensor (6). Following digitalization (11), the temporal signal is transformed (12) into a frequency spectrum . The modulus of the frequency spectrum is established (13) then the logarithmic derivative thereof is calculated (15) from whence there results a useful signal including a periodic signal onto which a linear signal is superposed. Following suppression (14, 161 of the linear signal the useful signal is applied (17) to a frequency window in order to eliminate the noise therefrom. The signal thus obtained is inversely transformed (18) to obtain a signal the amplitude of which is a function of the time and the envelope of this latter signal is determined (19). After having determined the temporal spread .DELTA.t separating the time origin from the envelope peak showing the greatest amplitude, the sought for distance e is calculated from the relation e = c?.DELTA.t/2, c representing the propagation velocity of the ultrasonic wave. The method may be employed for the non-invasive meas-urement of the thickness or diameter of an artery or indeed for the measurement of the thickness of thin layers.
ABSTRACT OF THE DISCLOSURE
This method for measuring the distance e separating two faces (3, 4) of an object (1) employs the ultrasonic wave diffused back by said faces and received by a sensor (6). Following digitalization (11), the temporal signal is transformed (12) into a frequency spectrum . The modulus of the frequency spectrum is established (13) then the logarithmic derivative thereof is calculated (15) from whence there results a useful signal including a periodic signal onto which a linear signal is superposed. Following suppression (14, 161 of the linear signal the useful signal is applied (17) to a frequency window in order to eliminate the noise therefrom. The signal thus obtained is inversely transformed (18) to obtain a signal the amplitude of which is a function of the time and the envelope of this latter signal is determined (19). After having determined the temporal spread .DELTA.t separating the time origin from the envelope peak showing the greatest amplitude, the sought for distance e is calculated from the relation e = c?.DELTA.t/2, c representing the propagation velocity of the ultrasonic wave. The method may be employed for the non-invasive meas-urement of the thickness or diameter of an artery or indeed for the measurement of the thickness of thin layers.
Description
2 ~
The present invention relates to a method for measuring the distance e separating two faces of an object, such method including the emitting of an ultrasonic wave propagating at the velocity c in the direction of the object and the reception by a sensor of the echoes generated by reflection of said wave from the aces of said object, said sensor delivering a composite electrical signal g(t) exhibiting an amplitude varying as a function of the time.
BACKG~OUND OF THE INVEN~ION
Several methods have already been proposed which employ the ultrasonic way in order to measure the distance separating two faces of an object. If the distance separating the two faces is substantial, there is in ~eneral no difficulty in interpreting the double signal received in the form of an echo since the faces generate two peaks which are clearly separated in time and which permit, in knowing the speed of propagation of the wave in the medium in which it is propagated, to calculate the distance which separates the faces in question. However, as soon as said faces are very close to one another and spaced a distance on the order of the wave length from one another, the echoes which they produce may appear merged.
Coating thicknesses applied to a substrate which thicknesses comprise between 5 and 20 ~m, have been apt to be measured by ultrasonic interferometry employing a frequency varying from 90 to 510 MHz. Such a method is described in the article entitled "Measurement of the Thickness of Thin Layers by Ultrasonic Interferometry"
appearing in J. Appl. Phys. 55 (1), January, 1984, pages 194-198. However, this method furnishes a satisfactory result only if it is applied to a sta~le object, the dimensions of which are not subject to variations during the measurement.
2~ 5~
The article entitled "Digital Signal Processing Method for Multilayered Media Thickness Measurement" which appeared in IEEE Ultrasonic Symposium 1988, pages 1027-1029, describes an ultrasonic measurement carried out on thin layers of paint applied in the automobile industry.
In this case the layers exhibit thicknesses between 20 and 110 ~m. The central frequency here is around 75 M~z with a band width of 100 MHz. The measurement method is based on the analysis of a power cepstrum, defined as being the inverse Fourier transformation of the logarithm of the power spectral density of the radio frequency signal.
Here however one is in the presence of superficial layers which permit working at a high frequency from whence derives a result enabling good resolution. This method would not suit for the measurement of the thickness of an object situated within a medium attenuating the ultrasonic wave and which would thus oblige working at lower frequencies (10 MHz) as will be the case in a possible employment of this invention.
The cepstrum method is likewise the object of a thesis presented 9th February, 1979 at the National Polytechnical Institute of Grenoble by J.C. Balluet, which thesis is entitled "Les opérateurs cepstres, applications à la séparation d'échos rapprochés". In this publication it will be noted however that the method employed for echoes which are very close to one another gives results which are difficult to interpret and that the observed peaks are located at the limit of visibility.
An ultrasonic echography method applied in vivo on living tissues is described in the article entitled "~ange Resolution Improvement by a Fast Deconvolution Method", which appeared in Ultrasonic Imaging 6, 1984, pages 435-451. Numerous problems of mathematical physics may be formulated by means of an equation of convolution of the form y = h * x, where x and y denote a system of 2 ~
input and outp~t respectively and where h is the pulse response of the system. The deconvolution is the inverse problem which consists in reconstructing the input on the basis of an output measured experimentally. Thanks to the method set forth in this article, it has been possible to study with precision the bottom of the eye and to follow the variations of the thickness of an arterial wall during the cardiac cycle. The method proposed however presents the disadvantage of separatinq only with great difficulty echoes which are very close in time.
Non-invasive sensors enabling the measurement of the interior diameter of an artery are known from patent US-A 4 370 985. The indicated method gives no solution for measuring the actual thickness of the artery.
SUM ARY OF THE INVENTION
From the contents of the literature briefly set forth hereinabove, it is apparent that the methods employed either are poorly adapted to the medium in which the measurement must be eff~cted, or give results which are difficult to interpret when it concerns measurement of very small thic~nesses. To overcome these difficulties, the present invention is characterized by the fact that the composite electrical signal received by the sensor is subjected to a processing which includes the following sequence of steps:
a) one transforms the composite electrical signal g(t) coming from the sensor into a fre~uency spectrum G(f) the amplitude of which is a function of the frequency, said spectrum exhibiting a real part RetG(f)3 and an imaginary part Im[G~f)~, b~ one calculates the order n modulus of the frequency spectrum ¦G~f)¦ = ~RetG~f~3 ~ Im~G(f)] } I , c) one calculates the logarithmic derivative ( ¦G~f)¦n) t¦G~f)¦n of the modulus calculated in step b) in order to separate the existing echoes from the composite electrical signal, said logarithmic derivative exhibiting a useful part centered substantially on the frequency of the ultrasonic wave and a part having noise located on either side of the useful part, said useful part comprising a periodic signal onto which a substantially linear signal is superposed, d) one applies a frequency window to the logarithmic derivative in order to isolate the useful part thereof and to eliminate the part having noise, e) one applies a transformation inverse to that of step a), at least to the periodic signal of said useful part, f) one determines the envelope of the inverse transformation, said envelope exhibiting a series of variable amplitude peaks spaced out over time, g) one determines the temporal spread ~t separating the time origin from the peak showing the greatest amplitude due to the contribution of the periodic signal from the useful part, the sought for distance e separating the two faces of the object being expressed by the relation e =
c ~t/2.
The invention will now be understood from reading the following description given by way of example and illustrated ~y the drawings.
BRIEF DESCRIPTION OF THE DRAWINGS
- Figure 1 is a schematic view of an example of the measurement arrangement showing an object for which it is desired to measure the thickness of the wall, and a sensor delivering a composite electrical signal to a computer, this signal undergoing the processing which is the obiect of this invention;
- figure 2 is a flow chart showing the state of the various steps o~ the method according to the invention - 6 - 2~
and according to one example of processing of the signal;
- figure 3 shows the composite signal delivered by the sensor when the thickness of the wall is on the order 0.75 mm, the wall being here that of a tube of plastic material;
- figure 4 shows the same signal for a thickness of 0.3 mm;
- figure 5 shows the same signal for a thickness of 0.1~ mm;
- figure 6 is a diagram showing the real part of the frequency spectrum obtained in transforming the signal issuing from the sensor for the signal of figure 5;
- figure 7 shows the imaginary part of the frequency spectrum obtained in transforming the signal issuing from the sensor for the signal of figure 5;
- figure 8 is a diagram showing the squared modulus of the frequency spectrum shown on figures 6 and 7;
- figure 9 is a diagram showing the logarithmic derivative of the squared modulus shown on figure 8;
- figure 10 is a diagram similar to that of figure 9 in which a curve adjustment has been carried out;
- fi~ure 11 is a diagram showing the signal of figure 10 after having applied thereto a frequency window;
- figure 12 is a representation of the inverse transformation of the signal of figure 11;
- figure 13 is a diagram where there has been represented the envelope of the signal of figure 12 with an expanded abscissa scale, which envelope enables calculating the thickness of the tube wall;
- figure 14 shows the envelope diagram resulting from a similar measurement applied in vitro onto the artery of an animal, and - figure 15 is a diagram obtained under the same conditions as those of figure 14, but for which the step shown on figure 10 has not been carried out.
DESCRIPTION OF THE PREFERRED EMBODIMENT
Figure 1 shows in cross-section an object 1 for which the thickness of the wall 2 is to be measured. The wall is bounded by an external face 3 and an internal face 4. In order to proceed with this measurement, there is emitted an ultrasonic wave S by means of an ultrasonic sensor emitter-receiver 6, each wave being propagated at velocity c in the direction of the object. In emission, the sensor is excited by an oscillator 7, the frequency of which is on the order of 10 MHz. In the example, sensor 6 likewise serves as a receiver of the echoes generated by reflection of the wave from the faces 3 and 4 of the object. This sensor delivers a composite electrical signal g(t) which exhibits an amplitude A varying as a function of time t. This electrical signal is subjected to processing according to a method which is the object of this invention.
All the processing may be carried out by means of a personal computer (PC) represented at 8 on figure 1. At each step of the method it is possible by means of display 9 to visualize the calculation effected. It will be here mentioned that the entire method according to the invention may be carried out by an I~M AT personal computer to which there is added a type ~087 mathematical coprocessor.
Other apparatus could be employed, those which have just been indicated constituting only one possible example.
The flow chart of figure 2 shows different blocks which each illustrate one operation of the method. The disposition here shown is that which is the most complete.
The description which is to follow will however indicate that certain blocks can be eliminated without straying however from the principal idea of the invention.
Oscillator 7 excites a piezoelectric sensor 6 which emits an ultrasonic wave having a central frequency of - 8 - 2 ~
about 10 MHz. This same sensor gathers in the signal g(t) diffused back by faces 3 and 4 of wall 2 of object 1.
This signal exhibits an amplitude A varying as a function of the time t and is presented in an analog form at the output of sensor 6. Following amplification in the amplifier 10 it is convenient, but not absolutely necessary, to digitalize the signal g(t) by means of an analog-digital converter 11. In the case taken here as example, the signal is sampled at 100 MHz over 8 bits, then stored in a memory including 2048 points.
The object employed here in order to explain the method of the invention is a tube of plastic material for which one proposes to measure the thickness of the wall.
As one can measure this thickness by other well known means as for example by means of a micrometer, one will be able to confront the results of the measurement of the two systems which will permit verifying the results obtained according to the method of the invention.
The diagrams of figures 3, 4 and 5 each show the signal g(t) following digitalization and for three particular cases of the thickness of the tube wall. In these diagrams, there has been indicated the amplitude of the signal on the ordinate (relative value) and the time on the abscissa (value given in microseconds).
~ igure 3 shows the signal g(t) for a wall thickness e = 0.75 mm. Here one clearly distinguishes the two echoes 20 and 21 sent back by the exterior face 3 and interior face 4 of the wall. If one measures the time ~t on the diagram separating the two main peaks of the signal, one finds ~t = 1 ~s. Since one furthermore knows the velocity c of wave propagation which is equal in the particular case taken as example, to 1.5 mm/ ~s, one may immediately calculate the thickness of the wall which is equal to e = c ~t¦2, the division by 2 taking into account the sending 9 2~
and return of the wave. By replacing c and ~t in the formula by 1.5 and 1 respectively, one will find for the thickness e the value of 0.75 mm. One determines that in the case where echoes are clearly separated, the measure of the thickness is immediate and that it is not necessary to appeal to a method according to the invention.
Figure 4 shows the signal g~t) for a wall thickness o~ 0.3 mm. ~ere the two echoes are closer together and if one proceeds with the same operation as previously, one will be able to obtain a result corresponding more or less to reality according to the peaks which one will have chosen.
Figure 5 shows the signal g(t) for ~ wall thickness of 0.15 mm. Here the two peaks to be chosen are difficult, indeed impossible to determine. The echoes appear as merged and it is the purpose of this invention to propose a calculation method capable of isolating the peaks to be taken into consideration in order to measure the thickness of the wall. This method is now going to be applied to the case of the tube giving as initial signal g(t) that represented on figure 5. For this we will refer back to figure 2 as well as to the series of diagrams shown in figures S to 13.
Exposition of the Method a1 ~he composite electrical signal g(t) present at the output of block 11 and represented on figure 5 initially undergoes an operation in block 12, this operation consistinq in transforminq said signal dependinq from time into a frequency spectrum, the amplitude A of which depends on the frequency f. This operation is generally prece~ed by an offset correction in the case where the signal g~t) is provided therewith. In a completely general manner, it may concern here a o 2 ~ 2 ~ 3 transformation in z introduced by Hurewiez for the study of servomechanisms. On this subject reference may be had to the work entitled "Theory of Servomechanisms" published by McGraw-Hill Book Co., Inc. New York, chapter 5, 1947.
In the case which here interests us, there will be used to advantage a special case of the z transformation which is the well known Fourier transformation. Following transformation, and at the output of block 12, the signal G(f) appears in the form of a real part ~e[G(f)] shown on figure 6 and an imaginary part Im[G(f)~ shown on figure 7.
b) The two parts, real and imaginary, of the frequency spectrum G(f) are then introduced into block 13 (figure 2), the purpose of which is to calculate the n order modulus of said spectrum according to the relation ¦G(f)¦n = {Re[G(f)]2 + Im[G(f)]2}n/2. In the special case, one will employ to advantage the squared modulus of the frequency spectrum by making n = 2. This squared modulus is represented on figure 8 which shows a very pronounced peak over the central emission frequency of the ultrasonic sensor (about 10 MHz).
c) The modulus obtained at the output of block 13 then introduced into block 15 (figure 2) the role of which is the calculation of the logarithmic derivative of the modulus defined by the quotient of the derivative of the modulus ( ¦G(f) In) over the actual modulus ¦G(f)¦nO There it concerns the most important step of the invention thanks to which it will be possible to separate the composite electrical signal from the two echoes which are sought after. At the output of block 15, the logarithmic derivative appears in the form represented on figure 9 which shows distinctly two parts: a useful part 25, substantially centered on the ultrasonic wave frequency, and a noise portion 26 situated on either side of the useful part. The useful part 25 which is the interesting part, comprises a periodic signal 27 onto which is superposed a substantially linear signal 28. These two components of the useful part are physically explained:
the starting signal g(t) includes in effect a double peak to be detected, this double peak being mixed with the signal emitted by the sensor which may be assimilated to a sinusoidal signal modulated in amplitude by a gaussian.
The Fourier transformation G(f) of the signal g(t) further shows a mixture and when one calculates therefrom the logarithmic derivative, the peaks are transformed into a periodic signal 27 and the signal emitted by the sensor into a substantially linear signal 28. It will here be noted that this separation is easily obtainable in the frequency domain, while not being possible in the time domain which explains the transformation carried out in step a).
In a preferred but not essential version of the invention, one then suppresses the linear signal 28 which is present on diaqram 9 in the form of a slope. The suppression of this slope gives an advantage which will appear subsequently. Several methods may be employed in order to suppress this slope or linear signal, these methods giving rise to a signal represented on figure 10 and referenced G1(f).
A first method consists in visualizing on a screen the signal represented on figure 9, then proceeding with the suppression of the linear signal in rectifying this signal by means of trials or successive approximations until one obtains a graphically satisfactory result of the type shown on figure 1Q. This will be the case when the slope 28 will have been merged with the frequency axis. Such operation is symbolized by block 16 on figure 2.
A second method employs a mathematical way. One 2 $ ~
begins by calculating from the expression of the squared modulus (n = 2) obtained following step b), the moment ml of order 1 and the moment m2 of order 2 of said modulus.
These moments are given by the following relation in making n = 1 for m1 and n = 2 for m2:
i=~ .
IG(fi)l2 fin m =
~, I G(fj)l2 i=--One next determines the straight line passing by f = m1 and having for slope p = k - m2, which straight line is expressed by the relationship p(~-m1), where f is equal to the frequency and k is a constant. Finally, one subtracts the straight line defined hereinabove from the logarithmic derivative obtained in step c) in a manner to obtain:
G1(f) = (¦G(f)¦2) /¦G(f)¦2 (f In the flow chart of figure 2, bloc~ 14 represents the calculation of moments m1 and m2 of the modulus obtained at block 13 and block 16 represents the actual suppression of the slope or linear signal which is brought about by the subtraction indicated hereinabove.
A third method employs the characteristics of the substantially linear part known a priori by a reference measurement d) To the logarithmic derivative obtained in step c), one applies a frequency window in order to isolate the useful part 25 and to eliminate the noise part 26.
2 ~
_ ~3 _ This operation takes place immediately following step c) represented by block 15 in the case where one does not proceed to correction of the slope suggested hereinabove.
In the case where the correction of the slope is effected (block 16), the frequency window may be applied indifferently before or after this correction. In the flow chart of figure 2 which provides a slope correction, one opted for application of the window (block 17 of figure 2) after said correction (block 16) and the diagram of figure 11 shows the signal following application of the window. This window is a band-pass filter which, as is shown on figure 11, causes the disappearance of the noise part which is located on either side of the useful part.
Advantageously, the window shows a bell form of the type of window called Hanning (see the work entitled "The Measurement of Power Spectra" by R.B. Blackman and J.W.
Tuckey, Dover Publications New York 1958). Other windows may be employed and for this one may consult the work of M. Kunt entitled "Traitement numérique des signaux"
published by Dunod, 1981.
e) The signal following the window is rendered usable if one returns now into the time domain which is effected by block 18 on figure 2, which applies to the signal a transformation inverse to that which has been described at step a). By analogy to what has been said hereinabove, this operation may be very generally effected by an inverse transformation in z or, in the case taken here by way of example, by an inverse Fourier transformation. The signal following transformation is shown on figure 12. If one desi~nates by TF the inverse Fourier transformation, by G1(f) the logarithmic derivative of the squared modulus following correction of the slope and by H(f) the frequency window, the signal of figure 12 may be expressed by the relation:
A(t) = TF 1lG1(f) ~(f~]
2 ~ f~
f) The following step of the method consists in determining the envelope of the inverse transformation brought about in step e). This operation is effected by block 19 of figure 2 and appears as shown on figure 13.
As is shown by this latter figure, the envelope exhibits a series of peaks 30, 31, 32, etc. of variable amplitude and spaced out in time. In order to arrive at this, one will advantageously employ the Hilbert transformation as described in the work "Probability, Random Variables and Stochastic Process" of A. Papoulis published by McGraw-Hill, 1984. If one designates by TH the Hilbert transformation, the envelope of figure 13 is expressed by:
{[A(t)]2 + [TH A(t)]2~1/2 the value of A(t) being that obtained hereinabove.
g) From the envelope obtained in the preceding step, one determines the time spread ~t separating the time origin from peak 30 exhibiting the greatest amplitude, which peak is due to the contribution from the periodic signal of the useful part defined hereinabove, in recalling that it is the periodic signal cf the useful part in which are found the echoes sought after. Here the time origin is a relative origin which corresponds to the first echo. ~he interval ~t as measured corresponds then very well to the interval separating the two echoes. This origin differs from that represented on figure S. Such comes about from the fact that in employing the squared modulus calculated in step b) one loses the information which contains the time separating the time origin shown on figure 5 from the first echo.
In the special case, since the slope of the useful signal has been suppressed, the peak having the greatest 1 5 2 ~
amplitude due to the contribution of the periodic signal is indeed the first encountered and referenced 30, the distance e between the two faces 3 and 4 of the wall 2 of tube 1 being expressed by the relation e = c at/2.
On the graph of figure 13, one measures ~t = 0.2 ~s. If c = 1.5 mm/~s as indicated hereinabove, one finds e =
1.5- 0.2/2 = 0.15 mm which corresponds well to the thickness of the tube which one knew from other means of measurement was equal to 0.15 mm.
Additional Considerations The plastic tube employed hereinabove in order to explain the method of the invention gives well defined signals particularly due to the rigidity of the tube subjected to measurement. In numerous cases, the object for which the thickness is to be measured has no rigidity;
furthermore, it is generally inaccessible and does not allow measurement by mechanical means. This is the case for instance of an artery in vivo for which one wishes to measure in a non-invasive manner the thickness of the wall.
Figure 14 presents the signal envelope obtained in the next to last step of the method according to the invention from a swine femoral artery in vitro. In this case, there has been effected a slope correction as indicated hereinabove in step c). The envelope shows several peaks and without any contestation it is the peak of the greatest amplitude noted at 0.56 ~s which must be taken into consideration.
The diagram of figure 15 is an envelope obtained from the same artery but where the slope correction has not been effected. One readily finds at 0.56 ~s a peak corresponding to that sought after but this peak is preceded by another peak 40 of greater amplitude, such peak originat-2 ~ 2 ~
ing from the contribution of the linear part (or from the slope) of the useful signal which h~re is preponderant.
This teaching militates in favour of the systematic employment of the slope correction in all cases where the object to be measured shows little stability and/or is not accessible.
It has already been mentioned that the method may be employed for the measurement of the thickness of an artery. The same method may be profitably employed for measuring the interior diameter of the artery and generally for the measurement of any object exhibiting a thin layer for which one is concerned with measuring the thickness.
The present invention relates to a method for measuring the distance e separating two faces of an object, such method including the emitting of an ultrasonic wave propagating at the velocity c in the direction of the object and the reception by a sensor of the echoes generated by reflection of said wave from the aces of said object, said sensor delivering a composite electrical signal g(t) exhibiting an amplitude varying as a function of the time.
BACKG~OUND OF THE INVEN~ION
Several methods have already been proposed which employ the ultrasonic way in order to measure the distance separating two faces of an object. If the distance separating the two faces is substantial, there is in ~eneral no difficulty in interpreting the double signal received in the form of an echo since the faces generate two peaks which are clearly separated in time and which permit, in knowing the speed of propagation of the wave in the medium in which it is propagated, to calculate the distance which separates the faces in question. However, as soon as said faces are very close to one another and spaced a distance on the order of the wave length from one another, the echoes which they produce may appear merged.
Coating thicknesses applied to a substrate which thicknesses comprise between 5 and 20 ~m, have been apt to be measured by ultrasonic interferometry employing a frequency varying from 90 to 510 MHz. Such a method is described in the article entitled "Measurement of the Thickness of Thin Layers by Ultrasonic Interferometry"
appearing in J. Appl. Phys. 55 (1), January, 1984, pages 194-198. However, this method furnishes a satisfactory result only if it is applied to a sta~le object, the dimensions of which are not subject to variations during the measurement.
2~ 5~
The article entitled "Digital Signal Processing Method for Multilayered Media Thickness Measurement" which appeared in IEEE Ultrasonic Symposium 1988, pages 1027-1029, describes an ultrasonic measurement carried out on thin layers of paint applied in the automobile industry.
In this case the layers exhibit thicknesses between 20 and 110 ~m. The central frequency here is around 75 M~z with a band width of 100 MHz. The measurement method is based on the analysis of a power cepstrum, defined as being the inverse Fourier transformation of the logarithm of the power spectral density of the radio frequency signal.
Here however one is in the presence of superficial layers which permit working at a high frequency from whence derives a result enabling good resolution. This method would not suit for the measurement of the thickness of an object situated within a medium attenuating the ultrasonic wave and which would thus oblige working at lower frequencies (10 MHz) as will be the case in a possible employment of this invention.
The cepstrum method is likewise the object of a thesis presented 9th February, 1979 at the National Polytechnical Institute of Grenoble by J.C. Balluet, which thesis is entitled "Les opérateurs cepstres, applications à la séparation d'échos rapprochés". In this publication it will be noted however that the method employed for echoes which are very close to one another gives results which are difficult to interpret and that the observed peaks are located at the limit of visibility.
An ultrasonic echography method applied in vivo on living tissues is described in the article entitled "~ange Resolution Improvement by a Fast Deconvolution Method", which appeared in Ultrasonic Imaging 6, 1984, pages 435-451. Numerous problems of mathematical physics may be formulated by means of an equation of convolution of the form y = h * x, where x and y denote a system of 2 ~
input and outp~t respectively and where h is the pulse response of the system. The deconvolution is the inverse problem which consists in reconstructing the input on the basis of an output measured experimentally. Thanks to the method set forth in this article, it has been possible to study with precision the bottom of the eye and to follow the variations of the thickness of an arterial wall during the cardiac cycle. The method proposed however presents the disadvantage of separatinq only with great difficulty echoes which are very close in time.
Non-invasive sensors enabling the measurement of the interior diameter of an artery are known from patent US-A 4 370 985. The indicated method gives no solution for measuring the actual thickness of the artery.
SUM ARY OF THE INVENTION
From the contents of the literature briefly set forth hereinabove, it is apparent that the methods employed either are poorly adapted to the medium in which the measurement must be eff~cted, or give results which are difficult to interpret when it concerns measurement of very small thic~nesses. To overcome these difficulties, the present invention is characterized by the fact that the composite electrical signal received by the sensor is subjected to a processing which includes the following sequence of steps:
a) one transforms the composite electrical signal g(t) coming from the sensor into a fre~uency spectrum G(f) the amplitude of which is a function of the frequency, said spectrum exhibiting a real part RetG(f)3 and an imaginary part Im[G~f)~, b~ one calculates the order n modulus of the frequency spectrum ¦G~f)¦ = ~RetG~f~3 ~ Im~G(f)] } I , c) one calculates the logarithmic derivative ( ¦G~f)¦n) t¦G~f)¦n of the modulus calculated in step b) in order to separate the existing echoes from the composite electrical signal, said logarithmic derivative exhibiting a useful part centered substantially on the frequency of the ultrasonic wave and a part having noise located on either side of the useful part, said useful part comprising a periodic signal onto which a substantially linear signal is superposed, d) one applies a frequency window to the logarithmic derivative in order to isolate the useful part thereof and to eliminate the part having noise, e) one applies a transformation inverse to that of step a), at least to the periodic signal of said useful part, f) one determines the envelope of the inverse transformation, said envelope exhibiting a series of variable amplitude peaks spaced out over time, g) one determines the temporal spread ~t separating the time origin from the peak showing the greatest amplitude due to the contribution of the periodic signal from the useful part, the sought for distance e separating the two faces of the object being expressed by the relation e =
c ~t/2.
The invention will now be understood from reading the following description given by way of example and illustrated ~y the drawings.
BRIEF DESCRIPTION OF THE DRAWINGS
- Figure 1 is a schematic view of an example of the measurement arrangement showing an object for which it is desired to measure the thickness of the wall, and a sensor delivering a composite electrical signal to a computer, this signal undergoing the processing which is the obiect of this invention;
- figure 2 is a flow chart showing the state of the various steps o~ the method according to the invention - 6 - 2~
and according to one example of processing of the signal;
- figure 3 shows the composite signal delivered by the sensor when the thickness of the wall is on the order 0.75 mm, the wall being here that of a tube of plastic material;
- figure 4 shows the same signal for a thickness of 0.3 mm;
- figure 5 shows the same signal for a thickness of 0.1~ mm;
- figure 6 is a diagram showing the real part of the frequency spectrum obtained in transforming the signal issuing from the sensor for the signal of figure 5;
- figure 7 shows the imaginary part of the frequency spectrum obtained in transforming the signal issuing from the sensor for the signal of figure 5;
- figure 8 is a diagram showing the squared modulus of the frequency spectrum shown on figures 6 and 7;
- figure 9 is a diagram showing the logarithmic derivative of the squared modulus shown on figure 8;
- figure 10 is a diagram similar to that of figure 9 in which a curve adjustment has been carried out;
- fi~ure 11 is a diagram showing the signal of figure 10 after having applied thereto a frequency window;
- figure 12 is a representation of the inverse transformation of the signal of figure 11;
- figure 13 is a diagram where there has been represented the envelope of the signal of figure 12 with an expanded abscissa scale, which envelope enables calculating the thickness of the tube wall;
- figure 14 shows the envelope diagram resulting from a similar measurement applied in vitro onto the artery of an animal, and - figure 15 is a diagram obtained under the same conditions as those of figure 14, but for which the step shown on figure 10 has not been carried out.
DESCRIPTION OF THE PREFERRED EMBODIMENT
Figure 1 shows in cross-section an object 1 for which the thickness of the wall 2 is to be measured. The wall is bounded by an external face 3 and an internal face 4. In order to proceed with this measurement, there is emitted an ultrasonic wave S by means of an ultrasonic sensor emitter-receiver 6, each wave being propagated at velocity c in the direction of the object. In emission, the sensor is excited by an oscillator 7, the frequency of which is on the order of 10 MHz. In the example, sensor 6 likewise serves as a receiver of the echoes generated by reflection of the wave from the faces 3 and 4 of the object. This sensor delivers a composite electrical signal g(t) which exhibits an amplitude A varying as a function of time t. This electrical signal is subjected to processing according to a method which is the object of this invention.
All the processing may be carried out by means of a personal computer (PC) represented at 8 on figure 1. At each step of the method it is possible by means of display 9 to visualize the calculation effected. It will be here mentioned that the entire method according to the invention may be carried out by an I~M AT personal computer to which there is added a type ~087 mathematical coprocessor.
Other apparatus could be employed, those which have just been indicated constituting only one possible example.
The flow chart of figure 2 shows different blocks which each illustrate one operation of the method. The disposition here shown is that which is the most complete.
The description which is to follow will however indicate that certain blocks can be eliminated without straying however from the principal idea of the invention.
Oscillator 7 excites a piezoelectric sensor 6 which emits an ultrasonic wave having a central frequency of - 8 - 2 ~
about 10 MHz. This same sensor gathers in the signal g(t) diffused back by faces 3 and 4 of wall 2 of object 1.
This signal exhibits an amplitude A varying as a function of the time t and is presented in an analog form at the output of sensor 6. Following amplification in the amplifier 10 it is convenient, but not absolutely necessary, to digitalize the signal g(t) by means of an analog-digital converter 11. In the case taken here as example, the signal is sampled at 100 MHz over 8 bits, then stored in a memory including 2048 points.
The object employed here in order to explain the method of the invention is a tube of plastic material for which one proposes to measure the thickness of the wall.
As one can measure this thickness by other well known means as for example by means of a micrometer, one will be able to confront the results of the measurement of the two systems which will permit verifying the results obtained according to the method of the invention.
The diagrams of figures 3, 4 and 5 each show the signal g(t) following digitalization and for three particular cases of the thickness of the tube wall. In these diagrams, there has been indicated the amplitude of the signal on the ordinate (relative value) and the time on the abscissa (value given in microseconds).
~ igure 3 shows the signal g(t) for a wall thickness e = 0.75 mm. Here one clearly distinguishes the two echoes 20 and 21 sent back by the exterior face 3 and interior face 4 of the wall. If one measures the time ~t on the diagram separating the two main peaks of the signal, one finds ~t = 1 ~s. Since one furthermore knows the velocity c of wave propagation which is equal in the particular case taken as example, to 1.5 mm/ ~s, one may immediately calculate the thickness of the wall which is equal to e = c ~t¦2, the division by 2 taking into account the sending 9 2~
and return of the wave. By replacing c and ~t in the formula by 1.5 and 1 respectively, one will find for the thickness e the value of 0.75 mm. One determines that in the case where echoes are clearly separated, the measure of the thickness is immediate and that it is not necessary to appeal to a method according to the invention.
Figure 4 shows the signal g~t) for a wall thickness o~ 0.3 mm. ~ere the two echoes are closer together and if one proceeds with the same operation as previously, one will be able to obtain a result corresponding more or less to reality according to the peaks which one will have chosen.
Figure 5 shows the signal g(t) for ~ wall thickness of 0.15 mm. Here the two peaks to be chosen are difficult, indeed impossible to determine. The echoes appear as merged and it is the purpose of this invention to propose a calculation method capable of isolating the peaks to be taken into consideration in order to measure the thickness of the wall. This method is now going to be applied to the case of the tube giving as initial signal g(t) that represented on figure 5. For this we will refer back to figure 2 as well as to the series of diagrams shown in figures S to 13.
Exposition of the Method a1 ~he composite electrical signal g(t) present at the output of block 11 and represented on figure 5 initially undergoes an operation in block 12, this operation consistinq in transforminq said signal dependinq from time into a frequency spectrum, the amplitude A of which depends on the frequency f. This operation is generally prece~ed by an offset correction in the case where the signal g~t) is provided therewith. In a completely general manner, it may concern here a o 2 ~ 2 ~ 3 transformation in z introduced by Hurewiez for the study of servomechanisms. On this subject reference may be had to the work entitled "Theory of Servomechanisms" published by McGraw-Hill Book Co., Inc. New York, chapter 5, 1947.
In the case which here interests us, there will be used to advantage a special case of the z transformation which is the well known Fourier transformation. Following transformation, and at the output of block 12, the signal G(f) appears in the form of a real part ~e[G(f)] shown on figure 6 and an imaginary part Im[G(f)~ shown on figure 7.
b) The two parts, real and imaginary, of the frequency spectrum G(f) are then introduced into block 13 (figure 2), the purpose of which is to calculate the n order modulus of said spectrum according to the relation ¦G(f)¦n = {Re[G(f)]2 + Im[G(f)]2}n/2. In the special case, one will employ to advantage the squared modulus of the frequency spectrum by making n = 2. This squared modulus is represented on figure 8 which shows a very pronounced peak over the central emission frequency of the ultrasonic sensor (about 10 MHz).
c) The modulus obtained at the output of block 13 then introduced into block 15 (figure 2) the role of which is the calculation of the logarithmic derivative of the modulus defined by the quotient of the derivative of the modulus ( ¦G(f) In) over the actual modulus ¦G(f)¦nO There it concerns the most important step of the invention thanks to which it will be possible to separate the composite electrical signal from the two echoes which are sought after. At the output of block 15, the logarithmic derivative appears in the form represented on figure 9 which shows distinctly two parts: a useful part 25, substantially centered on the ultrasonic wave frequency, and a noise portion 26 situated on either side of the useful part. The useful part 25 which is the interesting part, comprises a periodic signal 27 onto which is superposed a substantially linear signal 28. These two components of the useful part are physically explained:
the starting signal g(t) includes in effect a double peak to be detected, this double peak being mixed with the signal emitted by the sensor which may be assimilated to a sinusoidal signal modulated in amplitude by a gaussian.
The Fourier transformation G(f) of the signal g(t) further shows a mixture and when one calculates therefrom the logarithmic derivative, the peaks are transformed into a periodic signal 27 and the signal emitted by the sensor into a substantially linear signal 28. It will here be noted that this separation is easily obtainable in the frequency domain, while not being possible in the time domain which explains the transformation carried out in step a).
In a preferred but not essential version of the invention, one then suppresses the linear signal 28 which is present on diaqram 9 in the form of a slope. The suppression of this slope gives an advantage which will appear subsequently. Several methods may be employed in order to suppress this slope or linear signal, these methods giving rise to a signal represented on figure 10 and referenced G1(f).
A first method consists in visualizing on a screen the signal represented on figure 9, then proceeding with the suppression of the linear signal in rectifying this signal by means of trials or successive approximations until one obtains a graphically satisfactory result of the type shown on figure 1Q. This will be the case when the slope 28 will have been merged with the frequency axis. Such operation is symbolized by block 16 on figure 2.
A second method employs a mathematical way. One 2 $ ~
begins by calculating from the expression of the squared modulus (n = 2) obtained following step b), the moment ml of order 1 and the moment m2 of order 2 of said modulus.
These moments are given by the following relation in making n = 1 for m1 and n = 2 for m2:
i=~ .
IG(fi)l2 fin m =
~, I G(fj)l2 i=--One next determines the straight line passing by f = m1 and having for slope p = k - m2, which straight line is expressed by the relationship p(~-m1), where f is equal to the frequency and k is a constant. Finally, one subtracts the straight line defined hereinabove from the logarithmic derivative obtained in step c) in a manner to obtain:
G1(f) = (¦G(f)¦2) /¦G(f)¦2 (f In the flow chart of figure 2, bloc~ 14 represents the calculation of moments m1 and m2 of the modulus obtained at block 13 and block 16 represents the actual suppression of the slope or linear signal which is brought about by the subtraction indicated hereinabove.
A third method employs the characteristics of the substantially linear part known a priori by a reference measurement d) To the logarithmic derivative obtained in step c), one applies a frequency window in order to isolate the useful part 25 and to eliminate the noise part 26.
2 ~
_ ~3 _ This operation takes place immediately following step c) represented by block 15 in the case where one does not proceed to correction of the slope suggested hereinabove.
In the case where the correction of the slope is effected (block 16), the frequency window may be applied indifferently before or after this correction. In the flow chart of figure 2 which provides a slope correction, one opted for application of the window (block 17 of figure 2) after said correction (block 16) and the diagram of figure 11 shows the signal following application of the window. This window is a band-pass filter which, as is shown on figure 11, causes the disappearance of the noise part which is located on either side of the useful part.
Advantageously, the window shows a bell form of the type of window called Hanning (see the work entitled "The Measurement of Power Spectra" by R.B. Blackman and J.W.
Tuckey, Dover Publications New York 1958). Other windows may be employed and for this one may consult the work of M. Kunt entitled "Traitement numérique des signaux"
published by Dunod, 1981.
e) The signal following the window is rendered usable if one returns now into the time domain which is effected by block 18 on figure 2, which applies to the signal a transformation inverse to that which has been described at step a). By analogy to what has been said hereinabove, this operation may be very generally effected by an inverse transformation in z or, in the case taken here by way of example, by an inverse Fourier transformation. The signal following transformation is shown on figure 12. If one desi~nates by TF the inverse Fourier transformation, by G1(f) the logarithmic derivative of the squared modulus following correction of the slope and by H(f) the frequency window, the signal of figure 12 may be expressed by the relation:
A(t) = TF 1lG1(f) ~(f~]
2 ~ f~
f) The following step of the method consists in determining the envelope of the inverse transformation brought about in step e). This operation is effected by block 19 of figure 2 and appears as shown on figure 13.
As is shown by this latter figure, the envelope exhibits a series of peaks 30, 31, 32, etc. of variable amplitude and spaced out in time. In order to arrive at this, one will advantageously employ the Hilbert transformation as described in the work "Probability, Random Variables and Stochastic Process" of A. Papoulis published by McGraw-Hill, 1984. If one designates by TH the Hilbert transformation, the envelope of figure 13 is expressed by:
{[A(t)]2 + [TH A(t)]2~1/2 the value of A(t) being that obtained hereinabove.
g) From the envelope obtained in the preceding step, one determines the time spread ~t separating the time origin from peak 30 exhibiting the greatest amplitude, which peak is due to the contribution from the periodic signal of the useful part defined hereinabove, in recalling that it is the periodic signal cf the useful part in which are found the echoes sought after. Here the time origin is a relative origin which corresponds to the first echo. ~he interval ~t as measured corresponds then very well to the interval separating the two echoes. This origin differs from that represented on figure S. Such comes about from the fact that in employing the squared modulus calculated in step b) one loses the information which contains the time separating the time origin shown on figure 5 from the first echo.
In the special case, since the slope of the useful signal has been suppressed, the peak having the greatest 1 5 2 ~
amplitude due to the contribution of the periodic signal is indeed the first encountered and referenced 30, the distance e between the two faces 3 and 4 of the wall 2 of tube 1 being expressed by the relation e = c at/2.
On the graph of figure 13, one measures ~t = 0.2 ~s. If c = 1.5 mm/~s as indicated hereinabove, one finds e =
1.5- 0.2/2 = 0.15 mm which corresponds well to the thickness of the tube which one knew from other means of measurement was equal to 0.15 mm.
Additional Considerations The plastic tube employed hereinabove in order to explain the method of the invention gives well defined signals particularly due to the rigidity of the tube subjected to measurement. In numerous cases, the object for which the thickness is to be measured has no rigidity;
furthermore, it is generally inaccessible and does not allow measurement by mechanical means. This is the case for instance of an artery in vivo for which one wishes to measure in a non-invasive manner the thickness of the wall.
Figure 14 presents the signal envelope obtained in the next to last step of the method according to the invention from a swine femoral artery in vitro. In this case, there has been effected a slope correction as indicated hereinabove in step c). The envelope shows several peaks and without any contestation it is the peak of the greatest amplitude noted at 0.56 ~s which must be taken into consideration.
The diagram of figure 15 is an envelope obtained from the same artery but where the slope correction has not been effected. One readily finds at 0.56 ~s a peak corresponding to that sought after but this peak is preceded by another peak 40 of greater amplitude, such peak originat-2 ~ 2 ~
ing from the contribution of the linear part (or from the slope) of the useful signal which h~re is preponderant.
This teaching militates in favour of the systematic employment of the slope correction in all cases where the object to be measured shows little stability and/or is not accessible.
It has already been mentioned that the method may be employed for the measurement of the thickness of an artery. The same method may be profitably employed for measuring the interior diameter of the artery and generally for the measurement of any object exhibiting a thin layer for which one is concerned with measuring the thickness.
Claims (12)
1. A method for measuring the distance e separating two faces of an object which method includes the emission of an ultrasonic wave propagating at the velocity c toward the object and the reception by a sensor of the echoes generated by reflection of said wave from the faces of said object, said sensor delivering a composite electrical signal g(t) showing an amplitude varying as a function of time, said composite electrical signal being processed by the following sequence of steps:
a) transforming the composite electrical signal g(t) coming from the sensor into a frequency spectrum G(f) the amplitude of which is a function of the frequency, said spectrum exhibiting a real part Re[G(f)] and an imagin-ary part Im[G(f)], b) calculating the order n modulus of the frequency spectrum ¦G(f)¦n = {Re [G(f)]2 + Im[G(f)]2}n/2, c) calculating the logarithmic derivative (¦G(f)¦n)'/¦G(f)¦n of the modulus calculated in step b) in order to separate the existing echoes from the composite electrical signal, said logarithmic derivative exhibiting a useful part centered substantially on the frequency of the ultrasonic wave and a part having noise located on either side of the useful part, said useful part comprising a periodic signal onto which a substantially linear signal is superposed, d) applying a frequency window to the logarithmic derivative in order to isolate the useful part thereof and to eliminate the part having noise, e) applying a transformation inverse to that of step a) at least to the periodic signal of said useful part, f) determining the envelope of the inverse transform-ation, said envelope exhibiting a series of variable ampli-tude peaks spaced out over time, g) determining the temporal spread .DELTA.t separating the time origin from the peak showing the greatest amplitude due to the contribution of the periodic signal from the useful part, the sought for distance e separating the two faces of the object being expressed by the relation e =
c?.DELTA.t/2.
a) transforming the composite electrical signal g(t) coming from the sensor into a frequency spectrum G(f) the amplitude of which is a function of the frequency, said spectrum exhibiting a real part Re[G(f)] and an imagin-ary part Im[G(f)], b) calculating the order n modulus of the frequency spectrum ¦G(f)¦n = {Re [G(f)]2 + Im[G(f)]2}n/2, c) calculating the logarithmic derivative (¦G(f)¦n)'/¦G(f)¦n of the modulus calculated in step b) in order to separate the existing echoes from the composite electrical signal, said logarithmic derivative exhibiting a useful part centered substantially on the frequency of the ultrasonic wave and a part having noise located on either side of the useful part, said useful part comprising a periodic signal onto which a substantially linear signal is superposed, d) applying a frequency window to the logarithmic derivative in order to isolate the useful part thereof and to eliminate the part having noise, e) applying a transformation inverse to that of step a) at least to the periodic signal of said useful part, f) determining the envelope of the inverse transform-ation, said envelope exhibiting a series of variable ampli-tude peaks spaced out over time, g) determining the temporal spread .DELTA.t separating the time origin from the peak showing the greatest amplitude due to the contribution of the periodic signal from the useful part, the sought for distance e separating the two faces of the object being expressed by the relation e =
c?.DELTA.t/2.
2. A method as set forth in claim 1 wherein immediat-ely before or after step d) the substantially linear signal appearing at step c) is suppressed, the peak having the greatest amplitude as determined in step g) being then the greatest amplitude peak of the envelope obtained in step f).
3. A method as set forth in claim 1 or in claim 2 wherein step a) is preceded by a digitalization of the composite electrical signal, said digitalization being brought about by means of an analog-digital converter and the signal thus obtained being stored in a memory.
4. A method as set forth in claim 1 or in claim 2 wherein the transformation effected in step a) is a Four-ier transformation and the inverse transformation effected in step e) is an inverse Fourier transformation.
5. A method as set forth in claim 1 or in claim 2 wherein the envelope appearing in step f) employs a Hil-bert transformation of the signal appearing in step e).
6. A method as set forth in claim 1 or in claim 2 wherein the frequency window appearing in step d) is a Hanning window.
7. A method as set forth in claim 2 wherein the useful part of the signal appearing at step c) is visualized and one proceeds to suppress the substantially linear signal by means of successive trials until a graphically satisfac-tory signal is obtained.
8. A method as set forth in claim 2 wherein suppres-sion of the substantially linear signal appearing after step c) or step d) is obtained by the following sequence of steps:
h) calculating from the expression of the order n modulus obtained in step b), for which one places n =
2, the order 1 moment m1 and the order 2 moment m2, i) determining the line guided by f = m1 having as slope p = k?m2, said line being expressed by the relation p(f - m1) where f is equal to the frequency and k is a constant, and j) subtracting the line obtained in step i) from the logarithmic derivative obtained in step c).
h) calculating from the expression of the order n modulus obtained in step b), for which one places n =
2, the order 1 moment m1 and the order 2 moment m2, i) determining the line guided by f = m1 having as slope p = k?m2, said line being expressed by the relation p(f - m1) where f is equal to the frequency and k is a constant, and j) subtracting the line obtained in step i) from the logarithmic derivative obtained in step c).
9. A method as set forth in claim 2 wherein the substantially linear signal is suppressed by using the characteristics of said signal known a priori by a reference measurement.
10. Using the method as set forth in claim 1 for the non-invasive measurement of the wall thickness of an artery.
11. Using the method as set forth in claim 1 for the non-invasive measurement of the interior diameter of an artery.
12. Using the method as set forth in claim 1 for the measurement of the thickness of thin layers.
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR8909926A FR2650071B1 (en) | 1989-07-20 | 1989-07-20 | PROCESS FOR PROCESSING AN ELECTRICAL SIGNAL |
| FR8909926 | 1989-07-20 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CA2021509A1 true CA2021509A1 (en) | 1991-01-21 |
Family
ID=9384066
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CA002021509A Abandoned CA2021509A1 (en) | 1989-07-20 | 1990-07-19 | Method of processing an electrical signal |
Country Status (7)
| Country | Link |
|---|---|
| US (1) | US5197019A (en) |
| EP (1) | EP0409054B1 (en) |
| JP (1) | JPH0357907A (en) |
| CA (1) | CA2021509A1 (en) |
| DE (1) | DE69000081D1 (en) |
| FR (1) | FR2650071B1 (en) |
| NO (1) | NO903643L (en) |
Families Citing this family (47)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CH685420A5 (en) * | 1991-09-27 | 1995-07-14 | Asulab Sa | A method of ultrasonic measurement of the position of at least one movable wall and apparatus for implementing this method. |
| US5349859A (en) * | 1991-11-15 | 1994-09-27 | Scientific Engineering Instruments, Inc. | Method and apparatus for measuring acoustic wave velocity using impulse response |
| US5404754A (en) * | 1993-02-10 | 1995-04-11 | Shell Oil Company | Ultrasonic detection of high temperature hydrogen attack |
| US5723791A (en) * | 1993-09-28 | 1998-03-03 | Defelsko Corporation | High resolution ultrasonic coating thickness gauge |
| US5491668A (en) * | 1994-05-13 | 1996-02-13 | Western Atlas International, Inc. | Method for determining the thickness of a casing in a wellbore by signal processing pulse-echo data from an acoustic pulse-echo imaging tool |
| DE19612925C1 (en) * | 1996-04-01 | 1997-05-22 | Fraunhofer Ges Forschung | Coated workpiece material characteristics testing system |
| US5608165A (en) * | 1996-05-06 | 1997-03-04 | Ford Motor Company | Ultrasonic thickness gauge for multilayer plastic fuel tanks |
| ES2352340T3 (en) * | 1997-07-05 | 2011-02-17 | Hudson-Sharp Machine Company | APPLIANCE FOR THE APPLICATION OF RESELLABLE CLOSURES ON A FILM BAND. |
| JP2889568B1 (en) * | 1998-05-18 | 1999-05-10 | 正男 伊藤 | Vascular thickness measurement device and arteriosclerosis diagnosis device |
| EP1055133A1 (en) | 1998-12-15 | 2000-11-29 | Koninklijke Philips Electronics N.V. | Ultrasound process and apparatus for the determination of the location of a parietal surface in a tissue and of the absolute radius of an artery |
| IT1306874B1 (en) * | 1999-06-25 | 2001-10-11 | Consiglio Nazionale Ricerche | NON-INVASIVE ACOUSTIC METHOD FOR THE DETECTION AND MAPPING OF THE POST-PAINTINGS IN FRESCO PAINTINGS AND WALL PAINTS IN GENERAL, AND |
| US6397680B1 (en) * | 2000-07-24 | 2002-06-04 | National Research Council Of Canada | Ultrasonic spectroscopy apparatus for determining thickness and other properties of multilayer structures |
| US6494097B1 (en) | 2000-09-05 | 2002-12-17 | Elias Edmond Shihadeh | Method and apparatus for measuring thickness of a layer in a multi-layered object |
| JP3643562B2 (en) * | 2002-02-08 | 2005-04-27 | コーリンメディカルテクノロジー株式会社 | Pulse wave velocity measuring device |
| US20040127797A1 (en) * | 2002-06-07 | 2004-07-01 | Bill Barnard | System and method for measuring bladder wall thickness and presenting a bladder virtual image |
| US7520857B2 (en) * | 2002-06-07 | 2009-04-21 | Verathon Inc. | 3D ultrasound-based instrument for non-invasive measurement of amniotic fluid volume |
| US20090062644A1 (en) * | 2002-06-07 | 2009-03-05 | Mcmorrow Gerald | System and method for ultrasound harmonic imaging |
| US7819806B2 (en) * | 2002-06-07 | 2010-10-26 | Verathon Inc. | System and method to identify and measure organ wall boundaries |
| US8221322B2 (en) * | 2002-06-07 | 2012-07-17 | Verathon Inc. | Systems and methods to improve clarity in ultrasound images |
| US8221321B2 (en) | 2002-06-07 | 2012-07-17 | Verathon Inc. | Systems and methods for quantification and classification of fluids in human cavities in ultrasound images |
| US20100036252A1 (en) * | 2002-06-07 | 2010-02-11 | Vikram Chalana | Ultrasound system and method for measuring bladder wall thickness and mass |
| US20060025689A1 (en) * | 2002-06-07 | 2006-02-02 | Vikram Chalana | System and method to measure cardiac ejection fraction |
| US20090112089A1 (en) * | 2007-10-27 | 2009-04-30 | Bill Barnard | System and method for measuring bladder wall thickness and presenting a bladder virtual image |
| GB2391625A (en) * | 2002-08-09 | 2004-02-11 | Diagnostic Ultrasound Europ B | Instantaneous ultrasonic echo measurement of bladder urine volume with a limited number of ultrasound beams |
| US20080262356A1 (en) * | 2002-06-07 | 2008-10-23 | Vikram Chalana | Systems and methods for ultrasound imaging using an inertial reference unit |
| AU2003260288A1 (en) * | 2002-09-11 | 2004-04-30 | Beta Lasermike | Method and apparatus for determination of layer thickness in a multi-layer structure |
| US20040100282A1 (en) * | 2002-09-11 | 2004-05-27 | Henrik Christensen | Method and apparatus for determination of layer thickness in a multi-layer structure |
| US6745632B1 (en) * | 2003-06-03 | 2004-06-08 | Joseph Ernest Dryer | Method for measuring ultrasonic transit times |
| US20050038343A1 (en) * | 2003-07-10 | 2005-02-17 | Alfred E. Mann Institute For Biomedical Research At The University Of Southern California | Apparatus and method for locating a bifurcation in an artery |
| KR20050014051A (en) * | 2003-07-29 | 2005-02-07 | 안희태 | Distance Measuring Method and Device by Frequency Separation with Ultrasonic |
| CA2598335A1 (en) * | 2004-02-17 | 2005-09-01 | Verathon Inc. | System and method for measuring bladder wall thickness and mass |
| FR2886534A1 (en) * | 2005-06-03 | 2006-12-08 | Theraclion Soc Par Actions Sim | IMAGING AND PROCESSING HEAD OF LIVING ORGANS AND METHOD OF MANUFACTURING |
| FR2886551B1 (en) * | 2005-06-03 | 2007-09-07 | Theraclion Soc Par Actions Sim | METHOD FOR DETERMINING DISTANCE AND PROCESSING APPARATUS USING THE SAME |
| FR2886533B1 (en) | 2005-06-03 | 2007-09-14 | Theraclion Soc Par Actions Sim | IMAGING AND PROCESSING HEAD OF LIVING ORGANS AND METHOD OF MANUFACTURING |
| US7434988B1 (en) | 2006-05-17 | 2008-10-14 | Enertechnix, Inc. | Low pressure acoustic pyrometer signal generator |
| US7954378B2 (en) * | 2007-02-07 | 2011-06-07 | Pepperl + Fuchs Gmbh | Ultrasonic sensor and method for determining a separation of an object from an ultrasonic sensor |
| US8167803B2 (en) * | 2007-05-16 | 2012-05-01 | Verathon Inc. | System and method for bladder detection using harmonic imaging |
| NO330292B1 (en) * | 2007-09-12 | 2011-03-21 | Det Norske Veritas As | Acoustic thickness paints using gas as a coupling medium |
| JP2009103459A (en) * | 2007-10-19 | 2009-05-14 | Ricoh Elemex Corp | Ultrasonic plate thickness measuring device |
| JP5157610B2 (en) * | 2008-04-17 | 2013-03-06 | Ｊｆｅスチール株式会社 | Method for measuring the thickness of thin steel pipes |
| US8225998B2 (en) * | 2008-07-11 | 2012-07-24 | Es&S Innovations Llc | Secure ballot box |
| JP5658151B2 (en) * | 2008-08-07 | 2015-01-21 | ベラソン インコーポレイテッドＶｅｒａｔｈｏｎ Ｉｎｃ． | Apparatus, system and method for measuring the diameter of an abdominal aortic aneurysm |
| CN107635472A (en) | 2015-06-19 | 2018-01-26 | 神经系统分析公司 | Transcranial doppler detector |
| CA2996221C (en) * | 2015-08-24 | 2023-09-19 | Hifi Engineering Inc. | Method and system for determining the distance to an acoustically reflective object in a conduit |
| EP3399918A4 (en) | 2016-01-05 | 2019-08-21 | Neural Analytics, Inc. | Systems and methods for determining clinical indications |
| US11589836B2 (en) | 2016-01-05 | 2023-02-28 | Novasignal Corp. | Systems and methods for detecting neurological conditions |
| US10617388B2 (en) | 2016-01-05 | 2020-04-14 | Neural Analytics, Inc. | Integrated probe structure |
Family Cites Families (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2308909A1 (en) * | 1975-04-22 | 1976-11-19 | Commissariat Energie Atomique | Precise measurement of metallic tube dimensions - involves use of obstacles in ultrasonic path to determine coupling medium propagation velocity |
| BE840456A (en) * | 1975-04-22 | 1976-10-07 | DEVICE FOR PRECISE MEASUREMENT OF THE DIMENSIONS OF AN OBJECT BY ULTRA-SOUND | |
| US4063549A (en) * | 1975-12-22 | 1977-12-20 | Technicon Instruments Corporation | Ultrasonic method and apparatus for imaging and characterization of bodies |
| US4534220A (en) * | 1984-01-27 | 1985-08-13 | Kleeper John R | Apparatus for measuring the characteristics of an ultrasonic wave medium |
| US4648276A (en) * | 1984-01-27 | 1987-03-10 | Klepper John R | Apparatus for measuring the characteristics of an ultrasonic wave medium |
| US4817432A (en) * | 1985-09-27 | 1989-04-04 | Design Team Partners | Digital ultrasonic instrument for ophthalmic use |
| US4799177A (en) * | 1985-12-31 | 1989-01-17 | The Boeing Company | Ultrasonic instrumentation for examination of variable-thickness objects |
| US4855911A (en) * | 1987-11-16 | 1989-08-08 | Massachusetts Institute Of Technology | Ultrasonic tissue characterization |
| US4926870A (en) * | 1988-08-30 | 1990-05-22 | Osteo-Technology, Inc. | Method and apparatus for ultrasonic analysis of bone strength in vivo |
| FR2637400B1 (en) * | 1988-09-30 | 1990-11-09 | Labo Electronique Physique | DEVICE FOR IMPROVED PROCESSING OF AN ECHOGRAPHIC SIGNAL |
| FR2644913A1 (en) * | 1989-03-24 | 1990-09-28 | Labo Electronique Physique | ULTRASONIC ULTRASONIC IMAGING DEVICE USING AN IMPROVED ADAPTIVE FILTER |
| US5038615A (en) * | 1990-05-11 | 1991-08-13 | General Motors Corporation | Ultrasonic multilayer paint thickness measurement |
-
1989
- 1989-07-20 FR FR8909926A patent/FR2650071B1/en not_active Expired - Fee Related
-
1990
- 1990-07-11 DE DE9090113192T patent/DE69000081D1/en not_active Expired - Fee Related
- 1990-07-11 EP EP90113192A patent/EP0409054B1/en not_active Expired - Lifetime
- 1990-07-19 US US07/554,420 patent/US5197019A/en not_active Expired - Fee Related
- 1990-07-19 CA CA002021509A patent/CA2021509A1/en not_active Abandoned
- 1990-07-20 JP JP2190927A patent/JPH0357907A/en active Pending
- 1990-08-17 NO NO90903643A patent/NO903643L/en unknown
Also Published As
| Publication number | Publication date |
|---|---|
| FR2650071B1 (en) | 1991-09-27 |
| US5197019A (en) | 1993-03-23 |
| EP0409054A1 (en) | 1991-01-23 |
| NO903643L (en) | 1991-02-21 |
| DE69000081D1 (en) | 1992-06-04 |
| NO903643D0 (en) | 1990-08-17 |
| FR2650071A1 (en) | 1991-01-25 |
| EP0409054B1 (en) | 1992-04-29 |
| JPH0357907A (en) | 1991-03-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5197019A (en) | Method of measuring distance using ultrasonic waves | |
| US4803990A (en) | Examining moving objects by ultrasound echograpy | |
| He et al. | Acoustic dispersion and attenuation measurement using both transmitted and reflected pulses | |
| EP0864857A1 (en) | Method and apparatus for terahertz tomographic imaging | |
| US4928698A (en) | Device for measuring the speed of moving organs and blood flows by correlation | |
| JPS6377437A (en) | Movable matter examination apparatus by ultrasonic echography | |
| JPS59174152A (en) | Measuring system of ultrasonic medium characteristics | |
| US4844082A (en) | Ultrasonic examination apparatus | |
| NL8600444A (en) | DEVICE FOR ULTRA-SOUND DETECTION. | |
| JPH08193986A (en) | Nondestructive test device | |
| KR20040066795A (en) | Apparatus and method for indicating mechanical stiffness properties of body tissue technical field | |
| JP4795925B2 (en) | Ultrasonic thickness measurement method and apparatus | |
| US5109856A (en) | Device for measuring the speed of moving organs and blood flows by ultrasonic echography | |
| US5062430A (en) | Device for measuring the speed of blood flows by ultrasonic echography at an increased measuring speed | |
| JPH0715457B2 (en) | Ultrasonic echography inspection method and device | |
| WO2018208886A1 (en) | Ultrasonic inspection methods and systems | |
| Hanshaw et al. | An H∞ deconvolution filter and its application to ultrasonic nondestructive evaluation of materials | |
| KR20030086325A (en) | Ultra sonic method for measuring the thickness of a faintly reflective partial layer by addition of echo periods | |
| JP2818615B2 (en) | Ultrasonic measuring device | |
| JPH0676998B2 (en) | Precise measurement of ultrasonic round trip time by pulse reflection method | |
| EP0134470A2 (en) | Method and means for determining ultrasonic wave attenuation in tissue using phase locked loop | |
| JPH0542140A (en) | Probing device | |
| Renken et al. | Accuracy of resonant ultrasound spectroscopy: An experimental study | |
| Nickel et al. | Reduction of transducer reverberations by transducer movement-determination of the optimal transducer displacement | |
| JPH11108648A (en) | Measuring method of thickness and distance by using ultrasonic wave |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FZDE | Discontinued |