Wednesday, June 13, 2012

Natural Hazards - Tsunamis

This week's assigment looked at how GIS can be used in regards to Tsunamis.

Here is my map and process summary for this week:

Process Summary

Part I: Building a File Geodatabase
1. Copied the data from the R: drive and unzipped
2. Started ArcCatalog
3. Documented Metadata
4. Navigated to my S:\GIS4048-ApplicationsInGIS\W4-NaturalHazards-Tsunamis\Tsunami directory
5. Previewed the Roads shape file to make sure the Preview Tab works properly.
6. Created the Tsunami.gdb File Geodatabase
7. Created the Transportation Feature Dataset and set the Project Coordinate tp WGS1984
8. Added the Rails_UTM and the Roads_UTM shapefiles to the File Geodatabase
9. Started ArcMap
10. Now we need to convert the data in Excel into a spatial feature
11. First I added the Excel Sheet as a table
12. Next I right-Clicked it and chose “Display XY Data”
13. I chose to import the Coordinate system and chose my File Geodatabase and the Feature Dataset as the place to import from.
14. Added the Administrative and Damage_Assesment shapefiles and added the recommended shapefiles to each.
15. Made sure my File Geodatabase matched the lab instructions
16. Opened ArcToolbox.
17. Ran the Build Raster Attribute Table tool on each of the DEMs
18. Created the FukuDEM and SendDEM Raster Datasets
19. Load Data process completed
20. Calculated Statistics on both DEMs.  It is interesting to try to preview the data before running the statistics.   It appears solid black.   Post processing it appears normal.
21. Took a screenshot and prepared it as my first deliverable

Part II: Fukushima Radiation Exposure Zones
1. Started ArcMap
2. Opened Tsunami.mxd
3. Exported the Fukushima II Power Plant to its own Feature Class
4. Answered Q1
5. Created my Multiple Ring buffer.   Pretty cool tool !
6. Adjusted my Symbology
7. Clipped the Roads feature
8. Did a Select by Location query to answer Q2
9. Exported these cities out to Evac_zone_cities
10. Posted my screenshot for Q3
11. Answered questions 4 through 7 using that same tool and the statistics tool for the sum of population
12. Worked on my map deliverable

Thursday, June 7, 2012

SWOPPAT Announcements - New Tool

It has come to our attention that there has been a great demand for a tool that will allow SWOPPAT staff to add two additional columns to their shapefiles.   Acres and HectaAcres.    We at SWOPPAT GIS Programming have made this possible with the new SWOPPATAAHACT "The Save the World One Python Program at a Time Acre and HectaAcre Calculation Tool" (™ and © pending).

This tool can be accessed from your toolbar as follows:



This tool can save you countless hours.   Lets say for example that you create 3 shapefiles per week.  And lets say that calculating the Acres and HectaAcres tool you 3 minutes each time.    That would be 9 minutes a week.   When you look at this over the course of the year, you see that this tool can save you 468 minutes a year.    That is over 7 and a half hours a year, almost an entire work day !


Natural Hazards - Hurricanes

This week's assignment asked us to look at the damage caused by Hurricane Katrina in three counties in Mississippi.  There are 3 maps and a graph and a table that I have submitted here along with a Process Summary.
Deliverable #1- A map showing elevation/bathymetry of the Mississippi Coast counties with places, types of water, barrier islands, and hydrography.
Deliverable #2 - A map of flooded land of the Mississippi Coast after Hurricane Katrina.
Deliverable #3 - A map showing infrastructure and health facilities at risk from storm surge.

Deliverable #4 - A bar graph showing percentage of total flooded land by land-cover type.

Deliverable #5 - A table showing various land types that were flooded measured in acres and square miles.
Process Summary

