Difference between revisions of "Carneros Watershed: Map of Soil Erodibility Index"

From CCoWS Wiki
Jump to: navigation, search
 
(2 intermediate revisions by one user not shown)
Line 7: Line 7:
  
 
* Original source of data: SSURGO
 
* Original source of data: SSURGO
* Data analysis: [[Fred Watson and Kathy Pugh]]
+
* Data analysis: Kathy Pugh and [[Fred Watson]]
 
* Map produced by: Kathy Pugh
 
* Map produced by: Kathy Pugh
 
* Map produced using: [[ArcMap]]
 
* Map produced using: [[ArcMap]]
Line 15: Line 15:
 
== Data Analysis Steps Taken ==
 
== Data Analysis Steps Taken ==
  
#1. Download SSURGO data  
+
# Download SSURGO data  
#2. Receive zip file of zip files, each one like soil_ca069.zip.
+
# Receive zip file of zip files, each one like soil_ca069.zip.
#3. Extract, and open basic shapefile of map-units (MU): soilmu_a_ca069.shp in Arc
+
# Extract, and open basic shapefile of map-units (MU): soilmu_a_ca069.shp in Arc
#4. Can't read tables in this zip, because need special MS Access Template
+
# Can't read tables in this zip, because need special MS Access Template
#5. Download Access template from  http://soildatamart.nrcs.usda.gov/templates.aspx
+
# Download Access template from  http://soildatamart.nrcs.usda.gov/templates.aspx
#6. Open template in Access; Enable Macros; and it should prompt for location of your tables
+
# Open template in Access; Enable Macros; and it should prompt for location of your tables
#7. Wait a few minutes while Access imports your tables
+
# Wait a few minutes while Access imports your tables
#8. Create some custom queries in Access to extract some kind of Kf for each map unit
+
# Create some custom queries in Access to extract some kind of Kf for each map unit
##8a. Review the SSURGO data model: http://soildatamart.nrcs.usda.gov/documents/SSURGO_2_2_DataModel_Part_1.pdf
+
## Review the SSURGO data model: http://soildatamart.nrcs.usda.gov/documents/SSURGO_2_2_DataModel_Part_1.pdf
##8b. Each MapUnit contains several Components, each containing several Horizons.
+
## Each MapUnit contains several Components, each containing several Horizons.
##8c. First query selects the surface horizons:
+
## First query selects the surface horizons:
##8c1 - it joins tables muaggatt+component+chorizon,
+
## - it joins tables muaggatt+component+chorizon,
##8c2 - joining on keys mukey and cokey
+
## - joining on keys mukey and cokey
##8c3 - selecting only horizons whose hzdept_r is 0 (i.e. surface horizons)
+
## - selecting only horizons whose hzdept_r is 0 (i.e. surface horizons)
##8c4 - includes the kffact (Kf-factor) and comppct_r (Percentage of map-unit occupied by Component) fields that will be used below
+
## - includes the kffact (Kf-factor) and comppct_r (Percentage of map-unit occupied by Component) fields that will be used below
##8d. Second query takes first as input, and is a MakeTable query, making the table that will be used in Arc
+
## Second query takes first as input, and is a MakeTable query, making the table that will be used in Arc
##8d1 - computes an expression 'yournamekf: Sum([comppct_r]*[kffact]/100)'
+
## - computes an expression 'yournamekf: Sum([comppct_r]*[kffact]/100)'
##8d2 - it is a Summation query, and the Summation is done on the fredkf field (above)
+
## - it is a Summation query, and the Summation is done on the yournamekf field (above)
##8d3 - the other two fields are musym and mukey, and they are GroupBy fields
+
## - the other two fields are musym and mukey, and they are GroupBy fields
##8e. RUn the second query; it will make a table
+
## Run the second query; it will make a table
##8f. Export the table to DBF format (one way of getting it to Arc)
+
## Export the table to DBF format (one way of getting it to Arc)
#9. Back in Arc, add the DBF table as a layer
+
# Back in Arc, add the DBF table as a layer
#10. Join the soilmu_a_ca069.shp layer to the DBF layer, using the mukey
+
# Join the soilmu_a_ca069.shp layer to the DBF layer, using the mukey
#11. Change the soilmu_a_ca069.shp symbology to display quantitative based on the yournamekf field
+
# Change the soilmu_a_ca069.shp symbology to display quantitative based on the yournamekf field
  
 
== Alternate versions ==
 
== Alternate versions ==

Latest revision as of 07:44, 28 February 2011

Carneros Watershed: Soil Erodibility Index Map

Metadata

  • Original source of data: SSURGO
  • Data analysis: Kathy Pugh and Fred Watson
  • Map produced by: Kathy Pugh
  • Map produced using: ArcMap
  • Coordinate system: Geographic Coordinate System - GCS_North_American_1983
  • Datum: D_North_American_1983

Data Analysis Steps Taken

  1. Download SSURGO data
  2. Receive zip file of zip files, each one like soil_ca069.zip.
  3. Extract, and open basic shapefile of map-units (MU): soilmu_a_ca069.shp in Arc
  4. Can't read tables in this zip, because need special MS Access Template
  5. Download Access template from http://soildatamart.nrcs.usda.gov/templates.aspx
  6. Open template in Access; Enable Macros; and it should prompt for location of your tables
  7. Wait a few minutes while Access imports your tables
  8. Create some custom queries in Access to extract some kind of Kf for each map unit
    1. Review the SSURGO data model: http://soildatamart.nrcs.usda.gov/documents/SSURGO_2_2_DataModel_Part_1.pdf
    2. Each MapUnit contains several Components, each containing several Horizons.
    3. First query selects the surface horizons:
    4. - it joins tables muaggatt+component+chorizon,
    5. - joining on keys mukey and cokey
    6. - selecting only horizons whose hzdept_r is 0 (i.e. surface horizons)
    7. - includes the kffact (Kf-factor) and comppct_r (Percentage of map-unit occupied by Component) fields that will be used below
    8. Second query takes first as input, and is a MakeTable query, making the table that will be used in Arc
    9. - computes an expression 'yournamekf: Sum([comppct_r]*[kffact]/100)'
    10. - it is a Summation query, and the Summation is done on the yournamekf field (above)
    11. - the other two fields are musym and mukey, and they are GroupBy fields
    12. Run the second query; it will make a table
    13. Export the table to DBF format (one way of getting it to Arc)
  9. Back in Arc, add the DBF table as a layer
  10. Join the soilmu_a_ca069.shp layer to the DBF layer, using the mukey
  11. Change the soilmu_a_ca069.shp symbology to display quantitative based on the yournamekf field

Alternate versions

  • None

Links