System for analyzing and imaging respiratory tract sounds

SYSTEM FOR ANALYZING AND IMAGING RESPIRATORY TRACT SOUNDS

FIELD OF THE INVENTION

This invention relates to medical devices and methods, and more particularly to such devices and methods for analyzing body sounds.

BACKGROUND OF THE INVENTION

Body sounds are routinely used by physicians in the diagnosis of various disorders. A physician may place a stethoscope on a person’s chest or back and monitor the patient’s breathing in order to detect adventitious (i.e. abnormal or unexpected) lung sounds. The identification and classification of adventitious lung sounds often provides important information about pulmonary abnormalities.

It is also known to fix one or more microphones onto a subject’s chest or back and to record lung sounds. U.S. Patent No. 6,139,505 discloses a system in which a plurality of microphones are placed around a patient’s chest. The recordings of the microphones during inhalation and expiration are displayed on a screen, or printed on paper. The recordings are then visually examined by a physician in order to detect a pulmonary disorder in the patent. Kompis et al. (Chesfi 120(4), 2001) disclose a system in which M microphones are placed on a patient’s chest, and lung sounds are recorded. The recordings generate M linear equations that are solved using a least-squares fit. The solution of the system is used to determine the location in the lungs of the source of a sound detected in the recordings.

WO 03/057037

SUMMARY OF THE INVENTION

In the following description and set of claims, two explicitly described, calculable, or measurable variables are considered equivalent to each other when the two variables are proportional to one another.

)

5'

The present invention provides, in one of its embodiments, a system and method for recording and analyzing respiratory tract sounds produced in the respiratory tract. The system includes a plurality of N transducers (microphones) configured to be attached to an essentially planar region R of the individual’s back 10 or chest over the individual’s thorax. Positions in the region R are indicated by two-dimensional position vectors x=(x ,x) in a two-dimensional coordinate system defined in the planar region R. The ith transducer, for i=l to N, is fixed at a position x! in the region R and generates a signal, denoted herein by P(x״i), indicative of pressure waves in the body arriving at x,.

The transducers are typically embedded in a matrix that permits to affix them easily onto the individual’s skin. Such a matrix may typically be in the form of a vest or garment for easily placing over the individual’s thorax. As may be appreciated, different matrices may be used for differently sized individuals, for different ages, sexes, etc.

The N signals P(x!., t) are processed by signal processing circuitry. In accordance with the invention, the processing involves determining from the N signals an average acoustic energy, denoted herein by P(x, /1,/2), at least one position x in the region R over a time interval from ti to t2. The term acoustic energy at a location is used herein to refer to a parameter indicative of or 25 approximating the product of the pressure and the mass propagation velocity at that location.

In one embodiment, an average acoustic energy over a time interval from t! to /2 is obtained at a position of one of the microphones using the algebraic expression

WO 03/057037 ίζ P(x,,n,h)= Jp²(x_;,r)iZr (1) ii where x! is the position of the microphone.

In a more preferred embodiment, an average acoustic energy P a time interval from t/ to t₂ is obtained at a plurality of positions x!

of the microphones, for example using Equation (1), and then calculating P(x,t1,t2) at other locations x by interpolation of the ?(χ,,ί!,^) using any known interpolation method.

In a most preferred embodiment, the interpolation is performed to obtain an average acoustic energy P(x,t_i}t₂) at a position x = (x‘,x²) in the surface R using the algebraic expression:

P(x,t1,t2) = £P(x״r1,i2)g(x,x₍,a) (2) i»l where g(x, x,, σ) is a kernel satisfying (3) OCT

N g(x, x!,σ) is approximately equal to 1 (4) and where x, = (x],xf)is the position of the ith microphone and cr is a selectable parameter.

For example, the kernel may be used.

The system may optionally contain a display device for displaying the function P. The function P may be displayed on the display, for example using a gray level scale, as demonstrated in the examples below. A two dimensional graphical representation of the function P produces an image of the body region

WO 03/057037 that may be analyzed for the detection of a disorder in the body region similar to the analysis of images obtained by other imaging methods such as X-ray or ultrasound imaging.