1.      Set up my directory on S: Drive and copied over the necessary data from the R: drive
2.      Answered Questions 1-4 (Returned to them at several points during the lab)
3.      Opened ArcCatalog and previewed the data in both spatial and tabular form
4.      Examined the Metadata to complete the chart as part of Q5
5.      Organized my directory structure so that I have a S:\GIS4048-ApplicationsInGIS\W3-NaturalHazards-Hurricanes\Hurricanes\Project1\results path for my results
6.      Created a new map document in ArcGIS and called it coast1
7.      Created my Map Document Properties and set to relative paths.
8. Saved work
9.      Set my coordinate system by importing the coordinate system from the land cover feature
10.  Set up my Environment in Toolbox to handle current and scratch workspaces, Output Coordinate System, Raster Analysis settings and the mask to base the results on the 3 counties of interest
11.  Added the elevation raster.
12.  Answered Q6
13.  Symbolized the color ramp for the Elevation.   The tip about right clicking to see the names will come in handy at work and at my internship.
14.  Added the Counties and set up the labels
15.  Answered Q7
16.  Added the Places, symbolized, and answered Q8
17.  Added Islands to the map
18.  Added swamps and rivers and symbolized appropriately.
19.  This will be the base map that will give us a starting point for the rest of the assignment
20.  Answered Q9 and 10
21.  Saved my map (will com back to it later to get it ready to be a deliverable, as it is currently very busy)
22.  Created Coast2.mxd
23.  Documented the map, set the data frame properties, and set the environments just like the previous map
24.  Turned on Spatial Analyst
25.  Used the Math Tools in Spatial Analyst to convert meters to feet
26.  Used the logical toolset to get the lands where elevation is less than or equal to 15 feet and saved it as flooded_land
27.  Added Landcover
28.  Reclassified Landcover using the Spatial Analyst reclassify tool.
29.  Used Editor to label the reclassified data
30.  Saved the reclassified data as a layer file
31.  Added the counties and labeled them
32.  Obtained the Statistics on the Count field
33.  Answered Q12
34.  Created a new field called “Percent” as float
35.  Used the field calculator to set this new column to equal COUNT/1380797*100
36.  Verified that this totaled to 100% using the statistics tool.
37.  Saved my coast 2 document and graph.   I will come back to prep it for delivery later
38.  Created my coast3.mxd and set the description, coordinate system information, and environment
39.  Added my infrastructure data to the MXD and symbolized appropriately.
40.  Answered Q15 and Q16 ( I found a cool Spatial Analyst tool for this step called Spatial Analyst Tools | Extraction | "Extract Values to Points"
41.  Started my coast4.mxd document
42.  Documented the map, set the coordinate info, and set the environment settings.
43.  Prepared my deliverables and Blog Post Entries.

Wednesday, May 30, 2012

Natural Hazards / Earthquakes

This weeks' assignment really was a great way to shake off the cobwebs from the break and start making some maps again.   We concentrated on using some tools in ArcGIS that I hadn't used before and that will prove VERY helpful in my current job.

For part 1, we considered the impact on infrastructure in the event of an earthquake near Memphis, Tennessee.   On this map, I decided to make it 11" x 17" rather than 8.5" x 11" to try to best display all three scenarios (dams, Interstates, and railroads).

For Part 2, we looked at the peak ground velocity of the aftershocks.  

For Part 3, we examined aftershocks that were greater than magnitude 3.0

And for part 4, we looked at the intensity of the aftershocks.


We also created a couple of graphs along the way using the built in ESRI tools:


Here is my Process summary for the lab this week:

Process Sumary:
 Earthquakes Part I
1)      Copied all the data from the r: drive to the s: drive.
2)      Unzipped the files
3)      Opened up the NewMadrid.mxd file
4)      Used the identify tool to research the data
5)      It appears that Memphis would suffer the most damage in an earthquake in this region
6)      Performed a select by location to see the areas that would be affected
7)      Next we learned about the summary feature which I wasn’t aware of.   I can see this being VERY useful in my day to day work.
8)      We can see that there are 91 areas that would be greatly impacted by this event
9)      Next we calculated the population density of the counties.  Again, VERY useful lab so far and I have just started !
10)  Next we intersected the county data (including population density) with the New Madrid data
11)  Created the PopMMI table
12)   Using the summarize tool, we can see that 60,088,857 people live in the affected areas
13)  Created the graph as instructed
14)  Added my numeric values to the Modified Mercalli Scale data
15)  Created my clipped Istaterisk layer
16)  Next I did a SELECT BY ATTRIBUTES to determine all the areas where the intensity was 8 or greater
17)  Then I did a SELECT BY location to see which Interstates intersected these areas
18)  Clipped this data to the map to see my At Risk Interstates
19)  Did a similar process for Railroads in the are that had an index of 10 or higher
20)  Finished this section by doing the dams at risk as well
21)  For my map on this section, I wanted to try to think outside the box a little bit so I choose a different map size of 11x17 instead of the normal 8.5 x 11.   I hope this was OK for this map which will display all 3 results on it.
Earthquakes Part II
1)      Now on to the second part of the lab
2)      Opened the Northridge1.mxd map
3)      Made sure the Spatial Analyst extension was enabled (Under Customize, rather than Tools)
4)      Symbolized the dots for the building damage and arranged layers as instructed.
5)      Did a select by attribute to choose only the red and yellow tagged buildings
6)      Set up the Kernel Density tool as instructed
7)      Cleared the selection and turned off unused layers
8)      Set the lowest level of symbolization to “No Color”
9)      Used the effects toolbar to set visibility to 15%
10)  Saved the Damage Pattern as a layer file
11)  Added the layer back in and symbolized
12)  Opened the Geology layer
13)  Examined the Liquefaction layer
14)  Tuned on the Stations layer
15)  Opened the Spline Interpolation tool
16)  Ran this tool with the recommended settings
17)  Created the PGV Layer and set symbols/transparency/etc…
18)  Prepared Map
Earthquakes Part III
1)      Now on to the third part of the lab
2)      Opened up the Northridge2.mxd map
3)      Added the text file while changing the Projection to WGS
4)      Exported this data out to a shapefile
5)      Removed the CSV file from the project
6)      Selected the Main Shock by doing a select by attribute
7)      Created a new layer to represent the Main shock and symbolized as a red triangle
8)      Created my map of Aftershocks greater than 3 and symbolized appropriately.   I am a little worried that my dots are in a “grid” format rather than randomly scattered, but hopefully that is how it is supposed to look.   I actually reloaded the project from the R: drive to ensure that this should be fine.
9)      I didn’t see any of the ArcScene in this lab which it referred to
Earthquakes Part IV
1)      Now on to the fourth part of the lab
2)      Opened up the Northridge3.mxd map
3)      Created the Summary table similar to the process we did earlier in the lab
4)      Created my Graph and exported it out as a JPG file
5)      Posted my contents to dropbox and blog
I really enjoyed this weeks’ work even though it was very time consuming.   I learned a lot of new tips with ArcGIS that I can start applying to by job already.    I hope this is a sign of what kind of labs to expect for this class, because it really is helpful.   I tried to make my maps a little different from each other just to try to get back into the swing of making maps each week.




