This article has Open Peer Review reports available.

# Automated measurement of heterogeneity in CT images of healthy and diseased rat lungs using variogram analysis of an octree decomposition

- Richard E Jacob
^{1}Email author and - James P Carson
^{1}

**14**:1

https://doi.org/10.1186/1471-2342-14-1

© Jacob and Carson; licensee BioMed Central Ltd. 2014

**Received: **7 January 2013

**Accepted: **18 December 2013

**Published: **6 January 2014

## Abstract

### Background

Assessing heterogeneity in lung images can be an important diagnosis tool. We present a novel and objective method for assessing lung damage in a rat model of emphysema. We combined a three-dimensional (3D) computer graphics method–octree decomposition–with a geostatistics-based approach for assessing spatial relationships–the variogram–to evaluate disease in 3D computed tomography (CT) image volumes.

### Methods

Male, Sprague-Dawley rats were dosed intratracheally with saline (control), or with elastase dissolved in saline to either the whole lung (for mild, global disease) or a single lobe (for severe, local disease). Gated 3D micro-CT images were acquired on the lungs of all rats at end expiration. Images were masked, and octree decomposition was performed on the images to reduce the lungs to homogeneous blocks of 2 × 2 × 2, 4 × 4 × 4, and 8 × 8 × 8 voxels. To focus on lung parenchyma, small blocks were ignored because they primarily defined boundaries and vascular features, and the spatial variance between all pairs of the 8 × 8 × 8 blocks was calculated as the square of the difference of signal intensity. Variograms–graphs of distance vs. variance–were constructed, and results of a least-squares-fit were compared. The robustness of the approach was tested on images prepared with various filtering protocols. Statistical assessment of the similarity of the three control rats was made with a Kruskal-Wallis rank sum test. A Mann-Whitney-Wilcoxon rank sum test was used to measure statistical distinction between individuals. For comparison with the variogram results, the coefficient of variation and the emphysema index were also calculated for all rats.

### Results

Variogram analysis showed that the control rats were statistically indistinct (p = 0.12), but there were significant differences between control, mild global disease, and severe local disease groups (p < 0.0001). A heterogeneity index was calculated to describe the difference of an individual variogram from the control average. This metric also showed clear separation between dose groups. The coefficient of variation and the emphysema index, on the other hand, did not separate groups.

### Conclusion

These results suggest the octree decomposition and variogram analysis approach may be a rapid, non-subjective, and sensitive imaging-based biomarker for characterizing lung disease.

### Keywords

Lung imaging Disease detection COPD Emphysema Pulmonary Octree Variogram## Background

Emphysema is an obstructive pulmonary disease that results in airway expansion and tissue destruction. Early and accurate detection of emphysema is important for disease management and improved patient outcomes [1]. Emphysema typically results in heterogeneous air trapping and increased ventilation-perfusion inequality. Thus, signatures of tissue heterogeneity or regional air trapping may facilitate earlier and/or more accurate diagnoses of emphysema [2, 3].

It has been shown in many studies that CT is a reproducible and predictive modality for diagnosing and assessing emphysema [2, 4–7]. The most common way to quantify moderate to severe emphysema from CT images is to measure the emphysema index: the percentage of voxels below a preset Hounsfield Unit (HU) threshold. However, this index may not detect subtle or early onset emphysema [7], or, conversely, it may identify “false positive” regions that appear emphysematous in asymptomatic healthy young subjects with no smoking history [8]. In addition, threshold values vary between studies, typically within a range of-900 to-980 HU, with factors other than disease influencing CT densitometry [2, 4, 9, 10].

An approach to assessing heterogeneity that has been applied to myriad types of medical images is fractal analysis [11]. This approach exploits the scale-independent nature of fractals in systems, such as lung vasculature, whose variation in form or regularity is thought to be similar through different degrees of magnification [12]. Typically, a region (or regions) of interest (ROI) is selected, which is then evaluated for a mean value and/or repeatedly subdivided to relate differences in signal intensity across space and across ROI size. If the relationship is characterized by a power law, then the exponent, or fractal dimension, is taken as an indication of complexity, texture, or heterogeneity [13–15]. However, this method is based on an a priori assumption that power law relationship exists and that a fractal dimension can be determined. Furthermore, because this approach requires ROI selection, it is often subjective and can result in omission of large sections of the lung [13, 16].

