CA2527982A1 - Method and system for analyzing cardiovascular sounds - Google Patents
Method and system for analyzing cardiovascular sounds Download PDFInfo
- Publication number
- CA2527982A1 CA2527982A1 CA002527982A CA2527982A CA2527982A1 CA 2527982 A1 CA2527982 A1 CA 2527982A1 CA 002527982 A CA002527982 A CA 002527982A CA 2527982 A CA2527982 A CA 2527982A CA 2527982 A1 CA2527982 A1 CA 2527982A1
- Authority
- CA
- Canada
- Prior art keywords
- signals
- acoustic energy
- average acoustic
- determined
- cardiovascular
- 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
- A61B7/00—Instruments for auscultation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B7/00—Instruments for auscultation
- A61B7/02—Stethoscopes
- A61B7/04—Electric stethoscopes
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2562/00—Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
- A61B2562/02—Details of sensors specially adapted for in-vivo measurements
- A61B2562/0204—Acoustic sensors
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2562/00—Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
- A61B2562/04—Arrangements of multiple sensors of the same type
- A61B2562/046—Arrangements of multiple sensors of the same type in a matrix array
Abstract
A method and system for analyzing sounds originating in at least a portion of an individual~s cardiovascular system. N transducers, where N is an integer, are fixed on a surface of the individual over the thorax. The ith transducer is fixed at a location xi and generates an initial signal P(xi,t) indicative of pressure waves at the location xi; for i=1 to N. the signals P(xi,t)are processes so as to generate filtered signals in which at least one component of the signals P(xi,t)not arising from cardiovascular sounds has been removed.
The filtered signals may be used for generating an image of the at least portion of the cardiovascular system.
The filtered signals may be used for generating an image of the at least portion of the cardiovascular system.
Description
METHOD AND SYSTEM FOR ANALYZING
CARDIOVASCULAR 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
s ~ 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 or heartbeat in order to detect adventitious (i.e.
abnormal or unexpected) lung or heartsounds. 'The identification and classification of adventitious lung or heart sounds often provides important information about to pulmonary or cardiac 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 is 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.
(Chest, 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 2o used to determine the location in the lungs of the source of a sound detected in the-recordings.
_2_ U.S. Patent 5,285,788 discloses an ultrasound tissue imaging system having an acoustic transducer, and imaging means for producing an image of tissue.
The system also includes Doppler imaging means to produce a scanned acoustic image of moving tissue that is displayed superimposed on the ultrasound image.
s 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.
The present invention provides, in one of its embodiments, a system and 1o method for recording and analyzing cardiovascular sounds produced in the cardiovascular system. The system includes a , plurality of N transducers (microphones) configured to be attached to an essentially planar region R of the individual's back or chest over the individual's thorax. Positions in the region R are indicated by two-dimensional position vectors x=(xi,x2) in a two-dimensional is coordinate system defined in the planar region R. The ith transducer, for i=I to N;
is fixed at a position xZ in the region R and generates a signal, denoted herein by P(x~,t), 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 20 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(xi, t) are processed by signal processing circuitry. In accordance with the invention, the signals are filtered so as to remove one or more 2s components of the signals not arising from cardiovascular sounds (e.g:
respiratory tract signals). Cardiovascular sounds are typically in the range of 6 to 45 Hz, while respiratory tract sounds are typically in the range of 100 to 400 Hz. Thus, respiratory .sounds may be removed from the signals by filtering the signals, for example, with a band pass filter passing between 15 to 45 Hz.
CARDIOVASCULAR 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
s ~ 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 or heartbeat in order to detect adventitious (i.e.
abnormal or unexpected) lung or heartsounds. 'The identification and classification of adventitious lung or heart sounds often provides important information about to pulmonary or cardiac 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 is 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.
(Chest, 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 2o used to determine the location in the lungs of the source of a sound detected in the-recordings.
_2_ U.S. Patent 5,285,788 discloses an ultrasound tissue imaging system having an acoustic transducer, and imaging means for producing an image of tissue.
The system also includes Doppler imaging means to produce a scanned acoustic image of moving tissue that is displayed superimposed on the ultrasound image.
s 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.
The present invention provides, in one of its embodiments, a system and 1o method for recording and analyzing cardiovascular sounds produced in the cardiovascular system. The system includes a , plurality of N transducers (microphones) configured to be attached to an essentially planar region R of the individual's back or chest over the individual's thorax. Positions in the region R are indicated by two-dimensional position vectors x=(xi,x2) in a two-dimensional is coordinate system defined in the planar region R. The ith transducer, for i=I to N;
is fixed at a position xZ in the region R and generates a signal, denoted herein by P(x~,t), 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 20 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(xi, t) are processed by signal processing circuitry. In accordance with the invention, the signals are filtered so as to remove one or more 2s components of the signals not arising from cardiovascular sounds (e.g:
respiratory tract signals). Cardiovascular sounds are typically in the range of 6 to 45 Hz, while respiratory tract sounds are typically in the range of 100 to 400 Hz. Thus, respiratory .sounds may be removed from the signals by filtering the signals, for example, with a band pass filter passing between 15 to 45 Hz.
The N filtered signals (also indicated herein by P(x1, t)) may be processed in order to diagnose the state of the individual's cardiovascular system. This may be via an automatic differential diagnosis in which the results of the processing are compared to functions or parameters previously. stored in a database that are s known to be indicative of various disorders in the body region.
The filtered signals may also be processed to generate an image of the individual's cardiovascular system. The results of this processing are displayed on a display device, for example using a gray level scale, as demonstrated in the examples below. In the image, anatomic features of the heart such as the atria, to ventricles, septal walls, can be observed. The image may be visually or automatically analyzed for the detection of a disorder in the cardiovascular system similar to the analysis of images obtained by other imaging methods such as X-ray (scintigraphy) or ultrasound imaging (echocardiography).
A region or regions of the heart or cardiovascular system in a displayed Is 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 cardiovascular system. This includes murmurs and other hemodynamic 20 ' irregularities, cardiac effusion, narrowing of blood vessel, and other space containing lesions in the cardiovascular system, etc.
Additionally, a time interval can be divided into a plurality of sub intervals, and each subinterval processed separately. An image of the cardiovascular system for each of these subintervals may then be determined and displayed 2s sequentially on the display device. This generates a movie showing dynamic changes occurring in the cardiovascular system over the time interval. This allows viewing of the systoles and diastoles of the different parts of the heart during the heartbeat.
The filtered signals may also be processed to generate an image of the individual's cardiovascular system. The results of this processing are displayed on a display device, for example using a gray level scale, as demonstrated in the examples below. In the image, anatomic features of the heart such as the atria, to ventricles, septal walls, can be observed. The image may be visually or automatically analyzed for the detection of a disorder in the cardiovascular system similar to the analysis of images obtained by other imaging methods such as X-ray (scintigraphy) or ultrasound imaging (echocardiography).
A region or regions of the heart or cardiovascular system in a displayed Is 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 cardiovascular system. This includes murmurs and other hemodynamic 20 ' irregularities, cardiac effusion, narrowing of blood vessel, and other space containing lesions in the cardiovascular system, etc.
Additionally, a time interval can be divided into a plurality of sub intervals, and each subinterval processed separately. An image of the cardiovascular system for each of these subintervals may then be determined and displayed 2s sequentially on the display device. This generates a movie showing dynamic changes occurring in the cardiovascular system over the time interval. This allows viewing of the systoles and diastoles of the different parts of the heart during the heartbeat.
In a preferred embodiment, the processing involves determining from the N
signals an average acoustic energy arising from cardiovascular sounds, denoted herein by P (x, t~, tz) , at at least one position x in the region R over a time interval from t1 to t2. The erm "acoustic energy" at a location is used herein to refer to a s parameter indicative of or 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 t1 to t2 is obtained at a position of one of the microphones using the algebraic expression tZ
P(xt, ti, tz) = f P z (xr, t)dt ( 1 ) where xl is the position of the microphone.
In a more preferred embodiment, the processing involves obtaining an average acoustic energy P ( x1 , t; , t2 ) over a time interval from tl to t2 at a plurality of positions x1 of the microphones, for example using Equation (1), and then 1 s calculating P (x, t~, tz) at other locations x by interpolation of the P ( x; , t, , tz ) using any known interpolation method.
In a most preferred embodiment, the interpolation is performed to obtain an average acoustic energy P(x, t1 , tz ) at a position x = (x' , x2 ) in the surface R using the algebraic expression:
_ N _ P (x, t~> tz) = ~ P(x; , t~, tz)g(x, x, , ~) (2) where g(x, x; , o-) is a kernel satisfying °zg a6 ~ (3) N
g(x, x; , o-) is approximately equal to 1 (4) and where x; _ (x; , x? ) is the position of the ith microphone and o- is a selectable 2s parameter.
signals an average acoustic energy arising from cardiovascular sounds, denoted herein by P (x, t~, tz) , at at least one position x in the region R over a time interval from t1 to t2. The erm "acoustic energy" at a location is used herein to refer to a s parameter indicative of or 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 t1 to t2 is obtained at a position of one of the microphones using the algebraic expression tZ
P(xt, ti, tz) = f P z (xr, t)dt ( 1 ) where xl is the position of the microphone.
In a more preferred embodiment, the processing involves obtaining an average acoustic energy P ( x1 , t; , t2 ) over a time interval from tl to t2 at a plurality of positions x1 of the microphones, for example using Equation (1), and then 1 s calculating P (x, t~, tz) at other locations x by interpolation of the P ( x; , t, , tz ) using any known interpolation method.
In a most preferred embodiment, the interpolation is performed to obtain an average acoustic energy P(x, t1 , tz ) at a position x = (x' , x2 ) in the surface R using the algebraic expression:
_ N _ P (x, t~> tz) = ~ P(x; , t~, tz)g(x, x, , ~) (2) where g(x, x; , o-) is a kernel satisfying °zg a6 ~ (3) N
g(x, x; , o-) is approximately equal to 1 (4) and where x; _ (x; , x? ) is the position of the ith microphone and o- is a selectable 2s parameter.
For example, the kernel 1 _ I 2 2 _ 2 r')2 g~x~ xa ~ ~) - Exh - ~x x' ~) ' gxp - (.x x~ 'V (5) 26 ~ . 2~
may be used.
The system may optionally contain a display device for displaying the s 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 cardiovascular system. In the image anatomic features of the heart such as the atria, ventricles, septal walls, can be observed. The image may be analyzed for to the detection of a disorder in the cardiovascular system similar to the analysis of images obtained by other imaging methods such as X-ray (scintigraphy) or ultrasound imaging (echocardiography).
A region or regions of the heap or cardiovascular system in a displayed image that are suspected of including a pathological condition, may de identified in 1 s 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 cardiovascular system. This includes murmurs and other hemodynamic irregularities, cardiac effusion, narrowing of blood vessel, and other space 2o containing lesions in the cardiovascular system, 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 2s 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 or back and an average acoustic energy P determined in accordance with the invention for a plurality of sub intervals over one or more heartbeats. 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 heart over the heartbeat. This allows viewing of the systoles and diastoles of the different parts of s the hear during the heartbeat.
The signals P(x; , t) may also be subjected to additional analysis to detect abnormal heart sounds.
The present invention also provides a program storage device readable by machine, tangibly embodying a program of instructions executable by the machine to to perform method steps for determining for at least one time interval, an average acoustic energy function P wising from cardiovascular sounds using an algorithm involving at least one signal P(xi, t) indicative of pressure waves at a location x; on a body surface.
The present invention still further provides a computer program product Is comprising a computer useable medium having computer readable program code embodied therein analyzing sounds in at least a portion ' of an individual's cardiovascular system, 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 arising from the poi lion 20 of the cardiovascular system, P being determined in algorithm involving at least one signal P(xi, t) indicative of pressure waves at a location x; on a body surface.
The invention thus provides a system for analyzing sounds originating in at least a portion of an individual's cardiovascular system comprising:
(One) N transducers, where N is an integer, each transducer configured 2s to b~e fixed on a surface of the individual over the thorax, the ith transducer being fixed at a location x1 and generating an initial signal P(xa, t) indicative of pressure waves at the location x;; for i=1 to N; and (Two) a processor configured to receive the signals P(x; , t) and to filter the signals P(xt , t) so as to generate filtered signals in which at least _7-one component of the signals P(xt, t) not arising from cardiovascular sounds has been removed.
The invention thus further provides a system for analyzing sounds originating in at least a portion of an individual's cardiovascular system comprising:
(~.) N transducers, where N is an integer, 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 an initial signal P(x~, t) indicative of pressure waves at the location xt; for i=1 to N; and (Two) a processor configured to receive the signals P(x; , t) and to 1o generate therefrom an image of the at least portion of the cardiovascular system.
The invention thus further provides a method for analyzing sounds originating in at least a portion of an individual's cardiovascular system comprising:
(One) fixing on a surface of the individual over the thorax, N
Is transducers, where N is an integer, the ith transducer being fxed at a location x1 and generating an initial signal P(xr, t) indicative of pressure waves at the location x;; for i=1 to N; and (Two) processing the signals P(x; , t) so as to generate filtered signals in which at least one component of the signals P(x~, t) not arising from 2o cardiovascular sounds has been removed.
The invention thus further provides a method for analyzing sounds originating in at least a portion of an individual's cardiovascular system comprising:
(One) fixing on a surface of the individual over the thorax, N
transducers, where N is an integer, the ith transducer being fixed at a 2s. location xZ and generating an initial signal P(x~,t)indicative of pressure waves at the location x;; for i=1' to N; and (Two) processing the signals P(x; , t) so as to generate therefrom an image of the at least portion of the cardiovascular system.
The invention thus further provides a program storage device readable by machine, tangibly embodying a program of instructions executable by the machine to perform method steps for analyzing sounds originating in at least a portion of an individual's cardiovascular system, comprising: processing N initial signals s P(x; , t) , where N is an integer, the initial signals being indicative of pressure waves at a location x;; for i=1 to N, so as to generate filtered signals in which at least one component of the signals P(xr, t) not arising from cardiovascular sounds has been removed.
The invention thus further provides a computer program product comprising a computer useable medium having computer readable program code embodied therein for analyzing sounds originating in at least a portion of an individual's cardiovascular system, comprising: processing N initial signals P(x~ , t) , where N is an integer, the initial signals being indicative of pressure waves at a location x1; for i=1 to N, so as to generate filtered signals in which at least one component of the 1 s signals P(x~, t) not arising from cardiovascular sounds has been removed.
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 2o example only, with reference to the accompanying drawings, in which:
Fig. 1 shows a system for obtaining an analyzing cardiovascular sounds in accordance with one embodiment of the invention;
Fig. 2 shows a flow chart for carrying out a method of analyzing cardiovascular sounds in accordance with one embodiment of the invention;
2s Fig. 3 shows the locations of transducers on an individual's back for analyzing cardiovascular sounds;
Fig. 4 shows successive frames from a movie of the heart. of a healthy individual over one heart beat;
Fig. 5 shows successive frames from a movie of the heart and lungs of an individual over one respiratory cycle; and DETAILED DESCRIPTION OF THE INVENTION
Fig. 1 shows a system generally indicated by 100 for analyzing body sounds s 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, axe applied to a planar region of the chest or back skin of individual 110. The transducers 105 may be any type of sound transducer, such as a microphone or a Doppler shift detector. The transducers 105 may be applied to the to 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(x~, t) 125, represent the pressure wave at the location xt of the ith transducer (i=1 is 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 20 filtering so that only components in the signal within the range of cardiovascular sounds are analyzed. The signal may be divided into frequency bands, and each band analyzed separately.
An input device such as a computer keyboard I40 or mouse 145 is used to input relevant information relating to the examination such as personal details of 2s the individual 110. The input device 140 may also be used to input values of the times t1 and t2 . Alternatively, the times t1 and t2 may be determined automatically in a respiratory phase analysis of the signals P(xi , t) performed by the processor 135. The processor 135 determines an average acoustic energy P (x, t~, ta) over the - 1~ -time interval from t1 to tz at least one location x in the region R in a calculation involving at least one of the signals P(x~ , t) .
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 s 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 to invention in accordance with one embodiment. In step 200 the signals P(x~,t) are obtained from N transducers placed at predetermined locations xi for i from 1 to N
in a region R on the body surface. In step 205 values of ti and tz 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, t~, tz) is determined at least 15 one location x in the region R over the time interval t1 to tz _ 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 2o 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 The system and method of the invention were used to analyze cardiovascular sounds in an individual.
Example 1 Fig. 3 shows recording of signals over one heartbeat in an individual. A
two-dimensional coordinate system was defined on the individual's back. As shown in Fig. 3, 48 transducers were placed on the individual's back over the thorax, at the locations indicated by the circles 300. The curves 305 show the presumed contours s of the lungs, and the curve 306 shows the presumed contour of the heart. As can be seen, the transducers were arranged in a regular orthogonal lattice with spacing between the transducers in the horizontal and vertical directions of 2.5 cm.
The signals P(x~,t) from each transducer were then recorded over one heartbeat.
Each signal was filtered using a 6-45 Hz band pass filter, in order to remove respiratory to tract sounds. The heartbeat was divided into intervals of 0.1 sec duration, and for each interval, P(x, t~, tz) was obtained using Equations ( 1 ) and (2) above with the kernel g of Equation (5) with 6=36 pixels. Fig. 4 shows the images obtained by representing the obtained functions P(x; t~, tz) by gray level shading. The images may be displayed on the display device 150 in rapid succession so as to produce a is movie of the heart over a heartbeat. The movie can be analyzed to determine the values of basic parameters of heart function, such as left ventricular end diastolic (LVED) volume, left ventricular end systolic CLUES) volume, right ventricular end diastolic (RUED) volume, right ventricular end systolic (RVES), volume, left atrium end diastolic (LAED) diameter, right atrium end diastolic (LAES) diameter, 2o wall thickness of the inter-ventricular septum (systolic and diastolic), and parameters derivable from these parameters such as left ventricle stroke volume, left ventricular cardiac output, ejection fraction, left ventricular fractional shortening, inter-ventricular septal thickening. The, movie can also be analyzed in order to detect heart defects such as valve dysfunction and cardiac arrhythmia. .
2s Example 2 The signals P(x~, t) were obtained from each transducer as described in Example l, and~were then recorded over one respiratory cycle which includes about , S heartbeats. Each signal was divided into two sub-signals Pl(xl, t) and P2(xu t) of different frequency bands. The sub-signal Pl(xvt) was obtained by filtering the signal using a 6-40 Hz band pass filter. The sub-signal P2(x;, t) was obtained by filtering the signal using a 100-150 band pass filter. The sub-signal Pl(x;, t) consists primarily of heart sounds, while the sub-signal P2(x;, t) consists primarily of lung sounds. The PI(x;,t) sub-signal was analyzed by the method of the invention, and s the sub-signal P2(x;,t) was analyzed as disclosed in Applicant's co-pending U.S
Patent Application 10/33 x,742 filed on January 9, 2003. The signal P2(xv t) was divided into intervals of 0.25 sec duration, and the signal PI(x~, t) was divided into intervals of 0.1 sec duration. For each interval, functions . P (x, t~, t2) and P(x, t~, tz) were obtained from Pl(x;, t) and P2(xv t), respectively, using Equations (1 ) to and (2) above with the kernel g of Equation (5) with 6=36 pixels. The two functions are preferably displayed simultaneously on ~ a display device by intensity shading, using a different color for each function. Fig. 5 shows the images obtained by representing the obtained functions P(x, t~, tz) and P(x, t~, t2) simultaneously by gray level shading. The images may be displayed on the display device 150 in rapid ~s succession so as to produce a movie of the heart over a heartbeat. The movie can be analyzed to determine the values of parameters of heart function, such as cardiac output and blood ejection fraction. The movie can also be analyzed in order to detect hear defects such as valve dysfunction and cardiac arrhythmia.
may be used.
The system may optionally contain a display device for displaying the s 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 cardiovascular system. In the image anatomic features of the heart such as the atria, ventricles, septal walls, can be observed. The image may be analyzed for to the detection of a disorder in the cardiovascular system similar to the analysis of images obtained by other imaging methods such as X-ray (scintigraphy) or ultrasound imaging (echocardiography).
A region or regions of the heap or cardiovascular system in a displayed image that are suspected of including a pathological condition, may de identified in 1 s 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 cardiovascular system. This includes murmurs and other hemodynamic irregularities, cardiac effusion, narrowing of blood vessel, and other space 2o containing lesions in the cardiovascular system, 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 2s 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 or back and an average acoustic energy P determined in accordance with the invention for a plurality of sub intervals over one or more heartbeats. 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 heart over the heartbeat. This allows viewing of the systoles and diastoles of the different parts of s the hear during the heartbeat.
The signals P(x; , t) may also be subjected to additional analysis to detect abnormal heart sounds.
The present invention also provides a program storage device readable by machine, tangibly embodying a program of instructions executable by the machine to to perform method steps for determining for at least one time interval, an average acoustic energy function P wising from cardiovascular sounds using an algorithm involving at least one signal P(xi, t) indicative of pressure waves at a location x; on a body surface.
The present invention still further provides a computer program product Is comprising a computer useable medium having computer readable program code embodied therein analyzing sounds in at least a portion ' of an individual's cardiovascular system, 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 arising from the poi lion 20 of the cardiovascular system, P being determined in algorithm involving at least one signal P(xi, t) indicative of pressure waves at a location x; on a body surface.
The invention thus provides a system for analyzing sounds originating in at least a portion of an individual's cardiovascular system comprising:
(One) N transducers, where N is an integer, each transducer configured 2s to b~e fixed on a surface of the individual over the thorax, the ith transducer being fixed at a location x1 and generating an initial signal P(xa, t) indicative of pressure waves at the location x;; for i=1 to N; and (Two) a processor configured to receive the signals P(x; , t) and to filter the signals P(xt , t) so as to generate filtered signals in which at least _7-one component of the signals P(xt, t) not arising from cardiovascular sounds has been removed.
The invention thus further provides a system for analyzing sounds originating in at least a portion of an individual's cardiovascular system comprising:
(~.) N transducers, where N is an integer, 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 an initial signal P(x~, t) indicative of pressure waves at the location xt; for i=1 to N; and (Two) a processor configured to receive the signals P(x; , t) and to 1o generate therefrom an image of the at least portion of the cardiovascular system.
The invention thus further provides a method for analyzing sounds originating in at least a portion of an individual's cardiovascular system comprising:
(One) fixing on a surface of the individual over the thorax, N
Is transducers, where N is an integer, the ith transducer being fxed at a location x1 and generating an initial signal P(xr, t) indicative of pressure waves at the location x;; for i=1 to N; and (Two) processing the signals P(x; , t) so as to generate filtered signals in which at least one component of the signals P(x~, t) not arising from 2o cardiovascular sounds has been removed.
The invention thus further provides a method for analyzing sounds originating in at least a portion of an individual's cardiovascular system comprising:
(One) fixing on a surface of the individual over the thorax, N
transducers, where N is an integer, the ith transducer being fixed at a 2s. location xZ and generating an initial signal P(x~,t)indicative of pressure waves at the location x;; for i=1' to N; and (Two) processing the signals P(x; , t) so as to generate therefrom an image of the at least portion of the cardiovascular system.
The invention thus further provides a program storage device readable by machine, tangibly embodying a program of instructions executable by the machine to perform method steps for analyzing sounds originating in at least a portion of an individual's cardiovascular system, comprising: processing N initial signals s P(x; , t) , where N is an integer, the initial signals being indicative of pressure waves at a location x;; for i=1 to N, so as to generate filtered signals in which at least one component of the signals P(xr, t) not arising from cardiovascular sounds has been removed.
The invention thus further provides a computer program product comprising a computer useable medium having computer readable program code embodied therein for analyzing sounds originating in at least a portion of an individual's cardiovascular system, comprising: processing N initial signals P(x~ , t) , where N is an integer, the initial signals being indicative of pressure waves at a location x1; for i=1 to N, so as to generate filtered signals in which at least one component of the 1 s signals P(x~, t) not arising from cardiovascular sounds has been removed.
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 2o example only, with reference to the accompanying drawings, in which:
Fig. 1 shows a system for obtaining an analyzing cardiovascular sounds in accordance with one embodiment of the invention;
Fig. 2 shows a flow chart for carrying out a method of analyzing cardiovascular sounds in accordance with one embodiment of the invention;
2s Fig. 3 shows the locations of transducers on an individual's back for analyzing cardiovascular sounds;
Fig. 4 shows successive frames from a movie of the heart. of a healthy individual over one heart beat;
Fig. 5 shows successive frames from a movie of the heart and lungs of an individual over one respiratory cycle; and DETAILED DESCRIPTION OF THE INVENTION
Fig. 1 shows a system generally indicated by 100 for analyzing body sounds s 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, axe applied to a planar region of the chest or back skin of individual 110. The transducers 105 may be any type of sound transducer, such as a microphone or a Doppler shift detector. The transducers 105 may be applied to the to 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(x~, t) 125, represent the pressure wave at the location xt of the ith transducer (i=1 is 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 20 filtering so that only components in the signal within the range of cardiovascular sounds are analyzed. The signal may be divided into frequency bands, and each band analyzed separately.
An input device such as a computer keyboard I40 or mouse 145 is used to input relevant information relating to the examination such as personal details of 2s the individual 110. The input device 140 may also be used to input values of the times t1 and t2 . Alternatively, the times t1 and t2 may be determined automatically in a respiratory phase analysis of the signals P(xi , t) performed by the processor 135. The processor 135 determines an average acoustic energy P (x, t~, ta) over the - 1~ -time interval from t1 to tz at least one location x in the region R in a calculation involving at least one of the signals P(x~ , t) .
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 s 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 to invention in accordance with one embodiment. In step 200 the signals P(x~,t) are obtained from N transducers placed at predetermined locations xi for i from 1 to N
in a region R on the body surface. In step 205 values of ti and tz 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, t~, tz) is determined at least 15 one location x in the region R over the time interval t1 to tz _ 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 2o 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 The system and method of the invention were used to analyze cardiovascular sounds in an individual.
Example 1 Fig. 3 shows recording of signals over one heartbeat in an individual. A
two-dimensional coordinate system was defined on the individual's back. As shown in Fig. 3, 48 transducers were placed on the individual's back over the thorax, at the locations indicated by the circles 300. The curves 305 show the presumed contours s of the lungs, and the curve 306 shows the presumed contour of the heart. As can be seen, the transducers were arranged in a regular orthogonal lattice with spacing between the transducers in the horizontal and vertical directions of 2.5 cm.
The signals P(x~,t) from each transducer were then recorded over one heartbeat.
Each signal was filtered using a 6-45 Hz band pass filter, in order to remove respiratory to tract sounds. The heartbeat was divided into intervals of 0.1 sec duration, and for each interval, P(x, t~, tz) was obtained using Equations ( 1 ) and (2) above with the kernel g of Equation (5) with 6=36 pixels. Fig. 4 shows the images obtained by representing the obtained functions P(x; t~, tz) by gray level shading. The images may be displayed on the display device 150 in rapid succession so as to produce a is movie of the heart over a heartbeat. The movie can be analyzed to determine the values of basic parameters of heart function, such as left ventricular end diastolic (LVED) volume, left ventricular end systolic CLUES) volume, right ventricular end diastolic (RUED) volume, right ventricular end systolic (RVES), volume, left atrium end diastolic (LAED) diameter, right atrium end diastolic (LAES) diameter, 2o wall thickness of the inter-ventricular septum (systolic and diastolic), and parameters derivable from these parameters such as left ventricle stroke volume, left ventricular cardiac output, ejection fraction, left ventricular fractional shortening, inter-ventricular septal thickening. The, movie can also be analyzed in order to detect heart defects such as valve dysfunction and cardiac arrhythmia. .
2s Example 2 The signals P(x~, t) were obtained from each transducer as described in Example l, and~were then recorded over one respiratory cycle which includes about , S heartbeats. Each signal was divided into two sub-signals Pl(xl, t) and P2(xu t) of different frequency bands. The sub-signal Pl(xvt) was obtained by filtering the signal using a 6-40 Hz band pass filter. The sub-signal P2(x;, t) was obtained by filtering the signal using a 100-150 band pass filter. The sub-signal Pl(x;, t) consists primarily of heart sounds, while the sub-signal P2(x;, t) consists primarily of lung sounds. The PI(x;,t) sub-signal was analyzed by the method of the invention, and s the sub-signal P2(x;,t) was analyzed as disclosed in Applicant's co-pending U.S
Patent Application 10/33 x,742 filed on January 9, 2003. The signal P2(xv t) was divided into intervals of 0.25 sec duration, and the signal PI(x~, t) was divided into intervals of 0.1 sec duration. For each interval, functions . P (x, t~, t2) and P(x, t~, tz) were obtained from Pl(x;, t) and P2(xv t), respectively, using Equations (1 ) to and (2) above with the kernel g of Equation (5) with 6=36 pixels. The two functions are preferably displayed simultaneously on ~ a display device by intensity shading, using a different color for each function. Fig. 5 shows the images obtained by representing the obtained functions P(x, t~, tz) and P(x, t~, t2) simultaneously by gray level shading. The images may be displayed on the display device 150 in rapid ~s succession so as to produce a movie of the heart over a heartbeat. The movie can be analyzed to determine the values of parameters of heart function, such as cardiac output and blood ejection fraction. The movie can also be analyzed in order to detect hear defects such as valve dysfunction and cardiac arrhythmia.
Claims (48)
1. A system for analyzing sounds originating in at least a portion of an individual's cardiovascular system comprising:
(One) N transducers, where N is an integer, each transducer configured to be fixed on a surface of the individual over the thorax, the ith transducer being fixed at a location x i and generating an initial signal P(x i,t) indicative of pressure waves at the location x i; for i=1 to N; and (Two) a processor configured to receive the signals P(x i,t) and to filter the signals P(x i,t) so as to generate filtered signals in which at least one component of the signals P(x i,t) not arising from cardiovascular sounds has been removed.
(One) N transducers, where N is an integer, each transducer configured to be fixed on a surface of the individual over the thorax, the ith transducer being fixed at a location x i and generating an initial signal P(x i,t) indicative of pressure waves at the location x i; for i=1 to N; and (Two) a processor configured to receive the signals P(x i,t) and to filter the signals P(x i,t) so as to generate filtered signals in which at least one component of the signals P(x i,t) not arising from cardiovascular sounds has been removed.
2. The method according to Claim 1 wherein respiratory tract sounds have been removed.
3. The method according to Claim 1 wherein the initial signals are filtered by a band pass filter passing between 15 to 45 Hz.
4. The system according to claim 1 wherein the processor is further configured to generate an image of the at least portion of the cardiovascular system from at least one of the filtered signals.
5. The system according to Claim 1 wherein the processor is configured to generate an image of the at least portion of the cardiovascular system at a plurality of times or over a plurality of successive time intervals, and to display the images successively on a display device.
6. The system according to Claim 1 further comprising a two-dimensional display device.
7. The system according to Claim 4 wherein the processor is further configured to display an image of the at least portion of the cardiovascular system on a display device.
8. The system according to Claim 1 wherein the processor is further configured to determine an average acoustic energy ~(x,t1,t2) arising from the cardiovascular system at at least one position x over a time interval from a first time t1 to a second time t2, ~ being determined in an algorithm involving at least one of the processed signals.
9. The system according to Claim 8 wherein the processor is further configured to compare the average acoustic energy ~ to one or more predetermined functions ~ and determine a function ~0 from among the functions ~ most similar to ~.
10. The system according to Claim 9 wherein the processor is further configured to make a diagnosis based upon the determined function.
11. The system according to Claim 8 wherein the average acoustic energy ~ over a time interval from t1 to t2 is determined at a location x i of a transducer using the algebraic expression:
12. The system according to Claim 8 wherein the function ~ is determined at one or more locations x in an algorithm comprising:
(One) determining an average acoustic energy ~(x i,t1,t2) over a time interval from t1 to t2 at a plurality of locations x i of transducers; and (Two) determining an average acoustic energy ~(x,t1,t2) at at least one location x by interpolation of the determined ~(x i,t1,t2).
(One) determining an average acoustic energy ~(x i,t1,t2) over a time interval from t1 to t2 at a plurality of locations x i of transducers; and (Two) determining an average acoustic energy ~(x,t1,t2) at at least one location x by interpolation of the determined ~(x i,t1,t2).
13. The system according to Claim 12 wherein an average acoustic energy ~(x i.t1,t2) is determined over a time interval from t1 to t2 at a plurality of locations x i of transducers using the algebraic expression:
14. The system according to Claim 12 wherein an average acoustic energy is determined at at least one location x by interpolation of the determined ~(x i, t1, t2) using the algebraic expression:
where g(x,xi,.sigma.) is a kernel satisfying is approximately equal to 1. (4)
where g(x,xi,.sigma.) is a kernel satisfying is approximately equal to 1. (4)
15. The system according to Claim 14 wherein g(x,v i.sigma.) is the kernel
16. The system according to Claim 8 wherein the processor is further configured to display the average acoustic energy ~(x,t1,t2) on a display device.
17. The system according to Claim 1 wherein the processor is configured to determine an average acoustic energy over a plurality of successive time intervals, each average acoustic energy being determined using an algorithm involving at least one of the signals P(x i,t).
18. The system according to Claim 17 wherein the processor is configured to sequentially display on a display device a representation of each determined average acoustic energy.
19. The system according to Claim 1 wherein the processor is configured to determine an average acoustic energy over a plurality of successive time intervals, each average acoustic energy being determined using an algorithm involving at least one of the signals P(x i,t).
20. The system according to Claim 1 wherein the processor is configured to:
(One) for each of one or more frequency bands, (Twenty-seven) subject the signals (P,x i,t) to band pass filtering in the frequency band; and (ab) determine an average acoustic energy function for the frequency band based upon at least one of the filtered signals.
(One) for each of one or more frequency bands, (Twenty-seven) subject the signals (P,x i,t) to band pass filtering in the frequency band; and (ab) determine an average acoustic energy function for the frequency band based upon at least one of the filtered signals.
21. The system according to Claim 20 wherein the processor is configured to display one or more of the average acoustic energy functions determined for a frequency band on a display device.
22. A system for analyzing sounds originating in at least a portion of an individual's cardiovascular system comprising:
(a)10 N transducers, where N is an integer, each transducer configured to be fixed on a surface of the individual over the thorax, the ith transducer being fixed at a location x i and generating an initial signal P(x i,t)indicative of pressure waves at the location x i; for i=1 to N;
and (Two) a processor configured to receive the signals P(x i,t) and to generate therefrom an image of the at least portion of the cardiovascular system.
(a)10 N transducers, where N is an integer, each transducer configured to be fixed on a surface of the individual over the thorax, the ith transducer being fixed at a location x i and generating an initial signal P(x i,t)indicative of pressure waves at the location x i; for i=1 to N;
and (Two) a processor configured to receive the signals P(x i,t) and to generate therefrom an image of the at least portion of the cardiovascular system.
23. A method for analyzing sounds originating in at least a portion of an individual's cardiovascular system comprising:
(One) fixing on a surface of the individual over the thorax, N transducers, where N is an integer, the ith transducer being fixed at a location x i and generating an initial signal P(x i, t) indicative of pressure waves at the location x i; for i=1 to N; and (Two) processing the signals P(x i,t) so as to generate filtered signals in which at least one component of the signals P(x i,t) not arising from cardiovascular sounds has been removed.
(One) fixing on a surface of the individual over the thorax, N transducers, where N is an integer, the ith transducer being fixed at a location x i and generating an initial signal P(x i, t) indicative of pressure waves at the location x i; for i=1 to N; and (Two) processing the signals P(x i,t) so as to generate filtered signals in which at least one component of the signals P(x i,t) not arising from cardiovascular sounds has been removed.
24. The method according to Claim 23 wherein respiratory tract sounds are filtered from the initial signals.
25. The method according to Claim 23 wherein the initial signals are filtered by a band pass filter passing between 15 to 45 Hz.
26. The method according to claim 23 further comprising generating an image of the at least portion of the cardiovascular method from at least one of the filtered signals.
27. The method according to Claim 23 further comprising generating an image of the at least portion of the cardiovascular method at a plurality of times or over a plurality of successive time intervals, and displaying the images successively on a display device.
28. The method according to Claim 26 further comprising displaying an image of the at least portion of the cardiovascular method on a display device.
29. The method according to Claim 23 further comprising determining an average acoustic energy ~(x, t1, t2) arising from the cardiovascular method at at least one position x over a time interval from a first time t, to a second time t2, ~ being determined in an algorithm involving at least one of the processed signals.
30. The method according to Claim 29 further comprising comparing the average acoustic energy ~ to one or more predetermined functions F and determining a function ~o from among the functions ~ most similar to ~ .
31. The method according to Claim 30 wherein the processor is further configured to make a diagnosis based upon the determined function.
32. The method according to Claim 29 wherein the average acoustic energy ~ over a time interval from t1 to t2 is determined at a location x i of a transducer using the algebraic expression:
33. The method according to Claim 29 wherein the function ~ is determined at one or more locations x in an algorithm comprising:
(One) determining an average acoustic energy ~ (x t , t1, t2) over a time interval from t1 to t2 at a plurality of locations x i of transducers; and (Two) determining an average acoustic energy ~(x, t1, t2) at at least one Location x by interpolation of the determined ~(x t, t1, t2).
(One) determining an average acoustic energy ~ (x t , t1, t2) over a time interval from t1 to t2 at a plurality of locations x i of transducers; and (Two) determining an average acoustic energy ~(x, t1, t2) at at least one Location x by interpolation of the determined ~(x t, t1, t2).
34. The method according to Claim 33 wherein an average acoustic energy ~ (x t , t1, t2) is determined over a time interval from t1 to t2 at a plurality of locations x i of transducers using the algebraic expression:
35. The method according to Claim 33 wherein an average acoustic energy is determined at at least one location x by interpolation of the determined ~(x t, t1, t2) using the algebraic expression:
where g(x, x t, .sigma.) is a kernel satisfying
where g(x, x t, .sigma.) is a kernel satisfying
36. The method according to Claim 35 wherein g(x,v,.sigma.) is the kernel g(x, x t,.sigma.)=
37. The method according to Claim 29 further comprising displaying the average acoustic energy ~(x,t1,t2) on a display device.
38. The method according to Claim 23 further comprising determining an average acoustic energy over a plurality of successive time intervals, each average acoustic energy being determined using an algorithm involving at least one of the signals P(x i,t).
39. The method according to Claim 38 further comprising sequentially displaying on a display device a representation of each determined average acoustic energy.
40. The method according to Claim 23 further comprising determining an average acoustic energy over a plurality of successive time intervals, each average acoustic energy being determined using an algorithm involving at least one of the signals P(x i,t).
41. The method according to Claim 23 further comprising:
(c) for each of one or more frequency bands, (ca) subjecting the signals (P,x:p t) to band pass filtering in the frequency band; and (cb) determining an average acoustic energy function for the frequency to band based upon at least one of the filtered signals.
(c) for each of one or more frequency bands, (ca) subjecting the signals (P,x:p t) to band pass filtering in the frequency band; and (cb) determining an average acoustic energy function for the frequency to band based upon at least one of the filtered signals.
42. The method according to Claim 41 further comprising displaying one or more of the average acoustic energy functions determined for a frequency band on a display device.
43. A method for analyzing sounds originating in at least a portion of an is individual's cardiovascular system comprising:
(One) fixing on a surface of the individual over the thorax, N
transducers, where N is an integer, the ith transducer being fixed at a location x i and generating an initial signal P(x i, t) indicative of pressure waves at the location x i; for i=1 to N; and 2(0Two) processing the signals P(x i,t) so as to generate therefrom an image of the at least portion of the cardiovascular system.
(One) fixing on a surface of the individual over the thorax, N
transducers, where N is an integer, the ith transducer being fixed at a location x i and generating an initial signal P(x i, t) indicative of pressure waves at the location x i; for i=1 to N; and 2(0Two) processing the signals P(x i,t) so as to generate therefrom an image of the at least portion of the cardiovascular system.
44. A processor configured to receive the signals P(x i,t) and to generate therefrom an image of the at least portion of the cardiovascular system.
45. Use of the method of Claim 23 for diagnosing a cardiovascular disorder.
46. The use according to Claim 44 wherein the disorder is selected from the group comprising at least cardiac arrhythmia's and heart valve disorders.
47. A program storage device readable by machine, tangibly embodying a program of instructions executable by the machine to perform method steps for analyzing sounds originating in at least a portion of an individual's cardiovascular system, comprising: processing N initial signals P(x i, t), where N is an integer, the initial signals being indicative of pressure waves at a location x i; for i=1 to N, so as to generate filtered signals in which at least one component of the signals P(x i,t) not arising from cardiovascular sounds has been removed.
48. A computer program product comprising a computer useable medium having computer readable program code embodied therein for analyzing sounds originating in at least a portion of an individual's cardiovascular system, comprising: processing N initial signals P(x i,t), where N is an integer, the initial signals being indicative of pressure waves at a location x i; for i=1 to N, so as to generate filtered signals in which at least one component of the signals P(x i,t) not arising from cardiovascular sounds has been removed.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US47459503P | 2003-06-02 | 2003-06-02 | |
| US60/474,595 | 2003-06-02 | ||
| PCT/IL2004/000462 WO2004105612A1 (en) | 2003-06-02 | 2004-06-01 | Method and system for analyzing cardiovascular sounds |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CA2527982A1 true CA2527982A1 (en) | 2004-12-09 |
Family
ID=33490722
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CA002527982A Abandoned CA2527982A1 (en) | 2003-06-02 | 2004-06-01 | Method and system for analyzing cardiovascular sounds |
Country Status (11)
| Country | Link |
|---|---|
| US (1) | US7517319B2 (en) |
| EP (1) | EP1628576A1 (en) |
| JP (1) | JP2006526452A (en) |
| KR (1) | KR20060028682A (en) |
| CN (1) | CN100493458C (en) |
| AU (1) | AU2004243189A1 (en) |
| BR (1) | BRPI0410945A (en) |
| CA (1) | CA2527982A1 (en) |
| RU (1) | RU2346653C2 (en) |
| TW (1) | TW200509858A (en) |
| WO (1) | WO2004105612A1 (en) |
Families Citing this family (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8033996B2 (en) * | 2005-07-26 | 2011-10-11 | Adidas Ag | Computer interfaces including physiologically guided avatars |
| US7998091B2 (en) * | 2005-11-23 | 2011-08-16 | 3M Innovative Properties Company | Weighted bioacoustic sensor and method of using same |
| US8024974B2 (en) | 2005-11-23 | 2011-09-27 | 3M Innovative Properties Company | Cantilevered bioacoustic sensor and method using same |
| US8920343B2 (en) | 2006-03-23 | 2014-12-30 | Michael Edward Sabatino | Apparatus for acquiring and processing of physiological auditory signals |
| US20080139893A1 (en) * | 2006-12-08 | 2008-06-12 | Warren Lee | Apparatus And System For Sensing and Analyzing Body Sounds |
| US20100087746A1 (en) * | 2006-12-11 | 2010-04-08 | Naira Radzievsky | Method and system for analyzing body sounds |
| EP2249706B1 (en) * | 2008-03-04 | 2013-05-29 | Koninklijke Philips Electronics N.V. | Non invasive analysis of body sounds |
| US8218846B2 (en) | 2008-05-15 | 2012-07-10 | Superdimension, Ltd. | Automatic pathway and waypoint generation and navigation method |
| US8218847B2 (en) | 2008-06-06 | 2012-07-10 | Superdimension, Ltd. | Hybrid registration method |
| JP5603157B2 (en) * | 2010-07-15 | 2014-10-08 | 独立行政法人科学技術振興機構 | Emotion control device |
| US8491488B1 (en) * | 2010-10-01 | 2013-07-23 | Blaufuss Medical Multimedia Laboratories, LLC | Method and system for identifying cardiopulmonary findings by using a heart and lung sounds builder |
| EP3691532A2 (en) | 2017-10-04 | 2020-08-12 | Ausculsciences, Inc. | Auscultatory sound-or-vibration sensor |
| US11284827B2 (en) | 2017-10-21 | 2022-03-29 | Ausculsciences, Inc. | Medical decision support system |
| US10828009B2 (en) | 2017-12-20 | 2020-11-10 | International Business Machines Corporation | Monitoring body sounds and detecting health conditions |
| TWI775711B (en) * | 2022-02-22 | 2022-08-21 | 昌泰科醫股份有限公司 | Method and system of detecting specific physiological syndrome related to hyperactivity of liver-fire/heart-fire based on hemodynamic analysis and stringy pulse |
Family Cites Families (41)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3604852A (en) * | 1970-03-02 | 1971-09-14 | Howard Wise | Apparatus for the visual aesthetic display of sound |
| US3878832A (en) * | 1973-05-14 | 1975-04-22 | Palo Alto Medical Research Fou | Method and apparatus for detecting and quantifying cardiovascular murmurs and the like |
| US4220160A (en) * | 1978-07-05 | 1980-09-02 | Clinical Systems Associates, Inc. | Method and apparatus for discrimination and detection of heart sounds |
| US4289142A (en) * | 1978-11-24 | 1981-09-15 | Kearns Kenneth L | Physiological occurrence, such as apnea, monitor and X-ray triggering device |
| US4387722A (en) * | 1978-11-24 | 1983-06-14 | Kearns Kenneth L | Respiration monitor and x-ray triggering apparatus |
| SU993917A1 (en) | 1981-06-08 | 1983-02-07 | Научно-Исследовательский Институт Сельского Хозяйства Нечерноземной Зоны Усср | Device for auscalting animal udder |
| US4777961A (en) * | 1985-10-15 | 1988-10-18 | Bruce Saltzman | High sensitivity stethoscopic system and method |
| US4833625A (en) * | 1986-07-09 | 1989-05-23 | University Of Arizona | Image viewing station for picture archiving and communications systems (PACS) |
| US5213108A (en) * | 1988-02-04 | 1993-05-25 | Blood Line Technology, Inc. | Visual display stethoscope |
| US5010889A (en) * | 1988-02-04 | 1991-04-30 | Bloodline Technology | Intelligent stethoscope |
| US5218969A (en) * | 1988-02-04 | 1993-06-15 | Blood Line Technology, Inc. | Intelligent stethoscope |
| GB8825611D0 (en) | 1988-11-02 | 1988-12-07 | Stoneman S A T | Respiratory condition diagnosis & apparatus therefor |
| US5058600A (en) * | 1990-01-08 | 1991-10-22 | Center For Innovative Technology | Graphical readout of laryngotracheal spectra and airway monitor |
| US5016642A (en) * | 1990-04-11 | 1991-05-21 | Hewlett-Packard Company | Slow motion cardiac imaging |
| JPH0482538A (en) | 1990-07-25 | 1992-03-16 | Hitachi Ltd | Inhaling sound diagnosing apparatus |
| US5309922A (en) * | 1992-09-21 | 1994-05-10 | Center For Innovative Technology | Respiratory sound analyzer for use in high noise environments |
| US5285788A (en) * | 1992-10-16 | 1994-02-15 | Acuson Corporation | Ultrasonic tissue imaging method and apparatus with doppler velocity and acceleration processing |
| US5526442A (en) * | 1993-10-04 | 1996-06-11 | Hitachi Medical Corporation | X-ray radiography method and system |
| US5492125A (en) * | 1995-02-10 | 1996-02-20 | University Of Washington | Ultrasound signal processing apparatus |
| US6390977B1 (en) * | 1995-06-07 | 2002-05-21 | Alliance Pharmaceutical Corp. | System and methods for measuring oxygenation parameters |
| US5957866A (en) * | 1995-07-03 | 1999-09-28 | University Technology Corporation | Apparatus and methods for analyzing body sounds |
| US5774558A (en) * | 1995-10-30 | 1998-06-30 | Rsq, Llc | Sound imager |
| IL117146A0 (en) | 1996-02-15 | 1996-06-18 | Gull Medical Software Systems | Diagnosis of lung condition |
| DE69730224D1 (en) * | 1996-04-25 | 2004-09-16 | Per Bjorgaas | METHOD AND INSTRUMENT FOR EXAMINING THE HEART AND ARTERIES |
| US6168568B1 (en) * | 1996-10-04 | 2001-01-02 | Karmel Medical Acoustic Technologies Ltd. | Phonopneumograph system |
| US5844997A (en) * | 1996-10-10 | 1998-12-01 | Murphy, Jr.; Raymond L. H. | Method and apparatus for locating the origin of intrathoracic sounds |
| RU2127075C1 (en) * | 1996-12-11 | 1999-03-10 | Корженевский Александр Владимирович | Method for producing tomographic image of body and electrical-impedance tomographic scanner |
| US6140565A (en) * | 1998-06-08 | 2000-10-31 | Yamaha Corporation | Method of visualizing music system by combination of scenery picture and player icons |
| IL124900A0 (en) * | 1998-06-14 | 1999-01-26 | Tapuz Med Tech Ltd | Apron for performing ecg tests and additional examinations |
| US6135960A (en) * | 1998-08-31 | 2000-10-24 | Holmberg; Linda Jean | High-resolution, three-dimensional whole body ultrasound imaging system |
| US6139505A (en) * | 1998-10-14 | 2000-10-31 | Murphy; Raymond L. H. | Method and apparatus for displaying lung sounds and performing diagnosis based on lung sound analysis |
| US6396931B1 (en) * | 1999-03-08 | 2002-05-28 | Cicero H. Malilay | Electronic stethoscope with diagnostic capability |
| US7343195B2 (en) * | 1999-05-18 | 2008-03-11 | Mediguide Ltd. | Method and apparatus for real time quantitative three-dimensional image reconstruction of a moving organ and intra-body navigation |
| WO2001022883A1 (en) * | 1999-09-29 | 2001-04-05 | Siemens Corporate Research, Inc. | Multi-modal cardiac diagnostic decision support system and method |
| US6572560B1 (en) * | 1999-09-29 | 2003-06-03 | Zargis Medical Corp. | Multi-modal cardiac diagnostic decision support system and method |
| US6527729B1 (en) * | 1999-11-10 | 2003-03-04 | Pacesetter, Inc. | Method for monitoring patient using acoustic sensor |
| US6381351B1 (en) * | 1999-11-24 | 2002-04-30 | Direct Radiography Corp. | Weighted inverse topography method for digital x-ray image data processing |
| US6944330B2 (en) * | 2000-09-07 | 2005-09-13 | Siemens Corporate Research, Inc. | Interactive computer-aided diagnosis method and system for assisting diagnosis of lung nodules in digital volumetric medical images |
| KR100387201B1 (en) * | 2000-11-16 | 2003-06-12 | 이병훈 | Diaortic apparatus |
| GB0118728D0 (en) | 2001-07-31 | 2001-09-26 | Univ Belfast | Monitoring device |
| US20030130588A1 (en) * | 2002-01-10 | 2003-07-10 | Igal Kushnir | Method and system for analyzing respiratory tract sounds |
-
2004
- 2004-06-01 KR KR1020057023251A patent/KR20060028682A/en not_active Application Discontinuation
- 2004-06-01 US US10/857,451 patent/US7517319B2/en not_active Expired - Fee Related
- 2004-06-01 AU AU2004243189A patent/AU2004243189A1/en not_active Abandoned
- 2004-06-01 CA CA002527982A patent/CA2527982A1/en not_active Abandoned
- 2004-06-01 RU RU2005141556/14A patent/RU2346653C2/en not_active IP Right Cessation
- 2004-06-01 TW TW093115659A patent/TW200509858A/en unknown
- 2004-06-01 EP EP04735647A patent/EP1628576A1/en not_active Withdrawn
- 2004-06-01 CN CNB2004800194425A patent/CN100493458C/en not_active Expired - Fee Related
- 2004-06-01 JP JP2006508483A patent/JP2006526452A/en not_active Withdrawn
- 2004-06-01 WO PCT/IL2004/000462 patent/WO2004105612A1/en active Application Filing
- 2004-06-01 BR BRPI0410945-7A patent/BRPI0410945A/en not_active IP Right Cessation
Also Published As
| Publication number | Publication date |
|---|---|
| BRPI0410945A (en) | 2006-06-27 |
| TW200509858A (en) | 2005-03-16 |
| AU2004243189A1 (en) | 2004-12-09 |
| JP2006526452A (en) | 2006-11-24 |
| WO2004105612A1 (en) | 2004-12-09 |
| CN1819797A (en) | 2006-08-16 |
| RU2005141556A (en) | 2006-06-10 |
| US20040267149A1 (en) | 2004-12-30 |
| US7517319B2 (en) | 2009-04-14 |
| CN100493458C (en) | 2009-06-03 |
| RU2346653C2 (en) | 2009-02-20 |
| EP1628576A1 (en) | 2006-03-01 |
| KR20060028682A (en) | 2006-03-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6887208B2 (en) | Method and system for analyzing respiratory tract sounds | |
| EP2120720B1 (en) | Method and system for analyzing body sounds | |
| US7517319B2 (en) | Method and system for analyzing cardiovascular sounds | |
| JPH021217A (en) | Apparatus and method for detecting heart disease | |
| US20070197927A1 (en) | Method for detecting cardiovascular problems using micro or nano vibrations | |
| US20080221467A1 (en) | Method and systms for regional assessment of pulmonary function | |
| US20110130674A1 (en) | Method and System for Regional Assessment of Lung Physiology | |
| De Panfilis et al. | Multi-point accelerometric detection and principal component analysis of heart sounds | |
| MXPA05013006A (en) | Method and system for analyzing cardiovascular sounds | |
| JPH0690913A (en) | Living body diagnosing device | |
| IL162901A (en) | System for analyzing and imaging respiratory tract sounds |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FZDE | Discontinued |