DEVICE AND METHOD FOR ESTIMATION OF PULMONARY FUNCTION CHARACTERISTICS

Abstract

The present invention relates to the field of medical devices, and, more particularly, to a portable system for testing one or more lung functions, and novel techniques for noninvasive determination of one or more pulmonary function characteristics.

Claims

1. A device for determining at least one pulmonary function characteristic, the device comprising (i) at least one chamber comprising at least one gas inlet port, configured to receive gas flow; a distal element comprising at least one distal opening; at least one pressure transducer configured to detect pressure within the at least one chamber and to produce a signal upon detecting pressure at about a predetermined threshold; (ii) at least one shutter module connected to said distal element, said at least one shutter module comprising at least one shutter opening; and (iii) an actuator configured to displace, or induce displacement, of the at least one shutter module upon receiving the signal, or a signal derived from said signal, wherein the at least one shutter opening is configured to allow passage of gas flow from the at least one distal opening to the atmosphere through a first total effective area enclosed by a perimeter of the at least one shutter opening, when the shutter module is in a state of first end position; and wherein upon said displacement the first total effective area changes to a second total effective area, when the shutter module is in second end position.

2. The device of claim 1, wherein the chamber further comprises a laminar flow diffuser.

3. The device of claim 1, further comprising a mouthpiece connected to the at least one gas inlet port, and configured to enable delivery of exhaled gas flow to the at least one gas inlet port.

4. The device of claim 1, further comprising a handle connected to the at least one chamber.

5. The device of claim 3, wherein the mouthpiece further comprises an anti-bacteriological filter.

6. The device of claim 1, wherein the actuator is a rotary actuator and the shutter module is a shutter wheel having a circular form, and wherein the actuator is configured to displace, or to induce displacement of, the at least one shutter wheel by rotation.

7. The device of claim 6, further comprising a shutter pin and a hard stop comprising a groove, wherein the shutter pin is located within the groove, and wherein the shutter pin is configured to limit the range of displacement of the shutter module.

8. The device of claim 1, further comprising a control module, wherein the control module is configured to receive one or more signals from the at least one pressure transducer.

9. The device of claim 8, wherein the control module is configured to provide the one or more signals to the actuator.

10. The device of claim 8, wherein the control module is configured to transmit the one or more signals to an external device.

11. The device of claim 9, wherein the control module is configured to transmit the one or more signals to the external device in real time.

12. The device of claim 8, wherein the control module is configured to store the one or more signals.

13. The device of claim 6, wherein the shutter wheel comprises at least one shutter opening, wherein the at least one shutter opening allows passage of gas flow from the at least one distal opening through the area enclosed by the perimeter of the at least one shutter opening, or through any portion thereof.

14. The device of claim 6, wherein said shutter wheel and said distal element are detached from one another allowing passage of gas flow from the at least one distal opening through an area extending between the circumference of said shutter wheel and said distal element.

15. The device of claim 6, wherein the shutter wheel comprises at least one shutter edge opening, wherein the at least one shutter edge opening allows passage of gas flow from the at least one distal opening through the area enclosed by the perimeter of the at least one shutter edge opening, or through any portion thereof.

16. The device of claim 2, wherein the at least one pressure transducer is configured to detect pressure within the space between the laminar flow diffuser and the distal element.

17. A method for evaluating at least one pulmonary function characteristic in a subject in need thereof, comprising the steps of: (i) providing the device of claim 1; (ii) receiving gas exhalation from a subject in need thereof, through the gas inlet port; (iii) detecting pressure within the at least one chamber; (iv) producing pressure versus time curve and identifying a first peak pressure; (v) toggling the shutter module if the detected pressure post the first peak pressure is within the range of a predetermined pressure threshold; (vi) identifying a second peak pressure and an end point pressure in the pressure versus time curve; and (vii) deriving from said pressure versus time curve at least one value related to at least one pulmonary function characteristics of said subject in need thereof.

18. The method of claim 17, further comprising the step of identifying on the curve of step (vi) a pressure volume indicator point, wherein the pressure volume indicator point post the second peak pressure is a pressure within the range of a predetermined pressure threshold.

19. The method of claim 17 wherein said at least one value comprises any one or more of Thoracic Gas Volume (TGV) and Total Lung Capacity (TLC), and wherein the derivative of said pressure versus time curve includes at least one linear portion.

20. The method of claim 19, further comprising identifying on the curve of step (vi) a pressure volume indicator point, wherein the at least one linear portion is the derivative of the pressure during a time period encompassing the pressure volume indicator point.

21. The method of claim 17, wherein said at least one value comprises any one or more of TGV and TLC, and wherein said pressure versus time curve includes at least one exponential portion.

22. The method of claim 21, wherein the at least one exponential portion is during a time period extended between the second peak pressure and the end point pressure.

23. The method of claim 17, wherein said at least one value comprises any one or more of TGV and TLC, and wherein said pressure versus time curve includes at least one parabolic portion.

24. The method of claim 23, wherein the at least one parabolic portion is the pressure during a time period extended between the second peak pressure and the end point pressure.

25. The method of claim 18, further comprising the step of determining an airway resistance based on the predetermining pressure threshold and the second peak pressure.

26. The method of claim 25, further comprising the step of determining a lung compliance based on the airway resistance and the end point pressure.

27. The method of claim 18, further comprising determining a group of volume related indices; and producing a diagnostic event based on the group of volume related indices.

Description

BRIEF DESCRIPTION OF THE FIGURES

[0050] Some embodiments of the invention are described herein with reference to the accompanying figures. The description, together with the figures, makes apparent to a person having ordinary skill in the art how some embodiments may be practiced. The figures are for the purpose of illustrative description and no attempt is made to show structural details of an embodiment in more detail than is necessary for a fundamental understanding of the invention. For the sake of clarity, some objects depicted in the figures are not to scale.

[0051] In the Figures:

[0052] FIG. 1a constitutes a schematic diagram of a pulmonary function test device, in accordance with some embodiments, in a state of a first end position.

[0053] FIG. 1b constitutes a schematic diagram of a pulmonary function test device, in accordance with some embodiments, in a state of a second end position

[0054] FIG. 2 constitutes a view in perspective of a pulmonary function test device, in accordance with some embodiments.

[0055] FIG. 3 constitutes a cross-sectional view of the pulmonary function test device that is illustrated in FIG. 2.

[0056] FIG. 4 constitutes a partial cross-section view A that is marked in FIG. 3.

[0057] FIG. 5a constitutes a perspective front view of a shutter assembly, in accordance with some embodiments.

[0058] FIG. 5b constitutes a perspective back view of a shutter assembly, in accordance with some embodiments.

[0059] FIG. 6a constitutes a back view of a shutter assembly, in accordance with some embodiments, in a state of a first end position.

[0060] FIG. 6b constitutes a back view of a shutter assembly, in accordance with some embodiments, in a state of a second end position.

[0061] FIG. 7a constitutes a perspective front view of a shutter assembly, in accordance with further embodiments.

[0062] FIG. 7b constitutes a perspective back view of a shutter assembly, in accordance with further embodiments.

[0063] FIG. 8 shows an example of a graph of gauge pressure reading over time during a sequential forced gas exhalation event, in accordance with some embodiments.

[0064] FIG. 9 shows a flowchart of the system operation mode, in accordance with some embodiments.

[0065] FIG. 10 shows an example of a Bland-Altman plot of Total Lung Capacity values measured by a plethysmograph and a pulmonary function test device.

[0066] FIG. 11a shows a Total Lung Capacity histogram measured for a healthy subject, during a repeatability trial.

[0067] FIG. 11b shows a Total Lung Capacity histogram measured for a patient suffering from chronical asthma disease, during a repeatability trial.

[0068] FIG. 12a shows an example of the graph of gauge pressure reading over time during a sequential forced exhalation event, measured for healthy patients, in accordance with some embodiments.

[0069] FIG. 12b depicts an example of the graph of gauge pressure reading over time during a sequential forced exhalation event, measured for patients suffering from an obstructive lung disease, in accordance with some embodiments.

[0070] FIG. 12c depicts an example of the graph of gauge pressure reading over time during a sequential forced exhalation event, measured for patients suffering from a restrictive lung disease, in accordance with some embodiments.

DETAILED DESCRIPTION OF SOME EMBODIMENTS

[0071] In the following description, various aspects of the disclosure will be described. For the purpose of explanation, specific configurations and details are set forth in order to provide a thorough understanding of the different aspects of the disclosure. However, it will also be apparent to one skilled in the art that the disclosure may be practiced without specific details being presented herein. Furthermore, well-known features may be omitted or simplified in order not to obscure the disclosure. In the figures, like reference numerals refer to like parts throughout.