A promising new approach, detailed in recent work by Subramaniam et al. [17], demonstrates the use of analysis of CT image slices using a quadtree decomposition. This approach iteratively divides a 2D image into increasingly small squares. Division occurs whenever the range of intensities within the quadrant are above a pre-defined threshold. In this way, regions with homogenous intensity are characterized by larger squares. Their work measures heterogeneity by counting the number of squares-per-area in the lung or in a local region of the lung.

We propose to extend the quadtree decomposition to 3D by using octrees. In this way, the entire lung CT image volume is used. An octree iteratively divides a cubic volume into eight evenly sized cubes, or octants. The octree method is well-developed and is used, for example, in 3D mesh generation and computer graphics [18, 19], medical image registration [20], and finite element meshing in CT images [21]. The advantages of the octree here is that it facilitates a speed-up in 3D lung tissue analysis by non-subjectively subdividing the lungs into homogeneous regions, thereby allowing for focus on parenchyma by reducing partial volume averaging and eliminating edges [22].

To incorporate spatial information across the lung, we propose to couple the octree image decomposition approach with variogram analysis. Variograms are a spatio-statistical approach for measuring spatial variability by comparing sample value variances to the distance of separation [23, 24]. Variograms are well-established in geostatistics, but have also recently attracted interest in biomedical applications, such as characterizing magnetic resonance images of white matter [25]. In this pilot study using our approach, we are able to significantly differentiate from each other the three groups of elastase-dosed rats: control, distributed mild emphysematous disease, and region-specific severe disease.

## Methods

### Disease model

An elastase-induced model of emphysema was used in this study. Nine male Sprague-Dawley rats with an average weight of 212 ± 11 g were orally intubated and dosed intratracheally with: 250 U/kg elastase dissolved in 200 μL saline to the whole lung (n = 3), or 50 U/kg elastase in 200 μL saline to a single lobe (n = 3), or 200 μL saline as a control (n = 3). Dosing levels were based on our previous work in which emphysematous changes were detected using ^{3}He diffusion MRI and histology [26]. All animal use was approved by the Institutional Animal Care and Use Committee at Pacific Northwest National Laboratory.

### CT imaging

Three weeks following dosing, the rats were imaged using micro-CT. At this time, rats weighed 357 ± 10 g. The imaging procedure is described in detail in [27]. Briefly, rats were anesthetized, intubated, and mechanically ventilated at 1 Hz with 40% inhale and 60% exhale durations. Peak inhalation pressure was ~8 cmH_{2}O, and no peak end expiratory pressure was used so that images could be acquired at functional residual capacity. Anesthesia was maintained by providing 3-4% isoflurane in air (30% O_{2}, balance N_{2}). Sigh breaths were delivered periodically to maintain lung recruitment. A respiratory-gated GE eXplore 120 micro-CT scanner was employed with the following settings: 100 kV peak voltage, 50 mA tube current, 16 ms exposure time, and 360 projections with 1 degree angular steps. Images were reconstructed with supplied software to 150 μm isotropic resolution. Total imaging time was about 90 minutes due to the collection of multiple images throughout the breathing cycle; however, only the images acquired at full exhalation were used for the analyses herein. Post-mortem histology and tissue analysis were not possible due to *in situ* lung casting for other research purposes (i.e. supplying airway tree geometries for computational fluid dynamics models [28, 29]).

### Image preparation

A lung mask image was semi-automatically generated from the above-mentioned reconstructed images using ImageJ [30] and the ImageJ 3D Toolkit plug-in [31]. Starting from a seed-point inside the lung, a 3D connected threshold was applied with a threshold value empirically selected to exclude major vasculature and all external tissue, so that only lung tissue would be included in analysis. Mask boundaries were smoothed using the Region Dilate and Region Erode functions in succession. Applying the mask to the reconstructed images, all background was thus assigned an intensity value of 0, and all unmasked lung tissue retained original HU values. A five pixel diameter 3D median filter was then applied in order to reduce noise while contributing minimal blurring. Finally, the image canvas size was increased to 512 × 512 × 512 by zero-filling in each direction. We note that zero-filling to a power of 2 served to conveniently restrict all octants in the octree decomposition to isotropic cubes.

### Octree decomposition

There are multiple approaches to selecting the aforementioned user-defined threshold criteria. For example, in [17] boxes were subdivided if the HU range in the box exceeded 10% of the range of values in the entire image. Because our rats had different ranges of HU values depending on disease severity, we chose to use the standard deviation as the threshold criterion, as it does not require normalization. We describe how the threshold range was determined in Section Threshold range.

