Ranch Model for Tractor Treatment Areas

J W Ranch Study Area

Using ArcGIS 9.1 Model Builder

October 3, 2005

 

This model uses 3 input layers:

The model starts with the Slopes layer and then reclassifies and exports a layer with target slopes < 15°.  The reclassified layer is converted into features (Raster to Vector) and then it is clipped by the Rancher’s boundaries.  Next, the Soils layer is queried for good soil types and again clipped by the Rancher’s boundaries.  Finally, the good slopes (less than 15°) and the good soils (non-rocks) are intersected to produce a solution layer.

 

Model Flow Chart

 

Step 1 - Reclassify the Good Slopes

This procedure reclassifies slopes less than 15° as good (or having data) and all other slopes as null or having no-data. The input and reclassified layers cover the entire county and while processing this much data is slow, it is still faster than finding, mosaicking and processing individual quads.

 

Python Script:

Executing (Reclassify): Reclassify slopesfromdem Value "0 15 1;15 99 NODATA" D:\arcgis\1.dataCD\input_layers\reclass_slop1 DATA

Start Time: Mon Oct 03 18:21:39 2005

Executed (Reclassify) successfully.

End Time: Mon Oct 03 18:22:36 2005

(Elapsed Time: 57.00 seconds)

 

Step 2 – Convert Raster to Features

This step converts the reclassified raster file into a feature file.  Vectors or distinct features are necessary for the upcoming clip and intersect geoprocessing.

 

Python Script:

Executing (Raster to Polygon): RasterToPolygon D:\arcgis\1.dataCD\input_layers\reclass_slop1 D:\arcgis\1.dataCD\input_layers\RasterT_Reclass1.shp SIMPLIFY VALUE

Start Time: Mon Oct 03 18:22:37 2005

Executed (Raster to Polygon) successfully.

End Time: Mon Oct 03 18:23:30 2005

(Elapsed Time: 53.00 seconds)

 

Step 3 – Clip Slopes with Local Boundaries

This step clips the good slopes layer with the Rancher’s boundaries.  Clipping the countywide slopes with the study-area helps minimize the data required for future geoprocessing (Intersect routine in step 6).

 

Python Script:

Executing (Clip): Clip D:\arcgis\1.dataCD\input_layers\rastert_reclass1.shp bndry_jwranches D:\arcgis\1.dataCD\input_layers\RasterT_Reclass1_Clip.shp #

Start Time: Mon Oct 03 18:23:30 2005

Reading Features...

Cracking Features...

Assembling Features...

Executed (Clip) successfully.

End Time: Mon Oct 03 18:23:52 2005

(Elapsed Time: 22.00 seconds)

 

Step 4 – Query Non-Rocky Soils

This step uses a non-spatial SQL query to select good soil types from the NRCS soil dataset.  There are over 130 different soil types, so I’ve just eliminated soils with descriptions including rock, gravel, complex, water, etc.  Loam seems to be the surviving soil type in this selection.

 

 

Python Script:

Executing (Select): Select soilmu_nrcs D:\arcgis\1.dataCD\analysis\soilmu_nrcs_Select.shp ""musym" <> 'DBF' AND "musym" <> 'DCF' AND "musym" <> 'DDF' AND "musym" <> 'DFC' AND "musym" <> 'FAE' AND "musym" <> 'FBC' AND "musym" <> 'GBE' AND "musym" <> 'JAF' AND "musym" <> 'SFD' AND "musym" <> 'SHF' AND "musym" <> 'SKF' AND "musym" <> 'W' AND "musym" <> 'WC'"

Start Time: Mon Oct 03 18:23:52 2005

Executed (Select) successfully.

End Time: Mon Oct 03 18:23:57 2005

(Elapsed Time: 5.00 seconds)

 

Step 5 – Clip Soils with Local Boundaries

This step clips the good soils layer with the Rancher’s boundaries.

 

Python Script:

Executing (Clip (2)): Clip D:\arcgis\1.dataCD\analysis\soilmu_nrcs_Select.shp bndry_jwranches D:\arcgis\1.dataCD\analysis\soilmu_nrcs_Select_Clip.shp #

Start Time: Mon Oct 03 18:23:57 2005

Reading Features...

Cracking Features...

Assembling Features...

Executed (Clip (2)) successfully.

End Time: Mon Oct 03 18:23:58 2005

(Elapsed Time: 1.00 seconds)

 

Step 6 – Intersect the Good Slopes and Good Soils

This step intersects the good slopes with the good soils to produce solution areas for this suitability model. The intersection of soils and slopes isolates areas that are both less than 15° in transition and non-rocky soil types.

 

Python Script:

Executing (Intersect): Intersect "D:\arcgis\1.dataCD\input_layers\RasterT_Reclass1_Clip.shp '';D:\arcgis\1.dataCD\analysis\soilmu_nrcs_Select_Clip.shp ''" D:\arcgis\1.dataCD\analysis\Intersec_Slopes_Soils.shp ALL # INPUT

Start Time: Mon Oct 03 18:23:58 2005

Reading Features...

Cracking Features...

Assembling Features...

Executed (Intersect) successfully.

End Time: Mon Oct 03 18:23:59 2005

(Elapsed Time: 1.00 seconds)

 

 

Total Script Processing Time:    139 Seconds (approx. 2 ½ minutes)

 

Local Environment

CPU:               Pentium IV 2.5ghz

Memory:          512 mb

OS:                  XP

ArcGIS version:     9.1

Arc Extensions:     Spatial Analyst

 

Resulting Map & Data Layer: