Difference between revisions of "Great Lakes Bathymetry"
From TSAS Library
JStallings (talk | contribs) |
JStallings (talk | contribs) |
||
| Line 37: | Line 37: | ||
# Add GeoTIFF to model | # Add GeoTIFF to model | ||
# Analyze to DEM (Raster -> Analysis -> DEM [Terrain Models]) | # 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. | # If necessary, [http://www.qgistutorials.com/en/docs/raster_mosaicing_and_clipping.html merge the DEMs] | ||
* 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. | |||
<gallery mode="packed">File:Erie2bumps.png | |||
File:Eriesmooth.png</gallery> | |||
* 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> | ||
Revision as of 13:43, 11 June 2017
Ideas
- Process of creating a 3D map of the area
- Try this map with these instructions
- For OSM information, use Geofabrik.de and reformat if necessary as described here[1]
Resources
- Downloadable Great Lakes bathymetric maps from NOAA
- All bathymetric maps from NOAA
- The USGS Earth Explorer is easier to use than the above, but requires a login and download client
- Natural Earth provides shapefiles for cultural and physical elements
Create STL for CNC Router or 3D Printer
Proprietary Workflow (Aspire)
Video walkthrough here
- Data pulled from NOAA's ETOPO1 images, available directly here which offers a 1-minute resolution of worldwide bedrock
- Add raster image to QGIS
- Select extent desired (Raster -> Extraction -> Clipper)
- Open in Aspire, follow directions from video
Tool speeds used in above image:
| Print Pass | Tool Size | Feed Rate | Plunge Rate | Print Time |
|---|---|---|---|---|
| Rough Pass | .25" square bit | 15.0 mm/sec | 3.0 mm/sec | 1:30 |
| Finish Pass | .125" ball nose | 25.0 mm/sec | 6.0 mm/sec | 4:30 |
Open Source Workflow
- Download ArcGrid DEM from USGS National Map or Earth Explorer (3DEP, 1/3 arc-second DEM)
- Unzip file
- Open QGIS, add Vector, add layer file (.adf)
- Select extent (Raster -> Extraction -> Clipper)
- Save as GeoTIFF
- Add GeoTIFF to model
- Analyze to DEM (Raster -> Analysis -> DEM [Terrain Models])
- If necessary, merge the DEMs
- 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.
- Solution to bathymetric (negative) altitude values[2]: using Raster Calculator, add bottom value to all height values
- Create STL (Raster -> DEMto3D)[3]
Laser Cut, Low Resolution
Workflow
- Download shapefiles from NOAA. Superior bathymetry contours with unclear provenance is available here, with login
- Download and install QGIS and Inkscape[4]
- No need to extract files, open each in QGIS Desktop
- Add OpenStreetMap information if desired (the native support is not working at this time[5], but 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')[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
- Earth satellite imagery from USGS, esp. those over time
- ArcGIS contour map of the Great Lakes
- Live Access Servers hosted by NOAA/PMEL/TMAP provide information on air quality, ocean salinity, temperature, etc.
- LIDAR resources in case we want to machine it. I want to machine it.
- BBBike.org pulls from OSM, but extent is much too small.
- Using plywood (though very heavy) with the proper scale would approximate the contour lines of a topo map...
- Mosaicing and clipping rasters
- OceansMap, a well-aggregated set of realtime data about major bodies of water around the world, including the Great Lakes region
- Historical maps from the Army Corps of Engineers
- Portal for Great Lakes Observing System information from the Integrated Ocean Observing system
- The National Map combines many of the above resources
References
- ↑ Geographic Information Systems Stack Exchange. "Load .osm.pbf file in Qgis." here.
- ↑ Geographic Information Systems Stack Exchange. "Mars DTM has all negative elevations in QGIS. Can't build STL file for 3D print." here.
- ↑ Modified from Laguna Tools. "Max CNC Tutorial 3D Landscape." YouTube. here.
- ↑ Stack Exchange answer here
- ↑ "OpenStreetMap in QGIS." OpenStreetMap Wiki. here.
- ↑ Geographic Information Systems Stack Exchange. "QGIS choose contour line labels." Here