### Variogram analysis

The purpose of using variograms in this context was to determine the spatial relationships between regions of the lung that are relatively homogeneous. The octree decomposition isolated regions of relatively high homogeneity; therefore, each resulting box was treated as a homogeneous “voxel” using the mean value of the box as the voxel intensity. Furthermore, it was important to use octree boxes that were all the same size. The largest box sizes that resulted from the octree decompositions were mostly 8 × 8 × 8 with a few 16 × 16 × 16 boxes. In order to ensure a sufficient number of boxes that were distributed throughout the lung and to enhance statistical power [25], each 16 × 16 × 16 box was further broken down into 8 × 8 × 8 boxes, and all 8 × 8 × 8 boxes were used in the variogram analysis. This approach assumed that emphysematous disease was manifest on a scale greater than an 8 × 8 × 8 box, which for these images is 1.2 mm on a side–this is consistent with what we observed in this disease model in previous work [26]. Importantly, excluding the smaller boxes largely eliminated vasculature, conducting airways, boundaries, and other features that had high spatial variability [22].

*d*) of the differences in signal intensity

*I*is described in detail in many geostatistics texts, such as [23, 24], and is given by:

*d*is the distance between voxel centroids,

*x*is the voxel location in 3D, and

*N(d)*is the number of voxel pairs for a given

*d*. γ(

*d*) is calculated for all values of

*d*. For an image with

*n*voxels, the total number of voxel pairs (for all existing values of

*d*) is given by

We then calculated the average semi-variance at each distance. Results were plotted on a distance vs. variance graph, or variogram, which graphically represents the spatial dissimilarity within the image.

The purpose of the variograms was to characterize spatial relationships between boxes of relatively high homogeneity that were dispersed throughout the lung. However, as described by [25] in work with brain images, increased distances between voxels, or octants, lead to decreased likelihood that the voxels are related in any way. This is especially true in the lung, which is composed of largely independent lobes (five in the rat). Moreover, the spatial relationship between parenchymal signal intensity in different lobes is, in many contexts, most accurately described by a unique path up and back down the airway tree, potentially covering 20 or more orders of branching [33]. The spatial distances are analogous to the “regionalized variable” in the geological context of mineral distribution. Indeed, the regionalized variable, or the maximum distance *d*
_{
max
} over which the variogram is expected to be reliable, is estimated to be one half of the diameter of a region [25], while beyond *d*
_{
max
} any relationships between distance and variance are expected to be random. Therefore, we limited our variogram analysis to a region that represented half of the average dimension of the left lobe, which is typically the largest lobe in the rat. This assumes that the disease varies slowly compared to the box size but rapidly compared to the lobe size. Using the 3D images and a lung cast, *d*
_{
max
} was measured to be ≈ 8 mm (53 pixels) in the healthy rats, which was rounded up to the equivalent of seven face-bordering 8 × 8 × 8 boxes, or 56 pixels.

### Octree decomposition tests

We performed the following evaluations of the appropriateness and robustness of the octree and variogram approach.

#### Threshold range

The octree decomposition compared the standard deviation of the lung signal within a box to a threshold range to determine whether the box should be subdivided. The optimal threshold range was determined semi-empirically. As a starting point, we calculated the typical standard deviation of the lung tissue of the three control rats. This was accomplished by producing a histogram of each masked image and then fitting the main lung peak to a Gaussian curve. The mean of the control rats’ standard deviations (σ) was set as the initial threshold range. We then tested the effects of thresholds that were σ/3, 2σ/3, 4σ/3, and 2σ by rerunning the decomposition on the same images. We evaluated the effectiveness of each threshold range at grouping the control rats and at distinguishing the full-lung-dose rats from the control group in the variograms.

#### Image filtering

We examined the effects of image filtering by generating variograms from unfiltered images, images filtered with a five pixel diameter 3D median filter, images filtered with a 3D Gaussian blur of radius = 2, and images filtered with a 3D Gaussian blur of radius = 4. Filters were applied in ImageJ prior to applying the mask. The number of octree boxes and differences in the variograms were compared. Adjustments to the threshold range were made for each filter test based on the standard deviation of the filtered image.

#### Image translation