[0072] As used herein the terms subject and patient are interchangeable, and refer to the user of the devices and methods disclosed herein and include a healthy user, in the context of the measured values or a subject having, or being susceptible to have, a pulmonary disease or disorder. The subject may be a person or a mammal.

[0073] As used herein, the terms about or within the range of mean in the range of, roughly, or around. In general, the terms about or within the range of are used to modify a numerical value above and below the stated value by 20%. According to some embodiments, the term about or within the range of are used to modify a numerical value above and below the stated value by 15% thereof. According to some embodiments, the term about or within the range of are used to modify a numerical value above and below the stated value by 10% thereof.

[0074] Reference is now made to FIGS. 1a-1b. FIG. 1a constitutes a schematic diagram of a pulmonary function test device 100 in accordance with some embodiments, in a state of a first end position. FIG. 1b constitutes a schematic diagram of a pulmonary function test device 100 in accordance with some embodiments, in a state of a second end position. The device 100 comprises a main chamber 110 having an absolute pressure P.sub.M, and volume V.sub.M, a proximal end 114, a distal end 140, a pressure transducer 126 and shutter module 150. The pressure transducer 126 measures the pressure within main chamber 110. A mouthpiece 130 which may include an anti-bacteriological filter 132 is mounted on the gas inlet port 120 located at the proximal end 114 of the main chamber 110, in accordance with some embodiments. Shutter module 150 may either be in a state of first end position (see FIG. 1a) or in a state of second end position (see FIG. 1b). In FIG. 1a, in accordance with some embodiments, exhalation gas flow 190 flows from mouthpiece 130 through chamber 110 then through shutter module 150 to the atmosphere, having atmospheric pressure P.sub.A, via shutter opening 152, having a first total effective area A.sub.S1, when shutter module 150 is in first end position. In FIG. 1b, in accordance with some embodiments, exhalation gas flow 190 flows from mouthpiece 130 through chamber 110 then through shutter module 150 to the atmosphere via shutter opening 154, having a second total effective area A.sub.S2, when shutter module 150 is in second end position.

[0075] The terms flow, gas flow, exhalation gas flow and flow rate as used herein are interchangeable.

[0076] Reference is now made to FIGS. 2-4. FIG. 2 constitutes a view in perspective of a pulmonary function test device 200 in accordance with some embodiments. FIG. 3 constitutes a cross-sectional view of the pulmonary function test device 200 in accordance with some embodiments. FIG. 4 constitutes a partial cross-section view 300 that is marked in FIG. 3 in accordance with some embodiments. Mouthpiece 230 may be mounted on the gas inlet port 230 (see FIG. 3) of the main chamber 210, serving as an interface between a subject and the pulmonary function test device 200, in accordance with some embodiments. Mouthpiece 230 is disposable, in accordance with some embodiments. Mouthpiece 230 comprises an anti-bacteriological filter 232, in accordance with some embodiments. A handle 224 may be configured to manually holding the device 200 while being used, in accordance with some embodiments. A shutter assembly may comprise a distal element 240, a shutter module 250, an extension spring 280, a shutter break 270, also termed hereinafter a hard stop, and a shutter pin 276 (see FIG. 4).

[0077] In accordance with some embodiments, the shutter module 250 may be in the form of a shutter wheel (see FIG. 2).

[0078] In accordance with some embodiments, the device is portable. In accordance with some embodiments, the device is hand held.

[0079] An actuator holder 268 may be mounted on the distal element 240 of the shutter assembly, in accordance with some embodiments. The shutter assembly may be mounted on the main chamber 210, in accordance with some embodiments. An actuator 260 may be mounted on the actuator holder 268.

[0080] The term mounted on as used herein, includes, but is not limited to coupled, assembled, bound, and/or attached.

[0081] In accordance with some embodiments, the actuator holder 268 is a rotary actuator holder, and the actuator 260 is a rotary actuator, wherein it may rotate the shutter wheel 250 when actuated or toggled.

[0082] Pressure transducer 226 is installed on main chamber 210 and may be connected via an orifice 228 to the internal space of a laminar flow diffuser 234 (see FIG. 3 and FIG. 4), in accordance with some embodiments.

[0083] Without being bound by any theory or mechanism, the inclusion of a laminar flow diffuser 234 in the main chamber 210, may create a laminar exhalation flow. Laminar gas exhalation, also termed herein flow, may improve the accuracy of the forced vital capacity (FVC) calculation.

[0084] It may be advantageous to use a pressure transducer 226 with a pressure resolution of at least 2 Pa, configured to operate in a temperature range of 0 C. to 50 C. and with a response time of at least 0.001 second, such as the following absolute pressure sensors: SDX Series by Honeywell, MPS-301A by JL World, SCP1000 Series by VTI Technologies or XP-6000CA by Epson Toyocom.

[0085] In accordance with some embodiments, gauge pressure transducers may be applied. However, in such cases a special algorithm should be used for calculation of the absolute pressure level. Examples of gauge pressure transducers may include 40PC by Honeywell, PX40 Series by Omega, and MPXV7007GC64 or MPXV5010GC6T1 by Freescale Semiconductor.

[0086] In accordance with some embodiments, the pressure transducer 226 may be enclosed within the handle 224.

[0087] In accordance with some embodiments, an electronic control module 222 may be enclosed within handle 224 (see FIG. 3). The control module 222 may be configured to supply one or more signals to a rotary actuator 260, in accordance with some embodiments, wherein said one or more signals will enable switching of the shutter module 250 from a first end position, wherein gas flows from the main chamber 210 to the atmosphere via shutter opening 152 having a total effective area A.sub.S1, to a second end position wherein gas flows from the main chamber 210 to the atmosphere via shutter opening 154 having a total effective area A.sub.S2.

[0088] The terms total effective area A.sub.S1 and first total effective area A.sub.S1 as used herein are interchangeable.

[0089] The terms total effective area A.sub.S2 and second total effective area A.sub.S2 as used herein are interchangeable.

[0090] In accordance with some embodiments, the control module 222 is configured to receive the one or more signals of the pressure transducer 226 and either store said one or more signals, transmit it to any external device, either by hardline or wireless connection, or perform both storage and transmission, wherein data transmission may occur either in real-time or post storage.

[0091] In accordance with some embodiments, the control module 222 may process the one or more input signals of the pressure transducer 226, wherein said processed signals are either stored, transmitted to any external device, either by hardline or wireless connection, or both stored and transmitted, wherein data transmission may occur either in real-time or post storage.

[0092] In accordance with some embodiments, the handle 224 is designed in its size and shape so that the device may be handheld comfortably.

[0093] In accordance with some embodiments, the orifice 228 may be located distal to the laminar flow diffuser 234, such that it is measuring the pressure in a space caged between the laminar flow diffuser 234 and the distal element 240.

[0094] In accordance with some embodiments, main chamber 210 will not include a laminar flow diffuser 234, wherein pressure transducer 226 may measure pressure directly in the space caged within the main chamber 210.

[0095] The rotary actuator 260 may rotate the shutter wheel 250 via the coupling 236 and shaft 238. A hard stop 270 includes a groove, in which a shutter pin 276 is allowed to maneuver within its boundaries, and which is mounted on the shutter wheel 250. The shaft 238 may have a bearing unit 266, which is mounted on distal element 240 (see FIG. 4), in accordance with some embodiments. The bearing protection cover 264 protects the bearing unit 266 from being exposed to the exhaled gas flow.

[0096] In accordance with some embodiments, the rotary actuator 260 should overcome the torque M.sub.RA that counteracts the rotation of the shutter wheel 250, according to the following formula:

M.sub.RA=M.sub.FF+M.sub.RS+M.sub.AF

where M.sub.FF is the friction torque between the shutter wheel 250 and the distal element 240, M.sub.RS is the torque of the extension spring 280, and M.sub.AF is the moment of inertia of the shutter wheel 250.

[0097] It may be advantageous to use one of the following two types of the rotary actuators: stepper motor or rotary solenoid. An example of a suitable stepper motor is model NEMA 14HS10-0404S by STEPPERONLINE. Examples of suitable rotary solenoids are model M341-30-180-R by GEEPLUS, and model 3EVM by LEDEX.

[0098] Reference is now made to FIGS. 5a-b. FIG. 5a constitutes a perspective front view (facing the gas inlet port 220) of a shutter assembly 400, while FIG. 5b constitutes a perspective back view (facing the actuator 260) of same shutter assembly 400. The distal element 440 may have, in the example illustrated in FIG. 5a, six distal openings 446 in the form of straight slots. The shutter wheel 450 may have the same number (six) of shutter openings 456 shaped in the form of a T (see FIG. 5b).

