Characterization of Forested Landscapes from Remotely Sensed Data Using Fractals and Spatial Autocorrelation

Abstract

The characterization of forested landscapes is frequently required in civil engineering practice. In this study, some spatial analysis techniques are presented that might be employed with Landsat TM data to analyze forest structure characteristics. A case study is presented wherein fractal dimensions (FDs), along with a simple spatial autocorrelation technique (Moran's I ), were related to stand density parameters of the Oakmulgee National Forest located in the southeastern United States (Alabama). The results indicate that when smaller trees do not dominate the landscape (<50%), forested areas can be differentiated according to breast sizes and thus important flood plain characteristics such as ratio of obstructed area to total area can be estimated from remotely sensed data using the studied indices. This would facilitate the estimation of hydraulic roughness coefficients for computation of flood profiles needed for bridge design. FD and Moran's I remained fairly constant around the values of 2.7 and 0.9 (resp.) for samples with either greater than 50% saplings or less than 50% sawtimber and with ranges of 2.7-2.9 and 0.6-0.9 as the saplings decreased or the sawtimber increased. Those indices can also distinguish hardwood and softwood species facilitating forested landscapes mapping for preliminary environmental impact analysis.

(ProQuest: ... denotes non-US-ASCII text omitted.)

Mohammad Z. Al-Hamdan 1 and James F. Cruise 2 and Douglas L. Rickman 3 and Dale A. Quattrochi 3

Recommended by Kirk Hatfield

1, Universities Space Research Association at NASA Marshall Space Flight Center, National Space Science and Technology Center, NASA Global Hydrology and Climate Center, Huntsville, AL 35805, USA
2, Civil and Environmental Engineering Department, University of Alabama in Huntsville, Huntsville, AL 35899, USA
3, Earth Science Office at NASA Marshall Space Flight Center, National Space Science and Technology Center, NASA Global Hydrology and Climate Center, Huntsville, AL 35805, USA

Received 5 August 2011; Accepted 27 December 2011

1. Introduction

The characterization of forested landscapes is frequently required in civil engineering practice. Examples include estimation of quantities for clearing prior to construction projects, environmental impact analysis, and characterization of hydraulic roughness for flood plain studies. Increasingly, professionals are considering the use of remotely sensed data to aid in these estimations. In 2003, the Federal Highway Administration funded a project, conducted jointly by Mississippi State University and the University of Alabama in Huntsville, AL, USA. to evaluate the efficacy of using remotely sensed data in the planning and environmental impact analysis of transportation systems. Since transportation lines (both road and rails) frequently must traverse forested areas, one major focus of the project was the use of remotely sensed data to characterize these landscapes.

The specific transportation needs that directed the project were to characterize the hydraulic roughness of flood plains for bridge design and environmental impact analysis. Recent studies have confirmed the importance of rigid, unsubmerged vegetation in determining flow depths and velocities in shallow flow situations such as commonly encountered in wide flood plains [1-4]. The studies highlight the importance of stand density (ratio of obstructed area to total area) and trunk diameter. Musleh and Cruise [1] showed the relationship between these variables and classical hydraulic roughness parameters such as Manning's n and the Darcy f factor. Accurate land cover and species identification are also important aspects of environmental impact analysis, particularly in identification of wetlands under the forest canopy.

Several sources of remotely sensed data are currently available that might be useful for forest characterization purposes. The data can be from satellite or aircraft platforms and can be from either passive or active instruments. A large amount of research has been performed using active remotely sensed data, particularly airborne radar to estimate forest parameters [5-11]. Most passive remotely sensed data are much more easily accessible and cost effective than are active data. A significant amount of research has been performed on forest biomass estimation using passive instruments, particularly radiometric data [12-24]. However, studies that employ passive radiometric data (e.g., Landsat Thematic Mapper (TM), NOAA Advanced Very High-Resolution Radiometer (AVHRR), or the Moderate-Resolution Imaging Spectroradiometer (MODIS)) usually focus on the estimation of indirect measurement of biomass or canopy coverage such as the Leaf Area Index (LAI) or Normalized Difference Vegetation Index (NDVI). Thus, it would be of great benefit if passive radiometer data could be employed to characterize forest structure such as stand density and trunk size.

In this study, some spatial analysis techniques are presented that might be employed with Landsat TM data to analyze forest structure characteristics. The spectral characteristics (7 bands) and spatial resolution (30 m) of TM data make it very suitable for use in the analysis of even moderate-sized forested areas. In an earlier paper by Al-Hamdan et al. [25], the authors examined the impact of spatial and spectral resolution of remotely sensed data on the spatial indices that might be used for landscape characterization. That study concluded that Landsat TM might possess the ideal attributes for this purpose. A case study is now presented wherein fractal dimensions, along with a simple spatial correlation technique, were related to stand density parameters of the Oakmulgee National Forest located in the southeastern United States (Alabama).

2. Fractal Analysis and Spatial Autocorrelation Methods

One of the most essential issues in interpretation and analysis of remotely sensed data is the observation and measurement scale [26]. Scale is crucial to the characterization of geospatial data because many environmental processes and patterns are scale-dependent [27]. Recently developed spatial analysis approaches from a variety of science disciplines offer the possibility of highly efficient statistical characterization, analysis, and identification of spatial data in remotely sensed images. Spatial autocorrelation and fractal measurement are methods that have been developed to characterize the scaling property of spatial data [28]. Although the fractal technique and other textural analyses have been applied extensively [29-42], its use as a spatial technique for characterizing remote sensing images needs to be evaluated more in terms of its utility to help characterize forest growth characteristics such as stand sizes [25, 26, 43, 44], as demonstrated in this study.

2.1. Fractals

Classical geometry cannot provide tools for analysis of the forms of most spatial patterns of nature because they are so irregular and fragmented [44]. Fractal geometry was developed as a mathematical basis for characterizing complex natural patterns [45]. In classical geometry, a point has an integer topological dimension of zero, a line has one dimension, an area has two dimensions, and a volume has three dimensions [26, 44]. However, the fractal dimension (FD) is a noninteger value that exceeds the Euclidean topological dimension [44, 45]. The FD can vary between zero and one, between one and two, or between two and three; for a point pattern, a curve, or a surface, respectively [26, 44]. As the geometrical complexity of a perfectly flat two-dimensional surface ( FD=2.0 ) increases so that it begins to fill a volume, the FD values approach 3.0 [44].

The foundation for fractal analysis is self-similarity, which can be defined as a property of a curve or surface where each part is indistinguishable from the whole [26, 44-46]. In order to define the theoretical FD, the degree of self-similarity is used and expressed as a self-similarity ratio [26, 44, 45]. Thus, the FD of a curve can be defined as [26, 44] [figure omitted; refer to PDF] where _Ns is the number of similar copies and _rf is the scale reduction factor.

Measuring the length of the curve using various step sizes by a procedure called the walking-divider method [46] is a common way to estimate the FD value of a curve (e.g., a coastline) [26, 44]. In the case of irregular curves, the length increases as the measurement interval decreases, and a linear regression can capture such an inverse relationship between total line length and step size [26, 44, 47]: [figure omitted; refer to PDF] where L is the line length, _Ss is the step size, B is the slope of the regression, and C is a constant. For the case of a curve, FD can then be calculated by [figure omitted; refer to PDF]

For a raster-based remotely sensed image (surface), FD can be estimated in a similar approach using a method that is called the isarithm method [44], which was evolved from Goodchild [46], Shelberg et al. [48], and Lam and De Cola [49]. With the isarithm method, a mean FD from individual FD values of gray-scale contours is computed [26]. For each isarithm brightness value and each step size, the algorithm classifies each pixel below the isarithm value as white and each pixel above this value as black [26, 44]. Each neighboring pixel along the rows or columns is then compared to determine whether the pairs are both black or both white; if they are not of the same color, then an isarithm lies between the two neighboring pixels. The total number of boundary pixels is used to approximate the length of each isarithm line [26, 47]. A linear regression is performed using the logarithms of the total length of the boundary and the step size. The regression slope B is used to determine the FD of the isarithm line [26, 44], where [figure omitted; refer to PDF] Equation (4) differs from (3) because as a flat surface grows more complex, FD increases from a value of 2.0 and approaches 3.0 as the surface begins to fill a volume [26]. The final FD of the surface is the average of the FD values for those isarithms having a coefficient of determination " ^R2 " greater than or equal to 0.9 [26, 44, 47, 50].

It has been shown in many previous forestry research studies [51-60] that there is a significant positive relationship between crown width and the Diameter at Breast Height (DBH) for both hardwood and softwood species. Thus, an increase in the stand diameter implies an increase in the crown width and vice versa. As illustrated in Figure 1 [25], if continuous small crown trees are covering two adjacent remotely sensed pixels of a similar area, the result is two homogenous surfaces. Thus, the integration of the brightness levels within each pixel (i.e., pixel value) will be similar in magnitude. If the pixel values do not vary significantly, the result is less complexity in terms of fractals (smaller FD) and more homogeneity in terms of autocorrelation [25]. On the other hand, if the pixels are covered with large crown trees, the result is nonhomogenous adjacent pixels leading to more complexity in terms of fractals and less homogeneity in terms of autocorrelation [25].

Size class effect on remotely sensed data: (a) small crown trees and (b) large crown trees (courtesy of Al-Hamdan et al. [25]).

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

2.2. Spatial Autocorrelation

Another index that can be used to analyze the spatial autocorrelation of images is Moran's I . Moran's I [61] is an index of spatial autocorrelation which reflects the differing spatial structures of the smooth and rough surfaces [47]. Thus, it is a potential technique in characterizing and estimating measures of forest-related surface roughness based on crown diversity. Moran's I is a spatial correlation statistic and is calculated from the following formula [47]: [figure omitted; refer to PDF] where I(d) is Moran's spatial autocorrelation at distance d,_wij is weight at distance d , so that _wij is 1 if point j is within distance d of point i , otherwise _wij =0, _zi is deviation (i.e., _zi = _x-xmean for variable x ). W = the sum of all the weights where i...0;j .

Moran's I can vary from +1.0 for perfect positive autocorrelation (a clumped pattern) to -1.0 for perfect negative autocorrelation (a checker board pattern) [44]. Moran's I is different from the FD in that the FD is focused on the object shape, size, and the tortuosity of the edges of these objects [47]. Moran's I does not explicitly consider the shapes and sizes of objects once the weights, _wij , in (5) are determined [47].

3. Case Study

Oakmulgee National Forest is located in central Alabama and encompasses an area of 128,638 ha (317,861 acres) (Figure 2). According to the US Forest Service inventory analysis, there are three size classes present within the forest data sets, namely, sawtimber, poletimber, and saplings. This size classification was based on the diameter at breast height (DBH) of the tree trunk. The DBH values for those classes are greater than 22.86 cm (9.0 in), from 12.7 cm (5.0 in) to 22.6 cm (8.9 in), and from 2.54 cm (1.0 in) to 12.45 cm (4.9 in), respectively. Forest species includes softwood and hardwood trees. Longleaf-slash pine, shortleaf-loblolly, and cypress are examples of softwood trees (Figure 3(a)). White oak, red oak, hickory, sweetgum, ash, and yellow-poplar are examples of hardwood trees (Figure 3(b)). The forest database included digital maps linked to attribute tables containing the stand characteristics (i.e., tree types, sizes, and species) of the forest.

Figure 2: Location map of study area.

[figure omitted; refer to PDF]

Examples of (a) Softwood trees and (b) Hardwood trees.

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

Other spatial data were also collected to characterize the region for geoidentification and location purposes. These data included state and county maps, roads, streams, and topographic images (digital line graph data) that were useful in providing visual locations of points within the forest area. The sources of these data were the Environmental Systems Research Institute (ESRI) and the US Bureau of the Census Topologically Integrated Geographic Encoding and Referencing System (TIGER).

The remote sensing data consisted of Landsat TM images of the forest area taken in summer 1999. The TM instrument records energy reflected in seven bandwidths as shown in Table 1. As shown in the table, the spatial resolution of the Landsat TM images is 30 meters except for band 6 which is 120 meters. For consistency purposes, the data recorded in Band 6 were excluded from these analyses. The source for the Landsat TM images was the EROS Data Center of the United States Geological Survey (USGS). This study focused on a comparison of samples identified from the same Landsat TM scene. Thus, the authors believed that it was not necessary to atmospherically correct the satellite data because any atmospheric effects would prevail for all of the study sites; that is, all samples would have similar atmospheric effects. Additionally, previous work employing spatial statistical techniques in the Image Characterization and Modeling System (ICAMS), particularly fractal analysis of remote sensing data, has implied that there is little influence of atmospheric effects on the overall statistical analysis results [26, 43, 62].

Table 1: Landsat Thematic Mapper (TM).

Bands	Wavelength (micrometers)	Resolution (meters)
Band 1	0.45-0.52	30
Band 2	0.52-0.60	30
Band 3	0.63-0.69	30
Band 4	0.76-0.90	30
Band 5	1.55-1.75	30
Band 6	10.40-12.50	120
Band 7	2.08-2.35	30

4. Image Characterization and Modeling System (ICAMS)

ICAMS is a software module that was developed to measure, characterize, and model multiscale remotely sensed data [62, 63]. It contains a robust set of fractal measurement algorithms embedded in a GIS-type architecture. It contains a number of spatial measurement methods that can be used to address a variety of significant issues related to scale and fractal analysis.

The main functions of ICAMS are image input, image characterization, specialized functions, and image display and output. The image input subsystem includes basic image processing functions, such as file transformation, georeferencing, image registration, and a variety of image viewing capabilities [62]. The image characterization subsystem provides users with an array of spatial and nonspatial measures for characterizing image data. The non-spatial measures include basic descriptive statistics measures (i.e., mean, mode, median, and variance) and histograms. The spatial measures include fractal analysis [49], variogram analysis [64, 65], spatial autocorrelation indices [66, 67], and textural measures such as local variance [68]. The specialized subsystem of ICAMS provides aggregation routines for aggregating pixels to simulate multiscaled data for scale effect analyses [26]. The image output subsystem contains functions to output original images or derived products in two-dimensional or three-dimensional form. ICAMS provides the ability to calculate the FD of remotely sensed images using the isarithm method [49], variogram [65], and triangular prism methods [69]. Lam et al. [70] found that the isarithm method calculates the FD fairly accurately and more so than the other two methods, thus the isarithm method was used in this study.

5. Procedure

The procedure followed in the case study is illustrated in Figure 4. Raster and vector data were collected for the Oakmulgee National Forest as described previously. The raster data (TM images) were then imported into the ER Mapper remote sensing software package in order to trim, locate, and crop the images of the study area and to export the data to the ARCINFO 8.2 GIS package. The relevant vector data (state and county lines, roads, streams, DLG (topographic) data, and National Forest delineations) were also imported into ARCINFO and overlaid onto the raster images. Figure 5 illustrates some attributes of the data set. Samples were collected randomly covering of all parts of the forest. As shown in Figures 2 and 5, the area is bisected by the Cahaba River, so care was taken to obtain a representative number of samples from both sides of the river. In addition, the elevation of the forest ranged from 60 m to 170 m above mean sea level with an average of 130 m. So, again care was taken to obtain samples representative of all topographic conditions so as not to bias the results. Criteria for the selection of the window (number of pixels in a sample) were based on the resolution, minimum mapping unit size, and nature of the classes (size of the regions, characteristic scale, directionality, and spatial periodicity) to be identified [71]. It is understood that a smaller window size does not convey sufficient spatial or texture information to characterize land surfaces. On the other hand, if the window size is too large, too much information from other surfaces could be included and hence the algorithm might not be efficient. The sample size for this study was chosen to be 100 × 100 pixels based on the average sample size used in the research literature [28, 71]. As shown in Figure 5(e), a total of 52 samples were collected for Oakmulgee National Forest. It can be noted that the sample size in the study area was larger than 30 and hence the distribution of the sample means can be approximated reasonably well by a normal distribution for statistical purposes [72].

Figure 4: Study methodology.

[figure omitted; refer to PDF]

Overlaying and sampling process for Oakmulgee National Forest: (a) Landsat TM image, (b) counties, roads, and city locations (c) DLGs, (d) forest stands, and (e) overall layers and samples.

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

(c) [figure omitted; refer to PDF]

(d) [figure omitted; refer to PDF]

(e) [figure omitted; refer to PDF]

Information was then collected for each sample including size class (sawtimber, poletimber, saplings), species groups (hardwood and softwood), age, and elevation using the national forests vector GIS data obtained from the Forest Service and the digital elevations GIS data obtained from the Earth Resources Observation Systems (EROS) Data Center. For each sample, the percentage of each size class and species group was determined using the digital map of the national forest based on the area covered by each size class and species group in the GIS. Also, the average elevation for each sample was determined using the Digital Line Graph (DLG) data in the GIS. Table 2 shows the in situ forest data and the average elevations for all samples collected from all the study area. The raster (TM) data of the area were then imported to the GIS module, ICAMS, where the spatial analytical techniques of fractals and autocorrelation were conducted and the spatial indices of FD and Moran's I were found. Figure 6 shows an example of two subset images from TM Band 1 and their associated FD and Moran's I . The FD and Moran's I values were found for all bands of Landsat TM and for all 52 samples. The averages of FD and Moran's I for each sample were then calculated using the results of all Landsat TM bands except the thermal Band 6. Finally, the results were analyzed in terms of the relationships between the image complexity indices and the forest characteristics.

Table 2: In situ forest data of Oakmulgee National Forest, AL, USA.

Sample no.	*Size classes	Species	Evaluation (m)
Sawtimber (%)	poletimber (%)	Saplings (%)	Hardwood (%)	Softwood (%)
1	65	1	34	26	74	120
2	85	11	4	50	50	90
3	83	10	7	43	57	105
4	83	10	7	43	57	143
5	74	7	19	38	62	143
6	77	9	14	41	59	165
7	79	10	11	36	64	161
8	70	6	24	31	69	135
9	77	9	14	41	59	165
10	27	0	73	22	78	113
11	82	10	8	36	64	126
12	70	6	24	30	70	170
13	77	9	14	36	64	120
14	94	6	0	75	25	135
15	60	0	40	24	76	120
16	86	14	0	49	51	129
17	86	14	0	52	48	162
18	69	6	25	33	67	135
19	80	10	10	41	59	120
20	86	12	2	47	53	145
21	95	5	0	77	23	143
22	49	0	51	0	100	144
23	86	13	1	44	56	114
24	74	8	18	39	61	165
25	88	12	0	56	44	90
26	67	3	30	33	67	165
27	14	0	86	8	92	156
28	0	0	100	23	77	140
29	67	5	28	32	68	154
30	21	0	79	0	100	150
31	86	13	1	48	52	135
32	48	0	52	21	79	144
33	42	0	58	14	86	105
34	56	0	44	26	74	120
35	77	9	14	38	62	128
36	88	12	0	50	50	98
37	34	0	66	0	100	60
38	71	7	22	39	61	150
39	64	0	36	26	74	146
40	88	12	0	73	27	123
41	70	7	23	37	63	143
42	90	10	0	52	48	105
43	47	0	53	19	81	135
44	52	0	48	28	72	110
45	57	0	43	22	78	125
46	72	7	21	32	68	115
47	82	10	8	42	58	110
48	90	10	0	74	26	126
49	65	3	32	29	7	143
50	64	0	36	27	73	114
51	76	8	16	37	63	115
52	58	0	42	25	75	120

*Sawtimber: diameter at breast height (DBH) > 9 inch, poletimber: DBH = 5-8.9 inch, and saplings: DBH = 1-4.9 inch.

Two subset images from TM Band 1 and their associated FD and Moran's I .

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

6. Results

6.1. Analysis of Variance (ANOVA) Tests

First, the ability of the spatial complexity indices to distinguish between the different size classes and species groups was examined. ANOVA is a technique that is useful to test for the equality of several means simultaneously. It is a method for splitting the total variation of the data into meaningful components that measure different sources of variation [73]. One-way ANOVA is analysis of variance when there is one independent variable while two-way ANOVA is analysis of variance when there are two independent variables.

To study the effectiveness of remote sensing analytical techniques to distinguish between two size classes (i.e., two independent factors) and also between two species groups, a number of two-way ANOVA tests [73], were conducted using the average spatial analytical indices (i.e., FD, Moran's I ) as the dependent factor and the tree size classes or species groups as the independent factor. Rejecting the null hypothesis that the column effects are equal to zero means rejecting the hypothesis that there is no difference in the percentages of the two size classes or the two species groups. As discussed previously, the average values of FD and Moran's I were computed from the samples using the results of all TM bands except band 6 (thermal infrared).

Two-way ANOVA tests were conducted for two of the size classes at a time. The significance level was chosen to be 0.05 in order to be consistent with the research literature in the areas of remote sensing and land characterization [47, 51, 74, 75]. The first test rejected the null hypothesis of no difference between the percentages of sawtimber and poletimber classes at the 0.05 significance level using the average FD as the dependent variable. Similarly, ANOVA rejected the null hypothesis of the percentage of trees in the sawtimber and saplings size classes using the average FD. ANOVA also rejected the null hypothesis for the percentage of trees in the poletimber and saplings size classes using the average FD at the 0.05 level. Thus, it appears from these results that there is sufficient fidelity in the FD computed from the TM data that one can distinguish between timber size classes with a reasonable degree of statistical confidence.

The same battery of tests was also run using the average Moran's I as the dependent variable with similar results. The two-way ANOVA test results appear to show that Moran's I values computed from TM data can also be employed to distinguish between timber size classes with reasonable statistical confidence.

Two-way ANOVA tests were also conducted using different species (i.e., softwood and hardwood) as the independent factors. The tests rejected the null hypothesis of no difference in percentage of hardwood and softwood trees in the samples using both the FD and Moran's I as dependent variables at the 0.05 significance level.

6.2. Regression Analyses

The FD is plotted against the percentage of forest trees in each class category in Figures 7 through 9. First, the figures demonstrate that the FD has definite boundaries depending on the percentage of trees in the larger (sawtimber) and smaller (saplings) classes. The maximum FD is 2.89 for a sample with a saplings percentage of 0% and sawtimber percentage of 90%. Conversely, the minimum FD was 2.67 for a sample with a sapling percentage of 66% and sawtimber percentage of 34%. The figures also illustrate that the FD remains fairly constant around the value of 2.7 for samples with either greater than 50% saplings or less than 50% sawtimber. These observations are consistent with the reasoning outlined earlier, that is, trees with smaller crown widths will produce more homogeneous canopies (smaller FD) than will trees with larger crown widths. It appears from these results that when the percentage of smaller trees (saplings) within the forest becomes greater than 50%, then the canopy image becomes so homogeneous that the lower fractal bound is reached. Thus, the FD cannot be determinative for sapling percentages beyond 50%.

Figure 7: FD versus sawtimber.

[figure omitted; refer to PDF]

However, in the region between the boundaries, there appears to be a fairly linear relationship between the FD and the percentage of class category. Samples with greater than 50% saplings or less than 50% sawtimber were excluded from the data set (a total of 10 samples) and linear functions were fitted to the remaining data as shown in Figures 10-12. Both Figures 8 and 11 show that there is considerable more scatter in the poletimber data than is present in the other two classifications. This observation is also consistent with the general concept discussed previously in that the poletimber represents a forest structure that is dominated by neither small trees nor large trees and thus these samples would demonstrate a greater range of image complexity (and greater variation in FD) than would the sawtimber or saplings samples.

Figure 8: FD versus poletimber.

[figure omitted; refer to PDF]

Figure 9: FD versus saplings.

[figure omitted; refer to PDF]

Figure 10: FD versus sawtimber (bounds excluded).

[figure omitted; refer to PDF]

Figure 11: FD versus poletimber (bounds excluded).

[figure omitted; refer to PDF]

Figure 12: FD versus saplings (bounds excluded).

[figure omitted; refer to PDF]

The situation is similar in the case of Moran's I statistic (Figures 13-15). Again, the statistic appears to demonstrate definite bounds with a maximum of about 0.9 for samples with greater than 50% saplings and a minimum of about 0.6 for a sawtimber percentage of 94%. Again, as in the fractals case, there appears to be a fairly linear relationship between Moran's I and the class percentages between the boundaries (Figures 16-18). However, as in the case of fractals, Moran's I cannot be considered a definitive measure of stand density for samples with greater than 50% saplings.

Figure 13: Moran's I versus sawtimber.

[figure omitted; refer to PDF]

Figure 14: Moran's I versus poletimber.

[figure omitted; refer to PDF]

Figure 15: Moran's I versus saplings.

[figure omitted; refer to PDF]

Figure 16: Moran's I versus sawtimber (bounds excluded).

[figure omitted; refer to PDF]

Figure 17: Moran's I versus poletimber (bounds excluded).

[figure omitted; refer to PDF]

Figure 18: Moran's I versus saplings (bounds excluded).

[figure omitted; refer to PDF]

As for the species groups, there appears to be a fairly linear relationship between the FD and the species group percentages as well as between the Moran's I and the species group percentages (Figures 19-22). All the regressions showed an increase (positive slopes) in FD and a decrease (negative slopes) in Moran's I as the hardwood percentages increased while all the regressions showed a decrease (negative slopes) in FD and an increase (positive slopes) in Moran's I as the softwood percentages increased. The explanation is the same as given for the stand sizes case because softwood trees (e.g., Pine trees) are mostly with small crowns, while hardwood trees (e.g., Oak trees) likely have large crowns. As a matter of fact, the species group case had even stronger correlations with the average spatial indices than the DBH case. This can be due to the fact that remote sensing data do not measure DBH directly, but they measure crown reflectivity by satellite sensors. So, for a given species group, the reflectance value recorded by satellite sensors is a function of exposed projection area (canopy closure). The strong relationship between the spatial indices and both types of species groups hardwood and softwood is very important to the potential prediction process of the trunk diameter from the spatial indices, because it suggests that whether the tree is a softwood or a hardwood would not affect the relationship between its trunk size and the spatial indices. In other words, having different species groups would not skew the potential trunk size predictions. Thus, potential forest tree trunk size prediction models would be valid for any species groups.

Figure 19: FD versus softwood.

[figure omitted; refer to PDF]

Figure 20: FD versus hardwood.

[figure omitted; refer to PDF]

Figure 21: Moran's I versus softwood.

[figure omitted; refer to PDF]

Figure 22: Moran's I versus hardwood.

[figure omitted; refer to PDF]

7. Conclusions

The results of the case study presented herein appear to confirm the basic theory pertaining to the application of fractal and spatial analyses to forest canopies. The relationships between image complexity and forest stand characteristics appear to be confirmed up to a point. However, it appears that as the percentage of smaller diameter trees becomes greater, and particularly if it exceeds 50%, then the canopy image obtained from Landsat TM data becomes sufficiently homogeneous so that the spatial indices reach their lower limits and thus are no longer determinative.

It also appears, at least for the Oakmulgee forest, that the relationships between the spatial indices and forest class percentages within the boundaries can reasonably be considered linear. The linear relationship is much more pronounced in the sawtimber and saplings cases than in samples dominated by medium sized trees (poletimber). The large variation in the poletimber data is consistent with the theory of image complexity of forest surfaces as put forth in this work.

From a civil engineering applications perspective, the results indicate that up to a point where smaller trees dominate the landscape, forested areas can be differentiated according to breast sizes and thus important flood plain characteristics such as ratio of obstructed area to total area can be estimated from remotely sensed data. This would facilitate the estimation of hydraulic roughness coefficients (e.g., Manning's n ) for computation of flood profiles needed for bridge design. Even in the case where smaller trunks do dominate, then the knowledge of that fact alone would also facilitate the estimation of the required parameters. Likewise, remotely sensed indices such as fractals or Moran's I can distinguish hardwood and softwood species with good accuracy. Thus, forested landscapes can be mapped for preliminary environmental impact analysis.

These results appear to indicate that both fractal dimensions and spatial autocorrelation indices hold promise as means of estimating forest stand characteristics from remotely sensed images. However, additional work is needed to confirm that the boundaries identified for Oakmulgee National Forest and the linear nature of the relationship between image complexity indices and forest characteristics are generally evident in other forests. In addition, the effects of other parameters such as topographic relief and image distortion due to sun angle and cloud cover, for example, need to be examined.

Acknowledgments

This work was supported by a grant from the US Department of Transportation "Applications of remote sensing and related spatial technologies to environmental assessments in transportation." The assistance of the project manager, Dr. Roger King of Mississippi State University, and the USDOT program director, Dr. K. Thirumalai is gratefully acknowledged. The authors would also like to thank Dr. Nina Lam of Louisiana State University and Dr. Charles Emerson of Western Michigan University for their assistance.

References

[1] F. A. Musleh, J. F. Cruise, "Functional relationships of resistance in wide flood plains with rigid unsubmerged vegetation," Journal of Hydraulic Engineering , vol. 132, no. 2, pp. 163-171, 2006.

[2] D. S. Rhee, H. Woo, B. A. Kwon, H. K. Ahn, "Hydraulic resistance of some selected vegetation in open channel flows," River Research and Applications , vol. 24, no. 5, pp. 673-687, 2008.

[3] D. Poggi, C. Krug, G. G. Katul, "Hydraulic resistance of submerged rigid vegetation derived from first-order closure models," Water Resources Research , vol. 45, no. 10, 2009.

[4] E. Q. Hui, X. E. Hu, C. B. Jiang, F. K. Ma, Z. D. Zhu, "A study of drag coefficient related with vegetation based on the flume experiment," Journal of Hydrodynamics , vol. 22, no. 3, pp. 329-337, 2010.

[5] M. L. R. Sarker, J. E. Nichol, "Biomass estimation using texture measurements of dual polarization SAR data," in Proceedings of the Remote Sensing and Photogrammetry Society Conference, May, 2010.

[6] R. M. Lucas, A. A. Held, S. R. Phinn, S. Saatchi, S. L. Ustin, "Tropical forests," Remote Sensing for Natural Resource Management and Environmental Monitoring , pp. 239-315, John Wiley & Sons, Hoboken, NJ, USA, 2004.

[7] J. R. Santos, M. S. Pardi Lacruz, L. S. Araujo, M. Keil, "Savanna and tropical rainforest biomass estimation and spatialization using JERS-1 data," International Journal of Remote Sensing , vol. 23, no. 7, pp. 1217-1229, 2002.

[8] J. R. Santos, C. C. Freitas, L. S. Araujo, L. V. Dutra, J. C. Mura, F. F. Gama, L. S. Soler, S. J. S. Sant'Anna, "Airborne P-band SAR applied to the aboveground biomass studies in the Brazilian tropical rainforest," Remote Sensing of Environment , vol. 87, no. 4, pp. 482-493, 2003.

[9] J. E. S. Fransson, H. Israelsson, "Estimation of stem volume in boreal forests using ERS-1 C- and JERS-1 L-band SAR data," International Journal of Remote Sensing , vol. 20, no. 1, pp. 123-137, 1999.

[10] K. J. Ranson, G. Sun, "An evaluation of AIRSAR and SIR-C/X-SAR images for mapping northern forest attributes in Maine, USA," Remote Sensing of Environment , vol. 59, no. 2, pp. 203-222, 1997.

[11] P. A. Harrell, L. L. Bourgeau-Chavez, E. S. Kasischke, N. H. F. French, N. L. Christensen, "Sensitivity of ERS-1 and JERS-1 radar data to biomass and stand structure in Alaskan boreal forest," Remote Sensing of Environment , vol. 54, no. 3, pp. 247-260, 1995.

[12] J. E. Nichol, M. L. R. Sarker, W. Wagner, B. Székely, "Efficiency of texture measurement from two optical sensors from improved biomass estimation," in Proceedings of the 7th ISPRS TC Symposium--100 Years ISPRS, vol. 38, of IAPRS, Vienna, Austria, July, 2010.

[13] T. Kajisa, T. Murakami, N. Mizoue, N. Top, S. Yoshida, "Object-based forest biomass estimation using Landsat ETM+ in Kampong Thom Province, Cambodia," Journal of Forest Research , vol. 14, no. 4, pp. 203-211, 2009.

[14] D. Lu, "Aboveground biomass estimation using Landsat TM data in the Brazilian Amazon," International Journal of Remote Sensing , vol. 26, no. 12, pp. 2509-2525, 2005.

[15] A. Baccini, M. A. Friedl, C. E. Woodcock, R. Warbington, "Forest biomass estimation over regional scales using multisource data," Geophysical Research Letters , vol. 31, no. 10, 2004.

[16] O. N. Krankina, M. E. Harmon, W. B. Cohen, D. R. Oetter, Z. Olga, M. V. Duane, "Carbon stores, sinks, and sources in forests of northwestern Russia: can we reconcile forest inventories with remote sensing results?," Climatic Change , vol. 67, no. 2-3, pp. 257-272, 2004.

[17] D. P. Turner, M. Guzy, M. A. Lefsky, W. D. Ritts, S. van Tuyl, B. E. Law, "Monitoring forest carbon sequestration with remote sensing and carbon cycle modeling," Environmental Management , vol. 33, no. 4, pp. 457-466, 2004.

[18] G. M. Foody, D. S. Boyd, M. E. J. Cutler, "Predictive relations of tropical forest biomass from Landsat TM data and their transferability between regions," Remote Sensing of Environment , vol. 85, no. 4, pp. 463-474, 2003.

[19] Å Rosenqvist, A. Milne, R. Lucas, M. Imhoff, C. Dobson, "A review of remote sensing technology in support of the Kyoto Protocol," Environmental Science and Policy , vol. 6, no. 5, pp. 441-455, 2003.

[20] J. Dong, R. K. Kaufmann, R. B. Myneni, C. J. Tucker, P. E. Kauppi, J. Liski, W. Buermann, V. Alexeyev, M. K. Hughes, "Remote sensing estimates of boreal and temperate forest woody biomass: carbon pools, sources, and sinks," Remote Sensing of Environment , vol. 84, no. 3, pp. 393-410, 2003.

[21] R. F. Nelson, D. S. Kimes, W. A. Salas, M. Routhier, "Secondary forest age and tropical forest biomass estimation using thematic mapper imagery," BioScience , vol. 50, no. 5, pp. 419-431, 2000.

[22] T. Häme, A. Salli, K. Andersson, A. Lohi, "A new methodology for the estimation of biomass of conifer-dominated boreal forest using NOAA AVHRR data," International Journal of Remote Sensing , vol. 18, no. 15, pp. 3211-3243, 1997.

[23] G. L. Anderson, J. D. Hanson, R. J. Haas, "Evaluating landsat thematic mapper derived vegetation indices for estimating above-ground biomass on semiarid rangelands," Remote Sensing of Environment , vol. 45, no. 2, pp. 165-175, 1993.

[24] P. J. Curran, J. L. Dungan, H. L. Gholz, "Seasonal LAI in slash pine estimated with landsat TM," Remote Sensing of Environment , vol. 39, no. 1, pp. 3-13, 1992.

[25] M. Al-Hamdan, J. Cruise, D. Rickman, D. Quattrochi, "Effects of spatial and spectral resolutions on fractal dimensions in forested landscapes," Remote Sensing , vol. 2, no. 3, pp. 611-640, 2010.

[26] D. A. Quattrochi, C. W. Emerson, N. S.-N. Lam, H. L. Qiu, N. J. Tate, P. M. Atkinson, "Fractal characterization of multitemporal remote sensing data," Modeling Scale in Geographic Information Science , pp. 13-33, John Wiley & Sons, Hoboken, NJ, USA, 2001.

[27] M. Goodchild, R. Haining, S. Wise, "Integrating GIS and spatial data analysis: problems and possibilities," International Journal of Geographical Information Systems , vol. 6, no. 5, pp. 407-423, 1992.

[28] W. Zhao Multiscale analysis for characterization of remotely sensed images , Ph.D. dissertation, Louisiana State University, Baton Rouge, La, USA, 2001.

[29] S. Lovejoy, D. Schertzer, "Generalized scale invariance in the atmosphere and fractal models of rain," Water Resources Research , vol. 21, no. 8, pp. 1233-1250, 1985.

[30] J. R. Krummel, R. H. Gardner, G. Sugihara, R. V. O'Neill, P. R. Coleman, "Landscape patterns in a disturbed environment," Oikos , vol. 48, no. 3, pp. 321-324, 1987.

[31] M. J. MacLennan, P. J. Howarth, "The use of fractal geometry to identify ranges of scale variance in digital remotely sensed data," in Proceedings of the 21st International Symposium on Remote Sensing Environment, Ann Arbor, Mich, USA, October 1987.

[32] M. W. Palmer, "Fractal geometry: a tool for describing spatial patterns of plant communities," Vegetatio , vol. 75, no. 1-2, pp. 91-102, 1988.

[33] W. B. Cohen, T. A. Spies, G. A. Bradshaw, "Semivariograms of digital imagery for analysis of conifer canopy structure," Remote Sensing of Environment , vol. 34, no. 3, pp. 167-178, 1990.

[34] B. Zeide, "Fractal geometry in forestry applications," Forest Ecology and Management , vol. 46, no. 3-4, pp. 179-188, 1991.

[35] N. D. Lorimer, R. G. Haight, R. A. Leary, "The fractal forest: fractal geometry and applications in forest science," General Technical Report , no. NC-170, US Department of Agriculture, Forest Service, North Central Forest Experiment Station, St. Paul, Minn, USA, 1994.

[36] B. A. St-Onge, F. Cavayas, "Automated forest structure mapping from high resolution imagery based on directional semivariogram estimates," Remote Sensing of Environment , vol. 61, no. 1, pp. 82-95, 1997.

[37] B. Barbanis, H. Varvoglis, C. L. Vozikis, "Imperfect fractal repellers and irregular families of periodic orbits in a 3-D model potential," Astronomy and Astrophysics , vol. 344, no. 3, pp. 879-890, 1999.

[38] N. Coops, D. Culvenor, "Utilizing local variance of simulated high spatial resolution imagery to predict spatial pattern of forest stands," Remote Sensing of Environment , vol. 71, no. 3, pp. 248-260, 2000.

[39] J. F. Weishampel, D. L. Urban, H. H. Shugart, "Semivariograms from a forest transect gap model compared with remotely sensed data," Journal of Vegetation Science , vol. 3, no. 4, pp. 521-526, 1992.

[40] J. F. Weishampel, J. H. Sloan, J. C. Boutet, J. R. Godin, "Mesoscale changes in textural pattern of "intact" Peruvian rainforests (1970s-1980s)," International Journal of Remote Sensing , vol. 19, no. 5, pp. 1007-1014, 1998.

[41] J. F. Weishampel, J. R. Godin, G. M. Henebry, "Pantropical dynamics of "intact" forest canopy texture," Global Ecology and Biogeography , vol. 10, no. 4, pp. 389-397, 2001.

[42] J. B. Drake, J. F. Weishampel, "Simulating vertical and horizontal multifractal patterns of a longleaf pine savanna," Ecological Modelling , vol. 145, no. 2-3, pp. 129-142, 2001.

[43] C. W. Emerson, N. S.-N. Lam, D. A. Quattrochi, "A comparison of local variance, fractal dimension, and Moran's I as aids to multispectral image classification," International Journal of Remote Sensing , vol. 26, no. 8, pp. 1575-1588, 2005.

[44] N. S.-N. Lam, H. L. Qiu, D. A. Quattrochi, C. W. Emerson, "An evaluation of fractal methods for characterizing image complexity," Cartography and Geographic Information Science , vol. 29, no. 1, pp. 25-35, 2002.

[45] B. B. Mandelbort The Fractal Geometry of Nature , W. H. Freeman, New York, NY, USA, 1983.

[46] M. F. Goodchild, "Fractals and the accuracy of geographical measures," Mathematical Geology , vol. 12, no. 2, pp. 85-98, 1980.

[47] C. W. Emerson, N. S.-N. Lam, D. A. Quattrochi, "Multi-scale fractal analysis of image texture and pattern," Photogrammetric Engineering and Remote Sensing , vol. 65, no. 1, pp. 51-61, 1999.

[48] M. C. Shelberg, N. S.-N. Lam, H. Moellering, "Measuring the fractal dimensions of surfaces," in Proceedings of the 6th International Symposium on Automated Cartography, pp. 319-328, Ottawa, Canada, October 1983.

[49] N. S.-N. Lam, L. De Cola Fractals in Geography , Prentice Hall, Englewood Cliffs, NJ, USA, 1993.

[50] S. Jaggi, D. A. Quattrochi, N. S.-N. Lam, "Implementation and operation of three fractal measurement algorithms for analysis of remote-sensing data," Computers & Geosciences , vol. 19, no. 6, pp. 745-767, 1993.

[51] D. C. Bragg, "A local basal area adjustment for crown width prediction," Northern Journal of Applied Forestry , vol. 18, no. 1, pp. 22-28, 2001.

[52] D. Oladi Developing a framework and methodology for plantation assessment using remotely sensed data , Ph.D. thesis, University of New Brunswick, New Brunswick, Canada, 1996.

[53] J. Guavin, T. Hunt, R. Tardif, "Edmondstone Applied Forestry Technology Group," Documentation on Growth Phase Model, 1993

[54] W. A. Farr, D. J. Demars, J. E. Dealy, "Height and crown width related to diameter for open-grown western hemlock and Sitka spruce," Canadian Journal of Forest Research , vol. 19, no. 9, pp. 1203-1207, 1989.

[55] P. T. Sprinz, H. E. Burkhart, "Relationships between tree crown, stem, and stand characteristics in unthinned loblolly pine plantation," Canadian Journal of Forest Research , vol. 17, pp. 534-538, 1987.

[56] E. G. Roberts, R. D. Ross, "Crown area of free-growing loblolly pine and its apparent independence of age and site," Journal of Forestry , vol. 63, no. 6, pp. 462-463, 1965.

[57] B. C. Wile, "Crown size and stem diameter in red spruce and balsam fir," Forestry Pub. , no. 1056, Canada Department of Forestry. Forestry Research Branch, Ontario, Canada, 1964.

[58] J. H. Smith, G. R. Baily, "Influence of stocking and stand density on crown widths of douglas fir and lodgepole pine," Commonwealth Forestry Review , vol. 43, no. 3, pp. 243-246, 1964.

[59] P. E. Vezina, "More about the crown competition factor," Department of Forestry. Forest Research Branch, Contribution No. 505, 1963

[60] C. O. Minor, "Stem-crown diameter relations in southern pine," Journal of forestry , vol. 49, no. 7, pp. 490-493, 1951.

[61] A. D. Cliff, J. K. Ord Spatial Autocorrelation , Pion Limited, London, UK, 1973.

[62] N. S.-N. Lam, D. A. Quattarochi, H.-L. Qiu, W. Zhao, "Environmental assessment and monitoring with image characterization and modeling system using multi-scale remote sensing data," Applied Geography Studies , vol. 2, no. 2, pp. 77-93, 1998.

[63] D. A. Quattrochi, N. S.-N. Lam, H. L. Qiu, W. Zhao, D. A. Quattroch, M. F. Goodchild, "Image characterization and modeling system (ICAMS): a geographic information system for the characterization and modeling of multiscale remote sensing data," Scale in Remote Sensing and GIS , pp. 295-307, CRC/Lewis Publishers, Boca Raton, Fla, USA, 1997.

[64] P. A. Burrough, N. S.-N. Lam, L. De Cola, "Fractals and geostatistical methods in landscape studies," Fractals in Geography , pp. 87-121, Prentice Hall, Englewood Cliffs, NJ, USA, 1993.

[65] D. M. Mark, P. B. Aronson, "Scale-dependent fractal dimensions of topographic surfaces: an empirical investigation, with applications in geomorphology and computer mapping," Mathematical Geology , vol. 16, no. 7, pp. 671-683, 1984.

[66] M. F. Goodchild Spatial Autocorrelation, Concepts and Techniques in Modern Geography , of 47, Geo Books, Norwich, UK, 1986.

[67] N. S.-N. Lam, M. Fan, K. B. Liu, "Spatial-temporal spread of the AIDS epidemic, 1982-1990: a correlogram analysis of four regions of the United States," Geographical Analysis , vol. 28, no. 2, pp. 93-107, 1996.

[68] C. E. Woodcock, A. H. Strahler, "The factor of scale in remote sensing," Remote Sensing of Environment , vol. 21, no. 3, pp. 311-332, 1987.

[69] K. C. Clarke, "Computation of the fractal dimension of topographic surfaces using the triangular prism surface area method," Computers & Geosciences , vol. 12, no. 5, pp. 713-722, 1986.

[70] N. S.-N. Lam, H.-L. Qiu, D. Quattrochi, "An evaluation of fractal surface measurement methods using ICAMS (Image Characterization and Modeling System)," in Proceedings of the ACSM/ASPRS Annual Convention, Seattle, Wash, USA, April 1997.

[71] S. W. Myint Wavelet analysis and classification of urban environment using high-resolution multispectral image data , Ph.D. dissertation, Louisiana State University, Baton Rouge, La, USA, 2001.

[72] R. E. Walpole, R. H. Myers, S. L. Myers Probability and Statistics for Engineers and Scientists , Macmillan Publishing, New York, NY, USA, 1998.

[73] R. E. Walpole Introduction to Statistics , Macmillan Publishing, New York, NY, USA, 1982.

[74] M. A. Lefsky Lidar remote sensing of canopy height profiles: application to spatial and temporal trends in canopy structure , Ph.D. thesis, The University of Virginia, Charlottesville, Va, USA, 1997.

[75] S. E. Watts Determining forest productivity and carbon dyanamics in Southeastern Ohio from remotely-sensed data , Ph.D. dissertation, The Ohio State University, Columbus, Ohio, USA, 2001.