Sunday, May 20, 2012

GIS Programming Orientation


The programming language I am really looking forward to learning this semester is Python.   Python is an open-source project and it is freely usable and distributable, even for commercial use.   It strives to achieve “readability” which means it is relatively easy to read someone else’s code to determine the functionality.  It is "portable" which means that the code can be used on several platforms with very little modification.    It is "versatile" because it can be used as a scripting language but also supports fully object oriented programming as well.   ESRI has leveraged the strength of Python into its product line by implementing the ArcPy package which can “perform geographic data analysis, data conversion, data management, and map automation1    Python is the successor to VBA (Visual Basic for Applications) as the primary development language within the ESRI product line.   ModelBuilder can be used to generate Python code from a GUI (Graphic User Interface) view.  
An example of a program that could be generated using Python would be a drive time estimator.   You could click at any point on a map and Python could then compute the drive time to a particular location.  Python would leverage the geographic information that ArcGIS can provide to retrieve these results.

It is a great time to learn Python with the release of 10.1 coming out next month. 

Some pages I found to help fellow Argonauts include:


 Sources:
1 –“A quick tour of ArcPy” by ESRI Help http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//000v00000001000000

Saturday, May 19, 2012

Why I am looking forward to APPLICATIONS IN GIS....

The topic I am most interested in studying is “GIS for Location-Based Services”.     According to the ESRI Web site: “With these servers and services, you can build new Web and wireless applications faster or enhance existing ones with location, addresses, points of interest, dynamic maps, and routing directions.”    In the last few years we have seen a surge in the number of companies that take advantage of dynamic location based applications.   Augmented reality will continue this evolution of providing spatial information to the mobile user based on location and context.  Consider this concept video for “Google Glasses” : Project Glass     This concept shows some of the benefits that could be seen by using a Internet-Connected Heads-Up Display (HUD).  Spatial information such as current location, routing directions (both outside and inside of buildings), proximity messages to both people and locations, and location based services such as weather could be made available to a user.   Application Developers with GIS knowhow will be required to make these bold concepts appear in actual commercial projects.    That is why I am interested in both my classes this semester: APPLICATIONS IN GIS and GIS PROGRAMMING

Thursday, May 3, 2012

GIS 3015 Final - Cartographic Skills

Here is my submission for the Final Project for GIS3015.

This is a map that was created to show the participation rate by state as well as the mean scores by state for the ACT standardized test available to High School seniors nationwide.