Rat positioning during imaging can vary from animal to animal, and the octree decomposition should be independent of this. Therefore, we tested the effects of image translation on the variogram results. For the variogram analysis, the largest unit box retained after decomposition was 8 × 8 × 8; therefore, a shift in the image by 8 pixels in any (or all) dimensions should result in no change to the resulting variogram, since the original image had isotropic dimensions of 2^{n}. Conversely, maximum change would occur from a four-pixel shift. To test the significance of the effects of translation, the image of one rat was shifted in x, y, and z by four pixels and by eight pixels using ImageJ, and variogram results were compared to those of the original image.

#### Image rotation

Animal positioning can also have an apparent effect on image rotation, and octree decomposition should be insensitive to this. Using the same rat image as used in the translation test, we imitated an arbitrary 3D rotation by rotating the image π/4 radians about the x, y, and z axes using ImageJ. Results were compared to the original image and to the results of the translation test.

#### Image downsampling

To demonstrate the increased value of the combined octree and threshold approach for establishing the 8 × 8 × 8 boxes versus simply downsampling the image, we generated variograms using standard image downsampling as a comparison. For this, we applied a bilinear downsampling to the 512 × 512 × 512 image, creating a 64 × 64 × 64 image–each voxel the equivalent size of an 8 × 8 × 8 octree box. Then we generated variograms utilizing the intensity value of every voxel in the 64 × 64 × 64 image (excluding the background voxels, again defined as those voxels with intensities equal to zero).

### Emphysema index, coefficient of variation, and heterogeneity score

The percentage of lung below a HU threshold value, or emphysema index, was calculated for each rat to compare this conventional measurement of disease severity with the variogram results. Because there is no established HU threshold level in rat models of emphysema, we chose to count the percentage of voxels with HU values below two standard deviations from the control-group mean. This level was determined to be−717 HU. Calculations were made on the same masked images used for octree decomposition.

The CoV was calculated for each rat by fitting a histogram of the masked lung images to a Gaussian curve and taking the ratio of the standard deviation to the mean.

For comparison, we defined a new metric, the heterogeneity score (Δ), as the average difference between the mean variance of the control group and the spatial variance of each individual rat, in the range *d* ≤ *d*
_{
max
}.

### Statistical analysis

Comparisons of variograms of the three control group rats were made using the Kruskal-Wallis (KW) rank sum test with a null hypothesis of α = 0.05. This indicated whether there were significant differences within the group. For pairwise comparisons, a Mann-Whitney-Wilcoxon (MWW) rank sum test was employed, also with a null hypothesis of α = 0.05.

## Results

### CT images

### Octree decomposition tests

#### Threshold

The starting octree decomposition threshold level determined from the standard deviations of the control rats (Section Threshold range) was 60 ± 1 HU; therefore, 60 HU was used as the initial threshold level. Thus, the other threshold levels tested were 20, 40, 80, and 120 HU. Octree decompositions were performed using the different threshold levels on the three control rats, and variograms were compared using a KW rank sum test. The threshold level that showed the least difference among the control rats was 40 HU (p = 0.12), with 60 HU also showing no significant differences (p = 0.09). The 20, 80, and 120 HU thresholds did have significant differences among the controls (p = 0.02, p < 0.0001, and p < 0.0001, respectively). Based on this result, we chose to use the threshold level of 40 HU for all subsequent octree decompositions (unless otherwise noted).

#### Image filtering

Results of applying the different filters showed that the variograms from the unfiltered image and the image with the median filter were indistinct in an MWW test (p = 0.62). However, the relatively noisy unfiltered image resulted in about 20% fewer 8 × 8 × 8 boxes than the filtered image. This was in spite of a higher threshold that was used, 48 HU, based on the unfiltered image’s standard deviation. The radius = 2 and radius = 4 Gaussian filter results also did not differ significantly from the unfiltered image (p = 0.41 and p = 0.48, respectively) or from the median-filtered image (p = 0.71 and p = 0.84, respectively). A KW test showed that the radius = 2 filter resulted in control group variograms that were not distinct (p = 0.16), but the radius = 4 variograms were (p = 0.0001). However, because of the blurring caused by Gaussian filters, the number of 8 × 8 × 8 boxes that resulted from the Gaussian-filtered images was about 30% higher than the median filtered image, presumably because vascular structures and edges were not well preserved. Threshold levels used for the radius = 2 and radius = 4 images were 37 HU and 39 HU, respectively, based on the post-filtering standard deviations.

#### Effect of downsampling

#### Translation and rotation