A region or regions in a displayed image that are suspected of including a pathological condition, may de identified in the image, and this may be in a number of ways, for example, by different colors, by different patterns, by way of a written text, and many other ways. The term “pathological condition” refers to any deviation from the normal, healthy condition of the respiratory tract. This includes infection, inflammation, tumor, pleural effusion, pneumonia, narrowing of the airways, and other space containing lesions in the respiratory tract, etc.

Additionally, a time interval can be divided into a plurality of sub intervals, and an average acoustic energy P determined over the region R for two or more of the sub intervals. An image of P for each of these sub intervals may then be determined and displayed sequentially on the display device. This generates a movie showing dynamic changes occurring in the acoustic energy in the body region, over the time interval. For example, transducers may be placed on a person’s chest and an average acoustic energy P determined in accordance with the invention for a plurality of sub intervals over a breathing cycle. An image can be obtained for each of these sub intervals and displayed sequentially so as to generate a movie showing changes in the acoustic energy of the lungs over the breathing cycle.

The signals P(x״i)may also be subjected to band pass filtering before being analyzed by the method of the invention, so that an average acoustic energy is produced for one or more frequency bands of interest. The functions may be superimposed on the display device by representing each average acoustic energy function with a different color. Since respiratory sounds known as wheezes and crackles have different characteristic frequency ranges, band pass filtering can be used to identify these respiratory sounds. A region or regions in a displayed image

WO 03/057037 of wheezes or crackles may be identified in the image, for example, by a characteristic color, pattern, by way of a written text.

The present invention thus provides a system for analyzing sounds in at least a portion of an individual’s respiratory tract comprising:

(a) a plurality of N transducers, each transducer configured to be fixed on a surface of the individual over the thorax, the ith transducer being fixed at a location x,. and generating a signal P(x>,f) indicative of pressure waves at the location x₍.; for 1=1 to N; and (b) a processor configured to receive the signals P(x_l}f) and determine an average acoustic energy P(x, ti, ti) at at least one position x over a time interval from a first time t_x to a second time t₂, ?being determined in an algorithm involving at least one of the signals PM-

The present invention further provides a method for analyzing sounds in at least a portion of an individual’s thorax, comprising:

(a) obtaining N signals ?(x,r) for i=l to N, the signal P(x1,f) being indicative of pressure waves at the location x,; on a surface of the body over the thorax;