[0099] It is to be understood that the number of distal openings 446 on the distal element 440 or shutter openings 456 on the shutter wheel 450 may vary.

[0100] In accordance with some embodiments, the device comprises at least one distal opening 446.

[0101] In accordance with some embodiments, the device comprises at least one shutter opening 456.

[0102] In accordance with some embodiments, the number of distal openings 546 and the number of shutter opening 456 is identical.

[0103] In accordance with some embodiments, the number of distal openings 446 is different from the number of shutter opening 456.

[0104] A hard stop 470 includes a groove, in which a shutter pin 426 is allowed to maneuver within its boundaries, and which is mounted on the shutter wheel 450. The movement of shutter pin 426 is bounded between the groove's opposite sidewalls 472 and 474 (see FIG. 5b), which provide the limitation of the shutter wheel's 450 angle of rotation. In a first state of the shutter wheel, when the rotary actuator 260 is not toggled, the extension spring 480 exerts force on the shutter wheel 450, so as to rotate it counterclockwise until it reaches a stop position (first end position) due to the shutter pin's 426 inability to move beyond the hard stop's first edge 472.

[0105] When the shutter wheel 450 is positioned in a first end position, a first gas flow-flows from the distal openings 456 through the shutter openings 446 to the atmosphere.

[0106] Reference is now made to FIGS. 6a-b. FIG. 6a constitutes a back view (facing the actuator 260) of a shutter assembly 400 in a state of a first end position. FIG. 6b constitutes a back view of a shutter assembly 400 in a state of a second end position. In the first end position illustrated in FIG. 6a, six shaded areas 452, together forming a total effective area A.sub.S1, indicate the area through which passage of exhaled gas flow occurs in this state. When the rotary actuator 260 is actuated or toggled, it rotates the shaft 438 with the shutter wheel 450 until the shutter pin 426 is in contact with the hard stop's second edge 474, reaching a second end position illustrated in FIG. 6b. The six shaded areas 454 form together a total effective area A.sub.S2. The total effective area A.sub.S2 indicates the area through which a passage of exhaled air occurs and is also termed hereinafter a second gas flow. The first gas flow and the second gas flow, corresponding to effective areas A.sub.S1 and A.sub.S2, respectively, are affecting the device's resistance to gas flow exhalation.

[0107] In accordance with some embodiments, the distal element 440 and shutter wheel 450 may comprise distal openings 446 in the shape of straight slots and shutter openings 446 in the form of T shaped windows, respectively (see FIGS. 5a-b).

[0108] In accordance with some embodiments, distal openings 446 and shutter openings 456 may have any shape or size, wherein the resulting areas A.sub.S1 and A.sub.S2 are different. In a preferred embodiment, distal openings 446 and shutter openings 456 are configured so that A.sub.S1>A.sub.S2, that is to say that the flow resistance in the first end position, in which the total effective air channel area is A.sub.S1, is lower than the flow resistance in the second end position, in which the total effective air channel area is A.sub.S2.

[0109] Reference is now made to FIGS. 7a-b. FIG. 7a constitutes a perspective front view (facing the gas inlet port 220) of a shutter assembly 500, while FIG. 7b constitutes a perspective back view (facing the actuator 260) of same shutter assembly 500, in accordance with some embodiments. The distal element 540 may have, in the example illustrated in FIG. 7a, eight oval distal openings 546. The shutter wheel 550 may have the same number (eight) of the oval shutter openings 546 (see FIG. 7b), which are essentially identical, meaning that they are of the same size and shape.

[0110] It is to be understood that the number of distal openings 546 on the distal element 540 or shutter openings 556 on the shutter wheel 550 may vary.

[0111] In accordance with some embodiments, the device comprises at least one distal opening 546.

[0112] In accordance with some embodiments, the device comprises at least one shutter opening 556.

[0113] In a first state of the shutter wheel 550, when the rotary actuator 260 is not actuated nor toggled, the extension spring 580 exerts force on the shutter wheel 550, so as to rotate it clockwise until it reaches a stop position (first end position) due to the shutter pin's 576 inability to move beyond the hard stop's 558 second edge 574. In this position, distal openings 546 are positioned opposite to the shutter openings 556. In the first state the exhaled gas flow (first gas flow) flows through the shutter openings 556, thereby forming effective area A.sub.S1.

[0114] When the rotary actuator 260 is actuated or toggled, it rotates the shaft 538 with the shutter wheel 550 until the shutter pin 576 is in contact with the hard stop's first edge 572, reaching a second end position. In accordance with some embodiments, at least one of the shutter openings 556 is situated opposite to the smooth distal element wall 548 (see FIG. 7a). In the second end position the exhaled gas flow (second gas flow) flows through a portion of shutter openings 556 positioned opposite to distal openings 546, thereby forming effective area A.sub.S2.

[0115] In accordance with some embodiments, the device comprises an open gap (not shown) between the shutter wheel 550 and distal element 540, configured to allow gas flow from the distal opening 546 through the circumference of said gap, both at first end position and second end position. In accordance with some embodiments, the first state the exhaled gas flow (first gas flow) flows through the shutter openings 556, and through the open gap between shutter module 550 and the distal element 540, thereby forming effective area A.sub.S1.

[0116] In accordance with some embodiments, the number and position of distal openings 546 is configured so that in the second end position, some of the shutter openings 556 are situated opposite to the smooth distal element's wall 548, while the remaining shutter openings are positioned against opposing distal openings 546, in a manner that the perimeters of shutter openings556 through which passage of exhaled gas occurs in this state together with the perimeter of the shutter wheel 550 multiplied by the depth of the gap between shutter wheel 550 and distal element 540 (forming the total circumferential area of said gap), are forming together the total effective area A.sub.S2.

[0117] In accordance with some embodiments, all of the shutter openings 556 are situated opposite to the smooth distal element's wall 548 in the second end position, such that the effective area A.sub.S2 is only the area of the perimeter of shutter wheel 550 multiplied by the depth of the gap between shutter wheel 550 and distal element 540 (forming the total circumferential area of said open gap).

[0118] It is to be understood that at all states passage of airflow is available. During second end position, the passage is restricted in comparison to the passage during the first end position, such that only a portion of shutter openings is positioned opposite to distal openings. This relative restriction is also referred to as relative limitation.

[0119] In accordance with some embodiments, the device further comprises at least one shutter edge opening (not shown) located on the shutter module edge 520, wherein the at least one shutter edge opening is configured to allow gas flow from the distal opening 546 through the circumference of the at least one shutter edge opening, both at first end position and second end position.

[0120] It should be understood that A.sub.S2 is always a positive value, meaning that gas flow to the atmosphere must occur at second end position, either through at least one shutter opening 556, the gap between shutter wheel 550 and distal element 540, at least shutter edge opening, or any combination thereof.

[0121] In accordance with some embodiments, the area A.sub.S2 is in the range of 20% and 40% of the area A.sub.S1. It is to be understood that the relationship between A.sub.S2 and A.sub.S1 is not limited to the aforementioned ratios. In fact, the relationship between A.sub.S2 and A.sub.S1 may vary.

[0122] In accordance with some embodiments, as exemplified in FIGS. 5a-b and 6a-b, extension spring 480 is located such that it exerts force on the shutter wheel 450, which rotates it counterclockwise until shutter pin 426 is in contact with the hard stop first edge 472, designating such position as first end position, while the other position, in which the shutter pin 426 is in contact with the hard stop second edge 474 is designated as second end position.

[0123] In accordance with some embodiments, as exemplified in FIGS. 7a-b, extension spring 580 is located such that it exerts force on the shutter wheel 550, which rotates it clockwise until shutter pin 576 is in contact with the hard stop's second edge 574, designating such position as first end position, while the other position, in which the shutter pin 576 is in contact with the hard stop first edge 572 is designated as second end position.

[0124] It should be noted that additional technical solutions may be used by those skilled in the art to improve and extend some features of the device. For example different shutter module designs and displacement mechanisms, which may or may not be formed as circular elements, can be utilized to achieve different gas outlet, namely, different effective areas or resistances, between first end position and second end position of said shutter module.

[0125] The term displacement as used herein refer to a reposition, or change in position of the shutter module. The change is typically movement of the shutter module between two locations, where the distance between the position prior to displacement and the position post displacement has a range.

[0126] The terms displace, displacing, toggle and toggling are interchangeable and refer to the operation that induces, or causes, displacement.

