Difference between revisions of "Great Lakes Bathymetry"

From TSAS Library
Jump to: navigation, search
Line 48: Line 48:
* [http://cdm16021.contentdm.oclc.org/cdm/search/collection/p16021coll10/searchterm/great%20lakes/order/nosort Historical maps from the Army Corps of Engineers]
* [http://cdm16021.contentdm.oclc.org/cdm/search/collection/p16021coll10/searchterm/great%20lakes/order/nosort Historical maps from the Army Corps of Engineers]
* [http://www.glos.us/ Portal for Great Lakes Observing System information] from the Integrated Ocean Observing system
* [http://www.glos.us/ Portal for Great Lakes Observing System information] from the Integrated Ocean Observing system
* [https://nationalmap.gov/ The National Map] combines many of the above resources


==References==
==References==
<references />
<references />

Revision as of 08:01, 10 June 2017

Greatlakes.jpg

Ideas

Process of creating a 3D map of the area

Resources

CNC Carve

Workflow

  1. Download shapefiles from NOAA. Superior bathymetry contours with unclear provenance is available here, with login
  2. Download and install QGIS and Inkscape[1]
  3. No need to extract files, open each in QGIS Desktop
  4. Add OpenStreetMap information if desired (the native support is not working at this time[2], but the plugin does)
  5. Right-click shapefile in left layers panel, select Properties
  6. Choose Filter, then create expression, e.g. "DEPTH" IN ('0', '1', '25', '50', '60')[3]
It is helpful to use the Values box (make sure "Use unfiltered layer" is checked if changing an applied filter) to see the extent of the map's values

Workflow Take 2

3D Print

  1. Download ArcGrid DEM from USGS National Map or Earth Explorer (3DEP, 1/3 arc-second DEM)
  2. Unzip file
  3. Open QGIS, add Vector, add layer file (.adf)
  4. Select extent (Raster -> Extraction -> Clipper)
  5. Save as GeoTIFF
  6. Add GeoTIFF to model
  7. Analyze to DEM (Raster -> Analysis -> DEM [Terrain Models])
    This step creates an apparently inaccurate model (top). Using 3DEM creates a more accurate representation, but the extent was smaller on the second.
    Erie2bumps.png
    Eriesmooth.png
    • If necessary, merge the DEMs
    • Solution to bathymetric (negative) altitude values[5]: using Raster Calculator, add bottom value to all height values
  8. Create STL (Raster -> DEMto3D)[6]

A note on resolution.

Other Sources

References

  1. Stack Exchange answer here
  2. "OpenStreetMap in QGIS." OpenStreetMap Wiki. here.
  3. Geographic Information Systems Stack Exchange. "QGIS choose contour line labels." Here
  4. Geographic Information Systems Stack Exchange. "Load .osm.pbf file in Qgis." here.
  5. Geographic Information Systems Stack Exchange. "Mars DTM has all negative elevations in QGIS. Can't build STL file for 3D print." here.
  6. Modified from Laguna Tools. "Max CNC Tutorial 3D Landscape." YouTube. here.