(b) determining an average acoustic energy ?(x^^ijat at least one position x over a time interval from a first time t_x to a second time t₂, ?determined in an algorithm involving at least one of the signals.

The present invention also provides a program storage device readable by machine, tangibly embodying a program of instructions executable by the machine to perform method steps for determining for at least one time interval, an average acoustic energy function ?using an algorithm involving at least one signal P(xi,f) indicative of pressure waves at a location x, on a body surface.

WO 03/057037

PCT7IL03/00029

The present invention still further provides a computer program product comprising a computer useable medium having computer readable program code embodied therein analyzing sounds in at least a portion of an individual’s body, the computer program product comprising:

computer readable program code for causing the computer to determine, for at least one time interval, an acoustic energy function P, P being determined in algorithm involving at least one signal P(xi,t) indicative of pressure waves at a location η on a body surface.

BRIEF DESCRIPTION OF THE DRAWINGS

In order to understand the invention and to see how it may be carried out in practice, a preferred embodiment will now be described, by way of non-limiting example only, with reference to the accompanying drawings, in which:

Fig. 1 shows a system for obtaining an analyzing body sound in accordance with one embodiment of the invention;

Fig. 2 shows a flow chart for carrying out a method of obtaining analyzing body sounds in accordance with one embodiment of the invention;

Fig. 3 shows recording and analysis of signals over an inspiratory phase of a respiratory cycle; and

Fig. 4 shows recording and analysis of signals over an expiratory phase of a respiratory cycle.

Fig. 5a shows an image obtained on a healthy individual in accordance with the invention, and Fig. 5b shows a chest X-ray of the same individual;

Fig. 6 shows successive frames from a movie of the respiratory tract of a healthy individual;

Fig. 7a shows an image obtained on an individual with pleural effusion in accordance with the invention, and Fig. 7b shows a chest X-ray of the same individual;

WO 03/057037

Fig. 8 shows successive frames from a movie of the respiratory tract of an individual with pleural effusion;

Fig. 9a shows an image obtained on an individual with pneumonia in accordance with the invention, and Fig. 9b shows a chest X-ray of the same individual; and

Fig. 10 shows successive frames from a movie of the respiratory tract of an individual with pleural effusion.

DETAILED DESCRIPTION OF THE INVENTION

Fig. 1 shows a system generally indicated by 100 for analyzing body sounds in a three-dimensional region of an individual’s body in accordance with one embodiment of the invention. A plurality of N sound transducers 105, of which four are shown, are applied to a planar region of the chest or back skin of individual 110. The transducers 105 may be applied to the subject by any means known in the art, for example using an adhesive, suction, or fastening straps. Each transducer 105 produces an analog voltage signal 115 indicative of pressure waves arriving to the transducer. The analog signals 115 are digitized by a multichannel analog to digital converter 120. The digital data signals P(xi,t) 125, represent the pressure wave at the location xi of the ith transducer (i= 1 to N) at time t. The data signals 125 are input to a memory 130. Data input to the memory 130 are accessed by a processor 135 configured to process the data signals 125. The signals 125 may be denoised by filtering components having frequencies outside of the range of body sounds in the body region, for example, vibrations due to movement of the individual. Each signal 125 may also be subject to band pass filtering so that only components in the signal within a range of interest are analyzed.

An input device such as a computer keyboard 140 or mouse 145 is used to input relevant information relating to the examination such as personal details of the individual 110. The input device 140 may also be used to input values of the times and f₂. Alternatively, the times f! and A, may be determined automatically

WO 03/057037 in a respiratory phase analysis of the signals P(x_ift) performed by the processor

135. The processor 135 determines an average acoustic energy P(x,n,r2)over the time interval from r, to t₂ at least one location x in the region R in a calculation involving at least one of the signals T(x₍/).

The average acoustic energies are stored in the memory 130 and may be displayed on a display device 150 such as a CRT screen for diagnosis by a physician.

The processor 135 may also perform an automatic differential diagnosis by comparing the function P to functions stored in the memory and known to be indicative of various disorders in the body region.

Fig. 2 shows a flow chart diagram for carrying out the method of the invention in accordance with one embodiment. In step 200 the signals P(x1,f) are obtained from N transducers placed at predetermined locations x, for i from 1 to N in a region R on the body surface. In step 205 values of r, and t₂ are either input to the processor 135 using the input devices 140 or 145, or are determined by the processor. In step 210, an average acoustic energy P(x,t1,t2) is determined at least one location x in the region R over the time interval i, to t₂. In step 220 the average acoustic energy is displayed on the display 150 for at least one value of x. In step 230, it is determined whether a function P is to be determined over another time interval. If yes, the process returns to step 205. If not, the process terminates.

It will also be understood that the system according to the invention may be a suitably programmed computer. Likewise, the invention contemplates a computer program being readable by a computer for executing the method of the invention. The invention further contemplates a machine-readable memory tangibly embodying a program of instructions executable by the machine for executing the method of the invention.

Examples

WO 03/057037

The system and method of the invention were used to analyze lower respiratory tract sounds in an individual.

Fig. 3 shows recording and analysis of signals over an inspiratory phase of a respiratory cycle in an individual. A two-dimensional coordinate system was defined on the individual’s back. As shown in Fig. 3a, 48 transducers were placed on the individual’s back over the lungs at the locations indicated by the circles 300. The curves 305 show the presumed contours of the lungs. As can be seen, the transducers were arranged in a regular orthogonal lattice with a spacing between the transducers in the horizontal and vertical directions of 5 cm. The signals P(x1,t) were then recorded over one inspiratory phase of a breathing cycle (t! and t2 are the beginning and end respectively of the inspiratory phase). Each signal was filtered using a low-pass filter having a cut-off of 150Hz. The average value of each filtered function ?(»,/) over the inspiratory phase is indicated in Fig. 3a by means of gray level shading of each circle 300 with reference to the gray level scale 310. P(x,t\,t2) was obtained using Equations (1) and (2) above with the kernel g of Equation (5) with σ=36 pixels. Fig. 3b shows a 512 pixel x 512 pixel graphical representation of the function P(x, ti, t2) over the inspiratory phase also in reference to the gray level scale 310. In the graphical representation of the function P(x,t1,t2) shown in Fig. 3b, the contours of the lungs and heart are easily discernable.

Fig. 4 shows recording and analysis of signals over an expiratory phase of a respiratory cycle. As shown in Fig. 4a, 48 transducers were placed on an individual’s back at the same locations x, used in Fig. 3, as indicated by the circles 400. The curves 405 show the presumed contours of the individual’s lungs. The signals Ρ(χ!,ΐ) were then recorded over one expiratory phase of a breathing cycle (t! and r₂are the beginning and end respectively of the expiratory phase). Each signal was filtered using a low-pass filter having a cut-off of 150 Hz. The average value of each function P(x!,t) over the expiratory phase is indicated in Fig. 4a by means of gray level shading of each circle 400 with reference to the gray level scale

WO 03/057037

PCI7IL03/00029 . P(x, /1,/2) was obtained using Equations (1) and (2) above. Fig. 3b shows the function P(x, /1,/2) over the expiratory phase also in reference to the gray level scale 410. Comparison of Figs. 3b and 4b shows the volume change in acoustic energy between the inspiratory and expiratory phase of the respiratory cycle.

Movies showing changes in the lungs during a respiratory cycle were obtained using the method and system of the invention. Signals 125 were obtained and divided into time segments. Each segment was analyzed by the method of the invention and an image was generated. The images were displayed on the display device 150 in rapid succession so as to produce a movie of the respiratory tract over the respiratory cycle.

Fig. 5a shows an image of the respiratory tract of a healthy individual obtained over an entire respiratory cycle in accordance with the invention, and Fig. 5b shows a chest X-ray of the same individual. Fig. 6 shows 11 successive images obtained over successive 0.4 sec time intervals during a respiratory cycle of the individual. Each frame represents the processing of the recorded signals over a time interval of 0.4 sec. Frames 01 to 05 (obtained at times 0.0 to 1.6 sec) were obtained during the inspiratory phase of the respiratory cycle, while frames 06 to 11 (obtained at times 1.8 to 3.6 sec) were obtained during the expiratory phase. The sequence of images shown in Fig. 6 can be displayed in succession on a display device so as to create a movie of the respiratory tract over a respiratory cycle. The sequence of images shown in Fig. 6 shows complete filling and emptying of the lungs during the respiratory cycle, as would be expected in a healthy individual not having any space-filling lesions.

Fig. 7a shows an image of the respiratory tract of an individual with pleural effusion obtained over an entire respiratory cycle in accordance with the invention, and Fig. 7b shows a chest X-ray of the same individual. Fig. 8 shows 16 successive images obtained over successive 0.4 sec time intervals during a respiratory cycle of the individual. Each frame represents the processing of the recorded signals over a time interval of 0.4 sec. Frames 01 to 06 (obtained at times 0.0 to 2.0 sec) were

WO 03/057037 obtained during the inspiratory phase of the respiratory cycle, while frames 07 to 16 (obtained at times 2.4 to 4.0 sec) were obtained during the expiratoiy phase. The sequence of images shown in Fig. 8 can be displayed in succession on a display device so as to create a movie of the respiratory tract over a respiratory cycle. In the sequence of images shown in Fig. 8, lung tissue in the lower right lobe is not visualized indicating the absence of air flow in the lower right lung, as would be expected in an individual having a space-filling lesion as occurs in pleural effusion. Airflow in the upper portion of the right lung is also observed to be impaired.

Fig. 9a shows an image of the respiratory tract of an individual with pneumonia obtained over an entire respiratoiy cycle in accordance with the invention, and Fig. 9b shows a chest X-ray of the same individual. Fig. 10 shows 12 successive images obtained over successive 0.4 sec time intervals during a respiratory cycle of the individual. The sequence of images shown in Fig. 10 can be displayed in succession on a display device so as to create a movie of the respiratoiy tract over a respiratory cycle. Each frame represents the processing of the recorded signals over a time interval of 0.4 sec. Frames 01 to 06 (obtained at times 0.0 to 2.0 sec) were obtained during the inspiratory phase of the respiratory cycle, while frames 07 to 16 (obtained at times 2.4 to 4.0 sec) were obtained during the expiratoiy phase. In the sequence of images shown in Fig. 10, lung tissue in the lower .left lobe is not visualized indicating the absence of air flow in the lower left lung, as would be expected in an individual having a space-filling lesion as occurs in pneumonia. Airflow in the upper portion of the left lung is observed to be normal.