[0127] A control system (not shown) of the pulmonary function test device is configured to perform at least the following two main functions: data acquisition and data processing. The control system may be situated within the device, for example as part of the control module 222 (see FIG. 3), or as part of an external computer system, for example as a personal computer or laptop, wherein external computer system is in communication (either hardline or wireless) with the device. Both functions may be controlled by a computer system and a microprocessor firmware in communication with each other. During data acquisition, shutter module 250 is displaced due to one or more signals derived from patterns of acquired pressure signals, which in turn change the resistance to exhaled gas flow. An operator may initiate the operation process by providing an initiation command via software tool installed on a computer system, in accordance with some embodiments. A computer system may send commands, which may be related to initialization parameters, to an internal microprocessor through wireless (e.g. Wi-Fi, Bluetooth, near field communication, and the like) or hardline (e.g. serial, parallel, universal serial bus, and the like) communication means, which initiate data acquisition and data processing.

[0128] A signal from the pressure transducer may be digitized, for example, by 12 bits A/D converter with a predefined sample rate, and transmitted back to a computer system over, for example, USB protocol during a predefined time period. At the same time, the microprocessor processes the acquired pressure signal patterns and determines peak values and signal slope directions. When a signal's first peak followed by a signal's decay are determined, and the acquired values fall below a predefined threshold value, the microprocessor triggers an output signal to toggle shutter module displacement, changing the total effective area as cited above. A computer system may receive acquired data of the pressure signals in a continuous manner during the complete duration of the measurement process.

[0129] Data processing may be performed once all pressure data, composed of all acquired pressure signals during a measurement process, is transmitted to a computer system. The processing may contain, but is not limited to, the following steps: a digital filtration procedure, conversion of acquired pressure signals to absolute pressure values based on a pressure transducer calibration data, determination of the signal base flow, and determination of a start and finish points of a forced gas exhalation process. Additional steps of the processing may include, yet are not limited to, determination of the first and second pressure signal peaks, running algorithms and mathematical calculations, including such that are based on related values approximations and logical assumptions, as will be presented here forth, and a first and second peak of pressure signals.

[0130] The main clinical objectives of respiratory function tests are to provide diagnosis, assessment of severity, treatment monitoring and evaluation of prognosis. Spirometry measures vital capacity and force expiratory volume in one second. These parameters permit differentiation between restrictive and obstructive respiratory diseases. The spirogram flow-volume curves are used for diagnostics.

[0131] Reference is now made to FIG. 8, which shows an example of a graph of gauge pressure reading over time during a sequential forced gas exhalation event, in accordance with some embodiments. The time t.sub.1 represents the starting point of the forced gas exhalation, and time t.sub.5 represents the end of said forced gas exhalation. The pressure at t.sub.5 is P.sub.5, also termed end point pressure. In accordance with some embodiments, t.sub.5 at which the end point pressure is identified, is based on estimations as further detailed herein. Proper determination of the duration of the forced exhalation, i.e. the time between t.sub.1 and t.sub.5, may result in more reliable estimation of total lung capacity, as will be described here forth. In most conventional spirometry tests, the back extrapolation approach, which is well known to those skilled in the art, is used for determination of a start point t.sub.1. This approach may be used for the pulmonary function test device in accordance with some embodiments, wherein the start of the process is set to be at point (t.sub.1,P.sub.M1; see FIG. 8). A second approach for determination of the start point (t.sub.1,P.sub.M1) may be the pressure threshold method, according to which:

P.sub.M1=P.sub.A+P.sub.1

[0132] where P.sub.A is the atmospheric pressure, and P.sub.1 is a predetermined pressure threshold, which in accordance with some embodiments may be set as P.sub.1=50 Pa. Such a method for determining the start of the forced gas exhalation may be prone to errors associated with the use of an arbitrarily predetermined threshold. Nevertheless, the approach is known in the art and may be considered acceptable with a sufficient accuracy despite said errors. Other methods for determination of the start point (t.sub.1,P.sub.M1) may be used in the context of this disclosure.

[0133] The second point to be evaluated on the graph of pressure vs. time (see FIG. 8) is the point of first pressure peak (t.sub.M,P.sub.MM). The first pressure peak (t.sub.M,P.sub.MM) may be determined by to any known method in the art for determining a peak maximum in a graph, implemented on pressure data acquired during the operation of the device.

[0134] When pressure drops to (t.sub.2,P.sub.M2; see FIG. 8), an actuation signal, also referred to as a toggling signal, is transmitted to initiate the shutter module's displacement from first end position to second end position, wherein upon displacement the first flow, also referred to as flow rate, of exhaled gas, corresponding to the effective area A.sub.S1, changes to a second flow rate of exhaled gas, corresponding to the effective area A.sub.S2. The pressure point P.sub.M2 may depend on the values of effective areas A.sub.S1 and A.sub.S2, and may be roughly estimated, in accordance with some embodiments, by the following formulas:

P.sub.M2=P.sub.A+P.sub.2

where P.sub.M2 may be calculated as:

[00001]$.Math..Math.P_{M.Math..Math.2}=\frac{0.01.Math.P_A.Math.A_{S.Math..Math.2}}{A_{S.Math..Math.1}}$

[0135] The point (t.sub.3, P.sub.M3) represents a second pressure peak in the region between point (t.sub.2, P.sub.M2) and pressure volume indicator point (t.sub.4, P.sub.M4). Determination of this point may be performed by known methods in the art. This point may be used for calculating gas flow resistance and volumetric characteristics.

[0136] In pressure volume indicator point (t.sub.4, P.sub.M4) the pressure may be equal to the pressure in the point (t.sub.2, P.sub.M2), that is P.sub.M4=P.sub.M2. This point may be used for calculating volumetric characteristics.

[0137] The point (t.sub.5, P.sub.M5) represents the end of forced gas exhalation. This point may be determined, in accordance with some embodiments, by the deviation method, whereby forced exhalation is finished when the mouth pressure deviation is less than or equal to {dot over (P)}.sub.M5F, which may be estimated according to the following formula:

[00002]${\stackrel{.}{P}}_{M.Math..Math.5.Math.F}=\frac{0.1.Math.\left(P_{M.Math..Math.4}-P_{M.Math..Math.3}\right)}{t_4-t_3}$

[0138] In accordance with some embodiments, assuming that a maximal duration of forced or spontaneous gas exhalation for healthy subjects may be in the range of 4-6 seconds; and for subjects with airways obstruction or subject of older age, maximal duration of gas exhalation may typically be in the range of 6-8 seconds, it may be acceptable to assume a maximal recording time of pressure data as T.sub.REC=15-20 s. In accordance with some embodiments, the sampling rate of the pressure signals may be at least 1000 Hz (such that the sampling time is not higher than t.sub.S=0.001 s). In accordance with some preferred embodiments, sampling rate may be at least 5000 Hz (such that the sampling time is not higher than t.sub.S=0.0002 s), so that the minimal size of the memory module of the pressure transducer recording can be calculated according to the formula:

N.sub.M=T.sub.REC/t.sub.S

[0139] For embodiments where the sampling rate is 5000 Hz, minimal required size for the memory module is N.sub.M=T.sub.REC/t.sub.S=100,000 points.

[0140] Reference is now made to FIG. 9, which shows a flowchart of the system's operation mode, in accordance with some embodiments. First, the operating cycle of the device is started (step 910), which may be triggered, in accordance with some embodiments, by an operator of the system. From that moment, pressure P.sub.M within main chamber 210, measured by the pressure transducer 226, is stored (step 920) either in a memory module within the device, or in a memory module of a computer system, after being transmitted (either by hardwire or wireless transmission) thereto. The shutter module 250, from the beginning of the measurement process, and while not being actuated or toggled, is in a first end position, wherein exhaled gas flows to the atmosphere through first passage, corresponding to the effective area A.sub.S1 (step 930). Data processing identifies when first peak pressure P.sub.MM has been reached, after which the control system checks whether the pressure reached the level P.sub.M2 (step 940) or is within the range of that level. As long as the pressure has not reached the first peak P.sub.MM, or has not reached the level P.sub.M2 after reaching P.sub.MM, or has not reached pressures within the range of those levels, shutter module 250 remains in first end position. At the moment the pressure reaches P.sub.M2, or is within the range of that level, after reaching first peak pressure P.sub.MM, a signal may be sent to the actuator 260, such that, as a result the actuator 260 displaces, or induce displacement of, shutter module 250 to a second end position, wherein exhaled gas flows to the atmosphere through second passage, corresponding to the effective area A.sub.S2 (step 950). From that point, data processing identifies when second peak pressure P.sub.M3 has been reached, after which the control system checks whether the pressure is within the range of the level P.sub.M5 (step 960). As long as pressure has not reached the second peak P.sub.M3, or has not reached the level P.sub.M5 after reaching P.sub.M3, shutter module 250 remains in second end position. At the moment that the pressure reaches the level P.sub.M5, after reaching second peak pressure P.sub.M3, pressure data storage stops (step 970), and the operation cycle of the system ends (step 980), which may include the return of shutter module 250 to a state of first end position.

