<div dir="ltr"><div>Hi Roger,</div><div><br></div><div>To convert altitudes you just need to provide a Z... and to have the proper geoid grids. You are lucky, because SWEN17_RH2000 is available in PROJ-data</div><div>You can have a look at <a href="https://cdn.proj.org/">https://cdn.proj.org/</a> to see which grid files are available, and their value in a point.</div><div><br></div><div>To see the results, you can compare<br></div><div>echo 59.3 18.0 10 | PROJ_NETWORK=OFF cs2cs EPSG:4326+5773 EPSG:3006+5613<br>6577304.57      670817.92 10.00</div><div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><br></div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div>echo 59.3 18.0 10 | PROJ_NETWORK=ON cs2cs EPSG:4326+5773 EPSG:3006+5613<br>6577304.57   670817.92 10.12</div><div><br></div><div>If you do not want to use the network, you can include the gridfiles in your distribution. Set how PROJ_LIB works. There is an API to specify the paths in C as well.</div><div>For this particular case you need the gridfiles for EGM96 and  SWEN17_RH2000. You can see more info in</div><div>projinfo -s EPSG:4326+5773 -t EPSG:3006+5613</div><div><br></div><div>See that not every geoid file in the world is included in PROJ-data ;)<br></div><div><br></div><div>Cheers</div><div>Javier<br></div></div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature">.___ ._ ..._ .. . ._.  .___ .. __ . _. . __..  ... .... ._ .__<br>Entre dos pensamientos racionales <br>hay infinitos pensamientos irracionales.<br><br></div></div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, 25 Mar 2021 at 11:30, Roger Oberholtzer <<a href="mailto:roger.oberholtzer@gmail.com">roger.oberholtzer@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">We have usually been focused on latitude/longitude when using proj. We<br>
would now like to address altitude. I have looked in the proj docs,<br>
but I cannot say that I am much the wiser in this respect.<br>
<br>
We currently use the old C API (proj_api.h). But we are moving to the<br>
new C API. That's a work in progress. And perhaps the subject of other<br>
questions.... Here I am considering height.<br>
<br>
Our typical scenario is that we collect data from a receiver (high-end<br>
dual antenna - so good quality) in a moving vehicle. We collect<br>
LAT/LONG in WGS84, and then project it using proj. For example to<br>
SWEREF99. So it's basically EPSG:4326 ->  EPSG:3006, or the inverse.<br>
Of course we don't only use SWEREF99. It's just an example.<br>
<br>
The receivers are configured to collect altitude as the geoid in<br>
EGM96. I think it is in fact the default for our receivers. So I guess<br>
it is a common altitude.<br>
<br>
To date, we have just let the altitude go along for the ride when<br>
using proj. Our users can select if they want the altitude as supplied<br>
by the receiver, or as proj calculates. There is a difference in the<br>
values.<br>
<br>
I have never understood how altitude is described in proj. For<br>
example, using the example of EPSG:3006, the corresponding altitude<br>
could be SWEN17_RH2000 (See<br>
<a href="https://www.lantmateriet.se/en/maps-and-geographic-information/gps-geodesi-och-swepos/Om-geodesi/Transformationer/Geoidmodeller/" rel="noreferrer" target="_blank">https://www.lantmateriet.se/en/maps-and-geographic-information/gps-geodesi-och-swepos/Om-geodesi/Transformationer/Geoidmodeller/</a>).<br>
<br>
Can I use proj to get SWEN17_RH2000 from EGM96? Perhaps I just didn't<br>
recognize the information in the docs when I looked through it.<br>
<br>
Or should I be looking elsewhere?<br>
<br>
Sorry for being a noob here. But that's the state of my understanding<br>
with respect to altitude.<br>
<br>
-- <br>
Roger Oberholtzer<br>
_______________________________________________<br>
PROJ mailing list<br>
<a href="mailto:PROJ@lists.osgeo.org" target="_blank">PROJ@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/proj" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/proj</a><br>
</blockquote></div>