although many different imaging methods are suitable for use in the present invention , in a particularly preferred embodiment the present invention uses phase - contrast x - ray imaging ( pcxi ). pcxi exploits the phase change caused by x - ray refraction when passing between media of differing refractive indices to produce high contrast images of the lung . interference between transmitted and refracted x - rays produces high contrast images of the air / tissue boundaries compared with conventional x - ray absorption techniques . it is able to achieve this because the phase shift of the x - rays is generally more than three orders of magnitude greater than the absorption over the diagnostic x - ray energy range ( 20 kev - 90 kev ). an important consequence of this is that phase - contrast images can be recorded with significantly lower dose than conventional images , which is particularly important for both longitudinal studies and dynamic studies where repeated imaging is required . these benefits are particularly relevant to scientific use . for clinical use , reduced radiation dose is of value for reduction in cancer risk . since lung tissue motion is complex , dynamic and heterogeneous , the present invention is a novel application of the mathematical concept of divergence . in particular , even when the imaging used is of a two - dimensional nature , the divergence measure is highly correlated to changes in lung volume . fig1 ( a ) illustrates a plot of input pressure perturbations as measured with a pressure sensor at the airway opening . the input signal is composed of 9 distinct and different frequencies . fig1 ( b ) illustrates a plot of the vector divergence measured via a cross - correlation technique on phase contrast x - ray lung images . this is the global average of all locations across the lung (& gt ; 1000 locations ). each vector location produces its own response to the input pressure wave , correlating to the specific local lung properties in the region . fig2 ( a ) illustrates frequency response of the pressure input from fig1 , in the power spectrum . fig2 ( b ) illustrates frequency response of the measured vector divergence of fig1 . the 9 distinct frequencies can be obtained from the response of either pressure at the airway opening or the global average of vector divergence . fig3 illustrates lung maps showing lung expansion at specific frequencies . both the heart and the pressure inputs contribute to lung expansion at different frequencies , the heart &# 39 ; s first harmonic being at 3 . 6 hz . the pressure wave input frequencies are 4 , 6 , 10 and 14 hz . the heart harmonics can be seen and measured at 3 . 6 , 7 . 2 , 10 . 7 and 14 . 2 hz . above right is a larger version of the 4 hz oscillation to highlight the distribution of power of oscillations at that specific frequency . in the contour maps the large amount of information obtained from this technique is shown as each vector location (& gt ; 100 across the lung ) provides its own complete measure of local lung health . fig4 illustrates the u and v velocity components measured as ( a ) an fft of the global average of lung expansion or ( b ) an fft of each vector then globally averaged . the first method more clearly highlights the difference between u and v velocity components that are created by the heart . this can be used to identify the frequencies at which the heart has an effect as well as their relative magnitudes , thus allowing for very accurate filtering of the heart . fig5 illustrates amplitude of oscillation at 4 hz with the subject lying on its side . note that the bottom lung is supporting the weight of the other lung and the heart above it , and as a result has less measured oscillation amplitude . it appears that the more inflated the lung , the greater the oscillation amplitude . this technique is thereby suitable to measure regional affects for not only lung health but also lung mechanics and posture related changes . fig6 illustrates the power of input oscillations before ( blue ) and after the delivery of a dose of aerosol methacholine ( red ). there is a measurable decrease in the global average of expansion after the delivery of methacholine . fig7 ( a ) and fig7 ( b ) illustrates corresponding power maps of 4 hz oscillations before the delivery of methacholine ( fig7 ( a )) and post delivery of methacholine ( fig7 ( b )). fig8 ( a ) is a schematic diagram showing a suitable configuration for phase - contrast x - ray imaging . fig8 ( b ) and 8 ( c ) are examples of the possible increase in contrast for lung tissue obtained through phase - contrast x - ray imaging over absorption - based x - ray imaging . pxci exploits the phase change caused by x - ray refraction when passing between media of different refractive indices to produce high contrast images of the lung . interference between transmitted and refracted x - rays produces high contrast images of the air / tissue boundaries compared with conventional x - ray absorption techniques as illustrated by the images shown at fig8 ( b ) and 8 ( c ). fig9 ( a ) shows the 3d nature of x - ray illumination and velocimetric cross - correlation analysis . each 2d sampling region in the projection images represents a 3d volume for which a distribution of velocities may be present . the preferred parameter for the present invention is the modal velocity , which may significantly differ from the mean . fig9 ( b ) illustrates in vivo detection of lung tissue motion using instantaneous velocity of a healthy mouse lung ˜ 140 ms after the start of inspiration , shown as a vector field . vectors are reduced in number ( 293 of 2640 displayed ) for clarity . vectors are coloured according to magnitude ( from lowest ; blue , to highest ; red ) of velocity . the complete time sequence of inspiration consists of 70 instantaneous vector fields ( media ), one of which is shown . fig1 shows the empirical relationship between lung divergence and tissue expansion . fig1 ( a ) is a scatter - plot of divergence ( integrated throughout entire data series ) and lung volume ( measured by water plethysmography ) of a rabbit measured from its first breath . the solid line indicates a line of best fit , with an r 2 = 0 . 98 indicating excellent correlation between the data sets . fig1 ( b ) shows a time - series of lung volume ( measured by water plethysmography ) co - plotted with divergence ( integrated throughout the entire data series and normalised by the co - efficient determined by the fit in fig1 ( a ). fig1 shows physiological measures of lung pathology comparing the compliance for treated groups with controls ( statistically insignificant ) at 36 hours ( fig1 ( a )) and 6 days ( fig1 ( c )) after treatment . comparisons of the spontaneous tidal volumes ( v t ) at 36 hours ( fig1 ( b )) and 6 days after treatment ( fig1 ( d )). tidal volumes in controls are significantly lower than treated groups but are non - specific and global in nature . fig1 shows velocimetric measures of lung pathology comparing controls with groups 36 hours and 6 days after bleomycin exposure . frequency distribution of the divergence ( normalised to the average of controls ) is compared for treated groups ( n = 4 ) with controls ( n = 3 ). data are normalised by the average of the controls . as shown by fig1 ( a ), at 36 hours post treatment , treated mice have 24 % greater divergence on average and 14 % of treated lungs show difference 2 × the control average compared with less than 5 % for control lungs . as shown by fig1 ( b ), at 6 days post treatment , treated mice have 76 % greater divergence on average and 47 % of treated lungs show differences 2 × the control average compared with less than 4 % for control lungs . fig1 shows regional divergence within the lung and matching histology . fig1 ( a ) and 13 ( b ) are colour maps of regional divergence determined using x - ray velocimetry for typical ( a ) control , and ( b ) bleomycin - treated mice ( 6 days after exposure ). data are normalised by the average divergence across the control group and colour maps generated using the same colour scale ( see legend ). the mice treated with bleomycin ( b ) have dramatic regional alterations in the pattern of divergence . histological image ( c ) from lung imaged in ( a ) is typical of the control group . histological images ( fig1 ( d ) and 13 ( e )) from lung imaged in ( b ) are typical of the pathological group 6 days after bleomycin treatment . treated lungs show both regions of healthy tissue ( d ) and localised regions that are both hypercellular and endatemous ( e ). the following non - limiting example illustrates how the combination of pcxi and velocimetry can produce quantitative measures of regional lung motion , which can be used to differentiate between normal and abnormal lung tissue . furthermore the example illustrates that this technology is more sensitive and provides richer quantitative information for disease detection than other conventional measure such as global lung function tests and non - biased histological sampling . specifically , the present example illustrates the present invention when using single camera / 2d imaging . furthermore it illustrates the use of divergence as a measure of lung expansion . despite the two - dimensional nature of imaging , the divergence measure is highly correlated to changes in lung volume . this measure was evaluated in a bleomycin - induced lung injury model in immuno - deficient balb / c nude mice that are known to have a reduced inflammatory response to bleomycin compared to other strains . furthermore , mice are examined at 36 hours and 6 days after treatment to examine the early states of disease and determine whether disease progression can be detected . protocol : adolescent balb / c nude male mice were exposed to bleomycin ( 20 mg / kg in 20 ul saline ; n = 8 ; sigma - aldrich , australia ) or saline ( 20 ul ; n = 6 ) by intranasal instillation and lung function was tested daily using whole body plethysmography . mice were imaged at 36 h ( n = 4 ) or 6 days ( n = 4 ) after treatment . for imaging , mice were anaesthetized ( pentobarbital ; 15 mg / kg i . p . ), muscle relaxed ( pancuronium 1 mg / kg i . m . ), intubated and placed in a prewarmed ( 37 ° c .) water - filled plethysmograph . during imaging , mice were ventilated using a custom - designed ventilator at a peak inspiratory pressure of 20 cmh 2 o and end expiratory pressure of 2 cmh 2 o . inspiration and expiration times were 2 . 5 s and 1 . 5 s respectively . although this is significantly less than the normal ventilation rate for free - breathing mice , the imaging procedures only lasted for 5 breaths ( 20 s ) and was not expected to result in hypoxia or any changes in lung mechanics for these somnolent mice . following imaging , mice were killed ( pentobarbital ; 100 mg / kg i . p .) and the lungs fixed ( in 10 % formalin ) via the airways at a distending pressure of 20 cmh 2 o . paraffin - embedded sections ( 5 μm ) were stained with massons trichrome and used for histological analysis . 5 fields of view were chosen at random from at least 3 randomly selected sections per mouse to measure the relative volume density of abnormal parenchymal lung regions . then a subset analysis was performed to compare the relative tissue volume in normal and abnormal parenchymal regions using an unpaired t - test . phase - contrast x - ray imaging : studies were conducted in experimental hutch 3 of bl20b2 at the spring - 8 synchrotron in japan . the beamline consists of a bending magnet insertion device and si - 111 crystal monochromators , which generates a bright monochromatic x - ray beam . the x - ray beam transmits through the sample onto a scintillator , which converts the x - rays to visible light to be imaged by an optical detector system . imaging was conducted at 25 kev with a sample - to - detector distance of 2 m . images were acquired using an x - ray converter ( hamamatsu , bm5 ) and an emccd ( hamamatsu , c9100 - 02 ) camera ( fig8 ( a )) with an effective pixel size of 19 . 0 um . image acquisition occurred at 29 frames per second ( an exposure time of 20 ms with a 14 . 5 ms delay between exposures , corresponding to 34 . 5 ms between the start of frame acquisitions ) and was synchronized with ventilation to acquire 70 frames during the first part of inspiration and 30 frames during the first part of expiration for each breath . the mice were imaged in the upright position with all images acquired to obtain a frontal view of the entire thorax without the need for scanning or tiling . in all images ( fig9 ) the images are displayed without intensity inversion or laterally flipping and hence appear opposite in both regards in comparison to clinical x - ray images . velocimetry : the velocimetric analysis employed to measure lung motion is based on particle image velocimetry ( piv ); this is an established technique for measuring differential fluid velocities , including blood flow . piv determines the movement of particles from one image to the next , yielding information on both velocity and direction of particle movement . the basic concept is demonstrated in fig9 ( a ). images are paired and discretised into small sub - regions and cross - correlations are performed between the sub - regions in consecutive images . the position of the maximum of the cross - correlation function determines the most common ( modal ) inter - frame displacement of the structures within each sub - region . division of the displacement by the known inter - frame time yields the local modal velocity . x - ray velocimetry has been utilised for the measurement of flow within channels for blood flow and has recently been adapted to 3d analysis . the high contrast intensity patterns produced by pcxi of the lung can be used instead of having to introduce exogenous particles as is the practice in conventional piv . as a result , a comprehensive map displaying regional tissue velocities can be generated at all stages of the breathing cycle . whole animal motion was removed from image sequences by velocimetric analysis of upper vertebrae , followed by interpolation of images onto a static reference frame . lungs were isolated from images by band - pass filtering using the appropriate image frequencies and regions containing the lungs were identified and masked . following this pre - processing , velocimetric analysis was conducted for 5 consecutive inspirations using customised software and the date phase - average to produce data sets of 70 frames displaying the velocity vector fields throughout inspiration for each animal . divergence : at every time - point the spatial derivatives of the velocity fields can be evaluated and summed to form the two - dimensional divergence field . the spatial derivative distinguishes between bulk displacement of tissue and regional variations in tissue displacement , highlighting local differences in motion between regions . the local differences are directly related to local tissue expansion , and hence local variations in the divergence would be considered to be a measure of heterogeneity of tissue expansion and , by implication , tissue properties . for the sake of clarity of expression this projected divergence in motion will hereafter simply be referred to as the divergence . the total divergence over inspiration is the sum of the divergence between each pair of subsequent time points . as the data are integrated over the entire inspiration , total divergence is represented in a single map . statistics : unpaired one - tailed t - tests were used to compare mean tidal volume and histological parameters . two - way repeated measures anova was used to determine differences in frequency distributions of divergence and time of divergence . results were considered statistically significant at p & lt ; 0 . 05 . values are reported as mean +/− sem ( unless stated otherwise ). physiological analysis : plethysmography was used to measure tidal volume ( vt ) during both spontaneous breathing and mechanical ventilation . lung compliance was assessed only during mechanical ventilation , as airway pressure measurements were not available during spontaneous breathing . at 36 h after bleomycin treatment , the spontaneous vt was significantly increased from 1 . 4 +/− 0 . 1 ml / kg in saline - treated mice to 1 . 9 +/− 0 . 1 ml / kg in bleomycin - treated mice , whereas spontaneous breathing rates were reduced from 480 +/− 38 to 283 +/− 34 breaths / min . however , no significant differences in global lung compliance ( 0 . 68 +/− 0 . 04 vs 0 . 71 +/− 0 . 10 ml / cmh 2 o / kg ) could be detected during mechanical ventilation ( fig1 ). similarly , at 6 days after bleomycin treatment , the spontaneous vt remained significantly elevated ( 1 . 2 +/− 0 . 1 vs 1 . 6 +/− 0 . 2 ml / kg ) and the ventilation rate was reduced ( 440 +/− 13 vs 235 +/− 40 breaths / min ), but no significant affect on global lung compliance ( 0 . 5 +/− 0 . 04 vs 0 . 55 +/− 0 . 05 ml / cmh 2 o / kg ) was detected during mechanical ventilation ( fig1 ( c ) & amp ; 11 ( d )). velocimetry : determined using x - ray velocimetry , the velocity vectors define the timing and extent of regional lung motion throughout a breath . the vectors measured at mid - inspiration demonstrate that regional lung motion is very heterogeneous ( fig9 ) at this point in the breathing cycle . divergence — validation : to validate the divergence analysis , an existing , published , pcxi image data set was evaluated using the divergence analysis outlined herein . pcxi image were acquired as rabbit foetus lungs were slowly ventilated from their first breaths ( in situ ) inside a water plethysmograph — allowing concurrent measurements of lung air volume as it rises from zero . the pcxi images were analysed for velocimetry and divergence . the divergence data was then integrated throughout the entire series rather than just each breath and compared to the plethysmograph data . a scatter - plot of the integrated divergence ( fig1 ( a )) plotted against the lung volume as measured by water plethysmograph shows strong correlation between the two quantities ( r 2 = 0 . 98 ). the gradient of the line of best - fit could be considered as a combination of the image magnification ( pixel size ) and an equivalent thickness of the lung since this gradient represents the empirically - determined conversion vector to convert between the change in area ( in pixels ) and a change in volume ( ml ). to demonstrate the direct relationship between these quantities , the divergence data has been normalised by the empirical thickness and plotted ( fig1 ( b )) with the lung volume as measured by water plethysmography against time . despite some minor inaccuracies brought about by the imperfect sealing of the animal in the water plethysomgraph and the two - dimensional nature of the image data , the excellent agreement between these data demonstrates the direct link between divergence and tissue expansion and hence tissue mechanical properties . divergence — bleomycin treatment : divergence data for bleomycin treated mice and controls were normalised by the average of controls and accumulated into frequency distribution curves , demonstrating that major difference ( p & lt ; 0 . 001 ) in regional lung motion can be detected between saline - treated and bleomycin - treated mice both 36 h and 6 days after treatment . frequency distribution curves ( fig1 ) show that divergence , within individual regions of lung tissue , was on average 24 % greater in bleomycin - treated mice compared to saline - treated controls at 36 h after treatment , despite having the same global vt . moreover , in bleomycin treated mice at 36 h after treatment , 14 % of lung regions showed local differences twice the mean value in saline - treated controls ; less than 5 % of lung regions showed this degree of divergence in control mice , indicating a three - fold increase between bleomycin - treated and control mice . at 6 days after bleomycin treatment , the changes in lung motion were more enhanced than those detected at 36 h after treatment . indeed , despite no changes in global lung compliance , x - ray velocimetry detected a highly significant shift in the frequency distribution curve towards greater divergence ( fig1 ( b )). on average , the divergence of individual lung regions was 76 % greater in bleomycin - treated mice compared to saline - treated controls . furthermore 47 % of lung regions in treated mice show local differences in motion twice the mean value of saline - treated mice , while less than 4 % of lung regions showed this degree of divergence in saline - treated mice . hence there is nearly a 12 fold difference in divergence between bleomycin - treated mice and controls at 6 days after treatment . from these data , regional maps of divergence can be reconstructed and superimposed on the lung image to identify regions with abnormal motion ( fig1 ) . although ventilators , plethysmographs and spirometers can measure many characteristics of lung function , those measures reflect the integrated average of the entire lung . as a result these techniques have limited ability to detect regional lung disease until it is sufficiently widespread to influence total lung function . in contrast , the present invention uses the capabilities of x - ray velocimetry to non - invasively detect breath - by - breath alterations in regional lung motion that occur even during the early stages of lung disease . importantly , the velocimetric technique offers the advantages of detecting regional changes in lung function early , accurately and most importantly , in situ . superposition of the lung tissue velocity maps over the phase - contrast images acquire at mid - inspiration ( 140 ms after inspiration onset ) clearly demonstrates the heterogeneity of lung motion ( fig9 ). combining consecutive images reveals the changing dynamic of the velocity vector field during a breath and demonstrates that regional lung tissue motion is complex and non - linear . indeed , regional velocities are most likely influenced by local characteristics of regional compliance , the compliance and motion of nearby tissue as well as the proximity to structure such as the diaphragm , heart and chest wall . for example , lung tissue near the diaphragm displayed significantly more motion than tissue near the apex of the lung ( fig9 ( b )) which is likely due to differences in compliance as well as motion and activity of the chest wall that is immediately adjacent to the lung tissue . to accommodate the large differential in normal motion across the lung , a functional measure was derived from the velocity fields to identify regions with abnormal motion potentially caused by disease . specifically the local divergence was calculated , normalised to the average for all controls of reach treatment period ( eg 36 hours or 6 days ) so that differences could be detected . the measure of divergence was derived from integration of the velocity vector field within each region over an entire breath . it is well understood that the divergence of a velocity field relates to the local expansion or contraction of the object , which in this case is lung tissue . as such the divergence accounts for normal variations in tissue motion ( such as the increased motion near to the diaphragm compared to the apex ) and converts heterogeneous patterns of tissue motion ( displayed by normal tissue ) to a homogenous pattern of divergence . however it is anticipated that heterogeneous regions of tissue properties ( either resistance or compliance ) will result in local variations in divergence . however , lung tissue motion is three - dimensional and divergence measures are two - dimensional . therefore the most correct interpretation of the measure of divergence is that it directly relates to local heterogeneity of lung tissue motion cause by differences in expansion . if all ventilation parameters , such as inflation times , pressures and gas flows are kept constant , their local heterogeneity in motion reflects differences in the mechanical response of lung tissue across the lung . this altered response must be due to either changes in the tissue mechanical properties , or a constriction / dilation of the airways leading to local alteration in resistance or compliance . to test the ability of x - ray velocimetry and our subsequent analysis to detect abnormal lung motion , mice were exposed to bleomycin which resulted in progressive lung injury . inhaled bleomycin is well characterised and commonly used experimental model of pulmonary fibrosis that begins with the initiation of an inflammatory cascade . since balb / c nude mice ( an immuno - deficient strain ) were utilised , it is not surprising that the pulmonary fibrotic response was reduced in these mice compared with reports in other strains . this is likely because inflammatory responses are reduced in these mice , although a recent study has also observed a similar reduced response in conventional balb / c mice . in any event it is apparent that the induced pathological changes were insufficient to affect global lung compliance at the time points measured , as observed previously in conventional balb / c mice . the mice were deliberately examined at 36 hours and 6 days after treatment which encompasses the early stages of disease pathogenesis when the induced changes can be detected histologically but not physiologically . indeed a previous study has indicated that the maximum response to bleomycin occurs at 21 days for conventional balb / c mice which is the only time that functional changes can be detected . however , it is clear from the histological analysis ( fig1 ) that the bleomycin treatment was sufficient to cause observable focal lesions within lung tissue . the present velocimetric analysis in bleomycin - treated mice revealed highly significant changes in regional lung motion compared to saline - treated mice . despite having similar tidal volumes and inflation pressures during mechanical ventilation ( indicating no change in global compliance ), bleomycin treatment increased divergence across the lung by 24 % at 36 h and by 76 % at 6 days . this highly sensitive measure yielded a three fold difference between the groups after only 36 hours of treatment and was increased further after 6 days of treatment . it is apparent that it was not possible to detect a change in global lung compliance during mechanical ventilation at either 36 hours or 6 days after treatment , but that tidal volumes during spontaneous breathing were significantly increased by bleomycin treatment at both times ( fig1 ). the higher tidal volumes during spontaneous breathing were partially matched with lower ventilation rates resulting in minute ventilation ( ml / min / kg ) values that were slightly reduce by statistically similar in bleomycin and saline - treated mice . in summary pcxi combined with velocimetry can measure regional lung motion and define the regional velocity changes at each stage of the respiratory cycle . despite the large heterogeneity in normal motion across the lung , a detailed analysis of the velocity vectors can provide a very sensitive method for detecting abnormal motion caused by respiratory disease . a time - cycled pressure - limited ventilator was developed , the ventilator operating using labview &# 39 ; s virtual instrument ( vi ) controls to synchronise image acquisition with mechanical ventilation of small animals . a personal computer ( pc ) along with a data acquisition module ( ni usb - 6259 ) and national instruments lab view software were used to control the ventilator , as illustrated in fig1 . table 1 provides a detailed description of the ventilator components . the data acquisition system connects to the pc with a universal serial bus ( usb ) cable and can handle up to 4 analog outputs , 80 analog inputs and 48 digital input / output channels at 1 . 25 ms / s . labview can simultaneously control multiple devices and readily interface with external hardware since many hardware drivers are included in the programming library . the main advantage of the virtual interface is that all parameters ( eg air pressure and flow rate ) can be controlled remotely with real - time display . fig1 depicts the manner in which air was cycled around the ventilator and delivered to the lung . two pressure vessels , set at different pressures were used to control the pip and peep . at the start of inspiration the inspiratory solenoid valve opened whilst the expiratory solenoid remained closed ( fig1 and 15 ) for the entire set inspiratory time . air from the pip vessel flowed to the lung through the inspiratory solenoid via a variable restrictor valve ; this allowed the air to flow into the lungs until the airway pressure reached the pressure of the pip vessel . as the pip and peep vessels were of sufficiently large volume (˜ 1 l / box ), the volume change associated with opening and closing of the respiratory solenoids did not significantly influence the pressure within the vessels . this provided a high degree of stability with the applied pip and peep pressures . a variable restrictor valve allowed almost infinite variability of the rate of gas flow into the lung from the pressure vessel , which in turn was controlled remotely via the virtual interface . as a result the inspiratory pressure wave form could be varied and the length of an inspiratory pressure plateau ( as a proportion of inspiration time ) could be set by regulating the inspiratory airflow independently of this pip . upon expiration the states of the respiratory solenoids are simultaneously flipped , allowing the lungs to deflate to the lower peep level for a preset period . during expiration , airflow from the lung into the peep box could also be regulated via a restrictor valve , as shown in fig1 . both the inspiratory and expiratory times can be updated in software while the ventilator is operating . as a safety precaution the solenoids simultaneously closed when the ventilation sequence was terminated to prevent over - distension of the airways . a ventilator of this type is capable of performing high frequency ventilation up to at least 33 hz , which is well in excess of typical clinical usage of about 2 to 12 hz . while this invention has been described in connection with specific embodiments thereof , it will be understood that it is capable of further modification ( s ). this application is intended to cover any variations uses or adaptations of the invention following in general , the principles of the invention and including such departures from the present disclosure as come within known or customary practice within the art to which the invention pertains and as may be applied to the essential features hereinbefore set forth . as the present invention may be embodied in several forms without departing from the spirit of the essential characteristics of the invention , it should be understood that the above described embodiments are not to limit the present invention unless otherwise specified , but rather should be construed broadly within the spirit and scope of the invention as defined in the appended claims . the described embodiments are to be considered in all respects as illustrative only and not restrictive . various modifications and equivalent arrangements are intended to be included within the spirit and scope of the invention and appended claims . therefore , the specific embodiments are to be understood to be illustrative of the many ways in which the principles of the present invention may be practiced . in the following claims , means - plus - function clauses are intended to cover structures as performing the defined function and not only structural equivalents , but also equivalent structures . for example , although a nail and a screw may not be structural equivalents in that a nail employs a cylindrical surface to secure wooden parts together , whereas a screw employs a helical surface to secure wooden parts together , in the environment of fastening wooden parts , a nail and a screw are equivalent structures . it should also be noted that where a flowchart is used herein to demonstrate various aspects of the invention , it should not be construed to limit the present invention to any particular logic flow or logic implementation . the described logic may be partitioned into different logic blocks ( e . g ., programs , modules , functions , or subroutines ) without changing the overall results or otherwise departing from the true scope of the invention . often , logic elements may be added , modified , omitted , performed in a different order , or implemented using different logic constructs ( e . g ., logic gates , looping primitives , conditional logic , and other logic constructs ) without changing the overall results or otherwise departing from the true scope of the invention . various embodiments of the invention may be embodied in many different forms , including computer program logic for use with a processor ( e . g ., a microprocessor , microcontroller , digital signal processor , or general purpose computer and for that matter , any commercial processor may be used to implement the embodiments of the invention either as a single processor , serial or parallel set of processors in the system and , as such , examples of commercial processors include , but are not limited to merced ™, pentium ™, pentium ii ™, xeon ™, celeron ™, pentium pro ™, efficeon ™, athlon ™, amd ™ and the like ), programmable logic for use with a programmable logic device ( e . g ., a field programmable gate array ( fpga ) or other pld ), discrete components , integrated circuitry ( e . g ., an application specific integrated circuit ( asic )), or any other means including any combination thereof . in an exemplary embodiment of the present invention , predominantly all of the communication between users and the server is implemented as a set of computer program instructions that is converted into a computer executable form , stored as such in a computer readable medium , and executed by a microprocessor under the control of an operating system . computer program logic implementing all or part of the functionality where described herein may be embodied in various forms , including a source code form , a computer executable form , and various intermediate forms ( e . g ., forms generated by an assembler , compiler , linker , or locator ). source code may include a series of computer program instructions implemented in any of various programming languages ( e . g ., an object code , an assembly language , or a high - level language such as fortran , c , c ++, java , or html . moreover , there are hundreds of available computer languages that may be used to implement embodiments of the invention , among the more common being ada ; algol ; apl ; awk ; basic ; c ; c ++; conol ; delphi ; eiffel ; euphoria ; forth ; fortran ; html ; icon ; java ; javascript ; lisp ; logo ; mathematica ; matlab ; miranda ; modula - 2 ; oberon ; pascal ; perl ; pl / i ; prolog ; python ; rexx ; sas ; scheme ; sed ; simula ; smalltalk ; snobol ; sql ; visual basic ; visual c ++; linux and xml .) for use with various operating systems or operating environments . the source code may define and use various data structures and communication messages . the source code may be in a computer executable form ( e . g ., via an interpreter ), or the source code may be converted ( e . g ., via a translator , assembler , or compiler ) into a computer executable form . the computer program may be fixed in any form ( e . g ., source code form , computer executable form , or an intermediate form ) either permanently or transitorily in a tangible storage medium , such as a semiconductor memory device ( e . g ., a ram , rom , prom , eeprom , or flash - programmable ram ), a magnetic memory device ( e . g ., a diskette or fixed disk ), an optical memory device ( e . g ., a cd - rom or dvd - rom ), a pc card ( e . g ., pcmcia card ), or other memory device . the computer program may be fixed in any form in a signal that is transmittable to a computer using any of various communication technologies , including , but in no way limited to , analog technologies , digital technologies , optical technologies , wireless technologies ( e . g ., bluetooth ), networking technologies , and inter - networking technologies . the computer program may be distributed in any form as a removable storage medium with accompanying printed or electronic documentation ( e . g ., shrink wrapped software ), preloaded with a computer system ( e . g ., on system rom or fixed disk ), or distributed from a server or electronic bulletin board over the communication system ( e . g ., the internet or world wide web ). hardware logic ( including programmable logic for use with a programmable logic device ) implementing all or part of the functionality where described herein may be designed using traditional manual methods , or may be designed , captured , simulated , or documented electronically using various tools , such as computer aided design ( cad ), a hardware description language ( e . g ., vhdl or ahdl ), or a pld programming language ( e . g ., palasm , abel , or cupl ). hardware logic may also be incorporated into display screens for implementing embodiments of the invention and which may be segmented display screens , analogue display screens , digital display screens , crts , led screens , plasma screens , liquid crystal diode screen , and the like . programmable logic may be fixed either permanently or transitorily in a tangible storage medium , such as a semiconductor memory device ( e . g ., a ram , rom , prom , eeprom , or flash - programmable ram ), a magnetic memory device ( e . g ., a diskette or fixed disk ), an optical memory device ( e . g ., a cd - rom or dvd - rom ), or other memory device . the programmable logic may be fixed in a signal that is transmittable to a computer using any of various communication technologies , including , but in no way limited to , analog technologies , digital technologies , optical technologies , wireless technologies ( e . g ., bluetooth ), networking technologies , and internetworking technologies . the programmable logic may be distributed as a removable storage medium with accompanying printed or electronic documentation ( e . g ., shrink wrapped software ), preloaded with a computer system ( e . g ., on system rom or fixed disk ), or distributed from a server or electronic bulletin board over the communication system ( e . g ., the internet or world wide web ). “ comprises / comprising ” and “ includes / including ” when used in this specification is taken to specify the presence of stated features , integers , steps or components but does not preclude the presence or addition of one or more other features , integers , steps , components or groups thereof . thus , unless the context clearly requires otherwise , throughout the description and the claims , the words ‘ comprise ’, ‘ comprising ’, ‘ includes ’, ‘ including ’ and the like are to be construed in an inclusive sense as opposed to an exclusive or exhaustive sense ; that is to say , in the sense of “ including , but not limited to ”.
a method for dynamic investigation of a subject lung , the method comprising the steps of : imparting an oscillation to the lung at one or more forcing frequencies so as to elicit a lung response , sensing the response of the lung simultaneously with the imparting of the oscillation to elicit a lung response , choosing at least one parameter used in the sensing to define the lung motion associated with the lung response , comparing one of the chosen parameters at each forcing frequency with the response at the forcing frequency in at least one region of the lung , and recording the comparison of step .