[0141] The terms threshold and level as used herein are interchangeable.

[0142] Typical parameters determined by spirometry include, but are not limited to, any one or more of the following: [0143] 1) Spirometric Characteristics defined by forced expiratory flow parameters (spirogram). The main spirometric parameters are: Forced vital capacity (FVC); volume of gas exhaled in the first one second of exhalation (FEV1); forced expiratory flow between 25 and 75 percent of FVC (FEF25-75) i.e. the average expired flow over the middle half of FVC manoeuver; and peak expiratory flow (PEF). [0144] 2) Lung Volumetric Parameters: Total lung capacity (TLC); residual volume (RV) i.e. the volume of gas remaining in the lung after maximal exhalation; and thoracic gas volume (TGV) i.e., the absolute volume of gas in thorax at any point in time and any level of alveolar pressure. [0145] 3) Airway Resistance defined as the ratio of driving pressure to the rate of air flow. Resistance to flow in airways depends on whether the flow (laminar or turbulent), on the dimensions of an airway and on a viscosity of a gas. Total resistance to air flow includes three main components: (a) inertia of respiratory system (contributes negligibly to total resistance), (b) tissue resistance of lungs and chest walls (accounts for approximately 20% of total resistance), and (c) airway resistance (R.sub.AW) defined as the ratio of driving pressure to the rate of air flow (80% of total resistance). [0146] 4) Lung Compliance, or pulmonary compliance, refers to the extensibility of the lungs. It is expressed as a change in volume divided by a change in pressure. There are two types of lung compliance: static and dynamic. Static compliance of lungs is the change in volume for a given change in transpulmonary pressure with zero gas flow. Dynamic lungs compliance is compliance of lungs at any given time during actual movement of air.

[0147] The graph of pressure vs. time (see FIG. 8) is processes to derive at least one value related to at least one pulmonary function characteristics of the tested subject, such as Thoracic Gas Volume (TGV) and Total Lung Capacity (TLC), by applying the appropriate formulas. In accordance with some embodiments, a mathematical model of the pulmonary function test device considers two chambers: first camber represents is lung that has a varying volume TGV, and second chamber has constant volume V.sub.D, wherein V.sub.D is the sum of the following two volumes: the internal volume of the pulmonary function test device 100 V.sub.STS, and the anatomic dead-space volume V.sub.DA, wherein the volume V.sub.D may be assumed to be constant between the biological bronchial airway resistance (not shown) and shutter module 150. The anatomic volume V.sub.DA may be evaluated, for healthy individuals, as V.sub.DA=W.sub.P.Math.K.sub.V, where W.sub.P is patient weight in kg, and K.sub.V is proportional coefficient that is K.sub.V=2.2.10.sup.6 m.sup.3/kg. The internal volume of the pulmonary function test device may be defined as a sum of volumes V.sub.STS=V.sub.F+V.sub.MT+V.sub.AT+V.sub.M, wherein V.sub.MT is the volume of the mouth tube, V.sub.F is the biological filter internal volume, V.sub.AT is the volume of the adapter tube, and V.sub.M is the internal volume of the pulmonary function test device (see FIG. 1). In accordance with some embodiments, the mathematical model may further rely on the following assumptions: the exhaled air is the perfect gas; the first law of thermodynamics may be used for describing the pressure changing process in the pulmonary system and the pulmonary function test device; the pressure and temperature within the pulmonary function test device chambers and pulmonary system are homogeneous; and that the process may be assumed to be isothermal.

[0148] The differential equations of the mathematical model are:

