Hi every one! I am looking for an efficient way to generate equally spaced points between two coordinates on Earth’s Surface. Something like: range(point1, point2, length=10) I know there is the GMT wrapper, but it needs distances and it is a bit slow. Any ideas?

Interpolating along a small circle of a sphere should be simple. Just use normal range with the geographical coordinates and then compute the arc length for each interval knowing the radius of the circle at that latitude. For ellipsoidal computations I think the -G option of mapproject can be useful…

[image] marianoarnaiz: it is a bit slow. Using the example from the project module (2nd run) julia> tic(); D = project(origin=(-50,10), end_pt=(-10,30), step=10, km=true);toc() elapsed time: 0.0168785 seconds

Maybe Inwas using it wrong? I wish it could handle x,y,z too haha

What do you mean by x,y,z? One limitation of this method (using the proj module) is that computation is on the sphere, not the ellipsoid.

The module can produce points from a point to another: D = project(origin=(-50,10,0), end_pt=(-10,30,-100), step=10, km=true) But I wish it could: D = project(origin=(-50,10,0), end_pt=(-10,30,-100 [this depth in km]), step=10, km=true). The problem is that I am trying to find a good way to inte…

Have you looked into using Geodesy.jl ? It may not have exactly what you need, but maybe it would be easy enough to come up with a custom solution.

But rays are curved because velocity increases with depth. Interpolating along those rays is a kind of creating flow lines when velocity changes. Otherwise it’s just a linear interpolation between 2 3D points.

I have a few problems in my hands. The project output is great as an initial ray for ray tracing. But as the ray traces I need to interpolate the nodes I have in the ray (say 17 nodes and interpolate 10 times along the path of the ray). This is complicated because the interpolation libraries do not …

Is this a GMT plot? Cool. I’m still not sure that I fully understand the problem. Earth is not “round” on the vertical direction but due to velocity increase the rays are curved and any interpolation along the rays paths must take that into account. You can try also the GMT forum, there are seismo…

Yes, a simple plot to see it my ray tracer works well hehe. “roundness of the Earth”. By this I mean, if I interpolate from -67,10 to 10,20 It works But from -175, 0 to 175, 0 (via the minor arc) … I don’t know how to get this done well.

Let me give it a try! Thanks for the idea @joa-quim

Range in Earth's coordinates

New to Julia

PeterSimon April 1, 2021, 4:17pm 13

Also see here

Topic		Replies	Views
Linear interpolation to find a GPS location (latitude, longitude) between two given GPS locations (latitudes, longitudes) using Julia General Usage interpolations , geo	8	3328	June 12, 2019
Geodesy: how to calculate the straight line distance between two locations, which is represented by longitude and latitude? Geo diffeq , gmt , geodesy	29	8724	June 20, 2021
Use linear interpolation between two LLA coordinate system points General Usage interpolations , geo , geodesy	8	3122	September 8, 2019
Converting Longitude & Latitude to X Y on a map using Julia General Usage geo , gmt , geodesy	22	6089	November 4, 2022
Simple use of Geodesy.jl to get euclidean distances? General Usage geo , coordinate-transform , geodesy	23	596	August 29, 2024

Range in Earth's coordinates

Related topics