![\"fiximg1\"](\"http:\/\/blogs.oregonstate.edu\/geo599spatialstatistics\/files\/2013\/06\/fiximg1.gif\")
<\/a> ![\"fiximggood\"](\"http:\/\/blogs.oregonstate.edu\/geo599spatialstatistics\/files\/2013\/06\/fiximggood.gif\")
<\/a> \n<\/a><\/p>\n
Figure 1. Satellite image before and after georegistration correction (click images to see animation) The registration process was coded and automated using the R “Raster” package.\u00a0 The first step was to sample the image with points and extract small image chips or windows from around the point in the reference image and the misregistered image (see figure below).<\/p>\n <\/p>\n Figure 2. Sample reference and misregistered image and extract image chips from around the sample points.<\/p>\n Within the registration program, for each sample point, the misregistered image chip is “slid” over the reference image chip at all defined x and y offset distances and cross correlation coeffiecients are calculated for each shift based on the intersecting pixel values.\u00a0 Please click on the figure below to see an example of the moving window cross correlation.<\/p>\n <\/p>\n Figure 3.\u00a0 Animation example of moving window cross correlation analysis between reference and misregistered image chips for iterative xy offset shifts.<\/p>\n A cross correlation surface is produced for each sample point.\u00a0 The apex of the conical peaks represent the offset that is needed to correct the misregisted image.\u00a0 A 2nd order polynomial equation is created from all legitimate peaks and used to warp the misregistered image into its correct geographic position.<\/p>\n <\/p>\n Figure 4. Example cross correlation surfaces showing the best offset matches between the reference image and misregistered image chips<\/p>\n With all images in their correct geographic position, the “LandTrendr” change detection program was applied to the long time series of Landsat imagery to identify pixels that have been disturbed.\u00a0 A discriminant analysis of empirical variables related to spectral, shape, and topographic characteristics of identified disturbances was conducted to predict disturbance agent from a training set of labeled disturbances.\u00a0 The figure below depicts a mountain pine beetle outbreak identified in central Colorado.\u00a0 Click on the image to see the outbreak start and progression (colors represent magnitude change to forest: low-high\\blue-red)<\/p>\n <\/p>\n Figure 5: Mountain pine beetle outbreak spread as captured by annual Landsat satellite imagery<\/p>\n From the outbreak progression video above, you can see that the outbreak appears to move up slope.\u00a0 To find out if this is truly the case I extracted elevation values for all insect affected pixels per year and plotted the histogram of both elevation value density and frequency.\u00a0 Please click on the images below to see the shift in elevation and area affected as the outbreak progresses.<\/p>\n Figure 6. Animated histograms depicting the progression of mountain pine beetle progression up slope.<\/p>\n I was also curious about what the annual climate was doing before and during the outbreak.\u00a0 I extracted PRISM climate data from 1985 to 2008 in the region of the outbreak and plotted it with the count of insect-disturbed pixels.\u00a0 The figure shows that maximum and minimum annual temperature begin to increase 1 to 2 standard deviations from mean about 5 years before the outbreak really takes off.\u00a0 This corresponds with a 3-4 year drop in annual precipitation.\u00a0 These conditions could have drought stressed the forests and provided a highly productive climate for the beetle to reproduce multiple times in a season and avoid freeze-kill.<\/p>\n <\/p>\n Figure 7.\u00a0 Graph showing annual deviation from mean for PRISM climate variables and insect-disturbed pixels for 23 years.<\/p>\n In closing, I found the ArcMap spatial statistics unable to work with my raster format data, but was able to make a lot of progress in data preparation and analysis and exploration of satellite image detected insect outbreaks and corresponding environmental factors.<\/p>\n <\/p>\n <\/p>\n <\/p>\n","protected":false},"excerpt":{"rendered":" My goal for this class was to learn what tools available in ArcGIS’s Spatial Statistics toolbox might be useful in helping to answer my research questions regarding\u00a0 spatial-temporal relationships between mountain pine beetle outbreak\\spread and coinciding environmental characteristics such as topography, climate, and forest attributes.\u00a0 I quickly learned that the spatial statistics tools in ArcMap… Continue reading
\n<\/a><\/p>\n<\/a><\/p>\n<\/a>\u00a0 Click image to see animation<\/p>\n<\/a><\/p>\n<\/a> click image to see animation<\/p>\n<\/a><\/a>Click images to see animation
\n<\/a><\/p>\n<\/a><\/p>\n