Thursday, October 24, 2013

Collecting Data and Interpolating the battlefield



Introduction

While studying the Verdun battlefield in France, we collected large amounts of data that will be useful in analyzing and interpreting the battlefield. Evidence supporting a particular geographic research question is heavily dependent upon this data collection process. To provide empirical evidence to answer my research question, whether warfare may be considered a significant agent promoting ecological and pedological change, volumetric water content, soil temperature, and micro elevation data were collected on the Thiaumont ridge. This information will supplement the research conducted by Joseph Hupy and Thomas Koehler in explaining the change in vegetation and soil evolution. 

Methods

Volumetric water content, and temperatures were recorded on the cratered and non-cratered landscape. Moisture and temperature readings were collected in crater bottoms, crater sides, and crater tops (figure 1). Some areas of the landscape were so disturbed that crater tops were difficult to distinguish since they were often enveloped within one another. Coordinates of each record along with notes were organized using a Trimble Nomad global positioning system (Table 1). To minimize the amount of error involved with our equipment, volumetric water content was recorded five times for each crater position and the results were averaged together. The temperature probe was inserted 15 cm into the soil and left until the temperature stopped moving.

Figure 1: Collecting volumetric water content and temperature of the soil. Readings were taken on the crater bottoms, crater sides, and crater tops. A GPS was used to record the coordinates and characteristics of each site location.




Table 1: Each soil measurement was added to a table that is easily be imported into a GIS. The notes field consisted of various qualitative observations at that site location.


Along with soil moisture and temperature, information on the landscapes surface was also collected. Using a Total Station, we recorded accurate and precise micro elevation data to produce a three dimensional model of the surface. Elevation data was imported into an excel spreadsheet to include each points X,Y, and Z coordinates (Figure 2).

Figure 2: Using a Total Station, elevation data was recorded and compiled into an Excel file to be imported into a GIS. Each point represents the location where surface data was collected and contains an elevation (Z) field.

Once these fields were established, ArcMap was used to interpolate the surrounding elevations using complex algorithms to estimate values not physically recorded. The two methods of interpolation used were Kriging and Spline. Kriging uses a weighted algorithm dependent upon the distance between two points and their elevation values (Figure 3). Spline interpolation uses an algorithm that minimizes surface curvature and creates a model that passes exactly through each of the input values (ESRI) (Figure 4). Both interpolation methods provide a smooth and continuous representation of the area within each input point

Figure 3: A Kriging interpolation algorithm was used to estimate the elevation values between each recorded input. Kriging uses a weighted distance function to determine how output values are dependent upon the distance of multiple input values.

Figure 4: Spline interpolation estimates output values by altering the shape of a plane to pass exactly through each of the input points.
 
 

After interpolating the surrounding elevation values, the surface models were imported into ArcScene to be spatially examined in three dimensions. Both raster images were added to the program and processed to float on a custom surface as described in the previous post. This created a three dimensional model of the Thiaumont study area (Figures 5-6).

Figure 5: Kriging interpolation of Thiaumont ridge rendered in 3D.
Figure 6: Spline interpolation rendered in 3D. This method has a more pronounced representation of the surface because the model passes directly through each input value.

Figure 6: Spline interpolation rendered in 3D. This method has a more pronounced representation of the surface because the model passes directly through each input value.

Discussion

To fully examine the soils data, various processing tools will be used to determine the significance the cratered landscape has on the developing vegetation and soil. Soil moisture content and temperature will be compared to the surface data to determine if the artillery bombardment of the Verdun battlefield has significantly altered the landscape in a way that the recovered environment is set on a different pathway compared to the undisturbed environment. Hydrology tools will be used to examine how surface flow has been altered in heavily cratered areas.

Conclusion

The modeling of our study area using highly accurate and precise elevation data is important to be included with my soil measurements. To understand why this environment is behaving differently from its prewar environment, the data must be examined spatially. Within upcoming blog posts, I will begin quantitative analysis to provide evidence of this divergent pathway.

Wednesday, October 16, 2013

Collecting and Viewing DEM data in ArcScene

Introduction:

During the continuing government shutdown, a problem began to arise within the Geographic community. Reliance upon the United States Geological Survey (USGS) for gathering spatial information has begun to put a toll on GIS users within academia and the private sector alike. Creative methods of data collection have to be pursued to overcome such a dilemma. In this blog post I will cover the use of ASTER DEM data provided by Japan, to create a three dimensional model of the World War I Verdun battlefield.

Background Information

As mentioned above, relying on a single source for spatial data is considered poor practice within the geographic community. This reliance not only causes a dilemma in the wake of a government shutdown, but also exclude other valuable resources that should be considered.  Our search for different distributors of spatial data brought us to a collaborative project by the Ministry of Economy, Trade, and Industry of Japan (METI), and the United States National Aeronautics and Space Administration (NASA). METI and NASA are jointly developing a publicly accessible digital elevation model using the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) satellite.  ASTER provides land coverage between 83°N and 83°S with a resolution of 1 arc-second (30 meter) at a 95% confidence level (Group on Earth Observations).  This data offers a suitable alternative to the USGS in providing accurate and precise elevation data.

Methods

