Skip to main content

Geoprocessing

 This week we learned geoprocessing in ArcGIS Pro via the Modelbuilder and Python. 

First, I built a model in the ModelBuilder. My goal: 

  1.  Clips all soils to the extent of the basin 
  2.  Selects all soils that are classified as "Not prime farmland" under the [FARMLNDCL] attribute 
  3. Erases the "Not prime farmland" soil selection from the basin polygon 


First, I utilized the Clip tool to remove the Soil areas (input feature) from the Basin area (clip feature). [I think of the layers of the cookie dough sheet (clip feature) vs the cookie cutter (input feature).]

Second, I used the clipped data from step one to Select only the areas classified as ‘Not prime farmland,” utilizing the query Where FARMLNDCL is equal too ‘Not prime farmland.’ Here the attribute ‘Not prime farmland’ comes from the clipped clipped soil layer.

Lastly, I used the erased the selected data (Erase feature) from step 2 from the original soil shapefile (input feature).

The second part of the lab, I used python. Our desired outcomes were: 

  1. Adds XY coordinates to the hospitals shapefile 
  2. Creates 1000m buffer around the hospitals 
  3. Dissolves the hospital buffers into a single-feature layer 
  4. Prints geoprocessing messages after each tool runs 

To begin and add the XY coordinates, I based the script completely off of the ESRI help link. 

arcpy.management.AddXY(in_features)

I know I had to add an overwrite component and used ‘arcpy.overwriteOutput = True’ line based on our practice exercise. To see the desired history method for each script I used the print (arcpy.GetMessages()) syntax I learned in the exercise as well. I also manually created a custom print statement, stating what tool was running before the GetMessages statement. 

Also going off of the practice exercise, I used the following for the buffer script:

Arcpy. Buffer_managment (‘file’ , ‘file path’  , ‘buffer distance’)

The last step was the trickiest for me, dissolving the buffered file as I got a bit lost in the syntax.

arcpy.management.Dissolve(in_features, out_feature_class, {dissolve_field}, {statistics_fields}, {multi_part}, {unsplit_lines})

I used the above script, utilizing BUFF_DIST as the dissolve field, leaving the stat field blank (“”), using single rather than multi_part, and lastly ‘DISSOLVE_lines’. 

n

Labels:

Comments

Post a Comment

Popular posts from this blog

Positional Accuracy: NSSDA

 In this analysis, I compared the street and road intersect data collected for Alburquerque, NM by the City of Alburquerque and the application StreetMaps. I used an orthophoto base layer as the reference for this analysis, to compare and determine the accuracy of both the City and Streetmap layers using NSSDA procedures. The most difficult part of this analysis for me was how to determine what 20% per quadrant looks like. Because the reference map was divided into 208 quadrants, I had to determine how to subdivide all the quadrant's equality into 20%. After multiple trials and error, I decided to subdivide the entire area (208 sub-quadrants) into 4 equal-area subsections. In this way, I could do 5 random right intersection points per subsection or 20% per subsection.  Map 1: City of Albuquerque city map data.  Map 2: City of Alburquerque SteetMap data When selecting a random intersection to place the points within each quadrant, I choose a location that had data f...
Post a Comment
Read more

Isarithmic Mapping

  Map 1: Annual Precipitation, Washington State Map 1 is an Isarithmic map that follows the continuous phenomenon of rainfall in Washington state over a 30 year period. The data was created by the PRISM group at the Oregon State University in 2006, and then downloaded and amended by the U.S. Department of Agriculture, Natural Resources Conservation Service, National Geospatial Management Center in 2012. Eden Santiago Gomez, analyzed the data on 5/2/2021, to create the map above. Santiago Gomez created continuous tones for the data, also adding a hillshade effect. She then converted the floating raster data into Integer data via the geoprocessing tool Int (Spatial Analyst Tool) to bring out hypsometric tinting. Lastly, she added contours of the data via the Contour List tool.   How the precipitation data was derived and interpolated? The PRISM system has been continually developed over the past couple decades, utilizing physiographical maps and climate fingerprints as its ...
Post a Comment
Read more

Choropleth and Dot Mapping

 This week we explored choropleth and dot mapping. Choropleth is a thematic form of mapping that focuses on color units, whose color intensity is proportional to its corresponding data value. Dot mapping is also thematic. It uses either a proportional or graduated thematic symbol (like a circle), whose size increases due to its data value. Using ArcGIS pro, I analyzed the population densities of countries in Europe (person per square kilometer), as well as their wine consumption (liters per capita) to determine if there was a correlation between the two. In my choropleth map, I decided to use a natural breaks classification. I chose not to use Equal Interval because only 2 classes (with slight 3 rd class) were represented in the map, and it looked like almost just one color in the lower range. The standard deviation classification appeared to be more diverse at first glance but was actually skewed to the top ranges. I was then between Quantile and Natural Breaks. While both t...
Post a Comment
Read more