How to get 3D GPS p...
 
Notifications
Clear all

How to get 3D GPS positions from geographic latitude, longitude and altitude?

15 Posts
11 Users
0 Reactions
3 Views
(@roybar)
Posts: 1
Registered
Topic starter
 

I have a list of GPS measurements for a point for a period of about 2 years. They are latitude, longitude and altitude. I want to generate the time series for the east, north and up components. Do I just need to convert the geographic coordinates to UTM coordinates or should I do something else?

 
Posted : 21/04/2021 5:37 am
(@norman-oklahoma)
Posts: 7610
Registered
 

I gather that you want to convert your spherical coordinates (ie/lat-long) to grid coordinates. There are various projection systems designed to do that. UTM is one. Whether UTM is the most appropriate for your needs depends on where you are and what you are doing with the results.?ÿ

Note that conversion of the lat/long to grid is a separate function from that of converting the heights from ellipsoidal to orthometric.?ÿ?ÿ

 
Posted : 21/04/2021 11:43 am
(@dave-lindell)
Posts: 1683
 

If it's only one point I would suppose you could draw the actual movement to true scale. Just make the y-axis represent seconds of arc/altitude, not needing conversion. The x-axis could be time, and depending on the frequency of the observations (daily? weekly? monthly?) could be spaced at 0.2 inch/5 mm for weekly data to fit on an 11"x17" ledger sheet.

A graph with red (northing), blue (easting), and green (altitude/elevation), for example, would do nicely.

 
Posted : 21/04/2021 12:01 pm
(@geeoddmike)
Posts: 1556
Registered
 

The original poster wants to transform latitude, longitude and heights (LLh) to North, East and Up. What he wants are 3D values in order to generate a time series of differences.

I suggest reviewing http://www.mygeodesy.id.au/documents/LOCAL%20HORIZON%20COORDINATES.pdf

At this time, I cannot recollect any Internet-based tool that accepts as input LLh and returns ENU.

The paper linked above nicely describes the issue and shows how to do the rotations. Coding the problem is straightforward in Matlab or similar tools.

The paper is by Deakin. I highly recommend his site as a source for well-written papers on geodesy and related subjects.

 
Posted : 21/04/2021 12:15 pm
(@geeoddmike)
Posts: 1556
Registered
 

FWIW, the NGS Invers3D tool does provide results in both DXDYDX and DEDNDU. I also add a graphic from the Deakin paper referenced in my previous post.

9D08C87D 217F 4E3C 825B B534C8126C07
9CD24B68 84D2 4480 AF9C ED69BE335CAF

?ÿ

 
Posted : 21/04/2021 12:24 pm
(@john-hamilton)
Posts: 3347
Registered
 

Since he wants a time series, there is no need to convert ellipsoid height to orthometric. Both change by the same amount

 
Posted : 21/04/2021 1:26 pm
(@jonathan50)
Posts: 118
Registered
 

Isn't latitude, longitude and height already in a form of xyz coordinate?

 
Posted : 21/04/2021 10:36 pm
(@mightymoe)
Posts: 9920
Registered
 

@jonathan50

Seems like needless calculation when the answer is right there.

 
Posted : 22/04/2021 1:25 am
(@steven-metelsky)
Posts: 277
Registered
 

You have to be careful here. What is the basis of Lat and Long?

?ÿ

Geodetic or Geocentric? You have to know what you're dealing with so you can make the proper conversion.

?ÿ

 
Posted : 22/04/2021 4:45 am
(@john-hamilton)
Posts: 3347
Registered
 

@jonathan50

Latitude and Longitude are different lengths which vary by latitude (latitude much more slowly, but it does vary). If one wants to know how much something has moved, they are not as useful to plot directly.?ÿ

One can upload a file of lat/long to NCAT and get grid coordinates returned.

?ÿ

To answer the original question...just convert to grid (UTM, SPC, etc), nothing more is required.?ÿ

 
Posted : 22/04/2021 5:33 am
(@bill93)
Posts: 9834
 
Posted by: @steven-metelsky

Geodetic or Geocentric?

That isn't going to matter for most projects as long as you stay consistent.

 
Posted : 22/04/2021 5:41 am
(@paul-in-pa)
Posts: 6044
Registered
 

Depending on the Latitude, Longitude and Altitude? precisions he has, E, N and Up may just plot as an imprecise blob. We do need further explanation of what he calls Altitude.

Paul in PA

 
Posted : 22/04/2021 5:46 am
(@geeoddmike)
Posts: 1556
Registered
 

Pending clarification from the OP, I interpret his reference to ƒ??altitudeƒ? to mean height. Also, his reference to UTM coordinates as a potential solution makes he wonder.

If the OP has the output from a GPS-based survey expressed in latitude, longitude and ellipsoid height and wants to transform the results to a more understandable NEU, he should use the Deakin paper.?ÿ

IF he has the vector information from a reference/base station to his site/s, he can transform the DXDYDZ information to DNDEDU. I would prefer this approach.?ÿ

Some have argued that the DXDYDZ values are sufficient. As the OP wants differences in terms of North, East and Up an additional step is required. Take a look at the output of the NGS tool FORWARD3D:

95FE4C7A 943A 42BC 91B4 5D2EDFD8D060

This output is based on a forward azimuth, zenith and ellipsoid distances observed at the ƒ??fromƒ? station. You can calculate the NEU differences from the ƒ??fromƒ? station using simple trigonometry. ?ÿ?ÿ

There are lots of folks way smarter than me, if they can look at the DXDYDZ values alone and break the differences into NEU, more power to them. For me, the NEU differences are much more relevant when trying to understand what is happening to a site.?ÿ

As I do not know what the OP intends to monitor nor do I know if he has other data sources (e.g. differential leveling in NAVD88 or even NGVD29), I encourage him to consult with his NGS Regional Advisor see:?ÿ https://www.ngs.noaa.gov/ADVISORS/?ÿ

Awaiting correction or enlightenment,

?ÿ

DMM

 
Posted : 22/04/2021 9:48 am
(@loyal)
Posts: 3735
Registered
 

I have computed several ƒ??time seriesƒ? projects both for my private use as well as a couple of client projects. As with everything else, there is more than one way to skin this cat, and the choice of skinning technique (chili verses fur coat for SWMBO) depends A LOT on what you hope to determine.

The Lat/Lon/Height to N/E/U computations are the easy part!

The metadata supporting the input geographic coordinates (Latitude/Longitude/Height) of the subject station is the potential WILDCARD in the recipe (GIGO).

Without knowing the WHY, it's pretty hard to suggest a HOW.

 
Posted : 22/04/2021 11:15 am
(@steven-metelsky)
Posts: 277
Registered
 

@bill93 right

My line of thought was when receiving data from outside sources. You better know which type of Lat Long you're using.

Also, you need to be educated enough to understand what your collector is giving you vs. what you're plugging into it.

 
Posted : 23/04/2021 4:58 am