To begin collecting data I will direct you to the website where you will create a login account used to access the tiles of your study area. http://gdem.ersdac.jspacesystems.or.jp/login.jsp.  Once your account has been set up, navigate to the homepage by clicking ASTER GDEM at the top of the page. The tiles are accessed by selecting "Search" located on the side tab of the homepage. This should bring up an interactive map sectioned off by square tiles (Figure 1).

Figure 1: Interactive map used during the tile selection process. The four selection methods are listed at the top of the webpage under the search header. To select a method, use the tabs located above the map. 

 Use the map to zoom to your area of interest (AOI). For the purposes of my class our study area is a single tile located in Eastern France at roughly 5°W and 49°N. The interactive map has four methods of tile selection depending on the size of your study area. You can select the tiles of your AOI directly, select the tiles of your AOI by creating a polygon around it, use an existing shapefile containing the extents of your study area, or if you know the coordinates you can select the tile at that location. I used the select tiles directly method while paying attention to the coordinates of the cursor in the lower left hand corner of the map. 

Once you found the tile(s) within your AOI, click start located above the map to begin the selection process.  After selecting your tiles, click next below the map to bring up your selection queue and click next again.  If you are not already logged in, you will be prompted to do so otherwise you will be directed to the ASTER GDEM terms of use and policy agreement. It is very important that you select disaster as your purpose for the download as disaster relief/mitigation is the primary intended use of the data.  After selecting agree you can begin the downloading process. Navigate to your download and extract the files to your desired directory. This ends the data collection process for obtaining ASTER elevation data.

Now that you have collected your data, it can be imported into a GIS for manipulation and analysis. For our intended purposes we wanted to create a three dimensional model of the Verdun battlefield showing just how crucial it was for both combative forces to secure the north and south ridgelines east of the Meuse River. Since the tiles downloaded using the process above covered a larger area than necessary, the clip function was used to extract only my AOI. Using a separate raster, I designated the extent of the battlefield to create a new, smaller elevation raster. Figure 2 shows an overlay of my DEM clip on top of the ASTER tile.
Figure 2: A DEM clip (red to green diverging, dark) overlay on top of downloaded ASTER tile. The clip was created using a georeferenced raster of the Verdun battlefield

After creating the DEM clip of our AOI, it was imported into ArcScene to visualize the battlefield in three dimensions. When downloading data from the GDEM website, it comes preset to the WGS84 spatial reference. To ensure spatial integrity, the raster was re-defined and projected using the ETRS 1989 datum in UTM zone 31 North. Once the raster has been projected, it can be set to float on a custom surface within the layer properties, designating each pixel value from the DEM to represent change in elevation (Figure 3).

Figure 3: Setting base heights to float on a custom surface within the layer properties dialog box. This will create the three dimensional surface desired for the model. 


Once the raster has been set to float on a custom surface, it can be viewed in three dimensions. If for any reason your raster is not a three dimensional model, you may have to change the vertical exaggeration located in the general tab of the scene properties (Figure 4). If vertical exaggeration is set at zero, it will show no change in elevation which is obviously not an accurate depiction of Earth's surface. A vertical exaggeration set at "none" will maintain the elevational scale of 1:1; whereas, a vertical exaggeration set at "5" will embellish the changes in elevation five to one.

Figure 4: Setting the vertical exaggeration to emphasize changes in elevation. To view the model in 3D, make sure your models vertical exaggeration is not set to zero.

Results

Additional features within ArcScene may be used to enhance the quality of your model.  Changing the raster surface resolution can help in achieving either a "smoother" or more pronounced image (Figure 5 and Figure 6).
Figure 5: A smooth surface model using large X and Y cell sizes to show smooth changes in elevation
Figure 6: A rough surface model using a majority of the DEM's pixel values to depict changes in elevation.
The color scheme can also be changed to emphasize contrasts in elevation using warm and cool colors. The symbology can be changed from stretched to classified if you want to create color classes within a pre-defined elevation range. This may be useful when determining the suitability for a feature at a certain location (Figure 7 and Figure 8).  
Figure 7: Rendering the raster in a continuous, stretched color scheme to enhance the surface elevations. 

Figure 8: Rendering the raster in a classified color scheme to detect suitable locations of battlefield defensive locations.
Conclusion

Elevation data within a GIS has a number of practical uses. In future blog posts, I will continue to explore the Verdun battlefield by using various tools such as slope, aspect, and viewshed analysis.  Although I initially encountered the dilemma of losing the USGS as a data source, it was easily overcome by a quick search online. This brings me to my second point.

Collecting data from multiple sources should be considered best practice within the geospatial discipline. It offers a wide range of discretion to influence the integrity of your model.  Reliance on a single developer may be hazardous, acting as a limiting agent to your methods. Creditable alternatives exist beyond the conventional hosts such as USGS and the US Census Bureau. A quick google search came up with a number of sites distributing their data at no costs. Within the scope of my studies, these sources will be a valuable asset to my understanding of the Verdun battlefield.

ASTER information

Data sources

Monday, October 14, 2013



Creating an Online Story Map

An easy way to show spatial information highlighting a trip is to create an online story map. ESRI provides the online resources needed to do just that. In this blog post, I will walk you through the steps of creating an informative story map.

To begin you need to create an ESRI ArcGIS.com login account. ESRI offers a free 30 day trial on their website. To access the login page, visit this link: https://www.arcgis.com/home/signin.html 
If a longer account is needed, consult your schools geospatial facilitator.

Once your account has been created, you can begin accessing the software to create a story map.