|
|
||||||||
Departments of 1 Physiology and 2 Anesthesiology, Medical College of Wisconsin, Milwaukee 53226; 3 Department of Biomedical Engineering, Marquette University, Milwaukee 53201-1881; and 4 Research Service, Zablocki Veterans Affairs Medical Center, Milwaukee, Wisconsin 53295
| |
ABSTRACT |
|---|
|
|
|---|
The objective of this study was to develop an
X-ray computed tomographic method for measuring pulmonary arterial
dimensions and locations within the intact rat lung. Lungs were removed
from rats and their pulmonary arterial trees were filled with
perfluorooctyl bromide to enhance X-ray absorbance. The lungs were
rotated within the cone of the X-ray beam projected from a microfocal
X-ray source onto an image intensifier, and 360 images were obtained at
1° increments. The three-dimensional image volumes were reconstructed with isotropic resolution using a cone beam reconstruction algorithm. The vessel diameters were obtained by fitting a functional form to the
image of the vessel circular cross section. The functional form was
chosen to take into account the point spread function of the image
acquisition and reconstruction system. The diameter measurements
obtained over a range of vascular pressures were used to characterize
the distensibility of the rat pulmonary arteries. The distensibility
coefficient
[defined by D(P) = D(0)(1 +
P), where D(P) is
the diameter at intravascular pressure (P)] was ~2.8% mmHg and
independent of vessel diameter in the diameter range (about 100 to
2,000 mm) studied.
cone beam reconstruction; pulmonary arterial diameter; pulmonary blood flow distribution
| |
INTRODUCTION |
|---|
|
|
|---|
X-RAY MICROCOMPUTED TOMOGRAPHY (CT) promises to be a valuable source of vascular structure-function information and vascular phenotypes of small laboratory animals (7, 13, 25, 36). The rat has been widely used to study pulmonary vascular remodeling (12, 17, 29, 32, 38) and is being used increasingly for physiological genomics in general (14, 15, 23, 41, 46) and for physiological genomics of the pulmonary vasculature in particular (21, 26, 29, 42, 47). Thus the objective of this study was to develop an approach for X-ray micro-CT measurement of pulmonary arterial diameters and their three-dimensional locations within the intact rat lung. The method was applied to the measurement of rat pulmonary arterial distensibility, which is a potentially useful pulmonary vascular phenotype for identifying quantitative trait loci (15, 41) as well as an input to rat lung physiome models (4, 16, 18, 40) used to understand pulmonary vascular function.
| |
METHODS |
|---|
|
|
|---|
Animal Preparation
Lungs from four Sprague-Dawley rats (267 ± 57 g) were prepared for imaging as follows. Each rat was anesthetized with pentobarbital sodium (40 µg/g body wt ip), the trachea was clamped, and the chest was opened. Heparin (200 international units in 0.2 ml) was administered via injection into the right ventricle. The pulmonary artery was cannulated with a saline-filled catheter [polyethylene tubing 1.67 mm inner diameter (ID); 2.42 mm outer diameter] via the conus arteriosis and the heart was dissected away. The lungs were removed from the chest, and suspended from the trachea and pulmonary arterial cannula. The lungs were ventilated with a gas mixture containing 15% O2-6% CO2 in nitrogen at 40 breaths/min with end-inspiratory and end-expiratory tracheal pressures of 8 and 3 mmHg, respectively. This served to eliminate any atelectasis occurring during the dissection. The pulmonary artery cannula was connected to a perfusion system primed with a physiological salt solution containing 5% bovine serum albumin (3), and the lungs were perfused for about 5 min at a flow rate ranging from 5 to 40 ml/min to remove residual blood from the lung vessels. The maximum pulmonary artery pressure at 40 ml/min was ~10 mmHg, which was maintained for ~20 s of the 5 min. The perfusate exited via the severed pulmonary vein. Once cleared of blood, the lung, still suspended from the cannulas, was placed in a 41 mm ID plastic cylinder, or the right lung and cardiac lobar arteries and bronchi were ligated and the respective lobes removed so that the remaining left lung could be placed in a 24 mm ID cylinder. The cylinder axis was located at the center of a horizontal turntable so that the lungs could be rotated 360° around a vertical axis between the X-ray source and detector with no other significant X-ray absorbing objects passing through the beam. The airway pressure was set at 6 mmHg, and the salt solution in the reservoir that was connected to the arterial catheter was replaced by perfluorooctyl bromide (PFOB), which was allowed to fill the arterial tree at a pressure of ~20 mmHg. The PFOB provided high X-ray contrast for the vessel lumen, and the surface tension at the PFOB-aqueous interface prevented its entry into the capillary bed. Thus only the arterial vessels were filled. Then the arterial pressure was set at 30 mmHg relative to the horizontal plane through the center of the X-ray image of the lungs. The lungs were rotated continuously at one degree per second. The image acquisition sequence was 5 to 10 frames at 30 frames per second beginning at each 1° increment to acquire 360 X-ray image sets in ~6 min. The same procedure was repeated with the arterial pressure set at 21 mmHg, 12 mmHg, and 5.4 mmHg and, again, with the pressure returned to 30 mmHg. The actual intravascular pressure within each vessel, relative to atmospheric pressure at the level of the vessel, was obtained from the vessel's vertical distance from the pressure reference level at the central horizontal plane of the image and the PFOB density (1.94 g/ml).Imaging
The X-ray system included a Fein-Focus FXE-100.50 X-ray tube with 3-µm focal spot, a North American Imaging AI-5830-HP image intensifier set at either the 17.8 or 23 cm aperture, and a SMD 1M15 charge-coupled device (CCD) camera (Silicon Mountain Design, Colorado Springs, CO). The cylinder containing the rat lung was placed in the scanner so that its central axis was from 13 to 28 cm from the source. The source to image intensifier distance ranged from 50 to 60 cm, such that the geometric magnification was greater than ×4 and the half cone beam angle was <11.2°. Figure 1 is a schematic of the micro-CT imaging system demonstrating the magnification obtained by separating the object from the detector. The 5 or 10 consecutive frames comprising each image set were averaged to produce the stored image for each 1° of rotation. The same procedure was carried out to image various phantoms used to characterize and evaluate aspects of the imaging acquisition and analysis. After each rat lung or phantom imaging session, two additional images were obtained. One was a flood-field image with the lung removed from the beam. The flood-field image was used to correct for spatial variations in the X-ray beam and/or image-intensifier gain. The other image was of a phantom consisting of a uniform grid of 1-mm diameter stainless steel spheres (BBs) spaced at 1.5-cm intervals and embedded in a Plexiglas disk. The BB phantom was attached to the image-intensifier input surface and the acquired image was used to correct for spatial distortion (warping) due to the beam geometry and image intensifier.
|
Image Preprocessing
Each 8-bit planar image consisted of a 512 × 512 array of pixels ranging from 0 to 255 (minimum to maximum) X-ray intensity scale. Although the CCD camera had the capability of collecting 1,024 × 1,024 pixel images, we operated in the 512 × 512 mode for this study, because, with our present facilities, the larger image volume required impractical reconstruction time. The 512 × 512 reconstruction time on a Pentium-based 550 MHz personal computer was ~10 h. Before reconstructing the image volume from the 360 planar images, preprocessing of the individual images was performed in the following steps: 1) two-dimensional polynomial dewarping to correct the image-intensifier spatial distortion, 2) locating the axis of rotation and cropping the projection images to center on that axis, 3) flood-field division to correct for nonuniform illumination intensity, and 4) normalization of the intensity between projections to correct for any temporal drift.Dewarping.
The spatial distortion correction algorithm consisted of two steps.
First, the uncorrected individual bb center of mass coordinates (x',y') were determined from the BB phantom
image. These coordinates for each set of four neighboring BBs were
represented by a set of eight equations (Eq. 1 and
2 for each BB, where i is bb number, 1, 2, 3 or
4). These equations were solved simultaneously for the coefficients
a-h which map x',y' via
bilinear interpolation onto the known x,y
coordinates of the grid vertices.
|
(1) |
|
(2) |
|
|
Central axis location. The rotation axis was determined as the midpoint between the left- and right-most excursions of a high-contrast feature in a two-dimensional radon transform or sinogram using the same row of pixels extracted from each of the 360 projection images. Once this location was established, each planar image was cropped so that the axis of rotation occupied the central column of pixels. This resulted in a final preprocessed image size of 497 × 497 pixels.
Spatial normalization. The flood-field image acquired with the same X-ray voltage, current, and geometry as the lung projection images provided the unattenuated incident illumination intensity for each pixel. Thirty frames were averaged to suppress temporal noise (primarily quantum mottle). The spatially corrected, cropped projections were then divided by the flood-field image to provide pixel values proportional to X-ray absorption. The uniform background on the BB phantom image in Fig. 3d illustrates the result of the latter step.
Temporal normalization. Finally, three 20 × 20 pixel regions of interest (ROI) unobscured by the object in any projection were located within the volume. The mean intensity for the three ROI divided by the respective mean flood-field intensity was used as a scaling factor for each image. This scaling factor corrected for any temporal drift in the overall illumination intensity during planar image acquisition (24). This scaling step turned out to be unnecessary because there was no systematic temporal drift in the mean ROI intensities between projection images and the coefficient of variation (CV) in the mean ROI intensity for the 360 projection images comprising an image volume was only about 0.1% for these data sets.
Image Reconstruction
After preprocessing of the projection images, cone-beam reconstruction was performed on the projection data to yield an isotropic reconstruction matrix of 497 × 497 × 497 voxels representing volumes of ~2.5 × 2.5 × 2.5 cm3 to 4.0 × 4.0 × 4.0 cm3 (about 50-80 µm/side). The Feldkamp (22, 24) cone-beam algorithm was utilized to weight the projection data to account for the variation in the angle of the incident X-ray beam at each pixel. The reconstruction was obtained by first convolving the weighted projection data, perpendicular to the rotation axis, with the Shepp-Logan (43) filter kernel. The filtered data were then back-projected from every angle. The major advantage of the cone beam reconstruction was that it allowed for reconstruction of the entire volume from a single 360° scan. This minimized the scan time, which was important for maintaining tissue mechanical properties during the image acquisition period.The postreconstruction voxel gray scale values (I) were
equalized so that the range from air to PFOB was the same for all reconstructed volumes. First, a region of the volume known to contain
only air was selected and its I values were averaged.
i represents the average I value for air in a
particular image volume. The maximum I in the volume,
i, was then determined. Each voxel's intensity was then
adjusted according to Eq. 3
|
(3) |
To test the linearity of the resulting GSN with respect to object density following this equalization procedure, nine serial dilutions of iodinated contrast medium (0 to 0.282 g/ml of iodine in the form of diatrizoate meglumine) were loaded into lengths of 1.67-mm ID polyethylene tubing and imaged within a single volume. After reconstruction and equalization, the GSN value over the center of each cross section was determined. The coefficient of determination between GSN and iodine concentration was 0.996, indicating that the GSN was linearly proportional to X-ray absorbance.
The left side of Fig. 4 is a
three-dimensional rendering of the reconstructed image volume to
provide a sense of the data set obtained as a result of these
procedures. It was obtained after thresholding the CT volume to allow
appreciation of the gross structure. In what follows, the full,
unthresholded data set, exemplified by the two transaxial slices
depicted in Fig. 4, was used to make quantitative measurements.
|
| |
MEASUREMENTS |
|---|
|
|
|---|
Mapping
The first step required for determining vessel dimensions was to locate the vessel segment (the portion of a vessel between consecutive bifurcations) in three-dimensional space. This served two purposes. One was to obtain an orthogonal slice through the segment for subsequent diameter measurement. The other was to identify the vessel segment so that the same vessel segment could be found in a different reconstructed volume. The vessel segment was located by observing a vessel cross section while progressing through the sequence of transaxial slices, as depicted in Fig. 5. Figure 5, left, is a three-dimensional rendering of a portion of a rat pulmonary artery. Figure 5, right, is the corresponding set of transaxial slices demonstrating identification of two bifurcations along the artery. The single ellipsoid becomes bilobular and then splits into two. The location where there is visible separation between the two new ellipsoids represents the crotch of the bifurcation. A series of steps was performed to identify the vessel segment between subsequent bifurcations within the image volume. First, the central voxel coordinates (x1, y1, z1) of each of the new ellipsoids in the slice at which the single ellipsoid completely separates into two were recorded. Second, taking each new ellipsoid individually and proceeding transaxially through the image volume, the coordinates of the central voxel of the ellipsoid were recorded in the slice preceding the slice where lobulation reoccurred (x2, y2, z2). The central axis of a vessel segment was approximated as the line segment connecting (x1, y1, z1) and (x2, y2, z2). The mapping procedure and coordinate locations are illustrated in Fig. 6. A particular vessel segment could be identified by counting the number of upstream bifurcations along the pathway leading to that segment.
|
|
Locating Orthogonal Slices
The coordinates in three-dimensional space that defined the vessel central axis also provided an estimate of the vessel segment midpoint. Two vectors that passed through this midpoint, orthogonal to the central axis and to one another, described an orthogonal plane with the vessel segment midpoint as the origin. Cubic interpolation was performed to obtain the GSN values of the orthogonal slice. This procedure is also illustrated in Fig. 6. Figure 7 is an actual vessel segment cross section that is elliptical in the transaxial slice but circular in the orthogonal slice. Although elliptical orthogonal cross sections may be observed in some vascular beds, normal pulmonary arteries are nearly circular (49), as confirmed by observation of the many orthogonal slices in the present study.
|
Measuring Diameters
Approaches for estimating vessel cross-sectional dimensions from CT data sets include the brightness/area product (5, 31), the full-width half-maximum (FWHM) (9, 39, 48) and others (37). Each method is based on assumptions regarding 1) the form of the blurring function resulting during image acquisition and processing [the point spread function (PSF)], 2) the size of the vessel relative to the PSF, 3) the thickness of the CT slice, and 4) the obliqueness of the vessel segment axis relative to the orientation of the CT slice. To take full advantage of the resolution obtainable with the high contrast images collected in this study, concepts underlying these methods were extended to derive a general diameter measuring technique applicable for the entire range of vessel diameters [both larger and smaller than the FWHM of the PSF (FWPSF)] as follows.In a theoretical imaging system whose PSF is a delta function and
the CT slice has no thickness, the intensity of the line of pixels (a
line scan) across the minimum diameter of the elliptical cross section
of an obliquely cut vessel would appear as a rectangular function. In
reality, blurring is introduced in the image acquisition and
reconstruction and by finite pixel size and slice thickness. The latter
is eliminated by the use of only orthogonal slices, as indicated above.
The actual blurred image of the vessel is then essentially the PSF, or
blurring function, convolved with the theoretical vessel cross section.
The FWHM approach applied to the line scans recovers the vessel
diameter as long as the diameter is large enough compared with the
FWPSF, but it over estimates the diameter when the diameter approaches
the FWPSF. This can be appreciated by assuming a Gaussian distribution
for the PSF as in Eq. 4
|
(4) |
is the standard deviation. Half the maximum height
of the PSF occurs at x = ±1.175
. Thus the FWPSF is
~2.35
. As the actual vessel diameter decreases below FWPSF, the
FWHM of the blurred line scan approaches FWPSF. The smallest imaginable
diameter may be thought of as almost a delta function that when
convolved with the PSF would return the PSF, in which case the FWHM
method would return the FWPSF rather than the diameter.
The behavior of a line scan across a vessel image is further
demonstrated in Fig. 8, top,
illustrating three sets of rectangular functions representing line
scans of slices through the ideal (unblurred) image volumes containing
cylinders of one X-ray-absorbing medium embedded in three different
absorbing media (simulating different tissue densities). These
functions were convolved with the PSF of Eq. 4 to simulate
line scans of the reconstructed images (Fig. 8, bottom). The
graph in Fig. 9A is the FWHM
for the simulated line scans in Fig. 8, bottom. The FWHM
versus the diameter relationship is independent of the GSN of the
imbedding medium, and linearly proportional to diameter when the
diameter is larger than FWPSF, but it approaches the constant FWPSF
when the diameter is small. The GSN height (Gmax) of the
simulated line scans above that of the imbedding medium
(GSNm) is constant for large cylinders and then
becomes proportional to the diameter as the diameter decreases (Fig.
9B). Taken together, the quantity FWHM × (Gmax
GSNm), which approximates the area
between the line scan and GSNm is linearly proportional to
the diameter over the entire diameter range (Fig. 9C) as is
the actual area under the line scan (Fig. 9D). However, the
slope of the area versus diameter line depends on GSNm as
summarized by Eq. 5
|
(5) |
|
|
Assuming a Gaussian PSF, the functional form of a line scan through the
axis of a vessel can be approximated by Eq. 6, referred to
subsequently as the modified Gaussian function (MGF) and illustrated in
Fig. 10.
|
When the MGF is fit to the line scan, the integral of MGF
GSNm produces an area proportional to diameter. Thus a
method for measuring vessel diameters is achieved. However
|
|
(6) |
| |
RESULTS |
|---|
|
|
|---|
The MGS fit method was tested on a phantom consisting of
calibrated titanium alloy (90% titanium, 6% aluminum, and 4%
vanadium) wire segments of 50, 100, 200, 460, 500, and 1,000 µm
diameters. Two different surrounding absorbing media were obtained by
immersing a portion of the wire lengths in a silicone gel medium,
whereas the rest of the wire lengths were surrounded by only air.
Figure 11 shows the image data in the
same format as the simulations in Fig. 9.
|
To evaluate the effects of both the orientation and position of objects
within the cone beam field of view, another phantom was constructed
with 50 µm and 500 µm wires scattered in a medium of silicone gel
so that they were oriented from 0 to 45° relative to the axis of
rotation and 0 to 90% of the cone beam angle from the center to the
edge of the field of view. The CV in the estimated diameter was 14.5%
for the 50-µm wire and 9.4% for the 500-µm wire. Less than 4.1%
of the variation in estimated diameter could be attributed to the angle
relative to the axis of rotation and <7.4% to vertical distance above
or below the central plane perpendicular to the axis of rotation.
Although the wire phantom may not exactly reproduce the conditions of
an actual vessel lumen imbedded in its unique tissue medium, the
similarities of the average and CV of the standard deviations (
) of
the PSF obtained with wire phantoms [2.42 ± 0.41 pixels
(CV = 17%)], and the 783 vessel diameters measured for the
distensibility study indicated below [2.03 ± 0.36 pixels
(CV = 18%)] is one measure of consistency between the two.
Fig. 12 shows the results of the method
applied to a large (1,150 µm) and a small (43 µm) artery in an
image volume of a rat lung.
|
To investigate the distensibility of rat pulmonary arteries, 162 arterial vessel segments ranging from ~106 to 1,681 µm at the
lowest pressure and from about 151 to 2,743 µm at the highest pressure were measured at the five intravascular pressure settings. Over this range of pressures and vessel sizes, the diameter versus pressure relationship was nearly linear, as exemplified in Fig. 13 by a representative group of vessel
segments covering the range of diameters studied. Therefore, the
diameter versus pressure data were parameterized by a slope (
) and
intercept [D(0)] using Eq. 7 after
Yen et al. (49) and others (2)
|
(7) |
|
The values of
for all the measured vessel segments are plotted
versus D(0) in Fig.
14. Because the fractional differences in diameters over the pressure range studied are small, and because
and D(0) are highly correlated in the
regression analysis for Eq. 7, errors in the diameter
measurements can have a proportionately larger impact on the estimates
of
and D(0) than on the estimates of the
individual vessel diameters. This is reflected in the fact that the CV
in the paired measurements obtained at the high pressure was only
4.1%, whereas the average of the coefficients of variations in
and
D(0) (standard error of the estimate divided by
the estimated value) calculated for the individual vessels were 29 and
28%, respectively. Also the distributions of the coefficients of
variation were highly skewed as reflected by the fact that respective
medians were only 20 and 12%. To graphically represent the
distribution of the standard errors of the estimates in Fig. 14, the
gray level of the symbol representing
[D(0),
] for each vessel is inversely proportional to the standard error. The key observation from this weighting of the data symbols is that the outlying points tend to have
large standard errors. Thus the variability in the relationship between
and D(0) includes a component that is due
to the amplification of errors in the diameter measurements resulting
from the transformation resulting from fitting Eq. 7 to the
data. To evaluate the relationship between
and
D(0) for the data in Fig. 14, linear regression
was carried out assuming proportionate errors in both y and
x values, i.e., by minimizing the sum of the squared
distances orthogonal to the regression line (10). The
values were either weighted equally or they were weighted according to
the standard error of the estimate as indicated by the symbol gray
levels in Fig. 14. Also, the intercept was either a free parameter or
fixed at zero. The F-test was used to determine whether the
weighting significantly altered the fit and whether a nonzero intercept
was supported by the data. The results were essentially the same
whether the data were weighted or not, and there was no significant
improvement in the fits when the intercept was free. The latter implies
that the average vessel distensibility,
, defined as
/D(0) for an individual vessel, was
independent of vessel size with the average value approximated by the
slope of the Fig. 14 graph. The slope was 2.85 ± 0.19 for the
weighted and 2.82 ± 0.19%/mmHg (SE) for the unweighted fit. The
average of the values of
calculated for each individual vessel was
3.23 ± 3.56%/mmHg (SD) (CV = 110%), but the distribution
is quite skewed as expected for this transformation of the diameter
data (34). This is reflected by the fact that the median
of the individual values was 2.73%/mmHg, and thus closer to the
regressed values from Fig. 14.
|
Discussion
The key aspects of the methodology were generally discussed above as motivation for the various steps in the process. In summary, the approach taken was to use the geometric magnification obtainable with a small focal spot to obtain high-resolution CT images of the rat lung with the pulmonary arterial tree contrast enhanced by the brominated perfluorocarbon. The three-dimensional image volume was reconstructed with isotropic resolution using a cone beam reconstruction algorithm. The vessel diameters were then obtained by fitting a functional form to the image of the circular vessel cross section. The functional form was chosen to take into account the PSF of the image acquisition and reconstruction system. The diameter measurements obtained over a range of vascular pressures were used to characterize the distensibility of the rat pulmonary arteries.The vessel distensibility results indicate that over the vessel diameter range studied the distensibility of the rat pulmonary arteries (2.8%/mmHg) is within the range of values for pulmonary arteries of larger species studied previously. In summarizing the data available at the time, Al-Tinawi et al. (2) found that an average value of about 2%/mmHg was consistent with measurements on dog, cat, and human pulmonary arteries over the full range of arterial diameters. However, the variability from study to study tends to confound the issue of possible species differences. The distensibility of rat pulmonary arteries has been studied previously either by mounting isolated vessel segments on wires so that the force and distance between the wires could be measured (12, 28) or by cannulating isolated vessel segments such that the transmural pressure could be manipulated while the diameters were monitored via video microscopy (32). The distensibility calculated for the cannulated arteries in the 500- to 700-µm diameter range was also about 2.8%/mmHg (32). For the wire-mounted vessels segments in the 120-µm to 1-mm range studied, the average was about 3.7%/mmHg (12, 28). The somewhat larger latter values might have something to do with the different vessel preparation. The distensibility in the present study is defined by the intravascular pressure rather than transmural pressure. Transmural pressure is also affected by transpulmonary pressure and lung volume (1, 27, 45). Thus the distensibility values obtained in this study have to be considered to be specific for the transpulmonary pressure used. It is also possible that transmural pressure is a variable contributing to the range of previously reported values in other species. Such questions can be addressed in future studies using the methods described herein.
Although it has been suggested by some studies that pulmonary arterial distensibility may be vessel size dependent (33), the summary provided by Al-Tinawi et al. (2) suggested that the distensibility of the pulmonary arteries of the various species studied is nearly the same for vessels of any size, or that any systematic dependence on vessel size is small compared with the variability among vessels of a given diameter. In this regard, the present rat lung data are consistent with the results from the larger species. This diameter independence may be somewhat surprising given how vessel wall composition varies over this diameter range (17). However, there may be an adaptive advantage to the constancy of this mechanical property.
Figures 13 and 14 reflect two key relationships, namely, the nearly
linear diameter versus pressure relationship for the individual vessels
over the pressure range studied and the diameter independence of the
vessel distensibility over the diameter range studied. One functional
consequence of such relationships is that they tend to minimize the
sensitivity of the distribution of flow among the branches of a
heterogeneous, asymmetric, vascular tree to changes in total flow rate.
If the distensibility were diameter dependent, for example decreasing
or increasing in proportion to vessel diameter, increasing cardiac
output would tend toward a redistribution of flow to, or away from,
respectively, the smaller branch at an asymmetric bifurcation.
Likewise, because of the arterial-to-venous pressure drop through the
vascular network, a nonlinear diameter versus pressure relationship
might have a similar effect when the pressure drop is affected by
a change in total pulmonary flow. This can be appreciated by noting
that the inflow pressure (Pin) versus flow (
)
relationship for a distensible vessel segment having a cylindrical
configuration at D(0), length L,
outflow pressure Pout, viscosity µ, distensibility
,
defined by D(P) = D(0)(1 +
P), and Poiseuille flow (19, 30, 50) is
|
(8) |
is constant for all vessel
segments, regardless of D(0) and with
as a
function of D(0). With
constant, the relative distribution of flow among the branches and out the outlets is
independent of the total flow entering the tree. However, with
either directly or inversely proportional to
D(0) the flow distribution is dependent on the
total flow, as exemplified by the outlet flow distributions in Fig. 15.
|
|
Changes in total pulmonary blood flow have been accompanied by changes in the spatial distribution of pulmonary blood flow to lesser and greater extents in different experimental settings (6, 8, 11, 35, 44), and various factors may be involved. However, it might be assumed that the observed distensibility pattern helps to maintain the robustness of the flow distributing network in the face of changing flow, thereby helping to keep the pulmonary flow distribution within bounds appropriate for efficient gas exchange, regardless of the cardiac output.
One of the anticipated uses of this pulmonary vascular imaging approach is for the collection of more extensive pulmonary vascular morphometric data. In the past, pulmonary vascular morphometry has been carried out mainly on a small number of plastic casts of ostensibly normal lungs (20). Imaging methods such as these should be useful for expanding the morphometric database available for evaluating structure/function relationships in abnormal and remodeled pulmonary arteries as well. The concept behind microfocal X-ray imaging is that the small focal spot minimizes the penumbral contribution to the PSF, allowing for high resolution even when the object is close to the source. As magnification is inversely proportional to the distance between the X-ray source and object, a smaller object may be placed closer to the source and still be completely contained within the cone of the beam delimited by the source-to-detector distance and detector diameter. For example, the maximum resolution is greater for a mouse lung, for which the described procedures are also applicable, than for a rat lung. The rat was the subject of this study because of its history of use in studies of pulmonary vascular remodeling and because of its increasing importance in physiological genomics.
| |
ACKNOWLEDGEMENTS |
|---|
This study was supported by National Heart, Lung, and Blood Institute Grant HL-19298, The Whitaker Foundation, The W. M. Keck Foundation, The Falk Medical Trust, and the Department of Veterans Affairs.
| |
FOOTNOTES |
|---|
Address for reprint requests and other correspondence: C. A. Dawson, Research Service 151, Zablocki VAMC, 5000 West National Ave., Milwaukee, WI 53295 (E-mail: cdawson{at}mcw.edu).
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Received 23 March 2001; accepted in final form 1 May 2001.
| |
REFERENCES |
|---|
|
|
|---|
1.
Albert, RK,
Lamm WJE,
Rickaby DA,
Al-Tinwai A,
and
Dawson CA.
Lung inflation distends small arteries (<1 mm) in excised dog lungs.
J Appl Physiol
75:
2595-2601,
1993
2.
Al-Tinawi, A,
Madden JA,
Dawson CA,
Linehan JH,
Harder DR,
and
Rickaby DA.
Distensibility of small arteries of the dog lung.
J Appl Physiol
71:
1714-1722,
1991
3.
Audi, SH,
Bongard RD,
Okamoto Y,
Merker MP,
Roerig DL,
and
Dawson CA.
Pulmonary reduction of an intravascular redox polymer.
Am J Physiol Lung Cell Mol Physiol
280:
L1290-L1299,
2001
4.
Bassingthwaighte, JB.
Strategies for the physiome project.
Ann Biomed Eng
28:
1043-1058,
2000[Web of Science][Medline].
5.
Baxter, BS,
and
Sorenson JA.
Factors affecting the measurement of size and CT number in computed tomography.
Invest Radiol
16:
337-341,
1981[Web of Science][Medline].
6.
Beck, KC,
and
Rehder K.
Differences in regional vascular conductances in isolated dog lungs.
J Appl Physiol
61:
530-538,
1986
7.
Beighley, PE,
Thomas PJ,
and
Jorgensen ELR
3D architecture of myocardial microcirculation in intact rat heart: a study with micro-CT.
In: Analytical and Quantitative Cardiology, edited by Sideman S,
and Beyar R.. New York: Plenum, 1997, chapt 14, p. 165-176.
8.
Bernard, SL,
Glenny RW,
Erickson HH,
Fedde MR,
Polissar N,
Basaraba RJ,
and
Hlastala MP.
Minimal redistribution of pulmonary blood flow with exercise in racehorses.
J Appl Physiol
81:
1062-1070,
1996
9.
Block, M,
Liu Y,
Harris LD,
Robb RA,
and
Ritman EL.
Quantitative analysis of a vascular tree model with the dynamic spatial reconstructor.
J Comput Assist Tomogr
8:
390-400,
1984[Web of Science][Medline].
10.
Brace, RA.
Fitting straight lines to experimental data.
Am J Physiol Regulatory Integrative Comp Physiol
233:
R94-R99,
1977
11.
Caruthers, SD,
and
Harris TR.
Effects of pulmonary blood flow on the fractal nature of flow heterogeneity in sheep lungs.
J Appl Physiol
77:
1474-1479,
1994
12.
Coflesky, JT,
Jones RC,
Reid LM,
and
Evans JN.
Mechanical properties and structure of isolated pulmonary arteries remodeled by chronic hyperoxia.
Am Rev Respir Dis
136:
388-394,
1987[Web of Science][Medline].
13.
Corrigan, NM,
Chavez AE,
Wisner ER,
and
Boone JM.
A multiple detector array helical X-ray microtomography system for specimen imaging.
Med Phys
26:
1708-1713,
1999[Web of Science][Medline].
14.
Cowley, AW.
The emergence of physiological genomics.
J Vasc Res
36:
83-90,
1999[Web of Science][Medline].
15.
Cowley, AW,
Stoll M,
Greene AS,
Kaldunski ML,
Roman RJ,
Tonellato PF,
Schork NJ,
Dumas P,
and
Jacob HJ.
Genetically defined risk of salt sensitivity in an intercross of Brown Norway and Dahl S rats.
Physiol Genomics
2:
107-115,
2000
16.
Dao, N,
McCormick PJ,
and
Dewey CF, Jr.
The human physiome as an information environment.
Ann Biomed Eng
28:
1032-1042,
2000[Web of Science][Medline].
17.
Davies, P,
Maddalo F,
and
Reid L.
Effects of chronic hypoxia on structure and reactivity of rat lung microvessels.
J Appl Physiol
58:
795-801,
1985
18.
Dawson, CA,
and
Haworth ST.
A pulmonary circulation physiome (Abstract).
Ann Biomed Eng
28:
S-28,
2000.
19.
Dawson, CA,
Krenz GS,
Karau KL,
Haworth ST,
Hanger CC,
and
Linehan JH.
Structure-function relationships in the pulmonary arterial tree.
J Appl Physiol
86:
569-583,
1999
20.
Dawson, CD,
Krenz GS,
and
Linehan JH.
Complexity and structure-function relationships in the pulmonary arterial tree.
In: Lung Biology in Health and Disease, Complexity in Structure and Function of the Lung, edited by Hlastala MP,
and Robertson HT.. New York: Dekker, 1998, chapt. 13, p. 401-427.
21.
Dumas, P,
Granados M,
Dawson C,
and
Jacob H.
Genetic dissection of pulmonary hypertension in fawn-hooded rats by chromosomal substitution.
FASEB J
15:
A857,
2001.
22.
Feldkamp, LA,
Davis LC,
and
Kress JW.
Practical cone-beam algorithm.
J Opt Soc Am A
1:
612-619,
1984[Web of Science].
23.
Jacob, HJ.
Functional genomics and rat models.
Genome Res
9:
1013-1016,
1999
24.
Johnson, RH,
Hu H,
Haworth ST,
Cho PS,
Dawson CA,
and
Linehan JH.
Feldkamp and circle-and-line cone-beam reconstruction for 3D micro-CT of vascular networks.
Phys Med Biol
43:
929-940,
1998[Web of Science][Medline].
25.
Jorgensen, SM,
Demirkaya O,
and
Ritman EL.
Three-dimensional imaging of vasculature and parenchyma in intact rodent organs with X-ray micro-CT.
Am J Physiol Heart Circ Physiol
275:
H1103-H1114,
1998
26.
Karau, K,
Johnson R,
Molthen R,
Provost A,
Cowley A,
Jacob H,
and
Dawson C.
Pulmonary arterial remodeling in the fawn hooded rat revealed by microfocal X-ray tomography (CT) (Abstract).
FASEB J
14:
A711,
2000.
27.
Lai-Fook, SJ.
A continuum mechanics analysis of pulmonary vascular interdependence in isolated dog lobes.
J Appl Physiol
46:
419-429,
1979
28.
Leach, RM,
Twort CHC,
Cameron IR,
and
Ward JPT
A comparison of the pharmacological and mechanical properties in vitro of large and small pulmonary arteries of the rat.
Clin Sci (Lond)
82:
55-62,
1992[Medline].
29.
Lecras, TD,
Kim DH,
Gebb S,
Markham NE,
Shannon JM,
Tuder RM,
and
Abman SH.
Abnormal lung growth and the development of pulmonary hypertension in the fawn-hooded rat.
Am J Physiol Lung Cell Mol Physiol
277:
L709-L718,
1999
30.
Linehan, JH,
Haworth ST,
Nelin LD,
Krenz GS,
and
Dawson CA.
A simple, distensible-vessel model for interpreting pulmonary vascular pressure-flow curves.
J Appl Physiol
73:
987-994,
1992
31.
Liu, YH,
Hoffman EA,
and
Ritman EL.
Measurement of three-dimensional anatomy and function of pulmonary arteries with high-speed X-ray computed tomography.
Invest Radiol
22:
28-36,
1987[Web of Science][Medline].
32.
Madden, JA,
Keller PA,
Effros RM,
Seavitte C,
Choy JS,
and
Hacker AD.
Responses to pressure and vasoactive agents by isolated pulmonary arteries from monocrotaline-treated rats.
J Appl Physiol
76:
1589-1593,
1994
33.
Maloney, JC,
Roohalamini SA,
and
Wexler L.
Pressure-diameter relations of small blood vessels in isolated dog lung.
Microvasc Res
2:
1-12,
1970[Medline].
34.
Marsaglia, G.
Ratios of normal variables and ratios of normal means.
J Amer Stat Assoc
60:
193-204,
1965.
35.
Parker, JC,
Ardell JL,
Hamm CR,
Barman SA,
and
Coker PJ.
Regional pulmonary blood flow during rest, tilt, and exercise in unanesthetized dogs.
J Appl Physiol
78:
838-846,
1995
36.
Paulus, MJ,
Gleason SS,
Kennel SJ,
Hunsicker PR,
and
Johnson DK.
High resolution X-ray computed tomography: an emerging tool for small animal cancer research.
Neoplasia
2:
62-70,
2000[Web of Science][Medline].
37.
Pearlman, JD,
Leavitt M,
and
Newell JB.
A priori information in image analysis: assessment of intensity distribution for definition of shape and size of small vessels.
IEEE Trans Med Imaging
9:
461-465,
1990[Medline].
38.
Rabinovitch, M,
Konstam MA,
Gamble WJ,
Papanicolaou N,
Aronovitz MJ,
Treves S,
and
Reid L.
Changes in pulmonary blood flow affect vascular response to chronic hypoxia in rats.
Circ Res
52:
432-441,
1983
39.
Reinhardt, JM,
D'Souza ND,
and
Hoffman EA.
Accurate measurement of intrathoracic airways.
IEEE Trans Med Imaging
16:
820-827,
1997[Web of Science][Medline].
40.
Rudy, Y.
From genome to physiome: integrative models of cardiac excitation.
Ann Biomed Eng
28:
945-950,
2000[Web of Science][Medline].
41.
Saad, Y,
Garret MR,
and
Rapp JP.
Multiple blood pressure QTL on rat chromosome 1 defined by Dahl rat congenic strains.
Physiol Genomics
4:
201-214,
2001
42.
Sato, K,
Webb S,
Tucker A,
Rabinovitch M,
O'Brien RF,
McMurtry IF,
and
Stelzner TJ.
Factors influencing the idiopathic development of pulmonary hypertension in the fawn hooded rat.
Am Rev Respir Dis
145:
793-797,
1992[Web of Science][Medline].
43.
Shepp, LA,
and
Logan BF.
The Fourier reconstruction of a head section.
IEEE Trans Nuclear Science
NS-21:
21-42,
1979.
44.
Sinclair, SE,
McKinney S,
Glenny RW,
Bernard SL,
and
Hlastala MP.
Exercise alters fractal dimension and spatial correlation of pulmonary blood flow in the horse.
J Appl Physiol
88:
2269-2278,
2000
45.
Smith, JC,
and
Mitzner W.
Analysis of pulmonary vascular interdependence in excised dog lobes.
J Appl Physiol
48:
450-467,
1980
46.
Steen, RG,
Kwitek-Black AE,
Glenn C,
Gullings-Handley J,
Vanetten W,
Atkinson OS,
Appel D,
Twigger S,
Muir M,
Mull T,
Granados M,
Kissebah M,
Russo K,
Crane R,
Popp M,
Peden M,
Matise T,
Brown DM,
Lu J,
Kingsmore S,
Tonellato PF,
Rozen S,
Slonim D,
Young P,
and
Jacob HJ.
A high-density integrated genetic linkage and radiation hybrid map of the laboratory rat.
Genome Res
9:
AP1-8,
1999.
47.
Stelzner, T,
Hofmann TA,
Brown D,
Deng A,
and
Jacob HJ.
Genetic determinants of pulmonary hypertension in fawn-hooded rats.
Chest
111:
96S,
1997
48.
Wood, SA,
Zerhouni EA,
Hoford JD,
Hoffman EA,
and
Mitzner W.
Measurement of three-dimensional lung tree structures by using computed tomography.
J Appl Physiol
79:
1687-1697,
1995
49.
Yen, RT,
Fung YC,
and
Bingham N.
Elasticity of small pulmonary arteries in the cat.
J Biomed Eng
102:
170-177,
1980.
50.
Zhuang, FY,
Fung YC,
and
Yen RT.
Analysis of blood flow in cat's lung with detailed anatomical and elasticity data.
J Appl Physiol
55:
1341-1348,
1983
This article has been cited by other articles:
![]() |
H. D. Litzlbauer, C. Neuhaeuser, A. Moell, S. Greschus, A. Breithecker, F. E. Franke, W. Kummer, and W. S. Rau Three-dimensional imaging and morphometric analysis of alveolar tissue from microfocal X-ray-computed tomography Am J Physiol Lung Cell Mol Physiol, September 1, 2006; 291(3): L535 - L545. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Watz, A. Breithecker, W. S. Rau, and A. Kriete Micro-CT of the Human Lung: Imaging of Alveoli and Virtual Endoscopy of an Alveolar Duct in a Normal Lung and in a Lung with Centrilobular Emphysema--Initial Observations Radiology, September 1, 2005; 236(3): 1053 - 1058. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. S. Burrowes, P. J. Hunter, and M. H. Tawhai Anatomically based finite element models of the human pulmonary arterial and venous trees including supernumerary vessels J Appl Physiol, August 1, 2005; 99(2): 731 - 738. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. F. LaDisa Jr., Lars. E. Olson, R. C. Molthen, D. A. Hettrick, P. F. Pratt, M. D. Hardel, J. R. Kersten, D. C. Warltier, and P. S. Pagel Alterations in wall shear stress predict sites of neointimal hyperplasia after stent implantation in rabbit iliac arteries Am J Physiol Heart Circ Physiol, May 1, 2005; 288(5): H2465 - H2475. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. T. Reeves, J. H. Linehan, and K. R. Stenmark Distensibility of the normal human lung circulation during exercise Am J Physiol Lung Cell Mol Physiol, March 1, 2005; 288(3): L419 - L425. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. C. Molthen, K. L. Karau, and C. A. Dawson Quantitative models of the rat pulmonary arterial tree morphometry applied to hypoxia-induced arterial remodeling J Appl Physiol, December 1, 2004; 97(6): 2372 - 2384. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. C. Langheinrich, B. Leithauser, S. Greschus, S. von Gerlach, A. Breithecker, F. R. Matthias, W. S. Rau, and R. M. Bohle Acute Rat Lung Injury: Feasibility of Assessment with Micro-CT Radiology, October 1, 2004; 233(1): 165 - 171. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. S. Krenz and C. A. Dawson Flow and pressure distributions in vascular networks consisting of distensible vessels Am J Physiol Heart Circ Physiol, June 1, 2003; 284(6): H2192 - H2203. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |