Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.
Comment: Migrated to Confluence 5.3

...

  1. On the Desktopdouble-click the Computer icon > gisdata This PC > GISData (\\file-rnassmb.rdf.rice.edu\research\FondrenGDC) (RO:) >  > GDCTraining > 1_Short_Courses > > Data_Interpolation_and_Extraction.
  2. To create a personal copy of the tutorial data, drag the Interpolation folder onto the Desktop.
  3. Close all windows.

...

  1. At the top of the Contents pane, click the List By Data Source button.



  2. If you cannot see the entire file paths, hover over the right edge of the Contents pane until your cursor turns into a double-sided arrow.
  3. Click and drag the right edge of the Contents pane to the right until you can see the full file paths.

    Image Added

Now you can see the name of the Excel file from which each Sheet1$ originated, but to avoid further confusion, you will rename the Excel worksheet layers. This action will not rename the worksheets in the actual Excel file, but only how they are listed in the Contents pane in this particular Map.

...

  1. In the Contents pane, right-click the Sheet1$_Layer layer and select Data > Export Features.
  2. For ‘Output Feature Class:’, rename the feature class from “Sheet1Sheet1_Features_CopyFeatures” to CopyFeatures to Rainfall”.
  3. Ensure your ‘Copy Features’ window appears as shown below and click Run.

...

  1. At the top of the Contents pane, click the List by Drawing Order button.



  2. Click and drag the right edge of the Contents pane back to the left, so that it is only wide enough to display the short layer names and the List by... icons.

Interpolating Data

Interpolation allows you to estimate the values of new data points based on the values of existing data points. In this case, you would like to estimate the rainfall values at each of the house locations based on the rainfall values recorded at each of the rain gauge locations. The following section will teach you several methods of interpolating both point and line vector data into raster datasets, which provide estimated values of the data variables across the entire region. You will see how these different interpolation techniques lead to variations in the resulting raster datasets. It is important to select an interpolation technique best suited to your existing data and your desired analysis.

...

  1. At the top of the Geoprocessing pane, click the Back arrow button.
  2. At the top of the Geoprocessing pane, click the Toolboxes tab.
  3. Click the Spatial Analyst Tools toolbox > Interpolation toolset > IDW (Inverse Distance Weighted) tool.
  4. For ‘Input point features’, use the drop-down menu to select the Rainfall layer, which contains your rain gauge point locations.
  5. For ‘Z value field’, use the drop-down menu to select the Rainfall field that contains the rainfall measurements you wish to interpolate.
  6. For ‘Output raster’, rename the exported raster from “IdwIdw_Rainfall1” to Rainfall1 to IDW_Rain”.
  7.  Ensure that your 'IDW' tool Parameters tab appears as shown below and click Run.

...

  1. In the Contents pane, click the IDW_Rain layer to select it.
  2. In the ribbon on the top, click the Appearance tab
  3. In the Effects group, adjust the Layer Transparency slider to 40% 30%.

Notice that you can now see the basemap beneath the raster layer.

...

  1.  In the Contents pane, right-click the IDW_Rain layer and select Symbology.
  2. For ‘Symbology,’ use the drop-down menu to select Stretch.


    For ‘Color scheme’, use
  3. Click the drop-down Color scheme menu. Select the Show names check box at the bottom.
  4. Scroll down the color scheme list and select the Blues (Continuous) box to select a continuous blue color ramp.

You can now see the continuously interpolated raster containing rainfall values that were generated using the data collected at the specific rain gauge points. Notice how this IDW method results in distinct rings around each input data point, because it treats each data point as a local max or min.

...

  1. At the bottom of the Symbology pane, click the Geoprocessing tab.
  2. At the top of the Geoprocessing pane, click the Back arrow button.
  3. Within the Interpolation toolset, click the Natural Neighbor tool.
  4. Beginning with step 5, repeat all steps in the 'Inverse Distance Weighted Interpolation Method' section above, but rename the raster output “NN_Rain”.

...

  1. At the top of the Geoprocessing pane, click the Back arrow button.
  2. Collapse the Spatial Analyst Tools toolbox.
  3. Click the 3D Analyst Tools toolbox > Data Management toolset > TIN toolset > Create TIN tool.
  4. For ‘Output TIN’, typeTIN_Rain and click elsewhere to update the text field.
  5. For ‘Coordinate System’, use the drop-down menu to select the Rainfall feature class, which will import Houses layer, to ensure that the interpolated data will be in the same coordinate system used to map the original latitude and longitude data of the rain gauges.as the Houses layer, which contains the points at which we need to extract the interpolated values.

As you will see in the warning, it is recommended to only create a TIN in a projected coordinate system, rather than a geographic coordinates system. The process of projecting a layer from one coordinate system to another is outside the scope of this course, but is covered in the Introduction to Coordinate Systems and Projections course.

...

  1. At the top of the Geoprocessing pane, click the Back arrow button.
  2. Within the 3D Analyst Tools toolbox, collapse the Data Management toolset.
  3. Within the 3D Analyst Tools toolbox, click the Conversion toolset > From TIN toolset > TIN To Raster tool.
  4. For ‘Input TIN’, use the drop-down box to select the TIN_Rain TIN.
  5. For ‘Output raster’, rename the exported raster from “tintin_rain_tin” to tin to TINRast_Rain”.
  6.  Ensure that your 'TIN To Raster' tool Parameters tab appears as shown below and click Run.



  7. In the Contents pane, uncheck the TIN_Rain layer to view the new raster layer created from the TIN.
  8. Repeat the steps 11-16 above in within the 'Inverse Distance Weighted Interpolation Method' section on page 12 above to symbolize the TINRast_Rain layer as you did the other raster layers.

...

Now you have explored four methods for interpolating the rainfall data, so you will turn off the layers you have generated in this 'Point Interpolation' section.

  1. In the Contents pane, uncheck all layers, except for the Topographic basemap layer.

...

You will now try interpolating the Hurricane line feature class using the same TIN -based method you just used to interpolate the rainfall point data. Recall that the Hurricane feature class contains a Windspeed field in its attribute table, which tells you the peak average wind speed along each line at the time the hurricane landed.

  1. In the Contents pane, check the Hurricane layer to make it visible.
  2. At the top of the Geoprocessing pane, click the Back arrow button.
  3. Within the 3D Analyst Tools toolbox, collapse the Conversion toolset.
  4. Within the In ArcToolbox navigate to 3D Analyst Tools > toolbox, click the Data Management > TIN > Create TIN.
  5. For ‘Output TIN’, type “TIN_Wind” and click elsewhere to update the text field.
  6. toolset.
  7. The TIN toolset should already be expanded, but, if not, click the TIN toolset.
  8. Click the Create TIN tool.
  9. For ‘Output TIN’, typeTIN_Wind”.
  10. For ‘Coordinate System’, use the drop-down menu to select the Houses layer.
  11. For ‘Input Features’, use the drop-down menu Use the ‘Input Feature’ drop-down box to select the Hurricane feature class.
  12. Use the For ‘Height Field’ drop, use the drop-down box menu to select the Windspeed field as shown below.Use the ‘Coordinate System’ drop-down box to select the Houses layer.
  13.  Ensure that your 'Create TIN' tool Parameters tab appears as shown below and click Run.

Again, you need to convert the TIN to a raster, so that it is comparable to the other interpolation layers.

  1. Repeat steps 9-16 above 1-8 regarding TIN to raster conversion in the 'TIN Interpolation Method' section on pages 15 and 16above, but select TIN_Wind for the input TIN and name the output raster “TINRast_Wind”TINRast_Wind”and uncheck the TIN_Wind layer.
  2. In the Contents pane, click the arrow next to all of the raster and TIN layers to collapse their symbologies.

...

All of the raster interpolations you have created for wind and rain provide continuous value estimates over the entire hurricane area, but you are only interested in the data values at the particular locations of the individual houses filing insurance claims. You will extract point data for each of the house locations from each of the five raster layers you just created in a database-compatible format for the insurance company to use. Extracting data is one of the most important aspects of data manipulation, since interpolation is often useless if the results cannot be taken out of GIS and used for other procedures.

Extracting Multi Values to Points

  1. In the Toolboxes tab, click the 3D Analyst Tools toolbox to collapse it.

Before extracting data, you must define the projections of the TIN raster layers.

...

  1. At the top of the Geoprocessing pane, click the Back arrow button.
  2. Collapse the 3D Analyst Tools toolbox.
  3. Click the Spatial Analyst Tools toolbox > Extraction toolset > Extract Multi Values to Points tool.
  4. For ‘Input point features’, use the drop-down box to select the Houses layer.
  5. For ‘Input rasters’, use the drop-down box to select all five interpolation rasters you created earlier one at a time: IDW_Rain, Spline_Rain, NN_Rain, TINRast_Rain and TINRast_Wind.