[00003]$\{\begin{array}{c}{\stackrel{.}{P}}_L=\frac{1}{\mathrm{TGV}}.Math.\left[P_L.Math.Q_L-G_R.Math.R.Math.T\right]\\ {\stackrel{.}{P}}_M=\frac{R.Math.T}{V_D}.Math.\left(G_R-G_S\right)\end{array}$

where P.sub.L is the absolute pressure in the lung (alveoli), P.sub.M is the absolute pressure in the constant volume V.sub.D, Q.sub.L is lung flow rate, G.sub.R is airway resistance mass flow rate, G.sub.S is shutter module mass flow rate, R is gas constant, T is exhalation gas temperature.

[0149] The graph of pressure vs. time (see FIG. 8) may be further processes to estimate key indices employed in a typical pulmonary function test. The standard measured indices commonly measured by typical pulmonary function test devices are Forced Vital Capacity (FVC), Forced Expired Volume in 1 second (FEV1), Forced Expiratory Flow between 25 and 75 percent of FVC (FEF25-75%), and Peak forced Expiratory Flow (PEF).

[0150] FVC is the volume of gas, which is expired during a forced exhalation. FVC may be estimated as FVC=FVC.sub.1+FVC.sub.2, wherein FVC.sub.1 may be calculated in the sample interval of the pressure measurement between point t.sub.1 and t.sub.2 as a function of effective area A.sub.S1, and FVC.sub.2 may be calculated in the sample interval of the pressure measurement between points t.sub.2 and t.sub.5 as a function of effective area A.sub.S2.

[0151] FEV1 is the volume exhaled in first second of forced exhalation.

[0152] FEF25-75% is the average expired flow rate during the exhalation time period between 25% and 75% of the FVC. Estimation of FEF25-75% may be calculated according to the following formula:

[00004]$F.Math..Math.E.Math..Math.F_{25-75.Math.\%}=\frac{.Math..Math.t.Math.\underset{j=N.Math..Math.25}{\overset{j=N.Math..Math.75}{.Math.}}.Math.Q_j}{t_{75}-t_{25}}$

where t is sample time interval, Q is the flow rate, which may be estimated as a function of the atmospheric pressure, the absolute pressure in the main chamber 110, and outlet area A.sub.S, wherein:

[00005]$\{\begin{array}{cc}{A}_S=A_{S.Math..Math.1},& \mathrm{for}.Math..Math.t_{1}t{t}_2\\ A_S=A_{S.Math..Math.2},& \mathrm{for}.Math..Math.t_2<t{t}_5\end{array}$

where t.sub.25 is the time at which the forced vital capacity is 0.25.Math.FVC, and t.sub.75 is the time at which forced vital capacity is 0.75.Math.FVC, and where N25=t.sub.25/t, N75=t.sup.75/t.

[0153] PEF is the maximum value of the exhaled flow, which is estimated on the basis of the flow rate Q as cited above, at the maximum point (t.sub.M,P.sub.MM).

[0154] The shape of the pressure vs. time graph (see FIG. 8) may vary from subject to subject. Moreover, analysis of different indices of the curve may be indicative of the medical condition of a subject, as such curves may differ between healthy subjects, patients suffering from restrictive ling disease and patients suffering from obstructive lung disease. Several groups of indices, in accordance with some embodiments, may be calculated in order to test whether they can be used to produce a diagnostic event, by analyzing their associations with and diagnose the subjects' pulmonary conditions. These groups are: time related indices, volume related indices, and pressure, as well as pressure with volume, related indices.

[0155] The group of time related indices may include, but is not limited to, the following indices:

[00006]$T.Math..Math.I=\frac{t_4-t_3}{t_4-t_1}$

is time index

[00007]$T.Math..Math.C_1=\frac{t_3-t_2}{t_5-t_2}$

is the first time coefficient

[00008]$T.Math..Math.C_2=\frac{t_5-t_3}{t_5-t_1}$

is the second time coefficient

[00009]${\mathrm{TC}}_3=\frac{t_2-t_1}{t_5-t_2}$

is the third time coefficient

[0156] The group of volume related indices may include, but is not limited to, the following indices:

[00010]$\mathrm{VI}=\frac{\left({\mathrm{VC}}_4-{\mathrm{VC}}_2\right).Math.{\mathrm{TLC}}_{\mathrm{PR}}}{{\mathrm{VC}}_4.Math.\mathrm{FEV}.Math..Math.1}$

is volume index

[00011]${\mathrm{VCC}}_1=\frac{{\mathrm{VC}}_4-{\mathrm{VC}}_2}{{\mathrm{VC}}_4}$

is the first volume coefficient

[00012]${\mathrm{VCC}}_2=\frac{\mathrm{FEV}.Math..Math.1}{{\mathrm{VC}}_5}$

is the second volume coefficient

[00013]${\mathrm{VCC}}_3=\frac{{\mathrm{VCC}}_1}{\mathrm{VI}}$

is the third volume coefficient

[0157] The group of pressure, and pressure with volume, related indices may include, but is not limited to, the following indices:

[00014]$\mathrm{PI}=\frac{P_{M.Math..Math.3}-P_A}{P_{\mathrm{MM}}-P_A}$

is pressure index

TVI.sub.O=0.5.Math.(TI+TC2+VC3) is obstructive index

where VC.sub.2 is vital capacity at the time t=t.sub.2 [m.sup.3], VC.sub.4 is vital capacity at the time t=t.sub.4 [m.sup.3], VC.sub.5 is vital capacity at the time t=t.sub.5(VC.sub.5FVC) [m.sup.3], and TLC.sub.PR is predicted value of the TLC [m.sup.3], wherein TLC.sub.PR may be estimated by formulas that incorporate age, gender, height and weight.

[0158] The terms Vital Capacity or VC and Forced Vital Capacity or FVC are interchangeable.

[0159] The estimation of TGV and TLC, in accordance with some embodiments, may be performed by three algorithms. The main assumption of the first algorithm is that the derivative of the pressure P.sub.M along the graph of pressure vs. time varies linearly. In accordance with some embodiments, the derivative of P.sub.M may vary linearly between two points before and after t.sub.4.

[0160] Stated otherwise, the first algorithm is based on the assumption that the derivative of the pressure versus time curve is linear during a period of time that encompasses t.sub.4. This period includes t.sub.4 nested between two predetermined time periods, occurring before and after t.sub.4. The two time periods may be identical, similar or different from one another. According to the first algorithm the TGV and TLC may be estimated as:

[00015]${\mathrm{TLC}}_1=B_{11}.Math.\frac{{\stackrel{.}{P}}_{M.Math..Math.4}.Math.{\mathrm{VC}}_4-B_{12}.Math.A_{S.Math..Math.2}.Math.\sqrt{P_A.Math.\left(P_{M.Math..Math.4}-P_A\right)}}{{\stackrel{.}{P}}_{M.Math..Math.4}}$${\mathrm{TGV}}_1={\mathrm{TLC}}_1-{\mathrm{VC}}_4$

where {dot over (P)}.sub.M4 is pressure derivative in the constant volume V.sub.D at the time t=t.sub.4, B.sub.11 and B.sub.12 are constant coefficients.

[0161] The main assumption of the second algorithm is the pulmonary function test device's pressure between two points along the graph of pressure vs. time may be approximated to an exponential curve. In accordance with some embodiments, the pressure between the second pressure peak and the end of the forced gas exhalation may be approximated to an exponential curve. According to the second algorithm the TGV and TLC may be estimated as:

[00016]${\mathrm{TGV}}_2=B_{21}.Math.A_{S.Math..Math.2}.Math.\left(t_4-t_3\right).Math.\left[\ln \left(\frac{P_{M.Math..Math.3}-P_A}{P_{M.Math..Math.4}-P_A}\right)\right].Math.\sqrt{\frac{P_{M.Math..Math.4}-P_A}{P_{M.Math..Math.3}-P_A}}$${\mathrm{TLC}}_2={\mathrm{TGV}}_2+{\mathrm{VC}}_4$

where B.sub.21 is constant coefficient.

[0162] The main assumption of the third algorithm is that the pulmonary function test device's pressure between two points along the graph of pressure vs. time may be approximated to a parabolic curve. In accordance with some embodiments, the pressure between the second pressure peak and the end of the forced gas exhalation may be approximated to a parabolic curve. According to the third algorithm the TGV and TLC may be estimated as:

[00017]${\mathrm{TGV}}_3=\frac{B_{31}.Math.{\left(t_4-t_3\right)}^3}{3.Math.\ln \left(\frac{P_{M.Math..Math.4}}{P_{M.Math..Math.3}}\right)}+\frac{B_{32}.Math.{\left(t_4-t_3\right)}^2}{2.Math.\ln \left(\frac{P_{M.Math..Math.4}}{P_{M.Math..Math.3}}\right)}$${\mathrm{TLC}}_3={\mathrm{TGV}}_3+{\mathrm{VC}}_4$

where B.sub.31 and B.sub.32 is constant coefficients.

[0163] In accordance with some embodiments, TGV and TLC may be estimated as averages of the values calculated according to all three algorithms.

[0164] The estimation of airway resistance R.sub.AW, in accordance with some embodiments, may be performed with assumption that during the switching between the shutter openings with effective area A.sub.S1 and A.sub.S2 the pressure varies only in the constant volume V.sub.D, while in the lung volume the pressure remains constant or its changing is negligible. This assumption is based on the fact that the lung volume at least on the order of magnitude greater than the V.sub.D volume. Then the airway resistance may be evaluated as:

[00018]$R_{\mathrm{AW}}=\frac{\left(P_{L.Math..Math.3}-P_{M.Math..Math.3}\right)}{Q_{L.Math..Math.3}}$

where Q.sub.L3 is lung flow rate at the time t=t.sub.3, P.sub.L3 is the absolute pressure in the lung (alveoli) at the time t=t.sub.3, which may be estimated as:

[00019]$P_{L.Math..Math.3}=\frac{\left(P_{M.Math..Math.3}-P_{M.Math..Math.2}\right).Math.\left(1-2.Math.B_R\right)}{1-B_R},$

wherein B.sub.R is constant coefficient.

[0165] The estimation of lung compliance C.sub.LT, in accordance with some embodiments, may be performed using the time constant parameter. A time constant of a compartment is a function of the compartment's resistance and compliance. The lung is assumed to behave as the compartment in the trial. The time constant .sub.E expresses how quickly a compartment can react to pressure alteration and provides an indication of the filling or emptying velocity of a lung compartment. The lung consists of a large number of compartments with variable time constants. This heterogeneity is often exaggerated with lung disease, such as pneumonia or pulmonary fibrosis. Clearly, the more inhomogeneous the lung ventilation is, the wider the spectrum of regional time constants. This causes variation in the filling and emptying periods and the filling volumes for individual compartments. At a given pressure, a compartment with high resistance and good compliance fills slowly with a resulting large volume, which may be indicative of asthma, for example. Conversely, a compartment with poor compliance and low resistance fills quickly, resulting in a smaller volume, which may be indicative of pulmonary fibrosis, for example. Mathematically, the time constant is defined as compliance multiplied by the airway resistance, and the resulting value is given in units of seconds:

.sub.E=R.sub.AW.Math.C.sub.LT

where the time constant .sub.E=t.sub.E/3, wherein t.sub.E is the duration of exhalation.

[0166] The duration of exhalation may be estimated as the time of pressure change from its start level at point (t.sub.1,P.sub.M1) until it reaches finish level at point (t.sub.5,P.sub.M5) (see FIG. 8). The time constant for this assumptions is .sub.E=(t.sub.5t.sub.1)/3, and lung compliance is:

[00020]$C_{\mathrm{LT}}=\frac{t_5-t_1}{3.Math.R_{\mathrm{AW}}}$

[0167] It is appreciated that certain features of the invention, which are, for clarity, described in the context of separate embodiments, may also be provided in combination in a single embodiment. Conversely, various features of the invention, which are, for brevity, described in the context of a single embodiment, may also be provided separately or in any suitable sub-combination or as suitable in any other described embodiment of the invention. No feature described in the context of an embodiment is to be considered an essential feature of that embodiment, unless explicitly specified as such.

EXAMPLES

Example 1: Comparison of TLC Estimated by the Pulmonary Function Test Device and Plethysmograph

[0168] A trial including 45 subjects, measured by the pulmonary function test device and a full body plethysmograph, was performed to compare TLC results between both methods. The trial population included 17 healthy subjects, 17 patients with COPD and 11 patients with restrictive lung disease. The trial did not provide diagnosis, as all patients participating in the trial were clinically pre-diagnosed by conventional methods prior to participation.

[0169] The terms plethysmograph and full body plethysmograph as used herein are interchangeable.

[0170] A statistical analysis was performed, wherein the base indicator that reflects the correlation between two variables is a linear correlation coefficient R, given by the following formula:

[00021]$R=\frac{n.Math..Math.\left(S.Math.P\right)-\left(.Math.S\right).Math.\left(.Math.P\right)}{\sqrt{n.Math.\left(.Math.S^2\right)-{\left(.Math.S\right)}^2}.Math.\sqrt{n.Math.\left(.Math.P^2\right)-{\left(.Math.P\right)}^2}}$

where n is sample size, S indicates the data estimated by the pulmonary function test device, and P indicates the data measured by a Plethysmograph. The coefficient of determination R.sup.2 denotes the strength of the linear association between two variables.

[0171] The coefficient of variation (CV) represents the ratio of the standard deviation to the mean, given by the formula:

[00022]$\mathrm{CV}=\frac{}{}$

where is standard deviation, calculated as ={square root over ((S/P).sup.2/(n1))}, and is the mean, calculates as =(S/P)/n.

[0172] Table 1 summarizes trial results of TLC estimated either by the pulmonary function test device or by a Plethysmograph. The predicted value TLC.sub.PR is presented in Table 1 next to the value of TLC measured by a Plethysmograph, compared to TLC estimations according to the three algorithms in accordance with some embodiments, as well as the average of the three algorithms.

[0173] Table 2 presents statistical analysis of the results presented in Table 1. The correlation coefficient of the TLC estimated by first correlated algorithm is R.sub.10.964 and its coefficient of determination is R.sub.1.sup.293%. The correlation coefficient of the TLC determined by average of two correlated algorithms (values not shown) is R.sub.20.962 and its coefficient of determination is R.sub.2.sup.292.6%, and the correlation and determination coefficient for the averaged three correlated algorithms are R.sub.30.939 and R.sub.3.sup.288.2%, respectively. The coefficient variation of the TLC estimated by the average of all three algorithms is CV0.096. The statistical analysis of Table 2 presents a strong correlation with the results measured by a Full Body Plethysmograph, considered to serve as a gold standard for such measurements.

TABLE-US-00001 TABLE 1 Trial results of TLC as measured by a Plethysmograph and estimated by the pulmonary function test device Gender, Age (Y), Predicted Plethys. Pulmonary function test device [L] H (cm), W (kg) [L] [L] Algor. #1 Algor. #2 Algor. #3 Average Pre-diagnosis M, 61, 175, 80 6.85 6 5.89 6.03 5.79 5.90 Healthy F, 85, 149, 70 4.25 3.85 3.92 4.28 4 4.07 Healthy M, 71, 173, 92 6.74 5.88 5.5 5.49 5.37 5.45 Healthy M, 60, 171, 75 6.58 7.81 8.18 8.56 7.61 8.12 Healthy F, 35, 165, 85 5.2 4.61 5.04 4.99 4.86 4.96 Healthy M, 89, 160, 75 5.7 5.3 5.06 5.11 4.9 5.02 Healthy M, 49, 173, 88 6.66 5.82 5.97 5.93 5.67 5.86 Healthy M, 27, 191, 69 8.18 9.31 9.49 9.03 9.35 9.29 Healthy M, 48, 179, 70 7.22 7.15 7.37 7.06 7.06 7.16 Healthy F, 86, 157, 90 4.57 3.93 3.6 3.5 3.63 3.58 Healthy M, 66, 175, 62 6.9 5.16 5.3 5.17 5.22 5.23 Healthy M, 38, 174, 77 6.82 7.1 7.28 7.14 7 7.14 Healthy F, 71, 150, 68 4.11 4.09 3.54 3.77 3.48 3.60 Healthy F, 78, 155, 90 3.71 3.82 3.46 3.72 3.45 3.54 Healthy F, 78, 147, 60 3.91 3.73 3.36 3.79 3.36 3.50 Healthy M, 58, 166, 53 6.18 5.53 5.33 5.43 5.23 5.33 Healthy M, 37, 164, 80 5.9 4.9 5.33 4.91 5.2 5.15 Heavy Smoker F, 59, 159, 80 4.84 4.78 4.35 4.2 4.58 4.38 COPD (Chr. Asthma) F, 58, 170, 100 5.49 5.6 5.9 6.03 5.78 5.90 COPD (Asthma) F, 58, 171, 61 5.5 5.72 5.73 5.69 6.32 5.91 COPD M, 70Y, 170, 68 6.56 7.43 7.59 5.34 6.77 6.57 COPD (Emphisema) F, 66Y, 166, 84 5.26 4.68 4.78 4.16 4.75 4.56 COPD M, 70, 180, 75 7.28 7.18 7.48 6.22 7.45 7.05 COPD (Emphisema) F, 70, 158, 68 4.8 4.9 4.22 3.65 4.45 4.11 COPD (Asthma) M, 73, 162, 67 5.86 4.31 4.87 4.31 4.75 4.64 COPD M, 62, 175, 97 6.9 6.11 5.87 3.19 6.46 5.17 COPD M, 74, 168, 78 6.34 5.9 5.88 4.95 5.91 5.58 COPD M, 74, 158, 71 5.54 4.49 4.87 4.53 5.2 4.87 COPD F, 86, 155, 63 4.44 3.43 3.27 2.57 3.39 3.08 COPD (Asthma) M, 92, 164, 76 6.02 4.79 5.19 3.46 5.38 4.68 COPD (Asthma) F, 63, 155, 89 4.44 2.81 3.39 2.5 3.5 3.13 COPD M, 68Y, 174, 75 6.82 6.2 5.78 4.91 6.84 5.84 COPD M, 80, 163, 86 5.94 4.29 4.79 4.1 5 4.63 COPD M, 65, 165, 100 6.06 6.6 6.03 5.75 6.18 5.99 COPD (Asthma) F, 75, 149, 66 4.25 2.25 1.82 1.6 1.91 1.78 Restrictive (Fibrosis) M, 79, 163, 66 5.94 3.28 4 2.71 4.1 3.60 Restrictive (Fibrosis) M, 43, 176, 128 6.98 4.75 4.96 3.96 4.9 4.61 Restrictive M, 64, 167, 117 6.26 2.76 3.44 2.37 3.66 3.16 Restrictive M, 5, 177, 88 7.06 5.73 5.88 5.26 5.86 5.67 Restrictive M, 69, 168, 85 6.34 3.92 4.3 4.03 4.81 4.38 Restrictive F, 65, 154, 82 4.37 3.18 3.32 2.84 3.34 3.17 Restrictive F, 54, 161, 101 4.84 3.33 3.52 3.38 3.54 3.48 Restrictive M, 77, 175, 94 6.9 3.45 3.94 3.26 4.05 3.75 Rest. + Light COPD M, 70, 167, 86 6.26 3.85 4.52 3.93 4.63 4.36 Rest. + Light COPD M, 73, 164, 89 6.01 4.1 4.5 4.2 4.5 4.40 Restrictive (Fibrosis) * M denotes Male, F denotes Female, H denotes Height, W denotes Weight, Plethys. denotes results of plethysmography measurements, Algor. denotes Algorithm, Average refers to average of algorithms #1-#3.

TABLE-US-00002 TABLE 2 Trial statistical analysis results Average of Alg. Average of Alg. Algorithm #1 #1 & #2 #1, #2 & #3 Correlation 0.964 0.962 0.939 Coefficient R Coefficient of 0.930 0.926 0.882 Determination R.sup.2 Coefficient of 0.085 0.092 0.096 Variation CV

[0174] FIG. 10 shows a Bland-Altman plot of the difference between TLC values measured by Plethysmograph and the pulmonary function test device against their mean. Such a plot is considered to be more informative than the correlation coefficient R in determining how close the two methods of measurement are. The Bland-Altman method calculates the mean difference between two methods of measurement (the bias), and 95% limits of agreement as the mean difference (2, or more precisely 1.96). It is expected that the 95% limits include 95% of the differences between the two measurement methods. The plot in FIG. 10 is a presentation of the 95% limits of agreement, mainly intended for visual estimation of how well two methods of measurement agree. The smaller the range between these two limits the better the agreement is. Such a plot also provides an assessment of the possible relationship between the measurement error and the true value. In this case, the mean difference is defined as .sub.D=(PS)/n, and the standard deviation of the difference is .sub.D={square root over ((PS).sup.2/(n1))}. For the TLC data of the current trial, the mean difference is .sub.D=0.05 l and standard deviation of the difference is .sub.D=0.37 l. If the differences are normally distributed (Gaussian distribution), 95% of the differences will be included between .sub.D2.Math..sub.D and .sub.D+2.Math..sub.D, or more precisely, between .sub.D1.96.Math..sub.D and .sub.D+1.96.Math..sub.D, Assuming that the differences included within .sub.D2.Math..sub.D would not be clinically important, they may be referred to as the limits of agreement. For TLC data of the current trial, said limits of agreement are:

.sub.D2.Math..sub.D=0.052.Math.0.37=0.69 l

.sub.D+2.Math..sub.D=0.05+2.Math.0.37=0.79 l

[0175] Thus, TLC estimations based on measurements of by the pulmonary function test device may be 0.69 l below or 0.79 l above Plethysmograph measurements. Such values are known in the art to be acceptable for clinical purposes, as reported, for example, in R. M. Schwartzstein, and M. J. Parker, Respiratory Physiology: A Clinical Approach (integrated Physiology), (Lippincott Williams & Wilkins, 2005), which noted that the clinically acceptable range in the comparison of lung volume measurement between body plethysmography and helium dilution is 0.9 liter

Example 2: Repeatability Assessment

[0176] A trial involving two subjects was performed in order to test repeatability, which is defined by the ISO as the closeness of agreement between independent test results under repeatability conditions that are as constant as possible, where independent test results are obtained with the same methods, on identical test items, in the same laboratory, performed by the same operator, using the same equipment, within short intervals of time. The first subject was a healthy 62 years old male, and the second subject was a 60 years old female suffering from chronical asthma disease. The Coefficient of Repeatability (CR) can be calculated as 1.96 times the standard deviation of the differences between Plethysmograph (P) and pulmonary function test device (S) data:

[00023]$\mathrm{CR}=1.96.Math.\sqrt{\frac{.Math.{\left(P-S\right)}^2}{m}}$

where m is the number of observations for each subject. As known in the art, for example in S. A. McKenzie et al, Arch. Dis. Child; 87: pp. 248-251 (2002), the CR is the value below which the absolute differences between the two measurements, Full Body Plethysmograph and pulmonary function test device, would lie with 0.95 probabilities.

[0177] FIG. 11a shows trial results histogram for the first subject, and FIG. 11b shows trial results histogram for the second subject, wherein bins are of The number of observations was set to 50 for each subject.

[0178] Table 3 shows the results of the repeatability tests of both subjects, which may be interpreted as acceptable for medical practice.

TABLE-US-00003 TABLE 3 Repeatability test results Healthy subject Chronic Asthma patient Parameters (Plethysm. TLC = 6.0 L) (Plethysm. TLC = 4.78 L) Maximum (L) 6.19 5.05 Minimum (L) 5.66 4.42 Average (L) 5.88 4.73 CR (L) 0.356 0.302

Example 3: Airway Resistance Estimation

[0179] A pre-clinical trial was performed on 10 voluntary patients. Estimation of the airway resistance, which is the ratio between pressure and flow, was performed for the pulmonary function test device and a plethysmograph model ZAN 500, referred to as the gold standard.

[0180] Table 4 shows the results of the airway resistance assessment for both the pulmonary function test device and the plethysmograph. The results indicate that there is a good correlation between resistance estimates of both systems for all 10 subjects. The differences in the results between the plethysmograph and the pulmonary function test device may be attributed to two major reasons: one important difference between said device and the plethysmograph is that said device estimates the bronchial airway resistance value, while the plethysmograph enables measurement of the total airway resistance, and another difference between the two systems is that for said device, the resistance is calculated around the pressure P.sub.M=P.sub.M3 (see FIG. 8), and in the plethysmograph said resistance is measured at the time when flow rate in the range of 0.5-1 L/s.

Example 4: Lung Compliance Estimation

[0181] A pre-clinical trials were performed on 10 voluntary patients, measured by the pulmonary function test, was performed to validate method of lung compliance estimation. Table 5 shows estimation results, corresponding to lung compliance at point (t.sub.3, P.sub.M3).

TABLE-US-00004 TABLE 4 Airway resistant assessment Plethysmograph pulmonary function test Subject's Gender & Age (kPa .Math. s)/l device (kPa .Math. s)/l Male, 61 years 0.568 0.506 Male, 71 years 0.403 0.489 Male, 47 years 0.805 0.686 Female, 59 years 0.795 0.746 Female, 85 years 0.498 0.541 Male, 37 years 0.739 0.732 Female, 35 years 0.307 0.520 Male, 89 years 0.545 0.705 Female, 85 years 0.543 0.490 Male, 59 years 0.605 0.758

TABLE-US-00005 TABLE 5 Lung compliance Lung Compliance Subject's Gender & Age ([m.sup.3/Pa] 10.sup.6) Male, 61 years 1.153 Male, 71 years 1.48 Male, 47 years 1.413 Female, 59 years (Chron. Asthma) 0.936 Female, 85 years 0.962 Male, 37 years 1.387 Female, 35 years 1.168 Male, 89 years 2.108 Female, 85 years 1.544 Male, 59 years 1.493

Example 5: Diagnosis of Lung Diseases

[0182] The main clinical objectives of respiratory function measurements include diagnosis, severity assessment, treatment monitoring and evaluation of prognosis. Spirometry serves as a useful diagnostic device, which measures, amongst other parameters, vital capacity (VC) and force expiratory volume in 1 second (FEV1). These parameters permit differentiation between restrictive and obstructive respiratory diseases. Spirogram flow-volume curves are used for diagnosis. However, combined obstructive and restrictive diseases are not always distinguished as single ailments inform such curves.

[0183] An objective, in accordance with some embodiments, of the pulmonary function test device's time-pressure curves and diagnostic indices, such as TI, VI, PI and TVI.sub.O, is to provide a tool for lung disease diagnosis.

[0184] FIG. 12a-c shows typical patterns of the pulmonary function test device's time-pressure curves differences between healthy (FIG. 12a), restrictive (FIG. 12b) and obstructive (FIG. 12c) subjects, in accordance with some embodiments.

[0185] Table 6 shows different value ranges of the diagnostic indices, in accordance with some embodiments, that may allow discrimination between different types of lung diseases.

TABLE-US-00006 TABLE 6 Diagnostic indices for different conditions Subject Time Index Volume Index Pressure Index Obstructive Index Condition TI VI PI TVI.sub.0 Healthy 0.5 < TI 0.8 0.4 < VI 0.8 0.3 < PI 0.8 TVI.sub.0 < 2.0 Obstructive 0.8 < TI 1.2 1.5 < VI 4.5 1.4 < PI 4.0 2.0 TVI.sub.0 4.0 Restrictive 0.3 < TI 0.5 0.8 < VI 1.5 0.8 < PI 1.4 TVI.sub.0 < 2.0

[0186] A trial involving 12 subjects, including healthy subjects, as well as patients suffering from either obstructive or restrictive lung disease, was performed to assess the pulmonary function test device's diagnostic capabilities.

[0187] Table 7 shows the resulting indices and their interpretation for the pulmonary function test device's measurements, compared to results of traditional diagnostics. Results show that there is a good correlation between the diagnostic results of the pulmonary function test device and traditional diagnostics. Nevertheless, further statistical analysis of curve patterns and index values for a wider variety of lung diseases of larger patient population may provide more reliable and accurate diagnostic capabilities as interpreted for measurements made by the pulmonary function test device.

TABLE-US-00007 TABLE 7 Diagnostic results of the pulmonary function test device Time Volume Pressure Obstructive Subject's Gender & Index Index Index Index Traditional STS Age TI VI PI TVI.sub.0 Diagnostic Diagnostic Male, 61 years 0.655 0.671 0.507 1.02 Healthy Healthy Male, 71 years 0.552 0.44 0.3 0.82 Healthy Healthy Male, 47 years 0.512 0.405 0.312 0.77 Healthy Healthy Female, 35 years 0.585 0.612 0.796 0.937 Healthy Healthy Female, 59 years 0.824 2.81 2.87 2.37 Chr. Asthma Obstruct. Mail, 70 years 0.855 3.31 1.4 2.51 Emphysema Obstruct. Female, 58 years 0.835 2.91 3.27 2.3 COPD Obstruct. Female, 70 years 0.896 4.41 1.75 3.11 Asthma Obstruct. Female, 75 years 0.392 1.0 0.662 0.99 Fibrosis Restrict. Female, 54 years 0.548 0.81 1.17 1.0 Restrict. Restrict. Mail, 79 years 0.548 1.58 1.25 1.4 Fibrosis Restrict. Female, 65 years 0.502 1.05 0.88 1.09 Restrict. Restrict.

[0188] Although the invention is described in conjunction with specific embodiments thereof, it is evident that numerous alternatives, modifications and variations that are apparent to those skilled in the art may exist. It is to be understood that the invention is not necessarily limited in to its application to the details of construction and the arrangement of the components and/or methods set forth herein. Other embodiments may be practiced, and an embodiment may be carried out in various ways. Accordingly, the invention embraces all such alternatives, modifications and variations that fall within the scope of the appended claims.