AnswerBun.com

Measuring distances on custom ellipsoid

Geographic Information Systems Asked by Coel on August 20, 2020

I’m trying to set up a QGIS project for creating a map of a fictional planet smaller than Earth. See here for a helpful question that doesn’t quite touch on what I’m trying to accomplish (I don’t mind a non-cylindrical projection, and don’t have any shapefiles or rasters that the map has to "fit").

My primary goal is to be able to draw continents as if viewing their orthographic projection, while having the ability to check their ellipsoid measurements. After several hours of searching, I still haven’t found if that capability exists in QGIS, so I’ve contented myself with a lonlat world map. The problem is, I still can’t get the equatorial distance from -180 to 180 degrees to match the circumference I calculated.

Here is the proj4 string for the custom CRS:

+proj=lonlat +a=5733900 +lat_ts=0 +lon_0=0 +rf=700 +units=m +no_defs

(The fictional planet has an equatorial radius of 5733,9 km and a polar radius of 5725,7 km, and should have an equatorial circumference of 36027,1 km if I’ve calculated right.)

I have set Project Properties>General>Measurements to custom ellipsoid with the above mentioned equatorial and polar radii, and set the custom CRS as the project’s CRS.

I have created a vector layer with the custom CRS, a polyline that should show the bounds of the map. It is a rectangle from degree coordinates -180 to +180 x and -90 to +90 y. However, when I measure the equatorial distance across this rectangle (at y=0), I get wildly varying distances. Here is an example:
measurement issue from degree coordinates 0,0 to 0,+180 - half the circumference should be 18013,6
Even when I have defined the coordinate reference system as eqc, I still get an incorrect circumference at the equator.

What do I need to change in order to be able to get realistic measurements? Is there any way to work on all parts of the ellipsoid in one QGIS project (and have them have fairly orthographic dimensions when I zoom in on them), or should I be looking into some other software?

One Answer

Forget all kind of ellipsoidal calculations for non-Earth CRSes in QGIS.

But you can use the CRSes for the layers, and transform between geographic and projected non-Earth CRSes. If you export the layer to the custom eqc system, you can measure the Cartesian length of the geometries and they will be right.

Answered by Gabriel De Luca on August 20, 2020

Add your own answers!

Related Questions

Progress bar does not show progress during processing

1  Asked on December 9, 2020 by francisco-danubio-salas-rosett

   

Error trying to cast an ee.Image to a numpy array

1  Asked on December 9, 2020 by panchester

 

No repositories accesible

0  Asked on December 9, 2020 by user167525

 

Using subtype property of the esriaddinx for a AddIn button

0  Asked on December 9, 2020 by david-lewis

 

MBTiles on Leaflet map

0  Asked on December 9, 2020 by user1655130

     

Tilemill/CartoCSS: how to center label in a polygon

1  Asked on December 9, 2020 by vanuan

   

Filter Leaflet markers by input selectionin R Flexdashboard

1  Asked on December 9, 2020 by ryan-garnett

       

NDVI out of range?

1  Asked on December 8, 2020

   

QGIS Vector Toolbar Has No Tools

1  Asked on December 8, 2020 by charlie-white

     

Developing QGIS Action which loads a raster

2  Asked on December 8, 2020 by daniele-piccolo

       

Ask a Question

Get help from others!

© 2022 AnswerBun.com. All rights reserved. Sites we Love: PCI Database, MenuIva, UKBizDB, Menu Kuliner, Sharing RPP, SolveDir