Results of image translation and rotation were compared for one rat. For the eight-pixel translation, variogram results were exactly identical to that of the original image, as expected. The four-pixel shift was not statistically distinct from the original image as confirmed by a MWW test (p = 0.96). In addition, the result of the rotation showed no significant change from the original (p = 0.56) or from the four-pixel shift (p = 0.59). Thus, we confirm that shifts or rotations to the image (i.e. alternative positions of the lung during imaging) do not result in significant changes to the resulting variograms.

### Octree decomposition

### Variograms

*d*

_{ max }(see Section Variogram analysis). For

*d*≤

*d*

_{ max }a KW test showed that the control rats were statistically similar (p = 0.12), and a MWW test showed that the dose groups were each statistically distinct from the control group (p < 0.0001) and from one another (p < 0.0001).

### Emphysema index, coefficient of variation, and heterogeneity score

For comparison, we calculated the heterogeneity score Δ of each rat’s variance from that of the control group average; see Figure 7C. There is no overlap of Δ between dose groups (i.e. the highest Δ of the control group is lower than the lowest Δ of the full-lung group). This is not observed in Figure 7A or B, suggesting better sensitivity and specificity for Δ than the other metrics.

## Discussion

This work combined two different data analysis tools, octree decomposition and variograms, to study tissue heterogeneity in lung disease. We showed that this merged approach was better able to differentiate rats with mild emphysematous disease from the healthy control group than methods that relied on absolute HU values. The main criterion for octree decomposition was based on the standard deviation of HU values within an octree box. An advantage to this approach is that it avoids thresholding according to HU values, although sophisticated thresholding algorithms may be useful [35–37]; rather, it focuses only on heterogeneity-based signatures that may characterize disease [2]. We propose that a heterogeneity score Δ, the average distance of a rat’s variance from that of the control group average, may be useful to classify disease severity. Furthermore, to visualize the regions of the lung with the greatest heterogeneity, one could determine which boxes had the highest semi-variance within *d*
_{
max
} and map them back to the original image. This would provide 3D information about the spatial distribution of lung tissue heterogeneity and, potentially, disease distribution.

Another approach to a disease metric might be that of fitting the data to an established variogram model, most of which describe an asymptotic rise in variance (i.e. variance becomes independent of distance indicating that spatial relationships become random) [23]. This is seen to some degree in our data (see Figure 6). However, within the range of *d* ≤ *d*
_{
max
}, we found that a power law model generally fit our data better than other models. This model typically describes fractal behavior [38]. Our initial investigations into this model showed potential promise at using fit parameters to distinguish dose groups; however, results lacked statistical significance, and the full-lung dose group tended to not fit this model as well as the other two groups.

There was no single variogram model that satisfactorily fit the all the data over the entire range of distances, because the complex geometries found in the lung result in some problems for variograms. In particular, the direct linear path between two regions of the lung separated by large distances often crosses non-lung tissue, such as the heart, which are essentially treated as holes or voids in the geometry. Furthermore, neighboring lobes generally do not interact physiologically except through the vascular and airway trees, which may only connect regions through many orders of branching. As pointed out by Keil et al. [25], one could go to extraordinary measures to take into account structural distances (the physiological distance at which different regions interact) versus Euclidian (straight line) distances used herein. To limit the problem in this study, we constrained the distance of variogram analysis to approximately half the characteristic diameter of the largest lobe. However, to better understand the inter-and intra-lobe variance relationships, the lung could be segmented into lobes (if the image is of sufficient resolution to discern lobar boundaries), and the decomposition/variogram process repeated on the segmented images. Though, by ignoring inter-lobar variances, this approach would likely not capture information about disease that was confined to a single lobe, particularly if the entire lobe was affected homogeneously.

We employed octree decomposition prior to generating the variograms. In doing this, we assumed that the emphysematous disease is generally slowly varying over space, and that the disease causes changes to homogeneity on the order of or less than the lobar length scale but greater than the octree box length scale [4, 34]. This is consistent with what we previously observed using ^{3}He MRI in the same disease model [26]. Without this assumption, variograms would have to be made directly from the raw images. This is possible but impractical, because the semi-variance computation time (and resulting file size) is proportional to the number of voxel pairs, which rises approximately as *n*
^{2} (see Eq. 2). We measured the variogram computation time for 1078 octree-decomposed boxes from one rat to be 6.6 seconds (on a MacPro model 3.1), and we verified experimentally that the computation time indeed rose in proportion to the square of the number of voxels. Therefore, to create a variogram on the entire masked 3D image, which consisted of 1.48 × 10^{6} voxels (of lung tissue only), it would take us ≈ 1.25 × 10^{7} seconds, or about 5 months–with a resulting file size on the order of 20 GB. Therefore, octree decomposition dramatically reduces the computation time while focusing non-subjectively on regions of the lung that are of greatest interest. We note that the octree decomposition itself was performed in ~2 minutes.

