Difference between revisions of "Great Lakes Bathymetry"

From TSAS Library
Jump to: navigation, search
Line 11: Line 11:
* [http://www.naturalearthdata.com/downloads/ Natural Earth] provides shapefiles for cultural and physical elements
* [http://www.naturalearthdata.com/downloads/ Natural Earth] provides shapefiles for cultural and physical elements


==Laser Cut, Low Resolution==
==Create STL for CNC Router or 3D Printer==
===Workflow===
===Proprietary Workflow (Aspire)===
# Download shapefiles from [https://www.ngdc.noaa.gov/mgg/greatlakes/greatlakes.html NOAA]. Superior bathymetry contours with unclear provenance is available [https://databasin.org/datasets/88d305616c3e44cabbd465dde0d916c7 here], with login
Video walkthrough [https://www.youtube.com/watch?v=P3q5Fzb01ls here]
# Download and install [http://www.qgis.org/en/site/ QGIS] and [https://inkscape.org/en/download/ Inkscape]<ref>Stack Exchange answer [https://gis.stackexchange.com/questions/51980/how-to-convert-shapefile-to-vector-emf here]</ref>
# Data pulled from [https://maps.ngdc.noaa.gov/viewers/bathymetry/ NOAA's ETOPO1 images], available directly [https://www.ngdc.noaa.gov/mgg/global/global.html here] which offers a 1-minute resolution of worldwide bedrock
# No need to extract files, open each in QGIS Desktop
# Add raster image to QGIS
# Add OpenStreetMap information if desired (the native support is not working at this time<ref>"OpenStreetMap in QGIS." ''OpenStreetMap Wiki''. [http://wiki.openstreetmap.org/wiki/QGIS#OpenStreetMap_in_QGIS here].</ref>, but [https://plugins.qgis.org/plugins/OSMDownloader/ the plugin] does)
* Select extent desired (Raster -> Extraction -> Clipper)
# Right-click shapefile in left layers panel, select Properties
* Open in Aspire, follow directions from video
# Choose Filter, then create expression, e.g. "DEPTH" IN ('0', '1', '25', '50', '60')<ref>Geographic Information Systems Stack Exchange. "QGIS choose contour line labels." [https://gis.stackexchange.com/questions/121585/qgis-choose-contour-line-labels Here]</ref>
{{note|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}}


==Create STL for CNC Router or 3D Printer==
===Open Source Workflow===
===Open Source Workflow===
# Download ArcGrid DEM from USGS [https://viewer.nationalmap.gov/basic/#productSearch National Map] or [https://earthexplorer.usgs.gov/ Earth Explorer] (3DEP, 1/3 arc-second DEM)
# Download ArcGrid DEM from USGS [https://viewer.nationalmap.gov/basic/#productSearch National Map] or [https://earthexplorer.usgs.gov/ Earth Explorer] (3DEP, 1/3 arc-second DEM)
Line 29: Line 26:
# Save as GeoTIFF
# Save as GeoTIFF
# Add GeoTIFF to model
# Add GeoTIFF to model
# Analyze to DEM (Raster -> Analysis -> DEM [Terrain Models]){{note|A balance must be made between extent and resolution. 1:50,000 creates excellent models, but cannot contain an entire lake. Above that and the file crashes my system. A further [https://blogs.esri.com/esri/arcgis/2010/12/12/on-map-scale-and-raster-resolution/ note] on resolution.[[File:Erie2bumps.png|thumb|left]][[File:Eriesmooth.png|thumb|left]]}}
# Analyze to DEM (Raster -> Analysis -> DEM [Terrain Models])
#* A balance must be made between extent and resolution. 1:50,000 creates excellent models, but cannot contain an entire lake. Above that and the file crashes my system. A further [https://blogs.esri.com/esri/arcgis/2010/12/12/on-map-scale-and-raster-resolution/ note] on resolution.[[File:Erie2bumps.png|thumb|left]][[File:Eriesmooth.png|thumb|left]]
#* If necessary, [http://www.qgistutorials.com/en/docs/raster_mosaicing_and_clipping.html merge the DEMs]
#* If necessary, [http://www.qgistutorials.com/en/docs/raster_mosaicing_and_clipping.html merge the DEMs]
#* Solution to bathymetric (negative) altitude values<ref>Geographic Information Systems Stack Exchange. "Mars DTM has all negative elevations in QGIS. Can't build STL file for 3D print." [https://gis.stackexchange.com/questions/217203/mars-dtm-has-all-negative-elevations-in-qgis-cant-build-stl-file-for-3d-print here].</ref>: using Raster Calculator, add bottom value to all height values
#* Solution to bathymetric (negative) altitude values<ref>Geographic Information Systems Stack Exchange. "Mars DTM has all negative elevations in QGIS. Can't build STL file for 3D print." [https://gis.stackexchange.com/questions/217203/mars-dtm-has-all-negative-elevations-in-qgis-cant-build-stl-file-for-3d-print here].</ref>: using Raster Calculator, add bottom value to all height values
# Create STL (Raster -> DEMto3D)<ref>Modified from Laguna Tools. "Max CNC Tutorial 3D Landscape." ''YouTube''. [https://www.youtube.com/watch?v=Ef6f6dwQbCc here].</ref>
# Create STL (Raster -> DEMto3D)<ref>Modified from Laguna Tools. "Max CNC Tutorial 3D Landscape." ''YouTube''. [https://www.youtube.com/watch?v=Ef6f6dwQbCc here].</ref>


===Proprietary Workflow (Aspire)===
==Laser Cut, Low Resolution==
* Video walkthrough [https://www.youtube.com/watch?v=P3q5Fzb01ls here]
===Workflow===
* Imagery pulled from [https://maps.ngdc.noaa.gov/viewers/bathymetry/ NOAA's ETOPO1 images], available directly [https://www.ngdc.noaa.gov/mgg/global/global.html here] which offer a 1-minute resolution of worldwide bedrock
# Download shapefiles from [https://www.ngdc.noaa.gov/mgg/greatlakes/greatlakes.html NOAA]. Superior bathymetry contours with unclear provenance is available [https://databasin.org/datasets/88d305616c3e44cabbd465dde0d916c7 here], with login
* Select extent with QGIS as above
# Download and install [http://www.qgis.org/en/site/ QGIS] and [https://inkscape.org/en/download/ Inkscape]<ref>Stack Exchange answer [https://gis.stackexchange.com/questions/51980/how-to-convert-shapefile-to-vector-emf here]</ref>
* Open in Aspire, follow directions from video
# No need to extract files, open each in QGIS Desktop
# Add OpenStreetMap information if desired (the native support is not working at this time<ref>"OpenStreetMap in QGIS." ''OpenStreetMap Wiki''. [http://wiki.openstreetmap.org/wiki/QGIS#OpenStreetMap_in_QGIS here].</ref>, but [https://plugins.qgis.org/plugins/OSMDownloader/ the plugin] does)
# Right-click shapefile in left layers panel, select Properties
# Choose Filter, then create expression, e.g. "DEPTH" IN ('0', '1', '25', '50', '60')<ref>Geographic Information Systems Stack Exchange. "QGIS choose contour line labels." [https://gis.stackexchange.com/questions/121585/qgis-choose-contour-line-labels Here]</ref>
#* 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


==Other Sources==
==Other Sources==

Revision as of 13:30, 11 June 2017

Partial Erie bathymetry cut with MPCNC, map from NOAA, processed with QGIS and Aspire 8.5

Ideas

Resources

Create STL for CNC Router or 3D Printer

Proprietary Workflow (Aspire)

Video walkthrough here

  1. Data pulled from NOAA's ETOPO1 images, available directly here which offers a 1-minute resolution of worldwide bedrock
  2. Add raster image to QGIS
  • Select extent desired (Raster -> Extraction -> Clipper)
  • Open in Aspire, follow directions from video

Open Source Workflow

  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])
    • A balance must be made between extent and resolution. 1:50,000 creates excellent models, but cannot contain an entire lake. Above that and the file crashes my system. A further note on resolution.
      Erie2bumps.png
      Eriesmooth.png
    • If necessary, merge the DEMs
    • Solution to bathymetric (negative) altitude values[2]: using Raster Calculator, add bottom value to all height values
  8. Create STL (Raster -> DEMto3D)[3]

Laser Cut, Low Resolution

Workflow

  1. Download shapefiles from NOAA. Superior bathymetry contours with unclear provenance is available here, with login
  2. Download and install QGIS and Inkscape[4]
  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[5], 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')[6]
    • 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

Other Sources

References

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