0

So there is the cassini projection in Proj. It is possible to use it in the OSGeo4WShell. example: +proj=cass +lon_0=.... +lat_0=.... +x_0=0 +y_0=0 +ellps=bessel +units=m +no_defs

Is there a way to define such CRS in Python?

example for ecef system: ecef = pyproj.Proj(proj='geocent', ellps='WGS84', datum='WGS84')

But I want to use the cassini projection (to transform to ecef).

How do I define it?

Is CAssini implemented in pyproj?

1

pyproj is a python wrapper for PROJ. PROJ currently supports cassini, so pyproj should as well (https://proj.org/operations/projections/cass.html).

You likely want the CRS and Transformer classes:

import pyproj

cass = pyproj.CRS("+proj=cass +lon_0=.... +lat_0=.... +x_0=0 +y_0=0 +ellps=bessel +units=m")
ecef = pyproj.CRS(proj='geocent', ellps='WGS84', datum='WGS84')
transformer = pyproj.Transformer.from_crs(cass, ecef)

Also see: https://pyproj4.github.io/pyproj/stable/examples.html

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.