Image noise can confound octree decomposition and affect resulting variograms. Results of image filtering tests indicated that noise reduction using an edge-preserving filter resulted in more 8 × 8 × 8 octree blocks without significantly affecting variogram results. An alternative to octree decomposition is downsampling, which is a quick and straightforward approach to reducing image noise and size. However, we verified that the octree decomposition approach performed much better at separating dose groups than simply downsampling the image and then calculating the semi-variance using every voxel. The octree decomposition assures that only the homogeneous regions of the lungs are singled out for comparison, whereas downsampling the image blurs together proximal voxels, including vasculature, airways, and lung boundaries irrespective of signal intensity or tissue type. Thus, the downsampling approach apparently causes a loss of information. The radius = 4 Gaussian filter had a result similar to downsampling.

One limitation of this pilot study was the small number of animals in each group, which did not allow the statistical evaluation of specificity and sensitivity. Therefore, follow-on work will be required to validate these results and establish specificity and sensitivity [2]. This might be accomplished in conjunction with pulmonary function tests, conventional morphometric measurements [39], and histological techniques particularly sensitive to early emphysematous changes [40]. Additional future work should evaluate the performance of this method on clinical CT images as well as test the effectiveness for distinguishing different diseases and disease models.

## Conclusion

Results of this pre-clinical study of elastase-treated rats suggests that automated octree decomposition and variogram analysis based on image heterogeneity may provide a non-objective and sensitive metric for characterizing emphysematous lung disease, even in early disease stages. The method outperformed conventional approaches that utilize thresholding and absolute HU values. This approach may be applicable to human datasets and other diseases.

## Declarations

### Acknowledgements

The authors would like to thank T. Curry of PNNL for assistance with animal handling and B. Amidan of PNNL for helpful discussions. This project was supported by Award Number R01HL073598 from the National Heart, Lung, and Blood Institute and by PNNL through internal Laboratory Directed Research and Development LDRD DE-AC05-76RL01830. There was no involvement on the part of the sponsors in study design; in the collection, analysis, and interpretation of data; in the writing of the report; or in the decision to submit the paper for publication.

## Authors’ Affiliations

## References

