Note
Go to the end to download the full example code.
Cassini cylindrical projection
This cylindrical projection was developed in 1745 by César-François Cassini de Thury for the survey of France. It is occasionally called Cassini-Soldner since the latter provided the more accurate mathematical analysis that led to the development of the ellipsoidal formulae. The projection is neither conformal nor equal-area, and behaves as a compromise between the two end-members. The distortion is zero along the central meridian. It is best suited for mapping regions of north-south extent. The central meridian, each meridian 90° away, and equator are straight lines; all other meridians and parallels are complex curves.
clon0/lat0/scale or Clon0/lat0/width
The projection is set with c or C. The projection center is set by lon0/lat0, and the figure size is set with scale or width.
import pygmt
fig = pygmt.Figure()
# Use the ISO code for Madagascar (MG) and pad it by 2 degrees (+R2)
fig.coast(
region="MG+R2",
projection="C47/-19/12c",
frame="afg",
land="gray80",
water="steelblue",
)
fig.show()
Total running time of the script: (0 minutes 0.285 seconds)