fig1 shows a system generally indicated by 100 for performing regional assessment of the lungs in accordance with one embodiment of the invention . a plurality of n sound transducers 105 , of which four are shown , are applied to a planar region of the chest or back skin of individual 110 . the transducers 105 may be applied to the subject by any means known in the art , for example using an adhesive , suction , or fastening straps . each transducer 105 produces an analog voltage signal 115 indicative of pressure waves arriving to the transducer . the analog signals 115 are digitized by a multichannel analog to digital converter 120 . the digital data signals p ( x i , t ) 125 , represent the pressure wave at the location x i of the ith transducer ( i = 1 to n ) at time t . the data signals 125 are input to a memory 130 . data input to the memory 130 are accessed by a processor 135 configured to process the data signals 125 . the signals 125 may be denoised by filtering components such as components having frequencies outside of the range of lung sounds , for example , vibrations due to movement of the individual . each signal 125 may also be subject to band pass filtering so that only components in the signal within a range of interest are analyzed . an input device , such as a computer keyboard 140 or mouse 145 , is used to input relevant information relating to the examination such as personal details of the individual 110 . the input device 140 may also be used to input values of one or more times t 1 and t 2 that specify times at which the signals p ( x i , t ) are to be analyzed or that specify one or more time intervals over which no signals p ( x i , t ) are to be analyzed . the processor 135 calculates the value of a parameter at a plurality of locations over the lungs at the specified times or over the specified time intervals . in a preferred embodiment , the processor 135 is configured to calculate an average acoustic energy { tilde over ( p )}( x , t 1 , t 2 ) over a time interval from t 1 to t 2 at a plurality of locations x in the region r in a calculation involving at least one of the signals p ( x i , t ). the locations at which the parameter was calculated are divided into groups , where each group overlies a region of the lungs . the processor 135 is further configured to perform a regional assessment of the lungs . the regional assessment comprises for each of the groups determining the value of one or more regional parameters where each regional parameter is obtained in a calculation involving the parameter values calculated at the location in the region . for example , a regional parameter may be the sum of the parameters in the region , the maximum of the parameter value , the minimum or the average . the regional parameter values may be normalized by dividing the regional parameter by the sum of the regional parameter values . fig2 shows a flow chart diagram for carrying out the method of the invention in accordance with a preferred embodiment . in step 200 the signals p ( x i , t ) are obtained from n transducers placed at predetermined locations x i for i from 1 to n overlying the lungs . in step 205 values of one or more times are either input to the processor 135 using the input devices 140 or 145 , or are determined by the processor . in step 210 , a value of a parameter is determined at a plurality of locations x at the one or more input times or over one or more intervals . in step 220 a regional parameter is calculated in each of two or more predetermined lung regions . in step 225 , the total of the regional parameters is calculated . in step 230 , for each region , the regional parameters are normalized by dividing them by the calculated total to generate the regional assessment of the region . in step 240 , an image of the lungs is displayed on the display 150 in which the lungs are divided into the predetermined lung regions , and the normalized or non - normalized regional parameter for each region is displayed in the region in the image . the regional assessment is the total average acoustic energy in the region over the time interval or the total acoustic energy of the region divided by the total acoustic energy of the lungs . in a most preferred embodiment of the invention , an image of the lungs is generated from the average acoustic energies obtained over a time interval . u . s . pat . no . 6 , 887 , 208 to kushnir et al . discloses generating an image of the lungs from average acoustic energies calculated at a plurality of locations over the lungs . the image of the lungs is displayed on a display monitor with the lungs in the image being divided into the lung regions . it will also be understood that the system according to the invention may be a suitably programmed computer . likewise , the invention contemplates a computer program being readable by a computer for executing the method of the invention . the invention further contemplates a machine - readable memory tangibly embodying a program of instructions executable by the machine for executing the method of the invention . two subjects were subjected to regional assessment of lung function by vq scan and by the method of the invention . the first subject was a 35 year old male having a bmi ( body weight to height squared ) of 26 who never smoked . the second subject was a 71 year old male having a bmi of 30 who quite smoking five years prior to undergoing regional assessment of lung function . the second subject had a piy ( packs of cigarettes smoked per day times the number of years of to smoking ) of 150 for the regional assessment carried out by the method of the invention , a two - dimensional coordinate system was defined on the subject &# 39 ; s back . as shown in fig3 a , 48 transducers were placed on the individual &# 39 ; s back over the lungs at the locations indicated by the circles 300 . the curves 305 show the presumed contours of the lungs . as can be seen , the transducers were arranged in a regular orthogonal lattice with a spacing between the transducers in the horizontal and vertical directions of 5 cm . the signals p ( x i , t ) were then recorded . each signal was filtered using a low - pass filter having a cut - off of 150 hz . the average value of each filtered function p ( x i , t ) over the respiratory cycle is indicated in fig3 a by means of gray level shading of each circle 300 with reference to the gray level scale 310 . { tilde over ( p )}( x , t 1 , t 2 ) was obtained using equations ( 1 ) and ( 2 ) above with the kernel g of equation ( 5 ) with σ = 36 pixels . fig4 a shows an image 500 of the lungs obtained by the method of u . s . pat . no . 6 , 887 , 208 on the first subject . the image is a two - dimensional array of pixels x ( i , j ), where x ( i , j ) is the gray value or other intensity value at the pixel ( i , j ), where i and j are the column number and row number respectively of the pixel . the image 500 was divided into six regions using the algorithm shown in the flow chart diagram depicted in fig5 . in step 400 the intensity values in each column i are summed to yield column sums the graph 501 of the function a i is shown in fig4 a . the function a i has a local minimum 502 that identifies the boundary between the left lung 504 and the right lung 506 in the image 500 . in step 402 a vertical line 508 is introduced into the image at the boundary between the left and right lungs 504 and 506 , respectively . in step 404 the rows of the image in the right lung are summed to yield row sums the graph 511 of the function b j is shown in the image 500 adjacent to the right lung 506 . the top of the right lung is identified in step 406 as the highest row j for which b j exceeds a predetermined threshold value . a horizontal line 510 is then introduced into the image 500 at the top of the right lung in step 408 . the bottom of the right lung is identified in step 410 at the lowest row j for which b j exceeds a predetermined threshold value . a horizontal line 512 is then introduced into the image at the bottom of the right lung in step 410 . in step 412 , the rows of the image in the left lung are summed to yield row sums the graph 513 of the function c j is shown in the image 500 adjacent to the left lung 504 . the top of the left lung is identified in step 414 at the highest row j for which c j exceeds a predetermined threshold value . a horizontal line 514 is then introduced into the image at the top of the left lung in step 416 . the bottom of the left lung is identified in step 418 at the lowest row j for which c j exceeds a predetermined threshold value . a horizontal line 516 is then introduced into the image at the bottom of the left lung in step 420 . in step 422 the height of the right lung is calculated as the number of pixel rows in the image between the top and bottom of the right lung . in step 424 , the height of the right lung is divided by 3 and in step 426 , horizontal lines 520 and 522 are introduced into the image 500 so as to divide the right lung in the image into three regions , the right top rt , right middle rm and right bottom rb of equal height . in step 428 , the height of the left lung is calculated as the number of pixel rows in the image between the top and bottom in the left lung . in step 430 , the height of the left lung is divided by 3 , and in step 432 , horizontal lines 524 and 526 are introduced into the image so as to divide the left lung in the image into three regions the left top lt , left middle lm , and left bottom lb of equal height . now that the lungs in the image 500 have been divided into the six regions rt , rm , rb , lt , lm , and lb , the intensities of the pixels in each region are summed in step 434 . the sum for each region is a value of a regional assessment parameter for the region . in the case that the pixel intensities are calculated as disclosed in u . s . pat . no . 6 , 887 , 208 , the regional assessment that is obtained is indicative of the airflow in each region of the lungs . fig4 b shows the regional assessment of the same individual determined by vq scan . the image was divided into 6 regions and the fraction of radioactivity in each region was calculated , as is known in the art . the regional assessment of each region is shown in the region . fig6 a shows the regional assessment obtained on the second subject by the method of the invention , and fig6 b shows the regional assessment obtained on the second subject by vq scan .
the invention provides a system and method for regional assessment of lung physiology . the system includes a plurality of sound transducers configured to be fixed on a surface of the individual over the thorax . a processor is configured to receive signals generated by the transducers and to determine from the signals a value of a parameter in each of one or more regions of the lungs . the method of the invention includes obtaining signals indicative of pressure waves at locations over the thorax ; and determining from the signals a value of a parameter in each of the regions of the lungs .