7

I am creating a geodesic polyhedron using PostGIS. I am dividing the base icosahedron points (see Appendix A), into the first level of Class I subdivisions.

The setup

For simplicity let's take just the first triangle formed by these points:

CREATE TABLE pix (
    id bigserial
  , name text
  , geog geography(POINT,4326) -- use WGS84
);

insert into pix values (DEFAULT, 'China', st_point(122.3, 39.1));
insert into pix values (DEFAULT, 'Norway', st_point(10.53619898, 64.7));
insert into pix values (DEFAULT, 'Arabian sea', st_point(58.15770555, 10.44734504));

Which creates the base icosahedron triangle

enter image description here

The numbers of the triangle sides correspond to the row numbers shown in the table generated by this query

select p1.name
  , p2.name
  , st_distance(p1.geog, p2.geog) as dist
from pix p1
  cross join pix p2
where p1.id < p2.id

table of distances between 3 points

The points are roughly equidistant. There are small variances because the seed data above seems to have been generated assuming a sphere, not spheroid, which st_distance is returning - use_spheroid false in st_distance will yield much closer distances, but note that using false in all the presented queries does not change the problem described below.

First level of Class I subdivision

insert into pix (name, geog)
select p1.name || '-' || p2.name
  , st_project(p1.geog, st_distance(p1.geog, p2.geog) / 2, st_azimuth(p1.geog, p2.geog)) as geog
from pix p1
  cross join pix p2
where p1.id < p2.id

The above generates the mid-points of each of the pairs of points on the original triangle

enter image description here

And these have the following distances, generated by this query

select p1.name
  , p2.name
  , st_distance(p1.geog, p2.geog) as dist
from pix p1
  cross join pix p2
where p1.id < p2.id
  and p2.id > 3
order by dist

enter image description here

The problem

The distances of the segments 7, 8 and 9 (around 4000kms) are much longer than the distances of segments 1-6 (around 3500kms).

Why is this algorithm not creating something more close to equliateral triangles? Is my use of st_azimuth and st_project wrong? Or did I miss some more basic maths about projections of points onto spheres?

0
3
+50

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, when projected on a sphere. However, if we take our non-projected icosahedron and subdivide all triangles, the generated points will not be on that sphere. When these move outwards, they will be bigger, like the redtriangles in this image.

The first iteration

Also, these constructions are used in architecture. Using the Geodesic dome calculator, (image from their site), we find for an edge length B of 3.5, the corresponding length for A is 3.9:

enter image description here

So indeed it was the mathematics that you missed.

2
  • 1
    Since posting I've had growing suspicions that that was the problem. I had been seduced by beautiful pictures of seemingly regular DGGSs. The revelation from your answer is that regular DGGSs are not possible, and I am most surprised that this issue is not discussed more clearly in the literature I found. Or if I'm wrong, how can I go about making one? – poshest Nov 30 '20 at 13:35
  • There are exactly 5 ways to do such a thing, that correspond to the 5 platonic solids, with 4, 6, 8, 12 or 20 circles over the entire globe. For what you want, using the icosahedral geodesic is actually a good approach. I'm not whether it will become better as you increase the number of subdivisions... My intuition says not, because the distance from the icosahedron to sphere, and thus inflation when projecting, do not improve.. another thread – Gevaert Joep Dec 1 '20 at 0:07

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.