The order in which the rasters are listed is the same order in which the extracted value columns will appear in the output feature class attribute table.

  1. Check Bilinear interpolation of values at point locations to further interpolate the value, if the house point does not lie exactly at the center of a raster cell.
  2.  Ensure that your 'Extract Multi Values to Points' tool Parameters tab appears as shown below and click Run.

    Image Added

The Extract Multi Values to Points tool simply adds a field for the value of each raster to the end of your existing Houses attribute table, which is why you were not asked for an output location.

  1. In the Contents pane, right-click the Houses layer and select Attribute Table.

    Image Added

Notice the five new fields added to the end of the existing attribute table which contain the values of each of the five rasters at each house location. Notice also the differences in rainfall estimates resulting from the four different interpolation techniques. At many house locations, the value of the spline interpolation is the most different from the other interpolation method values.

You are now ready to export the data table into a format that can be used by the insurance company, without having to use ArcGIS software.

  1. At the top of the Geoprocessing pane, click the Back arrow button.
  2. Click the Conversion Tools toolbox > Excel toolset > Table To Excel tool.
  3. For ‘Input Table’, use the drop-down menu to select the Houses attribute table.

The output table location defaults to your project folder.

  1. For ‘Output Excel File’, rename the Excel file from Houses_TableToExcel to “Houses_Rain_Wind”.
  2. Check Use field alias as column header.
  3.  Ensure that your 'Table To Excel' tool Parameters tab appears as shown below and click Run.

...

The Extract Multi Values to Points tool simply adds a field for the value of each raster to the end of your existing Houses attribute table, which is why you were not asked for an output location.

...

Notice the five new fields added to the end of the existing attribute table which contain the values of each of the five rasters at each house location. Notice also the differences in rainfall estimates resulting from the four different interpolation techniques. At many house locations, the value of the Spline interpolation is the most different from the other interpolation method values.

You are now ready to export the data table into a format that can be used by the insurance company without having to use ArcGIS software.

  1. Click Conversion Tools > Excel > Table To Excel.
  2. Use the ‘Input Table’ drop-down box to select Houses.

The output table location defaults to your project folder.

  1. For ‘Output Excel File’, replace “Houses_TableToExcel” with “Houses_Rain_Wind”.
  2.  Check Use field alias as column header.
  3. Ensure your ‘Table To Excel’ pane looks like that shown below and click Run.

Now you have what the insurance company is looking for: a single text file containing the locations of the houses in question, as well as the corresponding rainfall and wind speed estimates to use in completing a damage quote for your customers.

...

Note that if you simply wanted a table containing the values of multiple rasters at each point and did not care about symbolizing or working with the points any further in the ArcGIS Pro environment, you could have used the Sample tool, which outputs a separate table, rather than adding the values to the end of your existing feature class attribute table.

...

For this example, you have learned that the rainfall data collected from the rain gauges had a common discrepancy. Each gauge misrepresented the rainfall amount by two inches uniformly across the state. There are several ways you could approach this problem, but you will add two inches uniformly across the interpolated rain raster to practice using the raster calculator.

...

interpolated rain raster to practice using the raster calculator.

  1. At the top of the Geoprocessing pane, click the Back arrow button.
  2. Within the Spatial Analyst Tools toolbox, click the Map Algebra toolset > Raster Calculator tool.
  3. Under ’Rasters’, double-click the IDW_Rain layer.
  4. Under ‘Tools’, double-click the + operator.
  5. Type "2".
  6. For ‘Output raster’, rename the exported raster from idw_rain_ras to “IDW_Rain_2”.
  7.  Ensure that your 'Raster Calculator' tool Parameters tab appears as shown below and click Run.

    Image Added

While the original IDW_Rain layer and the new IDW_Rain_2 layer will look the same visually, notice in the Contents pane, next to the continuous color ramps for the two layers, that the max and min values for the IDW_Rain_2 layer are now two higher than the max and min values for the original IDW_Rain layer.

...

You have now learned multiple methods of interpolating point and line data into raster data and extracting that raster data at corresponding point locations. You have also reviewed skills from other short courses, including mapping Mapping Locations Using XY dataCoordinates.