Skip to main content

Working with Rasters

 In this last lab assignment for GIS Programming, I learned how to work with rasters using the arcpy.sa (Spatial Analysis) module. Our goal was to create a single raster image that met the following goals through the original dataset: 

  • Forest landcover (classifications 41, 42, and 43) 
  • Slope between 5° and 20° •
  • Aspect between 150° and 270°
To achieve these goals, I created five temporary rasters files and then combined them into one raster layer. I came across a couple of issues while creating my script for this assignment. 

The first issue I came across was I had 1,2,3 as the second number component in the RemapValue statement. I had missed the part in the script template that stated to assign all the new values to 1.

My second major issue was also caused by skipping a line. I did not add the line to create the elevation raster. This gave me some very interesting issues. While the script ran through all the way, whenever I tried to map my final raster file I got 12 categories rather than just 2 or one single elevation sliding scale. Once I added in the elevation raster and associated it with the Slope and Aspect, my issues with mapping the data went away.

Below is a picture of the final output of my script, as well as the mapped final data. 

Image 1: Script final output 
Image 2: Mapped final data  


Module 6 Script FlowChart



Labels:

Comments

Post a Comment

Popular posts from this blog

Utilizing ERDAS Imagine to Analyze Map Features

 This week we learned how to utilize histograms and different bands to highlight different features in a map. On the following map that we worked on, dark bodies of water caused high peaks on the left of histograms while snow-peaked mountains were small blips on the far right. These simple distinctions help to quickly identify map features on a graph, that you can then utilize as a stepping stone to finding them on the image. I found it incredibly interesting how the different band layers highlighted different features on the map. Figure 1 below depicts three different features we found on the image.  Figure 1: Distinct features found on an image using ERDAS Imagine. Feature 1: Large body of water. Feature 2: Snow-capped mountains transitioning to thick vegetation. Feature 3: Shallow turbulent body of water near urbanized land, transitioning to deep calm body of water. 
Post a Comment
Read more

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

Scale Effect and Spatial Data Aggregation

 In this lab we first looked at how the scale at which we analyze data can affect geometric properties. Looking at data at three different scales (1:1200, 1:24000, 1:100000) I was surprised by my final calculations. I thought that there would be a lineal proportional relationship along the scales, but I did not find this. It makes sense that the greatest resolution [1:1200] would have the most detail of geometric properties, but I was surprised that 1:100000 had greater geometric properties than the medium resolution. I resampled the data using the Bilinear technique and found the following effect on DEM resolution.  As DEM resolution increases the average slope in degrees decreases. This makes sense because as you get closer, more in detail by getting closer, you are zooming in as actively seeing less features altogether.  In the last part of the lab, we looked into the Gerrymandering of districts in the USA. Gerrymandering essentially breaks up congressional distri...
Post a Comment
Read more