- West JB: Pulmonary Pathophysiology. 1998, Philadelphia: Lippincott Williams & Wilkins, 5Google Scholar
- Mets OM, De Jong PA, Van Ginneken B, Gietema HA, Lammers JW: Quantitative computed tomography in COPD: possibilities and limitations. Lung. 2012, 190 (2): 133-145. 10.1007/s00408-011-9353-9.View ArticlePubMedGoogle Scholar
- Robertson HT, Buxton RB: Imaging for lung physiology: what do we wish we could measure?. J Appl Physiol. 2012, 113 (2): 317-327. 10.1152/japplphysiol.00146.2012.View ArticlePubMedGoogle Scholar
- Litmanovich D, Boiselle PM, Bankier AA: CT of pulmonary emphysema–current status, challenges, and future directions. Eur Radiol. 2009, 19 (3): 537-551. 10.1007/s00330-008-1186-4.View ArticlePubMedGoogle Scholar
- Chong D, Brown MS, Kim HJ, Van Rikxoort EM, Guzman L, McNitt-Gray MF, Khatonabadi M, Galperin-Aizenberg M, Coy H, Yang K, et al: Reproducibility of volume and densitometric measures of emphysema on repeat computed tomography with an interval of 1 week. Eur Radiol. 2012, 22 (2): 287-294. 10.1007/s00330-011-2277-1.View ArticlePubMedGoogle Scholar
- Uppaluri R, Mitsa T, Sonka M, Hoffman EA, McLennan G: Quantification of pulmonary emphysema from lung computed tomography images. Am J Respir Crit Care Med. 1997, 156 (1): 248-254. 10.1164/ajrccm.156.1.9606093.View ArticlePubMedGoogle Scholar
- Besir FH, Mahmutyazicioglu K, Aydin L, Altin R, Asil K, Gundogdu S: The benefit of expiratory-phase quantitative CT densitometry in the early diagnosis of chronic obstructive pulmonary disease. Diagn Interv Radiol. 2012, 18 (3): 248-254.PubMedGoogle Scholar
- Irion KL, Marchiori E, Hochhegger B, Porto Nda S, Moreira Jda S, Anselmi CE, Holemans JA, Irion PO: CT quantification of emphysema in young subjects with no recognizable chest disease. AJR Am J Roentgenol. 2009, 192 (3): W90-W96.View ArticlePubMedGoogle Scholar
- Reske AW, Busse H, Amato MB, Jaekel M, Kahn T, Schwarzkopf P, Schreiter D, Gottschaldt U, Seiwerts M: Image reconstruction affects computer tomographic assessment of lung hyperinflation. Intensive Care Med. 2008, 34 (11): 2044-2053. 10.1007/s00134-008-1175-8.View ArticlePubMedGoogle Scholar
- Yuan R, Mayo JR, Hogg JC, Pare PD, McWilliams AM, Lam S, Coxson HO: The effects of radiation dose and CT manufacturer on measurements of lung densitometry. Chest. 2007, 132 (2): 617-623. 10.1378/chest.06-2325.View ArticlePubMedGoogle Scholar
- Chen CC, Daponte JS, Fox MD: Fractal feature analysis and classification in medical imaging. IEEE Trans Med Imaging. 1989, 8 (2): 133-142. 10.1109/42.24861.View ArticlePubMedGoogle Scholar
- Glenny RW, Robertson HT: Fractal properties of pulmonary blood flow: characterization of spatial heterogeneity. J Appl Physiol. 1990, 69 (2): 532-545.PubMedGoogle Scholar
- Copley SJ, Giannarou S, Schmid VJ, Hansell DM, Wells AU, Yang GZ: Effect of aging on lung structure in vivo: assessment with densitometric and fractal analysis of high-resolution computed tomography data. J Thorac Imaging. 2012, 27 (6): 366-371. 10.1097/RTI.0b013e31825148c9.View ArticlePubMedGoogle Scholar
- Kido S, Ikezoe J, Naito H, Tamura S, Machi S: Fractal analysis of interstitial lung abnormalities in chest radiography. Radiographics. 1995, 15 (6): 1457-1464. 10.1148/radiographics.15.6.8577968.View ArticlePubMedGoogle Scholar
- Uppaluri R, Mitsa T, Galvin JR: Fractal analysis of high-resolution CT images as a tool for quantification of lung disease. Medical Imaging 1995: Physiology and Function from Multidimensional Images. Edited by: Hoffman EA. 1995, Bellingham, WA: SPIE, 133-142. 2433View ArticleGoogle Scholar
- Kido S, Sasaki S: Fractal analysis for quantitative evaluation of diffuse lung abnormalities on chest radiographs: use of sub-ROIs. J Thorac Imaging. 2003, 18 (4): 237-241. 10.1097/00005382-200310000-00005.View ArticlePubMedGoogle Scholar
- Subramaniam K, Hoffman EA, Tawhai MH: Quantifying tissue heterogeneity using quadtree decomposition. Conf Proc IEEE Eng Med Biol Soc. 2012, 2012: 4079-4082.PubMedGoogle Scholar
- Shephard MS, Georges MK: Automatic 3-dimensional mesh generation by the finite octree technique. Int J Numer Meth Eng. 1991, 32 (4): 709-749. 10.1002/nme.1620320406.View ArticleGoogle Scholar
- Zhang JH, Owen CB: Octree-based animated geometry compression. Comput Graph-Uk. 2007, 31 (3): 463-479. 10.1016/j.cag.2006.12.002.View ArticleGoogle Scholar
- Szeliski R, Lavallee S: Matching 3-D anatomical surfaces with non-rigid deformations using octree-splines. Int J Comput Vision. 1996, 18 (2): 171-186. 10.1007/BF00055001.View ArticleGoogle Scholar
- Zhang YJ, Bajaj C, Sohn BS: 3D finite element meshing from imaging data. Comput Method Appl M. 2005, 194 (48–49): 5083-5106.View ArticleGoogle Scholar
- Dua S, Kandiraju N, Chowriappa P: Region quad-tree decomposition based edge detection for medical images. Open Med Inform J. 2010, 4: 50-57.View ArticlePubMedPubMed CentralGoogle Scholar
- Clark I, Harper WV: Practical Geostatistics 2000. 2000, Columbus, Ohio: Ecosse North America, LLCGoogle Scholar
- Gringarten E, Deutsch CV: Variogram interpretation and modeling. Math Geol. 2001, 33 (4): 507-534. 10.1023/A:1011093014141.View ArticleGoogle Scholar
- Keil F, Oros-Peusquens AM, Shah NJ: Investigation of the spatial correlation in human white matter and the influence of age using 3-dimensional variography applied to MP-RAGE data. Neuroimage. 2012, 63 (3): 1374-1383. 10.1016/j.neuroimage.2012.07.034.View ArticlePubMedGoogle Scholar
- Jacob RE, Minard KR, Laicher G, Timchalk C: 3D 3He diffusion MRI as a local in vivo morphometric tool to evaluate emphysematous rat lungs. J Appl Physiol. 2008, 105 (4): 1291-1300. 10.1152/japplphysiol.90375.2008.View ArticlePubMedPubMed CentralGoogle Scholar
- Jacob RE, Lamm WJ: Stable small animal ventilation for dynamic lung imaging to support computational fluid dynamics models. PLoS One. 2011, 6 (11): e27577-10.1371/journal.pone.0027577.View ArticlePubMedPubMed CentralGoogle Scholar
- Corley RA, Kabilan S, Kuprat AP, Carson JP, Minard KR, Jacob RE, Timchalk C, Glenny R, Pipavath S, Cox T, et al: Comparative computational modeling of airflows and vapor dosimetry in the respiratory tracts of rat, monkey, and human. Toxicol Sci. 2012, 128 (2): 500-516. 10.1093/toxsci/kfs168.View ArticlePubMedPubMed CentralGoogle Scholar
- Minard KR, Kuprat AP, Kabilan S, Jacob RE, Einstein DR, Carson JP, Corley RA: Phase-contrast MRI and CFD modeling of apparent (3)He gas flow in rat pulmonary airways. J Magn Reson. 2012, 221: 129-138.View ArticlePubMedPubMed CentralGoogle Scholar
- ImageJ. http://imagej.nih.gov/ij/,
- ImageJ Plugins: 3D Toolkit. http://ij-plugins.sourceforge.net/plugins/3d-toolkit/index.html,
- Jackins CL, Tanimoto SL: Oct-trees and their use in representing three dimensional objects. Computer Graphics and Image Processing. 1980, 14 (3): 249-270. 10.1016/0146-664X(80)90055-6.View ArticleGoogle Scholar
- Schulz H, Mühle H: Respiration. The Laboratory Rat. Edited by: Krinke GJ. 2000, San Diego: Academic Press, 323-336.View ArticleGoogle Scholar
- Spencer H: Pathology of the Lung. Volume Volume 1. 1985, New York: Pergamon Press, 557-594. 4Google Scholar
- El-Baz A, Gimel’farb G, Falk R, Holland T, Shaffer T: A new stochastic framework for accurate lung segmentation. Medical Image Computing and Computer-Assisted Intervention–MICCAI. Edited by: Metaxas D, Axel L, Szekely G. 2008, New York, NY: Springer, 322-330.Google Scholar
- Abdollahi B, Soliman A, Civelek AC, Li XF, Gimel’farb G, El-Baz A: A novel 3D joint MGRF framework for precise lung segmentation. Third International Workshop MLMI: 2012. 2012, Nice, France: Springer, 86-93.Google Scholar
- El-Baz A, Beache GM, Gimel’farb G, Suzuki K, Okada K, Elnakib A, Soliman A, Abdollahi B: Computer-aided diagnosis systems for lung cancer: challenges and methodologies. Int J Biomed Imaging. 2013, 2013: 942353-PubMedPubMed CentralGoogle Scholar
- Bohling GC: Introduction to Geostatistics. Kansas Geological Survey Open File Report no. 2007-26. 2007, 50-Google Scholar
- Hsia CC, Hyde DM, Ochs M, Weibel ER: An official research policy statement of the American Thoracic Society/European Respiratory Society: standards for quantitative assessment of lung structure. Am J Respir Crit Care Med. 2010, 181 (4): 394-418. 10.1164/rccm.200809-1522ST.View ArticlePubMedGoogle Scholar
- Jacob RE, Carson JP, Gideon KM, Amidan BG, Smith CL, Lee KM: Comparison of two quantitative methods of discerning airspace enlargement in smoke-exposed mice. PLoS One. 2009, 4 (8): e6670-10.1371/journal.pone.0006670.View ArticlePubMedPubMed CentralGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2342/14/1/prepub

### Pre-publication history

## Copyright

This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.