key: cord-0010495-aua7i3j9 authors: Gómez-Laberge, Camille; Arnold, John H.; Wolf, Gerhard K. title: A Unified Approach for EIT Imaging of Regional Overdistension and Atelectasis in Acute Lung Injury date: 2012-01-10 journal: IEEE Trans Med Imaging DOI: 10.1109/tmi.2012.2183641 sha: c40b8df3f113795275355b7b20952a2d7e013bd6 doc_id: 10495 cord_uid: aua7i3j9 Patients with acute lung injury or acute respiratory distress syndrome (ALI/ARDS) are vulnerable to ventilator-induced lung injury. Although this syndrome affects the lung heterogeneously, mechanical ventilation is not guided by regional indicators of potential lung injury. We used electrical impedance tomography (EIT) to estimate the extent of regional lung overdistension and atelectasis during mechanical ventilation. Techniques for tidal breath detection, lung identification, and regional compliance estimation were combined with the Graz consensus on EIT lung imaging (GREIT) algorithm. Nine ALI/ARDS patients were monitored during stepwise increases and decreases in airway pressure. Our method detected individual breaths with 96.0% sensitivity and 97.6% specificity. The duration and volume of tidal breaths erred on average by 0.2 s and 5%, respectively. Respiratory system compliance from EIT and ventilator measurements had a correlation coefficient of 0.80. Stepwise increases in pressure could reverse atelectasis in 17% of the lung. At the highest pressures, 73% of the lung became overdistended. During stepwise decreases in pressure, previously-atelectatic regions remained open at sub-baseline pressures. We recommend that the proposed approach be used in collaborative research of EIT-guided ventilation strategies for ALI/ARDS. M ECHANICAL ventilation is an essential part of caring for critically-ill patients with acute lung injury (ALI) or its more severe form acute respiratory distress syndrome (ARDS). Unfortunately, it has become clear after two decades of research that mechanical ventilation itself can exacerbate lung injury and, thus, increase the chance of death [1] , [2] . Multicenter trials attempting to reduce the mortality among ALI/ARDS patients by strategic selection of ventilation parameters have led to mixed results [3] - [6] . Their goal was to find a lung protective tidal volume ( ) and positive end-expiratory pressure (PEEP) according to measurements that characterize the lung as a whole. Advances in lung imaging research are now calling into question the presumption that whole-lung measurements are sufficient for identifying lung protective ventilation settings for ALI/ARDS. Computed tomography (CT) and, more recently, electrical impedance tomography (EIT) show that the mechanical properties of the ALI/ARDS lung are, in fact, very heterogeneous [7] , [8] . Furthermore, the percentage of collapsed, but potentially recruitable, lung throughout the patient population is highly variable [9] . In some patients, typical ventilator settings can lead to regional lung overdistension and atelectasis, which correlate with subsequent lung injury [10] , [11] . Consistent with these findings, a meta-analysis combining the data from multicenter trials has shown that only the subgroup of patients with ARDS showed a reduction in mortality to higher levels of PEEP [12] . Everything considered, ALI/ARDS affects lung mechanics in a heterogeneous manner. Therefore, a reduction in mortality among patients will most likely require a ventilation strategy based on regional indicators of lung mechanics. In this paper, we consolidate various techniques from recent electrical impedance tomography (EIT) research into a unified method for imaging regional lung overdistension and atelectasis during mechanical ventilation of the ALI/ARDS patient. Our objective is to encourage multicenter collaboration by proposing a comprehensive software solution for all data processing stages: from data acquisition to physiological interpretation. The design criteria we have selected for this purpose are: currency of the techniques used, compatibility with most EIT systems, transparency of the software implementation, physiological interpretability of output, and automaticity of processing where possible. We implement and demonstrate a software solution according to these criteria, which is available to all investigators, and amenable to updates and improvements. We believe that the output from an EIT system utilizing such software would be physiologically interpretable by clinical staff and directly comparable with output from other centers. The solution is implemented on the software platform of the Electrical Impedance and Diffuse Optical Reconstruction Software (EIDORS) project. 1 EIDORS is maintained on the Internet by EIT researchers worldwide; it is platform-independent and compatible with most EIT systems. The data acquisition and image reconstruction algorithms used here conform to the Graz consensus on EIT lung imaging (GREIT), as described by Adler et al. [13] . GREIT provides the broadest consensus on EIT lung imaging techniques for electrode placement, stimulation and measurement patterns, as well as image reconstruction algorithms. Upon this framework, we develop an automatic filtering technique to attenuate nonrespiratory signals and to detect peak inspiratory and expiratory time points for the measurement of regional tidal volume. The corresponding filtered images are then used for lung region classification based on the functional EIT (fEIT) method of Pulletz et al. [14] . We embed the fEIT method within an algorithm to automatically determine the threshold parameter for lung classification and to aggregate lung regions identified over time. The tidal volume images and ventilator pressure measurements are combined into maps of regional respiratory system compliance. 2 We represent temporal changes in respiratory system compliance with a pair of images: the first quantifies the maximum compliance at each pixel, and the second indicates the ventilatory pressure when it was achieved. These results are used to estimate the percentage and location of lung overdistension and atelectasis using the method of Costa et al. [15] . We evaluate the method's performance using EIT data acquired from patients with ALI/ARDS who underwent a lung recruitment maneuver followed by a PEEP titration strategy within our pediatric intensive care unit. The proposed methods were used to monitor lung mechanics in pediatric patients with ALI/ARDS participating in a clinical study of lung recruitment (ClinicalTrials identifier: NCT00830284). The study protocol was approved by the hospital's Institutional Review Board, and written informed consent was obtained from a parent or guardian prior to enrollment. At the time of enrollment, the patients met the ALI/ARDS criteria for less than 72 h [16] . Throughout the study, the patients were either deeply sedated demonstrating apnea or received neuromuscular blocking agents. The time course of the experimental protocol is divided into a baseline ventilation stage, a lung recruitment stage where airway pressure was increased sequentially, and a PEEP Titration stage where airway pressure was decreased sequentially (Fig. 1 ). During the baseline stage, patients were ventilated in a volume-controlled (VC) mode ( ), where PEEP was set according to clinical practice. We also measured the difference between the plateau pressure ( ) and PEEP at the airway opening (1) During the lung recruitment stage, patients were switched to pressure-controlled (PC) ventilation ( cm ), and PEEP was increased by 5 cm with the objective of meeting a predefined blood gas criterion of lung recruitment [17] (2) where and are the partial pressures of oxygen and carbon dioxide in arterial blood, respectively. We used this recruitment index, since it was shown to consistently coincide with anatomical recruitment measured by CT in patients with Fig. 1 . The experimental protocol is divided in three stages. Baseline PEEP was determined according to clinical practice. During the "baseline" and "PEEP titration" stages, patients were ventilated in VC mode with , and was measured by subtracting PEEP from the measured . During the "lung recruitment" stage, patients were ventilated in PC mode, and was set to 15 cm for each step. VC: volume-controlled; : plateau pressure; PC: pressure-controlled. early ALI/ARDS [17] . The lung recruitment stage was stopped before reaching the index if PEEP reached 35 cm , or if any of the following stopping criteria were met: arterial oxygen saturation , arterial pH , arterial lactate , mean arterial blood pressure decreases by 20%, or the required dosage of vasoactive agents increases by 50%. During the PEEP Titration stage, patients were switched back to VC mode ( ), and PEEP was decreased to the lowest possible setting that satisfied (2) . In our approach, PEEP was initially set at 20 cm and was decreased in steps of 2 cm until one of the following: the minimum PEEP satisfying (2) was found, cm , or if any of the aforementioned stopping criteria were met. EIT measurements were performed using either the Goe-MF II (CareFusion, San Diego, CA) or the Dräger EIT Evaluation KIT 2 (Dräger Medical, Lübeck, Germany), based upon availability of the devices in the intensive care unit. Both devices are designed to operate with similar specifications for the purpose of single frequency EIT. Stimulation current, electrode impedances and voltages are comparable in both devices. The image reconstruction and processing routines proposed here were applied to the raw voltage data stream from both systems. 1) Electrode Placement: Sixteen coplanar electrodes were placed equidistantly around the thorax at the level of the parasternal sixth intercostal space. The reference electrode was placed on the right side of the abdomen near the waistline. Electrodes #1 and #16 were symmetrically placed to the left and the right of the sternum, respectively, so that electrodes #8 and #9 straddled the spinal column. This configuration leads to transverse images in the radiological convention. Each image was reconstructed from voltage measurements acquired with adjacent current stimulation [13] (50 kHz, 5 mA rms), at a frame rate of 20 Hz. To maintain good signal quality, we ensured that electrode-skin resistance did not exceed 200 at any time. 2) Image Reconstruction: Image reconstruction was performed using the EIDORS implementation of the one-step linear Gauss-Newton solver with a regularization term for conductivity change and another accounting for electrode movement [13] . We calculated changes in regional impedance according to the normalized difference between electrode voltages at each time point and a reference , which was obtained by averaging voltage data corresponding to peak-expiration. We chose this reference because it would represent all impedance changes as positive quantities for all lung volumes. Thus, for each measurement, we have the nonnegative quantity . In the GREIT notation, changes in all of the voltage measurements and electrode positions are written in vector form . These can be calculated by solving the forward problem, where is obtained from a linear operation on the conductivity distribution over the regional elements representing the thorax. The linear model is (3) where is the Jacobian matrix for changes in conductivity and electrode position [18] , and is additive white noise. Each element of quantifies the change in regional conductivity. In later sections on physiological interpretation, we will use the reciprocal of , called the regional impedance , since it is more commonly used in the clinical literature on EIT. Furthermore, has been shown to correlate with changes in regional tidal volume [19] - [21] . The estimate of is obtained by solving the inverse problem as follows [22] : (4) where indicates the matrix transpose, and and are the regularization parameter and matrix, respectively. We tuned so that the amount of noise amplification remained constant between patient data sets. This was done for each data set by setting the noise figure to 0.60 [23] . The finite element models used for the forward and inverse problems are labeled d2t2 and d2t3, respectively, in the EI-DORS project. They are 2-D meshes with point electrodes that approximate the contour of the human thorax passing through the electrode plane. We reconstructed an image as described above for every time point during EIT scanning in order to obtain a image series during mechanical ventilation. The following methods pertain to calculations using the image series acquired during each step of the experimental protocol. 3) Intra-Tidal Impedance Changes: The image series was used to detect the time points of peak-inspiration and peak-expiration. An average time signal was computed by averaging all pixels within each image. The time signal was bandpass-filtered in order to isolate the respiratory component. The lower corner frequency (attenuating the signal by 3 dB) was . The choice for the upper corner frequency was determined automatically, since it depends on the respiratory rate. We used the magnitude of the Fourier-transformed signal to identify the frequency with the largest spectral peak greater than and set . We used the implementation of a 128-point zero-phase finite impulse response filter by calling the function filtfilt from the Matlab Signal Processing Toolbox. The rationale for using a zero-phase filtering technique was to preserve the accurate peak-inspiration and peak-expiration times by canceling the linear-phase shift associated with digital filters. These time points were used to estimate regional ventilation by calculating the average difference in between peak-inspiration and peak-expiration. Note that is insensitive to aerated but nonventilated lung regions. The lung region of interest (ROI) was determined by applying a threshold to the functional EIT image (fEIT) computed from an image series during tidal ventilation [14] . The pixel intensity of the fEIT image represents the standard deviation of during tidal breathing in the bandpass-filtered image series. Pixels are considered part of the lung ROI if they surpass a threshold value representing a percentage of the maximum standard deviation in the fEIT image. For each image series, we determined the value for automatically by searching for the ROI with the smallest number of inconsistently classified pixels when compared to ROI obtained for and . The ROI templates were then merged by the logical "or" operation, yielding the aggregate ROI, that we used for calculating regional respiratory system compliance, atelectasis, and overdistension. Changes in dynamic respiratory system compliance were estimated for each pixel in the EIT image by the ratio (5) We used changes in as a function of PEEP to estimate the magnitude of overdistension and atelectasis [15] . During each stage of the protocol, dynamic compliance was computed for each PEEP setting , and we determined each pixel's maximum compliance and the PEEP setting when it was achieved (6) (7) where returns the set of pressures that maximize the compliance signal for each pixel in the lung ROI. 6) Regional Overdistension and Reversible Atelectasis Estimation: The fraction of overdistended or atelectatic units in a given pixel at PEEP setting can be obtained in terms of and as follows [15] : This function is bounded between . For pressures above , a pixel within the lung ROI represents lung tissue that is overdistended. For pressures below , the corresponding lung tissue is atelectatic. We apply (8) to the EIT image pixel data acquired for every step of the protocol. We present our data in two formats: as images representing the regional distribution of lung overdistension and reversible atelectasis for each patient, and as graphs aggregating Corresponding intra-tidal difference image, a correlate of regional lung ventilation, obtained by subtracting the peak expiratory image from the peak inspiratory image. The colormap applies to all images and represents increases (cool colors) and decreases (warm colors) relative to an average peak-expiratory reference image. Images represent the transverse view in the radiological convention (ventral: image top; right side: image left). all patient data as a -weighted pixel average that is plotted versus throughout the entire lung ROI. Isolation of the respiratory signal with zero-phase filtering led to the accurate detection of peak inspiratory and expiratory time points. For each step of the protocol, 15 complete breaths were recorded on average. The total number of breaths recorded in all patients was 1087, of which 1044 (96.0%) were detected and 43 (4.0%) were missed. There were also 27 (2.4%) false detections, where we counted each inter-breath period as a true negative. Therefore, the sensitivity and specificity of the breath detection algorithm was 96.0% and 97.6%, respectively. Incorrect detections were either manually adjusted or discarded for further analysis. The filtered image series during a typical tidal breath is shown along with its corresponding intra-tidal difference image (Fig. 2) . The intra-tidal difference image estimates the regional tidal volume and is obtained by subtracting the peak expiratory image from the peak inspiratory image. We estimated the error in intra-tidal timing and regional volume by comparing the automatically detected time points against those Fig. 3 . Aggregation of the lung ROI template during the baseline (top row) and lung recruitment stages. The fEIT images from the filtered image series (left column) demonstrate how the signal locus shifts from ventral to dorsal lung regions. The corresponding ROI templates (right column) represent nonlung (gray) and lung (black) regions, where the standard deviation (SD) threshold was determined automatically. The aggregate ROI (bottom right) is obtained by applying the logical "or" operation on all ROI templates (black represents "true" and gray "false"). identified manually from the raw EIT signal. From each patient, we calculated the error for peak-inspiratory and peak-expiratory time points from five randomly-selected breaths. From these breaths, we also calculated the error in estimated tidal volume . The error for peak-inspiratory and peak-expiratory time point detection was and , respectively ( , ). Neither were significantly different from zero (paired -test, ). The larger error of the peak-expiratory time point was expected, since the expiratory phase of the signal is flat for at least 1/2 s; therefore, the time point of the minimum is susceptible to noise [ Fig. 2(a) ]. These timing errors led to the underestimation of tidal volume by ( , ). The intra-tidal impedance changes were detected after bandpass filtering. Throughout the protocol, in all patients, the respiratory component was adequately isolated. The mean and its standard error for the respiratory rate throughout the protocol was breaths/min ( ). The heart rate was beats/min ( ). The filter corner frequencies were and , respectively. Hence, on average, the respiratory signal was well centered in the passband, and the cardiac signal was attenuated by more than 3 dB. The fEIT image was computed on the bandpass-filtered image series in order to isolate the respiratory signal before lung ROI detection. A fEIT image for each step of the protocol, and the aggregate lung ROI was then obtained by performing the logical "OR" on all step-specific ROI (Fig. 3) . The parameter represents the threshold of standard deviation relative to the image maximum, above which pixels are included in the lung ROI. The recommended range for manual ROI selection is , [14] . The automatically-determined threshold had a mean of 0.24 and median 0.20 and, thus, was in agreement with manually selected lung ROI. Of the 73 set of data collected from all patients (one set per protocol step), 32 (44%) were below the manual range, 27 (37%) were within the range, and 14 (19%) were above the range. For each patient, we present the changes in regional respiratory system compliance separately for the lung recruitment and PEEP Titration stages. In order to display the maximum regional compliance at every pixel and the pressure at which this compliance was achieved, we use a pair of images labeled and whose pixel values correspond to the quantities satisfying (6) and (7), respectively (Fig. 4) . Each row of Fig. 4 shows the regional compliance data for one patient. The first pair of images shows the maximum regional compliance at each pixel (i.e., map) and the corresponding plateau pressure needed for each pixel to reach its maximum (i.e., map) during lung recruitment. The second pair of images shows the and maps in the same format for the PEEP Titration stage with the exception that the map quantifies PEEP. Because the maximum respiratory system compliance achieved by each patient was highly variable, ranging from 0.1-0.9 cm , the color scale was set according to the individual. Regions of high compliance are bright, while those of low compliance are dark. For each patient, the same color scale was used for both stages. The map uses dark colors to show regions that reached maximum compliance at low pressures, and bright colors for regions that reached maximum compliance at high pressures. Again, because each patient ended either stage at different pressures, the map color scale was set according to the individual. Patient #4 met safety-related stopping criteria during lung recruitment and, therefore, did not undergo PEEP titration. These results suggest that a considerable degree of heterogeneity exists throughout the lungs of each patient. All maps contain a patchy distribution of regions of low and high compliance. In most cases, compliance throughout the whole lung is not maximized at a single pressure as shown by the mosaic pattern of the maps. This heterogeneity of compliance does not change significantly during the protocol, since the maps obtained from the lung recruitment and PEEP Titration stages are very similar. There is virtually no difference in the lung ROI between stages, and the locations of high compliance also do not change remarkably. A comparison of the maps for each patient shows that the dorsal region compliance values reached during the lung recruitment stage are maintained at much lower pressures during PEEP Titration stage. Similar results exhibiting this form of lung volume hysteresis have been shown in spirometry studies of the ARDS lung [24] . In order to examine how the lung mechanics evolved during the lung recruitment and PEEP Titration stages, we pooled the pressure and compliance data from all patients for each step of the protocol (Table I) . The fourth column of Table I shows the average over the entire lung ROI. As a result of increasing plateau pressure by 23 cm , the respiratory system compliance dropped from a baseline value of 0.3 cm to 0.1 cm during the lung recruitment stage. Upon transition to the first step of the PEEP Titration stage ("T1" in Table I ), the effects of hysteresis were already apparent. Respiratory system compliance returned to the baseline level at 0.3 cm , effectively tripling the compliance measured at comparable pressure during step "R2" of lung recruitment. Reduction in pressure during PEEP Titration further increased compliance, which exceeded baseline and reached a maximum of 0.5 cm at step "T4." Reduction in pressure below values in step "T4" appeared to reduce compliance to levels smaller than baseline. The last column of Table I shows the percentage of pixels within the lung that met their maximum compliance during each step of the protocol. While pressure increased during the lung recruitment stage, the fraction of pixels at their maximum compliance fell from 52% to as low as 3% at step "R3." This effect was reversed during the PEEP Titration stage, where as much as 46% of pixels became maximally compliant at step "T4." This extent of maximal regional compliance was not sustained with the lower pressures in steps "T5-T8." Since our ring of EIT electrodes were closest to the basal lung regions, we asked whether our measurement of compliance is representative of the whole lung in ALI/ARDS patients. We compared the changes in respiratory system compliance measured by EIT to those measured by the ventilator by calculating the Pearson product-moment correlation coefficient between the time series of measurements from each instrument (Table II) . From our nine patients, we found that six had very strongly correlated signals that reached statistical significance ( , ). In two other patients, the correlations were strong but only showed a trend toward significance ( , , and , ). In contrast, the EIT data from patient #8 were weakly correlated ( , Fig. 4 . Maximum regional respiratory system compliance ( maps) and the corresponding ventilatory pressure ( maps) are shown for lung recruitment and PEEP Titration stages of each patient. The maps use a scale relative to each patient's maximum compliance achieved throughout the entire protocol. The maps represent the pressure at which each pixel's maximum compliance was achieved. For the recruitment stage, represents , and for the titration stage, represents PEEP. Bright colors represent maxima achieved at high pressures, while dark colors represent maxima achieved at low pressures. ). While this case suggests that EIT was rather insensitive to changes in the gross lung, we cannot exclude the possibility that patient manipulations during the protocol may have inadvertently interfered with the EIT measurements. The fraction of overdistended and atelectatic lung regions for each step of the protocol were obtained in terms of and according to (8) . The maps for regional overdistension and atelectasis are demonstrated for one patient in Fig. 5 . The color of each pixel within the lung corresponds to the value from (8) . Each image answers the question of "how much" and "where" overdistension (blue) and atelectasis (red) takes place during each step of the protocol. Plateau pressure and PEEP are shown above each image, and pixel-averaged overdistension and atelectasis are shown below. In this patient, overdistension took place in as much as 82% of the lung ROI during the lung recruitment stage. The most severe overdistension was in the ventral lung regions. Interestingly, during the fourth and final step of lung recruitment, overdistension fell back to 78%, suggesting the aeration of novel regions. The effect of this regional opening can be seen as a retreat of the interface between ventral overdistension (blue pixels) and adequate aeration (black pixels). Before commencing the PEEP Titration stage, the nondependent lung is almost completely overdistended. The benefit of these high pressures can be seen in the reversal of atelectasis (red pixels) in the dependent lung. A cm was necessary to reverse atelectasis in 5% of the lung down to 3%. At the end of the protocol, atelectasis was confined to 3% of the lung utilizing cm and cm . This was less than half of the pressure required to achieve the same degree of atelectasis reversal during the lung recruitment stage. The mean and standard error of pixel-averaged overdistension and atelectasis for all patients are shown in aggregate during lung recruitment and PEEP Titration in Fig. 6 . During the lung recruitment stage, the extent of lung overdistension increased monotonically from to with climbing from 30 to 50 cm [ Fig. 6(a) ]. This recruitment strategy resulted in the resolution of atelectasis in of the lung [ Fig. 6(a) ]. During the PEEP Titration stage, while PEEP was reduced from 20 to 14 cm , overdistension was maintained within approximately 20% of the lung, and atelectasis was maintained within 10% of the lung [ Fig. 6(b) ]. The In summary, stepwise increases in pressure could reverse atelectasis within 17% of the lung, appearing primarily in the dorsal regions. During the highest pressures applied, 73% of the lung became overdistended, appearing primarily in the ventral regions. Overdistension subsided during stepwise decreases in pressure, and previously-atelectatic regions remained open at pressures below those at baseline. Together, these results suggest that the described ventilation strategy can effectively recruit the lung, and that EIT imaging can be used to guide the reduction of ventilatory pressures to avoid re-collapse of the lung. Building upon a decade of experimental validation of EIT regional lung volume measurements [19] , [25] - [29] , there has been a recent surge in EIT research about mechanical ventilation of the ALI/ARDS lung [15] , [30] - [35] . These studies provide important insights into how the EIT data stream can be interpreted physiologically. Most recently, the following methods have been proposed: an intra-tidal respiratory system compliance index [15] , fEIT changes between PEEP settings [30] , a respiratory system compliance percentile change from a scaled baseline value [31] , a regional lung state classification index based on changes between PEEP settings [32] , the dorsoventral center-of-mass coordinate of a spectral power EIT image [33] , an intra-tidal global ventilation inhomogeneity index [34] , and ventilator-calibrated regional distributions of end-expiratory lung volume and potentially-recruitable lung volume [35] . In addition to this variety of proposed indexes of lung mechanics, each assumes that a reconstruction algorithm has already computed the images from the voltage measurements. However, it is difficult to compare results obtained from arbitrary reconstruction algorithms, since the ill-posed nature of the EIT conductivity problem can cause different algorithms to generate rather divergent images from the same data [36] . In this paper, we combined various EIT techniques to form a complete sequence for interactively quantifying regional lung overdistension and atelectasis from EIT voltage and airway pressure data. The methods were chosen according to criteria that we deemed essential for developing EIT-guided mechanical ventilation within the critical care setting. The lung imaging methods conform to the GREIT specifications [13] , which are recognized by a large collaborative network of scientific and clinical experts as being the state-of-the-art in EIT lung imaging. Image reconstruction was essentially automatic, since only the regularization parameter occasionally required adjustments when working with a new patient. The subsequent steps performed individual breath detection, lung ROI classification [14] , and regional lung mechanics analysis [15] . We reasoned that this approach would be helpful in the intensive care unit, since very few parameter adjustments are required, and the variables quantifying lung mechanics are familiar to the clinician. Also, since our estimates of regional tidal volume are obtained from the average difference between images at end-expiration and end-inspiration, we do not anticipate a bias due to changes in pulmonary perfusion. The EIT-based findings in our study of lung mechanics in ALI/ARDS are consistent with those from previous research using other modalities. First, lung recruitment studies using CT imaging have shown a heterogeneous distribution of lung density in ARDS [7] , where the dependent lung regions are predominantly collapsed. The respiratory system compliance data in our study also suggest that the ALI/ARDS lung in our patients consisted of a patchy distribution of high and low compliance regions, each requiring different airway pressures for maximizing regional compliance. The patients in our sample also exhibited a trend of progressive lung recruitment in the ventral-to-dorsal direction, supporting the proposition that dorsal regions are subject to the gravitational force from ventral regions and, therefore, require additional pressure to inflate. Second, hysteresis of the ALI/ARDS lung has been previously studied by comparing static pressure-volume data acquired during inflation and deflation [24] . We found that the regional and total respiratory system compliance data also indicated the presence of lung volume hysteresis, where improved compliance was achieved during the PEEP Titration stage at lower pressures than those used during baseline and the lung recruitment stage. Finally, the correlations between changes in respiratory system compliance measured by EIT and by the ventilator suggest that the anatomical landmarks for electrode placement are appropriate for monitoring the regional lung mechanics in patients with ALI/ARDS. Although lung CT imaging is the gold standard for anatomical recruitment, such measures were unjustified in our study. The risks associated with the transportation of critically-ill patients out of the intensive care unit as well as those from ionizing radiation exposure during childhood were unacceptable. The proposed method, of course, has important limitations. First, the method is limited to detecting reversible atelectasis, since EIT image reconstruction is currently best-suited to measure impedance changes rather than absolute measurements. Consequently, only lung regions that exhibit some airflow during the ventilation protocol can be identified. We also assume that the transverse area occupied by the lungs does not change appreciably during recruitment. Although CT studies of lung recruitment support this assumption [7] , [10] , [17] , its falsehood in a special case could affect the estimates of overdistension and atelectasis. Second, the processing steps following image reconstruction are relatively new and need further validation involving a broader range of patients with ALI/ARDS and different EIT systems. Making the programming sequence available to all centers is perhaps the most effective way to evaluate and optimize EIT-guided mechanical ventilation strategies for patients in respiratory distress. Although advanced alternatives to fEIT thresholding have been proposed for lung detection, e.g., least-squares regression [37] (supplement), [33] , principal component analysis [38] , and fuzzy logic pattern recognition [39] , we reasoned that standard fEIT thresholding would serve as a good benchmark because its limitations are well understood, and it will challenge new methods remain as simple, automatic and fast as possible. The calculations for regional lung mechanics have only been demonstrated in two case studies of ARDS in adults [15] and, therefore, will also require further evaluation in both adults and children. Third, the method will systematically underestimate the "true" static respiratory system compliance, which is best obtained from the quasi-static pressure-volume curve, where pressure and volume measurements are made during the best practical approximation of zero flow [24] . Quasi-static pressure-volume maneuvers are unfortunately difficult to perform and require ventilator disconnect, making the risk-to-benefit ratio for ALI/ARDS patients unattractive. It is clear from the epidemiological literature that a multicenter trial will be necessary to obtain a sufficient number patients to power a clinical trial of EIT-guided ventilation in ALI/ARDS. Therefore, the aforementioned challenges in data compatibility must be overcome before EIT-guided ventilation can be tested in terms of clinical outcomes, in particular, measures of morbidity and mortality. We recommend that the proposed approach be used as a framework for evaluating and improving each step of the sequence with the purpose of interactively guiding mechanical ventilation of the ALI/ARDS patient with EIT. Ventilator-induced lung injury: From the bench to the bedside Ventilator-induced lung injury: The anatomical and physiological framework Ventilation with lower tidal volumes as compared with traditional tidal volumes for acute lung injury and the acute respiratory distress syndrome Higher versus lower positive end-expiratory pressures in patients with the acute respiratory distress syndrome Ventilation strategy using low tidal volumes, recruitment maneuvers, and high positive endexpiratory pressure for acute lung injury and acute respiratory distress syndrome: A randomized conrolled trial Positive endexpiratory pressure setting in adults with acute lung injury and acute respiratory distress syndrome: A randomized controlled trial What has computed tomography taught us about the acute respiratory distress syndrome? Regional lung volume changes in children with acute respiratory distress syndrome during a derecruitment maneuver Lung recruitment in patients with the acute respiratory distress syndrome Tidal hyperinflation during low tidal volume ventilation in acute respiratory distress syndrome Lung opening and closing during ventilation of acute respiratory distress syndrome Higher vs lower positive end-expiratory pressure in patients with acute lung injury and acute respiratory distress syndrome: Systematic review and meta-analysis GREIT: A unified approach to 2D linear EIT reconstruction of lung images Comparison of different methods to define regions of interest for evaluation of regional lung ventilation by EIT Bedside estimation of recruitable alveolar collapse and hyperdistension by electrical impedance tomography The American-European consensus conference on ARDS. Definitions, mechanims, relevant outcomes, and clinical trial coordination Reversibility of lung collapse and hypoxemia in early acute respiratory distress syndrome Direct EIT Jacobian calculations for conductivity change and electrode movement Detection of local lung air content by electrical impedance tomography compared with electron beam CT Regional ventilation by electrical impedance tomography Electrical impedance tomography: A method for monitoring regional lung aeration and tidal volume distribution? Imaging of conductivity changes and electrode movement in EIT Electrical impedance tomography: Regularized imaging and contrast detection Pressure-volume curves of the respiratory system Imbalances in regional lung ventilation: A validation study on electrical impedance tomography Noninvasive assessment of lung volume: Respiratory inductance plethysmography and electrical impedance tomography Assessment of regional lung recruitment and derecruitment during a PEEP trial based on electrical impedance tomography Electrical impedance tomography compared with thoracic computed tomography during a slow inflation maneuver in experimental models of lung injury Regional lung volume changes during high-frequency oscillatory ventilation Bedside measurement of changes in lung impedance to monitor alveolar ventilation in dependent and non-dependent parts by electrical impedance tomography during a positive end-expiratory pressure trial in mechanically ventilated intensive care unit patients Regional tidal ventilation and compliance during a stepwise vital capacity manoeuvre Towards lung EIT image segmentation: Automatic classification of lung tissue state from analysis of EIT monitored recruitment manoeuvres Regional overdistension identified with electrical impedance tomography in the perflubron-treated lung Peep titration guided by ventilation homogeneity: A feasibility study using electrical impedance tomography A new non-radiological method to assess potential lung recruitability: A pilot study in ali patients Uses and abuses of EIDORS: An extensible software base for EIT Evaluation of an electrical impedance tomography-based global inhomogeneity index for pulmonary ventilation distribution Different approaches for quantifying ventilation distribution and lung tissue properties by functional EIT Data-driven classification of ventilated lung tissues using electrical impedance tomography The authors would like to thank A. Adler for helpful discussions and for having reviewed an earlier version of the manuscript.