Contour Tiles
cevely
Global Mapper UserTrusted User
in Bug Report
I have a large dataset of LAS tiles that I am trying to create contours from, with the result being a separate polyline shapefile for each LAS file. To reduce the memory usage to a manageable level, I have written a script to process a single tile at a time.
This script does the following for each of the LAS tiles:
- Loads the LAS points from the tile in question, as well as the points from the surrounding 8 tiles to cover the edges.
- Loads a breakline shapefile
- Generates an Elevation Grid using the points plus the breaklines as constraints for the bounding rectangle of all 9 tiles.
- Generates contours from the Elevation Grid, also to the bounds of all 9 tiles
- Exports a line shapefile from the contour layer, clipped to the bounds of the original tile.
This workflow works extremely well, with the one lone, huge exception: along the tile edge, the contours do not line up perfectly. The endpoints will typically be offset by around .01 feet or so.
I have tried several different combinations of parameters, including disabling smoothing, increasing the overedge and overriding the default grid size, but the resulting files still have the same issue.
I'm hoping there is simply a parameter I am overlooking. Any help you could give would be very much appreciated.
Thanks!
Chad Evely
This script does the following for each of the LAS tiles:
- Loads the LAS points from the tile in question, as well as the points from the surrounding 8 tiles to cover the edges.
- Loads a breakline shapefile
- Generates an Elevation Grid using the points plus the breaklines as constraints for the bounding rectangle of all 9 tiles.
- Generates contours from the Elevation Grid, also to the bounds of all 9 tiles
- Exports a line shapefile from the contour layer, clipped to the bounds of the original tile.
This workflow works extremely well, with the one lone, huge exception: along the tile edge, the contours do not line up perfectly. The endpoints will typically be offset by around .01 feet or so.
I have tried several different combinations of parameters, including disabling smoothing, increasing the overedge and overriding the default grid size, but the resulting files still have the same issue.
I'm hoping there is simply a parameter I am overlooking. Any help you could give would be very much appreciated.
Thanks!
Chad Evely
Comments
-
I recommend getting the latest release of Global Mapper to see if the problem has already been fixed. If the problem persists in the latest release, then your best bet is to send an e-mail describing your workflow, your sample data, and results to Blue Marble Geographics support (geohelp@bluemarblegeo.com). This will ensure that the support team sees the problem and can respond in a timely manner. Thank you.
Howdy, Stranger!
Categories
- 12.2K All Categories
- 5.3K Features Discussion
- 312 Downloading Imagery
- 1.2K Elevation Data
- 371 Georeferencing Imagery Discussion
- 580 GM Script Language
- 42 User Scripts
- 111 GPS Features
- 388 Projection Questions
- 788 Raster Data
- 1.3K Vector Data
- 6.4K Support
- 159 Announcement and News
- 881 Bug Report
- 554 SDK
- 1.1K Suggestion Box
- 3.6K Technical Support
- 520 Other Discussion
- 125 GIS Data Sources
- 26 Global Mapper Showcase
- 217 How I use Global Mapper
- 99 Global Mapper Forum Website