Difference between revisions of "Great Lakes Bathymetry"

From TSAS Library
Jump to: navigation, search
Line 9: Line 9:
* [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


==CNC Carve==
==Laser Cut, Low Resolution==
===Workflow===
===Workflow===
# 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
# 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
Line 18: Line 18:
# 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>
# 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}}
{{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}}
===Workflow Take 2===
* Try [http://apdrc.soest.hawaii.edu/las8/UI.vm#panelHeaderHidden=false;differences=false;autoContour=false;globalMin=-237;globalMax=1667;xCATID=C1354ADCD9A2E42563E67810557FE547;xDSID=etopo1;varid=b_bathy-etopo1;imageSize=auto;over=xy;compute=Nonetoken;catid=C1354ADCD9A2E42563E67810557FE547;dsid=etopo1;varid=b_bathy-etopo1;avarcount=0;xlo=-92.788119178879;xhi=-69.93315874592;ylo=41.247106546728;yhi=49.388349427462;operation_id=Plot_2D_XY_zoom;view=xy;ferret_image_format=gif;ferret_use_graticules=none;ferret_margins=false;ferret_deg_min_sec=false;ferret_contour_style=color_filled_plus_lines;ferret_contour_levels=20 this] map with [http://www.instructables.com/id/Lasercut-Bathymetric-Map-of-Arabian-Sea/ these] instructions
* For OSM information, use [http://download.geofabrik.de/north-america.html Geofabrik.de] and reformat if necessary as described [https://gis.stackexchange.com/questions/39112/load-osm-pbf-file-in-qgis here]<ref>Geographic Information Systems Stack Exchange. "Load .osm.pbf file in Qgis." [https://gis.stackexchange.com/questions/39112/load-osm-pbf-file-in-qgis here].</ref>


==Create STL for CNC Router or 3D Printer==
==Create STL for CNC Router or 3D Printer==

Revision as of 09:46, 10 June 2017

Greatlakes.jpg

Ideas

Process of creating a 3D map of the area

Resources

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[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

Create STL for CNC Router or 3D Printer

  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[4]: using Raster Calculator, add bottom value to all height values
  8. Create STL (Raster -> DEMto3D)[5]

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. "Mars DTM has all negative elevations in QGIS. Can't build STL file for 3D print." here.
  5. Modified from Laguna Tools. "Max CNC Tutorial 3D Landscape." YouTube. here.