[GIS] Aligning coordinate systems: gpx, LiDAR, shapefile

coordinatesgpxlidarshapefile

I am working on a project based on a lidar data set (for use in LidarViewer in Linux, a basic x,y,z viewer) that was downloaded from #1 below at Opentopo.org and looking to add onto this data set with a shapefile based layer and a gpx based layer.

In ArcMap 10.0 I am converting the shapefile (road perimeter around lake) into points, then Defining the Projection in Arc Toolbox's Projections and Transformation. Since Opentopo stated UTM Zone 10, I set that to the roads shapefile, exported, completed my .shp to .las(lidar file) conversion and sent it to LidarViewer. In LidarViewer my LiDAR data (only a section of the lake's perimeter) is floating the in the middle of the lake (with perspecitve to the road points).

So, I switched the Tahoe to the original NAD27 in ArcMap, went through the process again with no change in position. Do I have to have a shapefile based on NAD83 for this to work, or is there another way?

Lidar viewer imports .las or .txt/.csv files and only reads a x,y,z no coordinate system whatsoever.

The other related question is with the .gpx data that I am downloading from the #3 below, it is a .gpx file of a bike path around the perimeter of the lake. I am using the following link to convert .gpx into a text file (to process and import into Lidar Viewer) and the site provides the ability to output long/lat plus UTM coordinates (the output included both zone 10 and 11 South ) This .gpx file also does not project correctly. I was thinking that I should possible try another converter?

gpsvisualizer .com convert_format=gpx

Any input regarding the above two questions are appreciated, I can also further clarify any issues that may not seem clear.

Thanks!

1) opentopography-org -> Data -> South Lake Tahoe, Nevada

Horizontal Coordinates: UTM Zone 10 N NAD83 Meters [EPSG: 26910] – Vertical Coordinates: NAVD88

2) tahoe usgs gov/DLG UTM NAD27 zone 10 projection

3) bicycling. trimbleoutdoors. com/ViewTrip.aspx?tripID=153934

Best Answer

You are on the right track as far as trying to match up the coordinate systems of all of your datasets. The problem you are encountering is not so much a problem of the software projecting the data incorrectly, but of using the wrong operation to define the projection of your data.

This answer directly discusses the difference between the Define and Project commands in ArcGIS. It also has links to the relevant ESRI help files discussing the different commands.

Let's start with what factors you are able to control.

  1. You have a shapefile that defines your roads around Lake Tahoe, in UTM Zone 10 NAD 27, from the USGS. Generally, spatial data from there is going to have the projection already defined, meaning that when you look at the properties in ArcMap, it will list the Coordinate System. You can use ArcGIS to Reproject this layer as necessary, in this case, to the UTM Zone 10 NAD 83 Meters, projection.
  2. You have a .gpx file of GPS coordinates that defines a bike trip. You have the ability to convert this to a file of coordinates in either Lat/Lon or UTM. This file could then be brought into ArcGIS and reprojected as necessary.
  3. You know the projection of the Lidar data from opentopography.org. It is UTM Zone 10N NAD 83 Meters.

The one thing you are unable to change, because of the format, is the Lidar data, which you are going to read using the LidarViewer software. You thus, must ensure that the coordinate systems of all the other data layers, matche the Lidar, since you will be importing the other data into the Lidarviewer as well, to overlay.

Here is the process.

  1. Look at the Roads shapefile. In ArcMap, if you look at the properties, and shows that it has a coordinate system defined, you will need to use the Project command to Reproject the data from UTM 10 NAD 27, to UTM 10 NAD 83.
  2. a. Convert the .gpx data using the website, and save as Lat/Lon. This should give you a text file of coordinates in degrees.
    b. Add this using the Add X/Y Data command in ArcGIS, specifying the coordinate system as: Geographic: NAD 83.
    c. Use the Project tool on this new layer, and choose the UTM Zone 10, NAD 83 Meters coordinate system as your output.

What you should end up with is all 3 input layers being in the same coordinate system. You can now import them into the LidarViewer and they should overlay properly.

A note about the LidarViewer software: Even though it does not support coordinate systems directly, it will still overlay layers properly if they are all in the same coordinate system to begin with. This simply means that all of the projection has to be performed prior to the data import.

Related Question