Method for detecting blockage in a fluid flow vessel
11452460 · 2022-09-27
Assignee
Inventors
- Bret Kline (Columbus, OH, US)
- Peter Bakema (Denver, NC, US)
- Young Truong (Carrboro, NC, US)
- Richard Finlayson (Greenville, NC, US)
- Orville Day (Greenville, NC, US)
Cpc classification
A61B5/0285
HUMAN NECESSITIES
A61B5/7221
HUMAN NECESSITIES
A61B5/6844
HUMAN NECESSITIES
A61B2560/0431
HUMAN NECESSITIES
A61B2562/166
HUMAN NECESSITIES
A61B5/684
HUMAN NECESSITIES
A61B5/02007
HUMAN NECESSITIES
A61B5/6843
HUMAN NECESSITIES
International classification
A61B5/0285
HUMAN NECESSITIES
A61B5/00
HUMAN NECESSITIES
Abstract
A method for measuring sound from vortices in the carotid artery comprising: a first and second quality control provisions, wherein the quality control compares detected sounds to pre-determined sounds, and upon confirmation of the quality control procedures, detecting sounds generated by the heart and sounds from vortices in the carotid artery for at least 30 seconds. A method for determining stenosis of the carotid artery in a human patient consisting of a first step of placing a sensing device comprising an array and three sensing elements onto the patient, wherein a first sensing element is placed near the heart and the two remaining sensing elements are placed adjacent to the carotid arteries; the sensing elements then measure sounds from each of the three sensing elements, resulting in sound from three channels. The sound is measured in analog and modified to digital format and then each of the three channels are analyzed before a power spectral density analysis is performed. The power spectral density graph reveals peaks that are not due to noise, that are then analyzed to provide for a calculation of percent stenosis or complete occlusion of the carotid artery.
Claims
1. A method for measuring sound from vortices in the carotid artery comprising: a. performing a first quality control procedure on at least two sensing elements, wherein said quality control procedure is performed by playing a predetermined set of tones within a base unit, wherein said at least two sensing elements detect said set of tones and wherein said detected tones are compared to said predetermined set of tones, wherein the sensing elements are replaced if the comparison between said detected tones and said predetermined tones has a variance of more than 10% relative to a frequency; b. performing a second quality control procedure on at least two sensing elements, wherein said second quality control procedure is performed by detecting blood flow through the carotid artery and comparing said detected sounds to a predetermined sound signature, wherein the sensing element is repositioned if the detected sounds compared to the predetermined sound signature have a variance of more than 25% relative to the frequency, wherein if the variance in the second quality control procedure is more than 100% relative to the frequency, step (a) is repeated; and c. detecting sounds generated by the vortices in the carotid artery for at least 30 seconds.
2. The method of claim 1, wherein the sounds detected from the vortices in the carotid artery are between 40 Hz and 3000 Hz.
3. The method of claim 1, wherein a further step (d) comprises eliminating sounds from the carotid artery that are outside of the range of 40 Hz and 3000 Hz.
4. The method of claim 3, comprising a further step (e) comprising generating a power spectral density graph of the sounds from step (d).
5. The method of claim 1 comprising three sensor pods.
6. A method for measuring vortices produced in the carotid artery due to plaque accumulation in the artery comprising: a. performing a first quality control procedure on at least two sensing elements, wherein said quality control procedure is performed by playing a predetermined set of tones within a base unit, wherein said at least two sensing elements detect said predetermined set of tones forming detected tones, and wherein said detected tones are compared to said predetermined set of tones, wherein if said detected tones are within 10% of a frequency, the first quality control procedure is passed, and wherein if said detected tones are not within 10% of the frequency, the first quality control fails, then replace one or more sensing elements; b. performing a second quality control procedure on at least two sensing elements, wherein said second quality control procedure is performed by detecting sounds generated by blood flow through the carotid artery; wherein said at least two sensing elements detect said sounds generated by blood flow through the carotid artery, and said detected sounds are compared to a previously recorded sound signature, indicating an appropriate position for the one or more sensing elements if the detected sounds are within 25% of a frequency of the sound signature, or repositioning the one or more sensing elements if the detected sounds are greater than 25% of the frequency of the sound signal and wherein if the variance in the second quality control procedure is more than 100% relative to the frequency, step (a) is repeated; and c. detecting sounds generated by sounds from vortices in the carotid artery for at least 30 seconds.
7. The method of claim 6 comprising three sensor pods, wherein in step (c), detection of sounds generated by sounds from the vortices in the carotid artery are detected simultaneously by the sensor pods.
8. The method of claim 6, wherein the sounds detected in step (c) are between 20 and 3000 Hz.
9. The method of claim 6, further comprising: e. down sampling the detected sounds from step (b) from analog to digital at a sampling rate of 20 KHz; and f. removing sounds from the digital outside of the 40 Hz to 3000 Hz range.
10. The method of claim 9 comprising a further step (g) of generating a Power Spectral Density plot and detecting peaks in said plot.
11. The method of claim 10 comprising a further step (h) of determining percent stenosis from the peaks in said plot by calculating (1−f1/f2)×100.
12. The method of claims 1 or 6 performed by a device for detecting stenosis in the arterial circulatory system comprising a base and at least one sensor pod; said base comprising a processor and a speaker, capable of playing a predetermined sound through said speaker; said sensor pod comprising a circular piezo cap comprising a top and a bottom and an inner face and an outer face, with an opening between the top and bottom with the opening larger at the top than the opening at the bottom; a flange positioned on the inner face of the opening; a piezo having a top, a bottom, and a perimeter support; said piezo disposed of within said opening, with the bottom of the perimeter support engaged to and adhered to said flange; a printed circuit board having a ring shape and an outer diameter to fit within the opening and engaged to the bottom of said flange; and on said inner face one-half of an attachment means for securing said disposably piezo assembly to an assembly base.
13. A method for detecting stenosis of the arterial circulatory system comprising: performing a self-diagnosis quality control procedure on a sensor element by playing a predetermined sound signature from a speaker; detecting said predetermined sound signature with said sensor element; comparing said detected sound signature to said predetermined sound signature; proceeding to a second quality control procedure if said detected sound is within 25% of a frequency of the predetermined sound signature or replacing said sensor element if said detected sound is more than 25% from the frequency of the predetermined sound signature; placing said sensor element on an artery of interest; detecting a flow of fluid through said artery of interest wherein the sensing element is repositioned if the detected sounds compared to the predetermined sound signature have a variance of more than 25% relative to the frequency, wherein if the variance in the second quality control procedure is more than 100% relative to the frequency, the self-diagnosis quality control procedure is repeated; detecting a frequency of between 60 and 260 Hz to confirm proper position of said sensing element; moving said sensing element to a different position if a frequency between 60 and 260 Hz is not detected; upon detecting said frequency between 60 and 260 Hz, capturing data from said sensing element; plotting a Power Spectral Density Plot; calculating stenosis based on (1−f1/f2)×100.
14. The method of claim 13 further comprising performing a wavelet analysis after capturing data from said sensing element.
15. The method of claim 14 further comprising performing Burg's Method after the wavelet analysis.
16. The method of claim 15 further comprising performing Welch's method after performing Burg's Method.
17. The method of claim 13 wherein the calculation of stenosis is a binary calculation of greater than or less than 50%.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
(1)
(2)
(3)
(4)
(5)
(6)
(7)
(8)
(9)
(10)
(11)
(12)
(13)
(14)
(15)
(16)
(17)
(18)
(19)
(20)
(21)
(22)
(23)
(24)
(25)
(26)
(27)
(28)
(29)
(30)
(31)
(32)
(33)
(34)
(35)
(36)
(37)
(38)
(39)
(40)
(41)
(42)
(43)
(44)
(45)
(46)
(47)
(48)
(49)
(50)
(51)
(52)
(53)
(54)
(55)
(56)
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS
(57) The embodiments of the invention and the various features and advantages thereto are more fully explained with references to the non-limiting embodiments and examples that are described and set forth in the following descriptions of those examples. Descriptions of well-known components and techniques may be omitted to avoid obscuring the invention. The examples used herein are intended merely to facilitate an understanding of ways in which the invention may be practiced and to further enable those skilled in the art to practice the invention. Accordingly, the examples and embodiments set forth herein should not be construed as limiting the scope of the invention, which is defined by the appended claims.
(58) As used herein, terms such as “a,” “an,” and “the” include singular and plural referents unless the context clearly demands otherwise.
(59) All patents and publications cited herein are hereby fully incorporated by reference in their entirety. The citation of any publication is for its disclosure prior to the filing date and should not be construed as an admission that such publication is prior art or that the present invention is not entitled to antedate such publication by virtue of prior invention.
(60) The embodiments contemplate devices, systems, and methods for determining blockage in a fluid flow vessel. To reliably determine fluid flow, we need to determine that the components of the device are working properly, are clean and sanitary, are positioned in the correct locations for detection. Furthermore, the device needs to passively prevent ambient noise from entering the sensing device. However, active noise cancellation strategies can further eliminate ambient noise. Finally, processing strategies can be utilized to filter the collected data and to break it apart into useable packets of data for determination of occlusion in a fluid flow vessel.
(61) For many cases, fluid flow vessels include the arterial circulatory system, for example the carotid artery, but also the arteries of the heart, the coronary arteries. However, flow through industrial pipes can also be evaluated using the devices and methods described herein.
(62) Description of Ring Vortices being Detected
(63)
(64) The ring vortices are produced equidistant from each other at a distance between them equal to their diameter as they move downstream, as illustrated in
(65) In
(66) At RE less than 800 or greater than 2100, ring vortices do not form. The closer to 800 while still remaining below 800, the more string-like motions are seen, as seen in
(67)
(68) If no f1 appears in the PSD (between 60 and 260 Hz), there was insufficient energy in the flow emerging from the stenotic region for the vortices to reach Region II, in which the larger vortices appear, at the lower frequencies. This indicates the artery is heavily stenosed. If there is no f2, there is an insufficient amount of stenosis to create the smaller vortices (Region I) indicating a low level of stenosis (below 15%) as reported by Khalifa and Giddens [“Characterization and evolution of poststeotic flow disturbances,” Journal of Biomechanics 1981. Vol. 14, No. 5, pg292] who report that below 25% reduction in area due to stenosis (which corresponds to a reduction of 13% in diameter), no signal is picked up. If there is neither f1 nor f2, the indication is that there is a near blockage level of stenosis, as the vortices cannot be produced even when the velocity is sufficient to give RE between 800 and 2100.
(69) To measure the large ring vortices, we need to ensure that the device we are using contains properly sterile and functioning elements. Described herein are certain disposable components, methods for determining proper function of these elements, and methods for eliminating and reducing noise from the data sample in order to accurately and efficiently measure and quantify stenosis in the arterial circulatory system.
(70) Furthermore, these aspects and teachings can be applied into industrial structures. For example, these same perturbations that are present in industrial piping, such as fluid flow in gas an oil industries, production of fats, oils, and other consumer goods, chemical and biological production, and the like. Representative perturbations are depicted, for example in
(71) Replacement components provide for accurate and clean components that ensure greater chance of accuracy and reproducibility. Piezoelectric sensors have a variety of potential uses, but as described herein, they are being utilized as a contact microphone. The principle of operation of a piezoelectric sensor is that a physical dimension, transformed into a force, acts on two opposing faces of the sensing element. Detection of pressure variations in the form of sound is the most common sensor application, e.g. acting as a microphone, wherein the sound waves bend the piezoelectric material creating changing voltage. Accordingly, the piezo sensor can be placed on or near a sound to receive the sounds.
(72) Piezo sensors are especially used with high frequency sound in ultrasonic transducers for medical imaging and industrial nondestructive testing. However, piezo sensors are also frequently used for the detection and activation of a device, based on the ability to receive a signal and to then send an electronic signal, thereby acting as the actuator. In the embodiments herein, piezoelectric sensors (“Piezo”) are utilized for their ability to detect certain frequency sounds or vibrations caused by the distortion of a fluid flow vessel, specifically of the arterial circulatory system.
(73) Because of the sensitivity of these sensors, piezoelectric sensors can be somewhat fragile and can be broken from both normal use and misuse. Furthermore, as utilized in a medical device, there is the inherent need to ensure accuracy of each of the three piezoelectric sensors. Accordingly, any slight modification of the sensor may result in a modification of the input received and thus would result in erroneous data.
(74) Replacement components may be one of three different components as described herein. A first component may be a disposable piezo assembly, a second component may be a sensor pod, which comprises the disposable piezo assembly and a sensor base, and a third component may be a disposable array, comprising one or more sensor pods. In this manner, each component may be disposable to allow for easy replacement after use.
(75) Piezo sensors can include any number of materials. Typically, however, the sensor contains a portion of ceramic material and a metallic component. Piezo sensors may also use a polymer film configuration which exhibits a low acoustic impedance similar to that of human tissue, or made of metallic materials. These sensors, as used in the invention herein, are typically a circular shape with a diameter of about 3 inches. Typical piezos have a diameter from about 0.01 to about 6 inches for use in medical settings, with most typical sizes between about 0.5 to about 4 inches in diameter. For most applications, including industrial settings, a range of 0.01 inch to about 12.0 inches is preferred, wherein the size of the piezo is generally related to the diameter of the fluid flow vessel to be measured. In preferred embodiments, the fluid flow vessels are veins and arteries in the body, for which a 4.0 inch or smaller diameter piezo is preferable.
(76) There is no inherent frequency limit for a piezoelectric sensor. However, the limits of applications are usually determined by resonances associated with the shape and/or the size of the transducer design. The Piezo sensors utilized herein have a thickness of about 0.01 to 2.0 mm and are capable of detecting sounds between 10 Hz and 32 KHz and an amplitude of 0.0002 N/m2 to greater than 10 N/m2. In preferred embodiments, the piezos attached to a sensor pod detect sounds between about 20 to 3000 Hz, which are relevant towards measurements of fluid flow in the body. Typically, these sounds have an amplitude of between 0.002 N/m2 and 20 N/m2. While additional sounds are recorded, many of these sounds, i.e. the heart beat and extraneous noise, are removed from the data set through several filters.
(77)
(78) In the broadest sense, the piezo sensors are disposed of within a pod. On one side of the piezo is placed a sensor pad, for example those of 1, 2, 17 and 19. The sensor pad is then pressed against the skin or clothing of a patient to listen to the underlying circulatory system. The sensor pad allows for transmission of energy waves, sound and vibrations, which are received by the piezo element. Gel or other impedance matching substance may be applied to the skin facing surface of the pad.
(79) In view of
(80) There are several ways in which the piezoelectric elements can wear or be damaged including ordinary and standard use of the device. Ordinary wear may occur as the piezoelectric element wears from ordinary and standard use, and after about 10 to about 400 uses, the piezoelectric element breaks down so that the function and the electrical currents generated are different when comparing the first use to the 2.sup.nd, 5.sup.th, 10.sup.th, 25.sup.th, 50.sup.th, 75.sup.th, 100.sup.th, 200.sup.th, 300.sup.th, or 400.sup.th use and all numbers in between. Accordingly, to ensure that accurate results are received by each of the units, it is imperative to replace the unit that has worn to maintain consistent results.
(81) Additional wear or breakage can occur to the piezoelectric sensors by error or accident. For example, human error may lead to the array being dropped, or placed onto the base in a manner that breaks, bends, or otherwise damages the piezoelectric unit. Further damage may occur as clean sensor pads are attached and placed against the piezoelectric sensor for use on a patient.
(82) To ensure sanitary use of the device, the sensor pads are replaced between each use of the device. However, because the sensor pads are placed directly onto the piezoelectric unit, there is risk that human error may damage the piezoelectric sensor, either by too much force, or simply through improper pressure applied to the piezo when installing or removing a sensor pad.
(83) Ordinary wear or accidental damage is tested through routine quality control procedures performed in a self-diagnosis module. The sensor pods can be placed in a base or holding device that comprises a speaker embedded within the base which provides a predetermined sound that can be measured by each piezoelectric sensor. When the sensor device is activated for use, the sound, which can include both audible and inaudible sound waves, is played for between about 1 and about 20 seconds. During the time that the sound is playing, each of the piezoelectric sensors records the sound and a program then confirms that each of the three sensors is recording the appropriate sounds being played. If each of the three sensors detects the appropriate sounds, then the sensor device is ready for use. However, if one or more of the sensors detects sounds that do not match with the predicted sounds, the device will provide an alert, which may include lights, sounds, or other display elements, to alert the user of the device that one or more of the piezos needs to be replaced.
(84) An optional display screen attached to the base can further display the device and identify the sensor pod containing the piezo that failed the QC test. Another manner for identifying the failed sensor is to have lights that correspond to working or failed tests either on the base or on the sensor array itself. Once the failed piezo is identified, the user can then replace one or more of the components, as described herein, and then perform the QC test again to ensure that the device is now ready for use.
(85) Accordingly, in a preferred method, a piezo is replaced every 10 uses to ensure that there is no noticeable wear and tear on the piezo, and to prevent the possibility of erroneous data. Accordingly, the sensor device comprises a counter wherein the number of times that a test is run with each of the piezo is counted, so that the sensor device notifies a user that the piezo needs to be replaced, even if each of the piezos are working properly.
(86) In other embodiments, the piezos can be replaced every 1, 2, 5, 10, 25, 50, 75 uses, 100 uses, 125 uses, 150 uses, about every 200 uses, or about every 400 uses or a number in-between. The particular number of uses for each piezo will be determined through additional use of the devices in normal practices. However, to ensure sanitary and consistent results, it is preferred that the piezos are changed after no more than 100 uses.
(87) To facilitate easy changing of the disposable piezo assembly 85, the disposable piezo assembly 85 is able to easily attach to an underlying disposable sensor base 86, and to be replaced. For example, a simple threaded attachment mechanism allows the sensor pod to be removed from the sliding sensor pod base, which is attached to the sensor array. Alternatively quarter, or half-turn attachment means, magnetic attachment, and others as known to one of ordinary skill in the art are known.
(88)
(89) Near the vertex of the Y is a charging port, 820, and a PCB charging contact 131 disposed therein. This allows the array to be placed into a charging port and charge a central battery.
(90) Attached to the array is a sensor pod, made up of the components of a locking cap 125, a DBM 120, a PCB processor board 110, a PCB housing 115, a piezo cap 100, a piezo 90, and a disposable piezo assembly 85. These features are further detailed below. A disposable sensor pad 18 can be affixed to the piezo 90 via adhesives or by the natural adhesion of the pad material. For example, the piezo cap 100 can be attached to the PCB housing 115 in several ways, including as in
(91)
(92) In one embodiment, this disposable assembly 85 is the smallest disposable component, which allows for quick and easy replacement of the piezo without replacement of any further components (except for the disposable sensor pad 18, which is replaced for every use). The disposable assembly 85 comprises a quarter turn locking feature 101 that corresponds to a paired feature 116 on the PCB housing 115. This allows for a small turn of the disposable assembly 85 to remove the component and replace. Additional attachment mechanisms can be easily exchanged, for example magnetic, threaded engagement, or simply a threaded fastener or two that can be engaged for replacement. Finger capable fasteners can use a full, half, or quarter twist to secure a fastener between two components. A person of skill in the art will recognize that numerous options exist for attaching and detaching such components and that attaching means incorporates these listed and additional options not described in detail herein.
(93) The PCB housing contains a locking groove 117 that engages with and locks the elastomer DBM 120 to the PCB housing 115. In particular locking groove 117 engages locking key 121 between the locking cap 125 and the PCB housing 115. A locking cap 125 engages to a fastener 113 to secure the key 121. A second key 122, is also provided to lock the DBM 120 between the outer array housing 140 and the inner array housing 130. A further detail of these locking features are provided in
(94) While the disposable assembly 85 can be easily removed and replaced, it is also contemplated that the entire sensor pod can be removed and replaced easily. For example, removal of threaded fasteners 133 will allow for quick and easy replacement of the entirety of the pod, inclusive of the DBM 120. Furthermore, the DBM 120 can be held in place, and the locking cap 125 can reveal a threaded fastener 113 to replace the remaining components. In the Fig, the fastener 113 can be oriented in either direction to allow for quick replacement.
(95)
(96)
(97)
(98)
(99)
(100)
(101) The DBM 120 is a circular feature having an inner opening. At the outer edge of the DBM 120 is an outer flange 122. At the circumference of the inner opening, there is an inner flange 121. These flanges 122 and 121 are used to lock the DBM 120 into place between the array features 130 and 140, as well as between the locking cap 125 and the housing 115.
(102) Therefore, the DBM 120 is an elastomeric material, capable of allowing the attached piezo to flex in any direction, as well as move away from the surface to be compressed. This allows for a consistent pressure to be applied to the skin surface by the sensor pad 18, based on the rigidity of the membrane 120.
(103)
(104)
(105)
(106) In an ideal world, every patient would be the same shape and size, and modification of the structure would not be required. However, in practice, men, women, and children have significantly different shapes and sizes due to the amount of body mass, muscle, breast tissue, fat deposits, etc. Specifically, changes in body mass and shape between the neck and the torso create issues where the array must be modified to position one or more sensors in appropriate positions for acoustic sensing.
(107) Therefore, as used on human patients, a difficulty in such devices is that people come in all shapes and sizes and that the array must be easily modified to fit these different shapes and sizes. One option would be to utilize different sized, fixed position sensing elements, due to the fragile nature of the sensing elements. However, constant movement and replacement of the sensing elements from one device to another would likely result in more damage to the sensing elements and increase the risk for the need for frequent replacement of these elements. Therefore, an array with rails, both the neck and “Y” versions, provides the necessary stability and flexibility provides a great advantage in the array for use on patients.
(108) A particular feature of the sensor pods when affixed to an array is that they are adjustable and can be configured to account for the anatomical differences between individuals while remaining sufficiently rigid to support the sensing elements. Such flexibility can be seen in the depiction of
(109) The exploded view of
(110)
(111)
(112)
(113)
(114)
(115) A curved film piezo can be exchanged for any of the piezos in embodiments described herein. For example, the lower frame 604 may comprise a relevant attachment means, and further comprise a PCB contact point to allow for direct exchange with prior examples and figures.
(116)
(117) The sensor pods including both 85 and 86 components, are replaced, as necessary to allow for proper functioning of the piezo sensor. These replacements are performed as necessary, but at least every 10, 25, 50, 75, 100, 150, or 200 tests. When the sensor base 86 is replaced, the disposable piezo assembly 85 is also replaced. By contrast, in each test, sensor pads 18 are replaced.
(118) In certain preferred embodiments, the sensor pads 18 can be secured onto the piezoelectric unit via an adhesive, such as one of several common low tack adhesives for providing for a temporary securing of the sensor pad to the piezo element. Other embodiments may utilize a gel or other water or solvent based material that may secure the sensor pads without the need for an additional adhesive material. In further embodiments, the sensor pad fits into the sensor pod and secures onto the piezo without the need for any adhesive.
(119) A particular feature of the sensor pads described in the embodiments herein is the fact that the top face shape (that contacts the patient), and the bottom face shape (that contacts the piezo) are made so that when the top face contacts the patient and thus applies pressure to the sensor pad and through to the bottom face, the piezo does not flex when pressure is applied to the sensor pad. This is important to ensure consistency and accuracy of the piezo device. Therefore, the sensor pad, in certain embodiments, is designed such that the piezo does not flex when pressure is applied to the sensor pad. In a further preferred embodiment, the piezo flexes less than about 0.1%, 0.5%, 1.0%, 5.0%, 20%, and 25% and all percentages in between. Accordingly, in certain embodiments, the amount of flex is greater than zero (i.e. rigid and does not flex), but the amount of flex is minimized to maintain accuracy of the piezoelectric unit.
(120) It is also preferred that the sensor pads create a proper impedance matching with a patient. Accordingly, the sensor pad is designed to have a slight tackiness which ensures a proper impedance matching with the patient, which then successfully transfers sounds through to the piezo element so that the piezo can properly detect vibrations and noise signals from the patient.
(121) Therefore, in order to maintain both sterility of the medical device and proper function of the medical device, it is necessary to provide replaceable components. The entire device is a complex system comprising a display, a base unit, an array, a sensor base, a disposable piezo assembly, and a sensor pad. Each of the last four are disposable. The array itself can be disposed of after a number of uses, likely between 100-1000 uses. The array may lose elasticity to ensure proper fit on a patient, gain cracks, or simply lose stability. Each of these may increase variability and thus replacement is warranted.
(122) The sensor base as depicted in
(123) The disposable piezo assembly is intended for more frequent replacement than the base or the array, as the piezo is susceptible to wear or damage. Accordingly, frequent changes, such as between every use and every 10, 25, 50, or 100 uses is necessary for accurate results.
(124) The device is a complex system comprising multiple components, each working together to ensure that accurate results are obtained. Disposable components ensure that the system works properly, every time, and that it generates accurate and reliable data.
(125) A kit is envisioned with the system, wherein a plurality of sensor pads are provided, a plurality of disposable piezo assemblies are provided, at least two sensor base assemblies, and at least two arrays. Said kit can be used with a system comprising the base and a display, as well as necessary software and hardware for energizing and running the device through its necessary protocols.
(126) Quality Control Methodologies and Devices
(127) Now that we have a device that is clean and has readily replaceable components, we need to ensure that the device is properly functioning. Accordingly, we describe certain methods and embodiments that provide for self-diagnostic tests, active diagnostic tests, and guidance for properly positioning a sensor on a patient.
(128) The quality control protocols embodiments provide for a process or method for determining if a listening device, such as a piezoelectric device, or microphone, is properly functioning. This is a self-diagnostic quality control feature. A second test is an active quality control procedure, which is performed with sensors on a patient. The two tests can be used alone, each being sufficient to confirm that the sensor is working properly, or can be used together, to both ensure proper function and also proper placement of the sensors on a patient. When performed together, the tests are performed sequentially, first the self-diagnostic test and then the active diagnostic test on the patient.
(129) Accordingly, in preferred embodiments, methods exist for determining the proper function of the sensitive piezoelectric components.
(130)
(131) In one embodiment, disposed of within the base 300, and specifically adjacent to the cradle for each of the sensor pods 1, is a respective speaker 97. A computer is coupled to the base 300 for communication via a USB connection, Bluetooth, near field communication, RS-232, or the like. The computer couples to the speaker 97, and when the SDD (Stenosis Detection Device) is activated, a program is executed by the computer system so that it performs a diagnostic and quality control test on each of the sensor pods 1.
(132) The diagnostic and quality control procedure comprises a program that plays a known set of sounds generally corresponding to sounds that will be detected and recorded when measuring sounds on the body of a patient. These sounds include low and high frequency sounds, typically low amplitude. Once the sounds are played, the sensor pods 1 detect the sounds and convert the sound to a digital signal that is plotted and compared to a predetermined plot of the sounds that were played. Alternatively, an analog signal is generated and compared with the predetermined plot. Each of the sensor pods 1 is independently tested to determine if it meets an acceptable standard. In one embodiment, and error message is generated if the sensor pod output is not within 10 percent of the predetermined plot at a given data point. Other standards can be used to determine an error condition exists. A range of 1 to 50 percent at each data point can be used to determine if the sensor pod 1 is not functioning properly. Alternatively, the overall plot can be analyzed, instead of a point-by-point analysis, to determine if a sensor pod 1 is functioning properly. Typically, a sensor should be within 25% of a predetermined frequency.
(133) If any sensor pod is not detecting an appropriate sound, then the system will notify the user of an error. In most instances, the error means that a particular sensor pod has exceeded its useful lifetime and is due for replacement. These devices theoretically have a lifespan of several hundred uses under ideal conditions. However, in a medical office, the continuous placing of the array 5 on to a patient, and detecting and recording real sounds, may result in distortion after even a few uses. Accordingly, the system is able to determine whether the detected sounds are simply drift that is a slight change in the detected sounds, or whether there is an error or fault in one of the sensors. If there is only a slight drift, the system can calibrate each unit so that the measured noises from the system are consistent through use.
(134) If the measured sounds are greater than a tolerance of more than 10%, or more than 25% as defined for the occasion, the system notifies the user through images on a display, lights on the sensor pod, audible messages, or other manner to communicate the error, and identifies which sensor pod is faulty. A user can then quickly replace the faulty sensor pod or the disposable piezo assembly 85, and re-run the quality and calibration control program.
(135) After the sensor pod is replaced and the quality control program is re-run, and the replacement sensor pod is confirmed to be working properly, the system will alert that it is ready for placing on a patient. Each of the sensor pods can be appropriately placed onto the patient.
(136)
(137) The diagnostic and quality control procedure is depicted in a flow-chart of
(138) When performing the test in step 517, the sounds include low and high frequency sounds, typically at low amplitude corresponding to the range of sounds to be detected by the SDD device. Once the sounds are played, the sensor pods detect the sounds and convert the sound to digital 519. The criteria step 522 compares the digital sounds received to the actual sounds played. For example, a comparison can be made between amplitude and frequency, and overlaid to compare the two samples. Each of the sensor pods is independently determined to meet an acceptable standard, or tolerance for example within 50%, 25%, 10%, 5%, or within about 1% of the sounds based on the determined Hz and, optionally, the amplitude of the detected sounds. Simply comparison software can make these comparisons between the two sounds.
(139) If any sensor pod is not detecting an appropriate sound, then the system will notify the user of an error. In most instances, the error means that the particular sensor pod is due for replacement. While these devices may theoretically have a lifespan of several hundred uses under perfect conditions, the reality of a medical office and placing a device on or adjacent to a patient and detecting and recording real sounds may cause distortion after even a few uses. Accordingly, the system is able to detect and determine whether the sounds detected are simply drift that is a slight change in the detected sounds, or whether there is an error or fault in one of the sensors, thus requiring replacement. If there is only a slight drift, the system can calibrate each unit so that the measured noises from the system are consistent through use. An appropriate program on the system can make these changes to the data based on the actual versus detected sounds, through a simple calibration program. Accordingly, the played tones provide for the ability to both detect and calibrate the device before every use.
(140) If the measured sounds differ by more than the acceptable tolerance, the system engages the user through images on the display, lights on the sensor pod, audible messages, or other means for communicating error, and wherein the particular sensor pod that is faulty is identified. A user can then quickly replace the faulty sensor pod or disposable piezo assembly 85, and re-run the quality control program. An exploded view of a sensor pod is depicted in
(141) For example,
(142)
(143) When the disposable piezo assembly 85 is attached, it contacts the PCB Processor board 110, which assembles into a pocket in 115, and is captured by 85. In this manner, when a quality control test is performed, and a sensor is identified as faulty, the attachment means can be withdrawn and the disposable piezo assembly 85 can be removed and a new disposable piezo assembly 85 attached and the test re-run.
(144) In certain embodiments, it is advantageous to have the entire sensor pod replaced, not just the top disposable component. For example, the PCB board 110 may in some instances wear or be damaged. Alternatively, the diaphragm bellows membrane 120 may need replacement, or simply replacement is warranted because of contamination concerns. Accordingly, the entire piezo assembly can be replaced, by removing threaded fasteners 133 or by removing locking cap 125.
(145) The diaphragm bellows membrane 120 locks with certain features, to ensure that it can freely flex and compress to allow for the fit of the piezo against the body. The diaphragm bellows membrane 120 fits feature 121 into a locking groove 117, which traps locking feature 121 between locking cap 125 and the PCB housing 115. Locking feature 122 secures the diaphragm bellows membrane 120 between the inner array halve 130 and the outer array halve 140. This creates a flexible “drum head”.
(146) For each use of the piezo, a sensor pad 18 is also utilized for sanitary conditions and to ensure a quality sound contact to the piezo 90. The sensor pod 1 of
(147) After either replacement of the disposable component 85 or replacement of the entire sensor pod, the quality control program is re-run and the replacement sensor pod is confirmed to be working properly, the system will alert that it is ready for placing on a patient. Each of the sensor pods can be appropriately placed onto the patient, as depicted in
(148) As depicted in
(149) As with the quality control procedure on the base unit, once the sensor pods are placed on the patient, the operator can engage the device to begin detection and recording on the patient. Because the sounds that are being detected and recorded are known within a certain range of sounds, that is, the sounds are generally known to a certain frequency and amplitude, and a further quality control test is performed for a duration of between 1 and 30 seconds. This test provides a quality control diagnostic to ensure that the sensor pods are detecting proper sounds from the patient, and thus confirms two pieces of information: first the proper placement of the sensor pods on the patient; and second that the sensor has not failed in the time between initial quality control tests and placement on the patient.
(150) Since there are at least two and likely three sensor pods, each pod communicates with the computer identifying the detected sounds, which can be recorded by the system and compared in real time to a predicted sound. Accordingly, the sensor pod at the heart will predict a certain sound and the sensor pod(s) at the carotid arteries another sound. If one or more sensors does not detect the predicted sounds, a signal will engage to identify the sensor that is not properly detecting the predicted sound. This signal will alert the operator that the sensor pod needs to be adjusted to a different position to properly detect the sounds for the particular test.
(151)
(152) If the criteria is met, 513, then we proceed to start recording the data and processing the patient 516. However, if the criteria is not met, we need to first adjust the piezo on the patient 514. Adjustments can be just a few centimeters, or more as necessary, in order to get the piezo closer to the artery of interest. After adjustment the device again receives sounds from the patient 511 and compares the sounds to the expected sounds 512 to determine if the criteria is met.
(153) In certain instances, after movement and adjustment of the device, the piezo is still not finding the proper sounds. This can be due to continued improper placement or failure. Accordingly, it is best to replace the piezo 515 and start another quality control procedure as outlined above on the base.
(154) The embodiments of the system utilize variations of quality control programs for initial setup testing of the sensor pods and then for quality control testing of the proper position on the patient. A variety of alarms, indicators, or signals can be utilized in each of the quality control programs to ensure that the issue is detected and corrected.
(155) For the initial quality control program, when the sensor pods are still in the base unit cradle, it is appropriate to indicate a fault with a computer Graphical User Interface (GUI) as depicted in
(156) In other embodiments, a colored light system, such as a green or red light based on green being good, and red signaling an error with the sensor pod can be directly placed on the sensor pods (see
(157)
(158) The lights of
(159) In certain embodiments, a button on the device or on the base is pressed to perform the active diagnostic phase. However, in preferred embodiments, once the self-diagnostic test is complete, the active diagnostic phase immediately starts. The active diagnostic phase will continue, until either all sensors indicate green or one indicates red. Typically, this will last up to 30 seconds, at which time a red light will indicate to re-start the test, or to replace a sensor.
(160) If one sensor remains yellow or yellow with green/red, during the active diagnostic step, the lights, visual, and or audible alarms can further assist in positioning the device properly on a patient. For example, the light remaining yellow will turn to yellow and green, if the signal is better, or from yellow to yellow and red, if the signal is worse. Accordingly, the sensor can be moved in a proper direction towards the yellow/green until just a green light is indicated. Furthermore the GUI can be utilized in the same manner, with an indicator on the screen suggesting the direction to move the sensor. Ultimately, if a sensor pod does not detect the proper sounds from the patient, then one or more alarms will register and the operator will know that one or more sensor pods need to be replaced on the patient. In certain embodiments, the visual screen, a visual identifier will flash to aid the operator in placing the sensor pod in the proper location.
(161) In further embodiments, where a sensor pod is identifying an improper sound or not detecting a sound, a visual alarm may be generated, such as a red light, which indicates improper position or a sensor failure. The SDD can detect and compare the sounds in real-time, so the operator can then slowly move the sensor pod to a different location and wait a few seconds to see if the light turns from red to green, indicating a proper position. The operator can continue to move the sensor pod on the patient until it is indicated on either the sensor pod, on the array, or on the SDD device display that the position is correct.
(162) If the operator is unable to determine a proper location on the patient after 30 seconds, the SDD will alarm with a visual or audio signal to perform a base unit quality control procedure again to ensure that the sensor pods are all functioning correctly, or to simply replace the sensor that indicated failure. After replacement or if the sensor pods are determined to be functioning correctly, the operator can again restart the process of placing the sensor pods on the patient.
(163) Accordingly, a preferred embodiment for determining proper placement of sensor pods on a patient comprises a stenosis detection system comprising a base unit having a cradle, at least two sensor pods, a display and at least one alarm mechanism; wherein while the sensor pods are engaged in the base unit cradle a self-diagnostic quality control procedure is performed to confirm that the sensor pods are properly functioning. After confirmation of the proper function of each of the sensor pods, the devices can be placed onto a patient wherein an active quality control procedure is performed. The active quality control program is run for between 1 and 30 seconds wherein each sensor pod is communicating with the compute of the detection system in real-time to ensure that each of the sensor pods is measuring the appropriate sounds. Wherein the system provides for an audio or visual notification that the active quality control program is met, or wherein the system identifies one or more sensor pods that are improperly placed. Wherein the system then provides an alarm to any sensor pod that is not properly placed. Wherein a visual or audio mechanism is provided to provide real-time feedback as to the proper position for each sensor pod, and wherein one example provides for a red light for improper position and green light for a proper position. Certain embodiments utilize a yellow light to indicate that one or more of the self-diagnostic test or active diagnostic test are proceeding.
(164) Other audio or visual alarms or mechanism may be further included in the system so as to aid in the placement of the sensor pods on a patient.
(165) In preferred embodiments, the active quality control step on the patient provides for immediate real-time feedback to the correct placement of each sensor pod to ensure fast and reliable positioning of the sensor pods, and also to confirm fast, precise, and accurate detection and determination of stenosis on the patient.
(166) The method comprises: Performing a first base unit quality control test; confirming that each of the sensor pods is properly functioning; placing sensor pods on a patient; performing a second quality control test, wherein the sensor pods detect sound in real-time and compare said sound to a predicted sound; and indicating with an alarm whether the sensor pod is properly placed on the patient by comparing the detected sound in real-time to a predicted sound based on historical data.
(167) In a preferred embodiment the system uses a computer to run software to implement the features as described in the embodiments herein. Accordingly, the computer is connected to the array and/or to the sensor pods via a connection means either wired or wireless, as is known to one of ordinary skill in the art. The software comprises the various quality control procedures, as well as appropriate code to provide alarms and to notify of the need for replacement or modification. Further features include the ability to calibrate the system in view of a quality control test.
(168) Therefore, preferred embodiments of the disclosure comprise a method of confirming the proper position of a medical device upon a patient comprising: performing a first quality control procedure to ensure functioning of the sensor pods, comprising playing a predetermined set of sounds and comparing the predetermined sounds to the detected sounds; performing a second quality control procedure while detecting sounds from a patient wherein the test compares the detected sounds to sounds that are ordinarily present in detection of the particular artery or vessel of interest; and triggering an alarm wherein the detected sound does not meet the predicted sound, or triggering an approval if the detected sound confirms with the predicted sound.
(169) Noise Attenuating Strategies
(170) A major hurdle in creating a device that conforms to the necessary levels of accuracy is to ensure that the data received for each test is of the highest quality. By performing the prior quality control procedures, the devices are known to be functioning properly. However, it is necessary to now utilize passive and active noise attenuating strategies, as well as computer implemented de-noising strategies to generate clean and clear data. Accordingly, we need to eliminate noise from the data sample in any number of ways, so that the resulting data is clean and clear for quantification of stenosis.
(171) The noises that we are particularly measuring are subtle large ring vortexes. These vortexes are created as wall pressure fluctuations distal to a constriction (stenosis) in rigid or elastic pipes, or in arteries, reveal the presence of low-frequency maxima. These fluctuations are found to be associated with large-scale, medium-scale, or small-scale vortices (also called “eddies” if small), that are strong in the region distal to the constriction (called “stenosis” when in an artery).
(172) Normal blood flow in a heathy patient causes certain sounds which are detectable by our device. Patients which have stenosis in the carotid arteries will often have another 2 or 3 additional sounds that can be picked up by our device. Depending on the amount of stenosis and how many stenosed areas the sound will change. The carotid artery has a branch which feeds two main areas in the head. One main branch going to the brain and the other branch going to the face. The area that we test for is where the carotid artery branches into these two areas. Thus depending if there is stenosis in one branch or two can lead to multiple sounds being picked up. Because these sounds/vibrations are at such a low level it is vital to make sure as much external noise is eliminated as possible. Even small noises in the 20-3000 Hz range can overwhelm the noises we are looking for making noise elimination critical.
(173) With regard to flow and the noises created therein, some of the fluid-flow energy enters into the vortex motions distal to a constriction, which then results in an increase in the wall pressure amplitude, above that of turbulence alone, at the lower frequency end of the wall pressure power spectrum. These maxima are nearly Gaussian-shaped bell curves situated atop a broad, nearly flat spectrum at low frequencies that is due to turbulence within the pipe or artery. The maxima are always found at lower frequencies than the so-called “break” frequency characteristic of the turbulence spectrum where the latter changes quite abruptly from nearly flat to steep declining in intensity (when the logarithm of signal intensity is plotted versus a logarithmic frequency scale).
(174) Interestingly, measuring these maxima and plotting the power spectrum provides for a visual image of stenosis in an artery. Indeed, we have determined that by plotting the power spectrum on the y axis and amplitude on the x-axis, we can effectively determine the percentage of stenosis in the carotid arteries of a patient.
(175) These maxima (generally two in number) are the main features in the frequency power spectrum at low frequencies generated by the wall pressure fluctuations when there is a constriction as compared to the situation of no constriction yet fully developed turbulence. In order to analyze this data, we have developed devices and invented several methodologies and processes that reduce or eliminate extraneous noise from the data samples, to enable further spectrum analysis downfield.
(176) The device eliminates noise in several ways. One by using sound barriers/dampening material to eliminate external noise as much as possible as well as noise caused by the patient moving; i.e. passive noise cancelling. We also eliminate or cancel ambient noise with active noise cancelling strategies, whether generating opposing waves or subtracting ambient noise; finally, we de-noise the received data by methodologies related to data processing using Wavelet, Welch's and Burg's methods. Ultimately, we plot peaks on a PSD and calculate stenosis of an area of interest in the arterial circulatory system through comparing these peaks on the PSD.
(177) Passive Noise Cancelling Strategies and Methodologies
(178) A first set of strategies includes mechanical strategies to eliminate or reduce noise. We can also consider these strategies to be passive noise cancelling strategies.
(179) For example, in preferred embodiments, the yoke 5, as depicted in
(180)
(181) In further embodiments, it is advantageous to utilize gel on the skin of a patient that assists in forming a temporary seal between the pad and the skin of the patient. Certain oil and water based gels or liquids are useful in assisting with the seal.
(182)
(183) The sensor pod itself, therefore, must also attenuate and block out some of the ambient noise. This can be achieved through several features that are depicted in detail above in
(184) A second adhesive 92 connects to the Printed Circuit Board 105, and several PCB contacts 106 contact the spring pins 111 on the PCB processor board 110 to make electronic connections. A processing unit 112 is defined on the bottom of the PCB processor board and comprises a battery, memory, and a processor. Alternatively, a battery may be centrally located, and the processing unit may be centrally located. The Piezo cap 100 contains a groove 101 to receive a quarter-turn locking feature 116 that is located on the PCB housing 115. This housing, like the PCB cap 100 attenuates and reduces ambient noise penetration to the piezo 90. A screw 113 secures the PCB housing 115 to a diaphragm bellows membrane 120, which allows movement of the entire sensor pod in directions in the lateral and longitudinal axis. Accordingly, when a device is placed against a surface, the sensor pod will be able to move away from the surface, or laterally to create a better fit towards the skin of the patient. Furthermore, this diaphragm bellows membrane 120, being non-rigid, will reduce the transfer of vibration and movement from a person holding a device containing the sensor pod, such as an array. A locking mechanism 121 secures the inner portion of the diaphragm bellows membrane 120 between the locking groove 117 and the locking cap 125.
(185) Accordingly, an embodiment of the disclosure comprises passive noise cancellation strategies comprising a sensor pod (features 85 and 86 together) comprising a disposable piezo cap 85, having a piezo 90, a Piezo cap 100 having noise attenuating properties, and a PCB house assembly 86 having a PCB board 110, a diaphragm bellows membrane 120, and a PCB housing 115. A locking feature on the PCB housing 115 connects to the Piezo cap 100 to secure them together. The rear of the PCB house assembly 86 comprises a diaphragm bellows membrane 120 that allows for movement of the components to isolate them from ambient noise and vibrations. The device may further comprise a noise attenuating material 219 disposed of around the sensor pad 18 to passively waves from the piezo sensor 90.
(186)
(187) Active Noise Cancelling Strategies and Methodologies
(188) In addition to the passive noise cancelling features of the sensor pod assembly, a further strategy for reducing noise to the piezo includes active cancellation of noise, such as found in the frequency chart of
(189)
(190)
(191)
(192)
(193)
(194) Cancellation and subtraction of sound can be accomplished in two ways. First, the sounds from the second piezo can be inversed and literally subtracted from the first piezo. Second, the sounds can be eliminated in analog by sending in a negative background signal which eliminates the sound. The prior art details several noise cancelling headphones, which use an external microphone to detect sound. This sound is then processed by a computing system with the device, and identifies and generates an out of phase sound, being out of phase by 180 degrees. This, when combined with the external sound, effectively cancels out the sounds that are received. Either method is functional, though the subtraction method may be preferable in certain embodiments.
(195)
(196)
(197)
(198) A particular method comprises a method of reducing noise to a sensor comprising: engaging a first sensor to a patient and a second sensor to ambient air, adjacent to said first sensor; detecting noises from said patient and simultaneously detecting noises from ambient air with said second sensor; subtracting the noise from said second sensor from the data from said first sensor, which will remove the ambient noise from the data from the first sensor.
(199) A particular method utilizes a phase change detected from a sensor to modify the sounds received at an adjacent sensor; a first sensor placed on a patient to detect sounds from the patient; a second sensor placed adjacent to said first sensor but shielded from the sounds of the patient; performing a phase change on the sounds received in said second sensor and emitting a proportional sound in said phase change.
(200) Analysis Based Noise Filtration Methods
(201) Active and passive cancellation can provide for a dramatic reduction in the amount of noise that ends up in a set of collected data. However, even with these background strategies to reduce and eliminate noise, detection of low frequency sounds can often be understood as looking at sounds that are “in the weeds.” Accordingly, further processing may be necessary, in certain embodiments, to collect data, amplify the data and perform certain analysis using a computer to clarify the data for best analysis.
(202) Spectrum analysis, also referred to as frequency domain analysis or Power Spectral Density (“PSD”) estimation, is the technical process of decomposing a complex signal into simpler parts. As described above, many physical processes are best described as a sum of many individual frequency components. Any process that quantifies the various amounts (e.g. amplitudes, powers, intensities, or phases), versus frequency can be called spectrum analysis.
(203) Spectrum analysis can be performed on the entire signal. Alternatively, a signal can be broken into short segments (sometimes called frames), and spectrum analysis may be applied to these individual segments. Periodic functions (such as sin(t) are particularly well-suited for this sub-division when t (time) includes several cycles. General mathematical techniques for analyzing non-periodic functions fall into the category of Fourier analysis.
(204) The Fourier transform of a function produces a frequency spectrum which contains all of the information about the original signal, but in a different form. This means that the original function can be completely reconstructed (synthesized) by an inverse Fourier transform. For perfect reconstruction, the spectrum analyzer must preserve both the amplitude and phase of each frequency component. These two pieces of information can be represented as a 2-dimensional vector, as a complex number, or as magnitude (amplitude) and phase in polar coordinates (i.e., as a phasor). A common technique in signal processing is to consider the squared amplitude, or power. In this case the resulting plot is referred to as a power spectrum.
(205) In practice, nearly all software and electronic devices that generate frequency spectra apply a Fast Fourier Transform (“FFT”), which is a specific mathematical approximation to the full integral solution. Formally stated, the FFT is a method for computing the discrete Fourier transform of a sampled signal.
(206) Because of reversibility, the FFT is called a representation of the function, in terms of frequency instead of time; thus, it is a frequency domain representation. Linear operations that could be performed in the time domain have counterparts that can often be performed more easily in the frequency domain. Frequency analysis also simplifies the understanding and interpretation of the effects of various time-domain operations, both linear and non-linear. For instance, only non-linear or time-variant operations can create new frequencies in the frequency spectrum.
(207) The Fourier transform of a stochastic (random) waveform (noise) is also random. Some kind of averaging is required in order to create a clear picture of the underlying frequency content (frequency distribution). Typically, the data is divided into time-segments of a chosen duration, where time is long enough to include several cycles of typical frequencies, and transforms are performed on each one. Then the magnitude or (usually) squared-magnitude components of the transforms are summed into an average transform. This is a very common operation performed on digitally sampled time-domain data, using the discrete Fourier transform. This type of processing is called Welch's method or Entropy Maximum (Burg) method. These methods are known and understood by a person of ordinary skill in the art. When the result is flat, it is commonly referred to as white noise. However, such processing techniques often reveal spectral content even among data which appear noisy in the time domain.
(208) Accordingly, by taking a piezoelectric unit, capable of measuring sounds and vibrations at low amplitude and within a particular frequency range, we can measure the wall pressure fluctuations due to stenosis. Accordingly, the sensitive piezoelectric devices combined with amplifiers are placed onto the skin above the carotid artery and the piezoelectric device detects these sounds. The detected sounds are then passed through analog to digital converters before reaching a computer in which further amplification and an analysis of the signal occurs.
(209) In the case of the arterial circulatory system, the piezo is placed on the skin above the artery in the region of a suspected stenosis. In the case of a carotid artery the placement would be on the neck, slightly below the ear. The particular placement of the piezo and the location of the stenosis is suggested by Fredberg and Borisyuk. Indeed, in an artery, between the stenosis and the region where turbulence has significantly decayed, the intensities can be rather large, where the wall can be subjected to large fluctuating stresses imposed by the turbulent blood flow. [Fredberg 1974] The distance over which this occurs is estimated to be about 12D downstream, where D is the normal diameter of the carotid artery. Borisyuk [2010]. For a typical internal carotid D of 0.5 cm, that distance would be of the order of several cm.
(210) Detection of vortices generated due to flow in the carotid artery produce low intensity sounds that are related to development of stenosis in an artery. These low intensity sounds are sometimes difficult to detect and to pull out of the mass of noise being generated by the body. Accordingly, a highly specialized detection device using piezo devices for arteries that are near the surface. In the relevant frequency range of 20 Hz to about 3000 Hz the wavelengths are long compare to other lengths, such as artery length or thickness of tissue between the artery and the skin. In this case the surface is still within the “near field” of a wave (much closer than one wavelength), in which case the tissue acts as an incompressible medium. The energy in the near field of a wave is attached to the source and cannot propagate away. Thus there is no net energy flux out from the source. Because near-field pressure fluctuations cannot propagate away, they are generally called “pseudo-sound”.
(211) Borisyuk [2010] has been able to relate the shape of the power spectrum at the surface to the vortex structures in the blood flow distal to a constriction. He divides the region distal to a constriction into three: Region I. The flow separation region, in which a jet flow of higher velocity, in the center, acts separately from the slower flow outside the jet. Region II. The flow reattachment region. The two regions, I and II, constitute the “most disturbed flow region”. The length of the first two regions, I plus II, based upon extensive calculations, Borisyuk estimates to be less than 7D, where D is the normal diameter of the artery. Here, stenosis may be detected in several different arteries in the arterial circulatory system. For example, detection may be directed towards detecting stenosis in the Internal Carotid Artery (ICA) in an adult, in which D is approximately 0.7 cm but the internal carotid is typically 0.5 cm. Therefore, the total length of the regions spoken of, I and II, would be at most about 3.5 cm. Region III is the region of flow stabilization where flow develops into the less turbulent flow farther upstream. This region extends from at most, 7D to 12D, or in the case of the ICA, at most from about 3.5 cm to about 6 cm.
(212) Conservation of fluid requires that v=V (D/d){circumflex over ( )}2. Let lower case v be the flow velocity inside the constriction and capital V the flow velocity past the constriction. Let d be the diameter of the flow inside the constriction. Borisyuk suggests estimates of two characteristic ring vortex frequencies. The first, f1, of vortices inside the jet, with typical size d; the second, f2, of vortices between the jet and the outer wall, with typical size, D.
(213) Accordingly, Borisyuk provides for a broad disclosure that certain structures in the blood generate flow patterns. Based on these flow patterns, and separated into three regions, Borisyuk estimates characteristics of vortex frequencies. However, these estimations provide only a rough estimate as to a vortex structure.
(214) Accordingly, our method for determining stenosis consists in connecting the frequencies associated with largest intensities in the spectral domain to three frequencies, f1 thru f2 in order to obtain estimates of percentage stenosis of the artery, (1−d/D)×100.
(215) The method has been implemented in a computer language we convert to binary, encrypted to be packaged as one whole product, software and hardware. The particular software used to run the data analysis can be determined by a person of ordinary skill in the art.
(216) A particular embodiment comprises the following steps: A sensor device is placed on a patient and data is sampled from the patient and the sound/vibrations are converted from analog to digital. The data is streamed from the device with both of the sensors in one data stream. We break the data stream down into two streams, one for the left sensor and one for the right. We then begin the Wavelet analysis which takes out noise. After the Wavelet removes the noise a power spectral density analysis is done and we are given a power spectral density (PSD). This tells us what frequency noise is found within the data and how strong/powerful the noise is. Because the PSD gives transient noise smoothing the PSD must be done to correctly identify the strongest peaks within the data. After smoothing is done peaks are determined and based on the where the peaks are will determine the amount of stenosis or whether no stenosis is present. If there is one peak, No stenosis is present. If there are two or more peaks the patient has some stenosis present.
(217) Wavelets have been frequently used in digital signal processing and are often known as small waves. A wavelet is a real-valued integral function ψ: R.fwdarw.R satisfying Zψ(t) dt=0. For practical applications, it has n vanishing moments: Z t pψ(t) dt=0, p=0, 1, . . . , n−1. Consider the following family of dilations and translations of the wavelet function ψ defined by ψjk(t)=2−j/2ψ(2−j t−k), j, k=0, ±1, ±2. The terms j and 2j are called the octave and the scale, respectively. By construction, this family consists of orthogonal basis functions in the sense that for a given time series or observed signal or simply data y(t), it can be written as the sum of these basis functions in a unique way: y(t)=X j X k djkψjk(t), where djk is the discrete wavelet transform (DWT) of y(t) given by djk=Z y(t)ψjk(t) dt, j, k=0, ±1, ±2. In practice, data is decomposed into its rough approximation at the chosen resolution level J (signal of interest) and details on a finite number of resolution levels j(≤J). The latter will be considered as noise. Denoising is equivalent to removing the details to allow for improved fit and prediction of peaks in a PSD plot.
(218) An Example of the Process for Calculation:
(219)
(220) We next perform a wavelet analysis 71, to de-noise the data by removing low-frequency components 1-60 or 1-70 Hz. After the wavelet analysis we generate a Power Spectral Density (PSD) 73 using the denoised data, in combination with Welch and/or Burg's method. From this PSD plot, we detect a first spike, typically between 75-250 Hz, (74) though it can go as low as 60 Hz. Where lower peaks are present, the Wavelet is re-run to remove a lower set of data, so that the first peak is not obfuscated.
(221) If a first spike is present between 75 and 250 Hz, we continue data acquisition (74). In certain embodiments, if there is no spike in this range, the sensor is adjusted (72) and the data acquisition process is re-started. Using this embodiment, we effectively build in a mechanism to ensure proper placement of the sensor, to make sure we have good quality data. However, other sounds may be utilized as a predetermined sound for ensuring proper placement in other embodiments.
(222) Once we have a first spike between 75 and 250 Hz, a second spike is analyzed (75), as different from the first and less than 3000 Hz. (feature 75). If the second spike is not found in this range, we declare stenosis at less than 25%. If the second spike is in this range, then we can calculate stenosis by peak comparison using the formula. We use the formula (1−f1/f2)×1100%, where f1 is the base frequency for the ring vortices in the artery (between 60 and 260 Hz) and f2 is the frequency from the restricted ring vortices (below 3000 Hz). If f1 is not present, the artery is too stenosed to show a base ring vortex and therefore we conclude there is a very high level of stenosis. If f2 is not present then we conclude that there is insufficient stenosis to create a restricted ring vortex and thus we say there is a very low level of stenosis. If neither f1 nor f2 are present, the patient is stenosed to the point where ring vortices can no longer form. This patient has extremely high stenosis and needs to see a specialist as soon as possible.
(223) Example of Data Analysis
(224) Read in data and look for extraordinary features. The step is important for reviewing if the device has followed protocol or not, and whether the subject has complied with the data acquisition procedures.
(225) The function CVRData provides a pop-up menu asking a user to select data, followed with a graph plotting channels, selected from Left—channel 1, Right—channel 2, or center—channel 3. One or all channels can be selected.
(226) The data of
(227) To select a channel to analyze, we look at the following aspects:
(228) Ch=1; note that Left or Ch=1, Right or Ch=2, and center or Ch-3.
(229) Setup of basic parameters for data analysis. Variable x is one of the channels in the following formula x=y(ch:3:length(y));
(230) Fs is the sampling rate, wherein Fs=20,000;
(231) One second record: the variable t is used for data visualization by plotting the first Fs or one second record of the channel values. Accordingly we can use the data:
(232) t=(0:Fs)/Fs; subplot(111), plot(x(1:10*Fs)), title (‘Ten second channel plot’)
(233) The resulting channel plot is depicted in
(234) A periodogram is generated. In general, one way of estimating the PSD of a process is to simply find the discrete-time Fourier transform of the samples of the process (usually done on a grid with an FFT) and appropriately scale the magnitude squared of the result. This estimate is called the periodogram.
(235) Periodogram(x, hamming(length(x)), length(x), Fs); xlabel(‘Frequency (Hz)’).
(236)
(237) The number of frequencies plotted is 1+half of length (x) and the unit is Hertz (Hz).
(238) Welch's Method can be used as an improved estimator of the PSD. Welch's Method, as known to a person of ordinary skill in the art, consists of dividing the time series data into (possibly overlapping) segments, computing a modified periodogram of each segment, and then averaging the PSD estimates. The result is Welch's PSD estimate.
(239) The averaging of modified periodograms tends to decrease the variance of the estimate relative to a single periodogram estimate of the entire data record. Although overlap between segments introduces redundant information, this effect is diminished by the use of a nonrectangular window, which reduces the importance or weight given to the end samples of segments (the samples that overlap).
(240) However, as mentioned above, the combined use of short data records and nonrectangular windows results in reduced resolution of the estimator. In summary, there is a tradeoff between variance reduction and resolution. Once can manipulate the parameters in Welch's method to obtain improved estimates relative to the periodogram, especially when the SNR is low. This is illustrated in the following example:
(241) A signal such as x consisting of the left channel data pwelch(x); which is graphically represented in
(242) The graph of
(243) Parameters to be specified with the Welch's method must be considered. The first parameter is the segment length. Default length is (x)/8. In code we use SGM=100,000. The next parameter is percent of overlaps: novoerpals=50,000.
(244) Through these elections we obtain Welch's overlapped segment averaging PSD estimate of the preceding signal. Use a segment length of 100,000 samples with 50 overlapped samples. Use 1+length(x)/2 DFT points so that 100 Hz falls directly on a DFT bin. Input the sample rate to output a vector of frequencies in Hz. We can plot the result.
(245) Example: [Pxx,F]=pwelch(x, sgm, noverlaps, [ ], Fs); plot (f, 10*log 10(Pxx)). The result is the plot of
(246) We can further estimate PSD through autoregressive PSD estimate through use of Burg's Method. Burg's Method is a parametric method for estimating PSD. Below returns a frequency vector, F, in cycles per unit time. The sampling frequency, Fs, is the number of sample per unit time. If the unit of time is seconds, then F is in cycles/second (Hz). For real-valued signals, F spans the interval [0,fs/2] when nfft is even and [0,fs/2] when nfft is odd.
(247) The following formula assumes an AR(50) model to the data.
(248) [Pxx,F]=pburg(x, 50, [ ], Fs); plot(F,10*log 10(Pxx)). The result is plotted in
(249) We use AR(50) because we tested model orders starting from 5 through 50 and determined that AR(50) provided the cleanest data result.
(250) Reflection Coefficients for Model Order Determination
(251) The reflection coefficients are the partial autocorrelation coefficients scaled by −1. The reflection coefficients indicate the time dependence between y(n) and y(n−k) after subtracting the prediction based on the intervening k−1 time steps.
(252) Use of arburg to determine the reflection coefficients. Use the reflecting coefficients to determine an appropriate AR model order for the process and obtain an estimate of the process PSD. We use the following formula:
[a,e,k]=arburg(x,50);
(253) Stem(k, ‘filled’); title(‘Reflection Coeficients’); xlabel(“model Order’)
(254)
(255) To find frequencies, we zoom into the data. Bf=0.1000/129:3876
(256) Plot(0:1000/129:3876, 10*log 10(Pxx(1:51)))
(257) Legend (‘pburg PSD Estimate’); x label (‘Frequency (Hz)’); y label (‘Power/frequency (dB/Hz)’); title (‘PSD before denoising’). The result is the data of
(258) We can then experiment with several choices of parameters in the Welch's PSD estimate, for example with 20 percent overlaps. Sgm=10,000; noverlaps=2000; [Pxx,F]=pwelch9x, sgm, noverlaps, [ ], Fs); plot(F, 10*log 10(Pxx)). This results in the plot of
(259) We can also test PSD by Welch with no overlaps:
(260) Sgm=10000; noverlaps=0; [Pxx,F]=pwelch(x, sgm, noverlaps, [ ], Fs);
(261) Plot(F,10*log 10(Pxx)); xlabel(‘Frequency (Hz)’); ylabel(‘Magnitude (dB)’); title (‘PSD before nenoising’). This results in the plot of
(262) If we zoon in the range of 2K Hz, with:
(263) Uf=2000; plot (F1:uf), 10*log 10(Pxx)1:uf))
(264) xlabel(‘Frequency (Hz)’); ylabel(‘Magnitude (dB)’); title (‘PSD before nenoising’). This results in the plot of
(265) Finally, we can output with frequencies, for peak analysis with [Pxx, F]=pburg(D1, 50, [ ], Fs0′ and zoom to within 2000 Hz (though 3000 would be good as well).
(266) Plot (0; 1000/129:1938, 10*log 10(Pxx(1:26))) grid on;
(267) Legend (‘pburg PSD estimate’)
(268) xlabel(‘Frequency (Hz)’); ylabel(‘Magnitude (dB/Hz)’); title (‘Parametric PSD after denoising). This results in the plot of
(269) We then allow the software to define the peaks. Once identified, the peaks can be used to calculate stenosis by (1−d/D)×100.
(270) Accordingly, we know that ambient noise is present in any data set and we know some of the sounds that are always present. Furthermore, we know the sounds that we are trying to detect and have determined that these sounds are at range 20-3000 Hz. We can remove other sounds introduced through these sensitive machines and concept is to provide a claim that covers the external and internal steps being applied to generate clean data.
(271) In certain embodiments, we determine stenosis based upon a class of stenosis. For example a first class may be less than 25% stenosis. A second class may be less than 50% stenosis, less than 70% stenosis, less than 90% stenosis. Accordingly, a method may be to calculate a binary response of less than or more than 25% stenosis. Another method may be to calculate a binary response of less than or more than 50% stenosis. Another method may be to calculate a binary response of less than 70% or less than 90% stenosis.
(272) Calculation of stenosis in such binary decision charts allows for a broad and quick determination of risk to a patient. Furthermore, certain procedures may be medically recommended at a certain stenosis percentage. Accordingly, for example, when testing the coronary artery, it may be necessary only to determine a binary decision of more or less than 50% stenosis, as procedures are recommended for surgical action once stenosis reaches such threshold.
(273) Utilizing the devices, systems, and methods as described above, the present components can be utilized in a system to identify large ring vortices from a fluid flow vessel. We can then analyze the signal utilizing low frequency (Spectral) methods and assess the range of stenosis, occlusion.
(274) In preparing for a test, the system first goes through a series of calibration steps, ensuring correct receipt of the signals, correlating the signals from the two carotid arteries and the heart, and identifying the systolic time, the period of most rapid fluid flow. Once the signal is recorded, the system prepares the data for processing the digital signal to conduct a spectral analysis. Using the signal features, a statistical analysis is performed against multiple parameters to render a classification of degree of stenosis, occlusion or aneurysm within each fluid flow vessel. For stenosis of the carotid artery, the output renders a report indicating a range of blockage against the defined Nascet categories with a representation of the systolic events.
(275) In accordance with one embodiment, the sensor array one or more sensors, which are positioned proximate the fluid flow vessel to be examined. In some instances the sensors are placed onto an array for determination of stenosis of the carotid artery. An array comprises two branches, which are biased inward and can be bent/flexed outward to the proper position. To accommodate bodies of differing heights, additional modifications can be made to allow for the adjustment of the lower sensor with respect to the upper sensors (e.g., providing a telescoping or otherwise extendable portion or arrangement in the lower branch and/or the upper two branches).
(276) A particular feature of the array is that it is adjustable and can be configured to account for the anatomical differences between individuals, while remaining sufficiently rigid to support the sensing elements. Furthermore, the shape and design of the array is particular important to assist with orienting sensing elements to each portion of the array, wherein sensing elements can easily be positioned adjacent to the neck for appropriate positioning to sense the carotid artery. At the same time, the materials and the angles utilized in the array provide appropriate resistance and a gentle force to compress the sensing element to the side of the neck for sensing. The shape and material thus provide an important feature to gently, but securely assist in positioning of the sensing elements and for testing patients for stenosis of the carotid artery.
(277) The array is adjustably designed to fit a majority of adults and to be held by the patient or a third person when performing a carotid artery test. In a preferred embodiment, the array, when placed on the patient, imparts sufficient pressure on the patient so as to achieve a measurement of sufficient quality to accurately determine stenosis, while limiting the pressure applied to the carotid artery. The goal is for there to be sufficient pressure to assist in positioning the sensing elements, and maintaining their position for about 2-3 minutes during a test, but not such pressure as to significantly impact the shape and size of the carotid artery being assessed. Indeed, as a whole, the array and the sensing elements are designed to be a passive test that is non-emitting, non-invasive, and is configured so that anyone can conduct the test without requiring certification.
(278) In accordance with one embodiment, the sensor elements in collaboration with the software or application running on a PC or main computing unit, takes three readings simultaneously from the right and left carotid arteries in the neck and from the heart just below the sternum, calibrates the sound signature, filters and then digitizes data for analysis. A shielded cable transmits the signals to the main computing unit. In further embodiments, signals and data can be transmitted via other transmission means, including wireless, Bluetooth, or other suitable data transmission mechanisms.
(279) Therefore, a method for determining stenosis of the carotid artery in a human patient consists of a first step of placing a sensing device onto the patient, wherein a first sensing element is placed adjacent to the carotid arteries; the sensing elements then measure sounds from the carotid artery. The sound is measured in analog and modified to digital format and then analyzed before a power spectral density analysis is performed. The power spectral density graph reveals peaks that are then analyzed to provide for a calculation of percent stenosis or occlusion of the carotid artery.