11
votes
Accepted
Why is SQL Server 2019's STDistance Very Slightly Different Than The Vincenty Formula? (Same Up To ~0.0001km)
SQL Server's STDistance function does not compute geodesic distances. Its algorithm essentially generates straight lines in 3D, densifies them and maps them on the ellipsoid in a manner described more ...
- 3,786
7
votes
Get lat/long given current point, distance and bearing
What follows is a development of the formulas found in the Aviation Formulary by Ed Williams and in Movable Type Scripts by Chris Veness, from the formulas of spherical trigonometry, to solve the ...
- 13.4k
5
votes
calculating bounding box from known centre coordinate and zoom
Assuming you are using the Web Mercator projection (Google Maps, MapBox etc), the key lies in these two equations:
where λ is the longitude in radians and φ is geodetic latitude in radians. The x and ...
- 283
5
votes
Accepted
Euclidean and Geodesic Buffering in R
For most purposes it will be accurate enough to use a spherical model of the earth--and the coding will be easier and the calculations much faster.
Following suggestions by M. Kennedy in a comment, ...
- 69.3k
4
votes
Definition of "most northwesterly point" on a landmass?
Draw the landmass in the Mercator projection. Find the most NW point,
e.g., by scaning a NE-SW line across the map. That's it! The Mercator
projection is the one to use because it's the only one ...
- 3,201
4
votes
Accepted
Why are my icosahedron triangle subdivisions not equilateral (PostGIS)?
Yes, this is due to the tiled icosahedron to-sphere projections. I think of this process as:
Take a normal icosahedron, all corners lie on the same sphere, and therefore, the distance is also the same,...
- 1,111
4
votes
Accepted
How is it possible to calculate azimuth from a single lat lon cordinate with x, y offsets in meters?
After couple of days of searching and testing I finally managed to get this to work.
The first question was whether it was possible to calculate azimuth using an image's x,y pixel coordinate.
The ...
- 161
3
votes
Accepted
Complexify geography polygons in PostGIS
I think you can use ST_Segmentize to do this (especially if you have Geography objects).
Returns a modified geometry having no segment longer than the given max_segment_length. Distance computation ...
- 79.9k
3
votes
Figure out location of city direction sign in image - friday fun project
I am the author of this picture. So I can tell for sure the sign post is located in Thessaloniki, Greece.
- 31
3
votes
Accepted
Would nearest point using Geodesic distance and nearest point using Haversine distance be the same point?
For sure, "closest" will return different points, in general, for
great-circle (what you call haversine) and geodesic distances. For a
specific example, consider the set of points {A, B} where the ...
- 3,201
3
votes
Why is SQL Server 2019's STDistance Very Slightly Different Than The Vincenty Formula? (Same Up To ~0.0001km)
Without having access to the source code of SQL Server it may be impossible to say. As a comparison, PostGIS ST_Distance query from the geography example of https://postgis.net/docs/ST_Distance.html ...
- 61.9k
3
votes
Accepted
ArcPy Insert cursor losing vertices when writing Polygon features?
Try adding sr to arcpy.Polygon(polygon, sr)Here is the reason why:
I sometimes see significant coordinate movement when I do not supply a Spatial Reference to an Arcpy Geometry object such as a ...
- 7,455
3
votes
Accepted
What's wrong with using latitude and longitude as X and Y on a Cartesian space?
Everything really. Distances will become meaningless. Areas will become meaningless. Hence all outputs such as slope will become problematic. As you move in latitude then your "grids" will appear to ...
2
votes
Moving dead ahead: rhumb line, great circle or none
I know the question has answered already, but allow me to show you something cool for future use:
A live Google map where you can move the edge markers affecting the great circle and rhumb line same ...
- 21
2
votes
Buffer not providing accurate geodesic buffers in ArcGIS Desktop 10.0?
I just tested this in ArcMap 10.3 and it has produced accurate buffers based on the field 'Dist' in my point feature class: arcpy.Buffer_analysis(in_features="//pmgfile01/Users/jb/Documents/ArcGIS/...
- 7,275
2
votes
JTS Vividsolutions API: compute First (direct) geodetic problem
I would much rather suggest doing it manually since the math is simple. This will be much faster than using so many potentially complex methods.
public static Point createPointInDistanceAtAngle(Point ...
- 10.2k
2
votes
Creating a geodetic line that goes around the whole ellipsoid with pyproj
pyproj hooks up with the geodesic routines (I am the author of
these) in proj.
The geod_directline and geod_inverseline functions both initialize
objects which allow you to compute points on a ...
- 3,201
2
votes
Figure out location of city direction sign in image - friday fun project
If you want to measure (real world)- distance, don't draw circles / create buffers on a flat map canvas. Using EPSG:3857 / Web Mercator is especially a bad choice for measuring distances or areas, as ...
- 63.1k
2
votes
What's the name of the area between two geographic coordinates?
A quadrangle?
Yes, it is named: spheroidal quadrangle.
- 13.4k
2
votes
Accepted
When calcuating distance between points on Earth why are my Haversine vs. Geodesic calculations wildy diverging?
I appreciate your code because it contributes to geodesy, by allowing to make some interesting analyzes in a practical way, regarding spherical approximations for the calculation of length of geodesic ...
- 13.4k
2
votes
When calcuating distance between points on Earth why are my Haversine vs. Geodesic calculations wildy diverging?
There was a matrix algebra error in the Haversine formula. I updated the code in the question. I am getting much better agreement between Haversine and geodesic now:
On my actual dataset:
- 193
2
votes
Accepted
Figure out location of city direction sign in image - friday fun project
Just eyeballing this on a globe view with a measure tool like: https://map.terria.io/
Makes me think the location is in Southern Greece.
Using an online app like https://www.distancefromto.net we can ...
- 13.2k
2
votes
Choosing projection for running polygon intersections at global scale (i.e, geodesic intersections)?
You should not use a projection (except for making maps; or when you must use an algorithm that really cannot deal with angular data). For area computations in R, you should use the lon/lat data. That ...
- 10.2k
2
votes
Is order among distances preserved when using Geodesic and Cartesian distances?
No! And it's pretty easy to find counterexamples (where the cartesian distances
are ranked opposite to the geodesic distances). For example, for the
WGS84 ellipsoid, pick 3 points with positions
...
- 3,201
1
vote
Accepted
What's the name of the area between two geographic coordinates?
More informally than a spherioidal quadrangle, a much more common and useful term is "bounding box", sometimes shortened to "bbox".
- 8,598
1
vote
Leaflet.Geodesic.js line not drawn
var Geodesic = L.geodesic([[[52.5, 13.35], [33.82, -118.38]]], {
opacity: 1,
color: 'red',
steps: 50
}).addTo(map);
Instantiates a polyline object given an array [point_1, point_2, point_3] ...
- 29.9k
1
vote
Converge origins on destination leaflet geodesic
Simply instantiate your L.Geodesics inside the for (var i in sheet) loop, like
var fg = L.featureGroup().addTo(map);
var orleans = new L.LatLng(30.053206, -89.934873);
var origin;
Tabletop.init({
...
- 10.1k
1
vote
Calculate perimeter of a large polygon using geodesic calculations
You can achieve that using web mercator projection and addGeometryAtributes form ArcGIS
1
vote
Accepted
Calculate perimeter of a large polygon using geodesic calculations
In QGIS, it seems that if you have the project coordinate reference system set to an ellipsoid and not a projection, then calculating the perimeter via the field calculator will automatically use the ...
- 3,480
1
vote
Accepted
Understanding geodesic envelope?
It looks like you've figured things out OK. In case it helps, I give here the MATLAB/Octave code that I used to generate this Wikipedia figure showing the envelope:
The envelope is formed by the ...
- 3,201
Only top scored, non community-wiki answers of a minimum length are eligible