PHPCoord\CoordinateOperation\CoordinateOperationMethods

Complexity

Total Complexity

155

Size/Duplication

Total Lines	2497
Duplicated Lines	0 %

Test Coverage

Coverage

37.86%

Importance

Changes	1
Bugs	0	Features	0

146 Methods

Rating	Name	Duplication	Size	Complexity
A	EPSG_MADRID_TO_ED50_POLYNOMIAL()	0	2	1
A	EPSG_DEGREE_REPRESENTATION_CONVERSION_HDM_TO_DMSH()	0	2	1
A	EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_OSGM02_IRE()	0	2	1
A	EPSG_GEOGRAPHIC3D_TO_2D_CONVERSION()	0	2	1
A	EPSG_TIME_SPECIFIC_POSITION_VECTOR_TRANSFORM_GEOCEN()	0	2	1
A	EPSG_MERCATOR_VARIANT_B()	0	2	1
A	EPSG_NTV1()	0	2	1
A	EPSG_GEOCENTRIC_TRANSLATIONS_GEOG2D_DOMAIN()	0	2	1
A	EPSG_LAMBERT_CONIC_CONFORMAL_1SP()	0	2	1
A	EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_IGN1997()	0	2	1
A	EPSG_DEGREE_REPRESENTATION_CONVERSION_HDEG_TO_DMSH()	0	2	1
A	EPSG_VERTICAL_OFFSET_BY_GRID_INTERPOLATION_BEV_AT()	0	2	1
A	EPSG_COMPLEX_POLYNOMIAL_OF_DEGREE_3()	0	2	1
A	geocentricPositionVectorTransformation()	0	24	1
A	EPSG_LAMBERT_CONIC_CONFORMAL_2SP_MICHIGAN()	0	2	1
A	EPSG_LONGITUDE_ROTATION()	0	2	1
A	EPSG_VERTICAL_OFFSET_BY_GRID_INTERPOLATION_NZLVD()	0	2	1
A	EPSG_ALBERS_EQUAL_AREA()	0	2	1
A	EPSG_VERTICAL_OFFSET_BY_GRID_INTERPOLATION_ASC()	0	2	1
A	EPSG_TIME_DEPENDENT_COORDINATE_FRAME_ROTATION_GEOCEN()	0	2	1
A	EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_GRAVSOFT()	0	2	1
A	geographicAbridgedMolodensky()	0	35	3
A	EPSG_NTV2()	0	2	1
A	EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_EGM2008()	0	2	1
A	EPSG_LAMBERT_AZIMUTHAL_EQUAL_AREA()	0	2	1
A	EPSG_NADCON()	0	2	1
A	EPSG_BONNE()	0	2	1
A	EPSG_LABORDE_OBLIQUE_MERCATOR()	0	2	1
A	operationByEPSGCode()	0	18	1
A	EPSG_VERTICAL_PERSPECTIVE_ORTHOGRAPHIC_CASE()	0	2	1
A	EPSG_KROVAK_NORTH_ORIENTATED()	0	2	1
A	EPSG_EQUAL_EARTH()	0	2	1
A	EPSG_MERCATOR_SPHERICAL()	0	2	1
A	EPSG_VERTICAL_OFFSET_BY_GRID_INTERPOLATION_GTX()	0	2	1
A	EPSG_DEGREE_REPRESENTATION_CONVERSION_DMH_TO_DMSH()	0	2	1
A	EPSG_GEOCENTRIC_TOPOCENTRIC_CONVERSIONS()	0	2	1
A	EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_GTX()	0	2	1
A	EPSG_CHANGE_OF_VERTICAL_UNIT()	0	2	1
A	geographicPositionVectorMolodenskyBadekas()	0	38	2
A	EPSG_NEW_ZEALAND_MAP_GRID()	0	2	1
A	EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_OSGM_GB()	0	2	1
A	EPSG_LAMBERT_AZIMUTHAL_EQUAL_AREA_SPHERICAL()	0	2	1
A	EPSG_GEOG3D_TO_GEOG2D_PLUS_GRAVITYRELATEDHEIGHT_US_GTX()	0	2	1
A	EPSG_GEOG3D_TO_GEOG2D_PLUS_VERTICAL_AUSGEOID_V2()	0	2	1
A	geocentricCoordinateFrameMolodenskyBadekas()	0	33	1
A	EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_CI()	0	2	1
A	EPSG_REVERSIBLE_POLYNOMIAL_OF_DEGREE_2()	0	2	1
A	EPSG_POINT_MOTION_GEOCENTRIC_CARTESIAN()	0	2	1
A	geographicCoordinateFrameMolodenskyBadekas()	0	38	2
A	EPSG_TIME_DEPENDENT_POSITION_VECTOR_TFM_GEOG3D()	0	2	1
A	EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_AUSGEOID_V2()	0	2	1
A	EPSG_REVERSIBLE_POLYNOMIAL_OF_DEGREE_13()	0	2	1
A	EPSG_TIME_DEPENDENT_COORDINATE_FRAME_ROTATION_GEOG3D()	0	2	1
A	EPSG_BONNE_SOUTH_ORIENTATED()	0	2	1
A	EPSG_GUAM_PROJECTION()	0	2	1
A	EPSG_EQUIDISTANT_CYLINDRICAL_SPHERICAL()	0	2	1
A	EPSG_DEGREE_REPRESENTATION_CONVERSION_DMS_TO_DMSH()	0	2	1
A	EPSG_ORDNANCE_SURVEY_NATIONAL_TRANSFORMATION()	0	2	1
A	EPSG_POINT_MOTION_ELLIPSOIDAL()	0	2	1
A	EPSG_TIME_DEPENDENT_POSITION_VECTOR_TFM_GEOCENTRIC()	0	2	1
A	EPSG_VERTICAL_OFFSET()	0	2	1
A	EPSG_CARTESIAN_GRID_OFFSETS_FROM_FORM_FUNCTION()	0	2	1
A	EPSG_VERTICAL_PERSPECTIVE()	0	2	1
A	EPSG_MERCATOR_VARIANT_A()	0	2	1
A	EPSG_AFFINE_GEOMETRIC_TRANSFORMATION()	0	2	1
A	EPSG_HYPERBOLIC_CASSINI_SOLDNER()	0	2	1
A	EPSG_GENERAL_POLYNOMIAL_OF_DEGREE_4()	0	2	1
A	EPSG_LAMBERT_CONIC_CONFORMAL_2SP_BELGIUM()	0	2	1
A	geocentricCoordinateFrameRotation()	0	24	1
A	EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_EGM()	0	2	1
A	EPSG_GEOCENTRIC_TRANSLATION_BY_GRID_INTERPOLATION_IGN()	0	2	1
A	EPSG_AFFINE_PARAMETRIC_TRANSFORMATION()	0	2	1
A	EPSG_POINT_MOTION_GEOCEN_BY_GRID_INADEFORM()	0	2	1
A	EPSG_LAMBERT_CYLINDRICAL_EQUAL_AREA()	0	2	1
A	EPSG_GNTRANS()	0	2	1
A	EPSG_LAMBERT_CONIC_CONFORMAL_WEST_ORIENTATED()	0	2	1
A	EPSG_DEGREE_REPRESENTATION_CONVERSION_DM_TO_DMSH()	0	2	1
A	geographicCoordinateFrameRotation()	0	24	1
A	EPSG_HEIGHT_DEPTH_REVERSAL()	0	2	1
A	EPSG_NORWAY_OFFSHORE_INTERPOLATION()	0	2	1
A	EPSG_POLAR_STEREOGRAPHIC_VARIANT_B()	0	2	1
A	EPSG_NADCON5_2D()	0	2	1
A	EPSG_TIME_DEPENDENT_COORDINATE_FRAME_ROTATION_GEOG2D()	0	2	1
A	EPSG_PSEUDO_PLATE_CARREE()	0	2	1
A	EPSG_SWISS_OBLIQUE_CYLINDRICAL()	0	2	1
A	EPSG_COMPLEX_POLYNOMIAL_OF_DEGREE_4()	0	2	1
A	EPSG_GEOCENTRIC_TRANSLATIONS_GEOCENTRIC_DOMAIN()	0	2	1
A	geographicGeocentricConversion()	0	41	4
A	EPSG_GEOCENTRIC_TRANSLATIONS_GEOG3D_DOMAIN()	0	2	1
A	EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_CGG2013()	0	2	1
A	EPSG_GEOGRAPHIC3D_OFFSETS()	0	2	1
A	EPSG_POLAR_STEREOGRAPHIC_VARIANT_A()	0	2	1
A	EPSG_TRANSVERSE_MERCATOR()	0	2	1
A	EPSG_GEOGRAPHIC2D_WITH_HEIGHT_OFFSETS()	0	2	1
A	EPSG_GEOGRAPHIC2D_OFFSETS()	0	2	1
A	EPSG_CASSINI_SOLDNER()	0	2	1
A	EPSG_TIME_SPECIFIC_COORDINATE_FRAME_ROTATION_GEOCEN()	0	2	1
A	geographicMolodensky()	0	36	3
A	EPSG_POINT_MOTION_BY_GRID_CANADA_NTV2_VEL()	0	2	1
A	EPSG_GENERAL_POLYNOMIAL_OF_DEGREE_2()	0	2	1
A	EPSG_SIMILARITY_TRANSFORMATION()	0	2	1
A	EPSG_GENERAL_POLYNOMIAL_OF_DEGREE_3()	0	2	1
A	EPSG_DEGREE_REPRESENTATION_CONVERSION_HDMS_TO_DMSH()	0	2	1
A	EPSG_KROVAK_MODIFIED_NORTH_ORIENTATED()	0	2	1
A	EPSG_DEGREE_REPRESENTATION_CONVERSION_DEG_TO_DMSH()	0	2	1
A	EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_IGN2009()	0	2	1
A	EPSG_AXIS_ORDER_REVERSAL_GEOGRAPHIC3D_HORIZONTAL()	0	2	1
A	EPSG_LAMBERT_CONIC_NEAR_CONFORMAL()	0	2	1
A	EPSG_LAMBERT_CONIC_CONFORMAL_2SP()	0	2	1
A	EPSG_DEGREE_REPRESENTATION_CONVERSION_DEGH_TO_DMSH()	0	2	1
A	EPSG_POPULAR_VISUALISATION_PSEUDO_MERCATOR()	0	2	1
A	EPSG_TRANSVERSE_MERCATOR_SOUTH_ORIENTATED()	0	2	1
A	EPSG_KROVAK_MODIFIED()	0	2	1
A	EPSG_NEW_ZEALAND_DEFORMATION_MODEL()	0	2	1
A	EPSG_REVERSIBLE_POLYNOMIAL_OF_DEGREE_4()	0	2	1
A	EPSG_GEOGRAPHIC_TOPOCENTRIC_CONVERSIONS()	0	2	1
A	EPSG_VERTICAL_OFFSET_AND_SLOPE()	0	2	1
A	EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_AUSGEOID98()	0	2	1
A	EPSG_AXIS_ORDER_REVERSAL_2D()	0	2	1
A	EPSG_VERTICAL_OFFSET_BY_GRID_INTERPOLATION_VERTCON()	0	2	1
A	EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_NZGEOID()	0	2	1
A	EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_PNG()	0	2	1
A	EPSG_HOTINE_OBLIQUE_MERCATOR_VARIANT_A()	0	2	1
A	EPSG_POLAR_STEREOGRAPHIC_VARIANT_C()	0	2	1
A	EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_SA_2010()	0	2	1
A	EPSG_TIME_DEPENDENT_POSITION_VECTOR_TFM_GEOG2D()	0	2	1
A	EPSG_OBLIQUE_STEREOGRAPHIC()	0	2	1
A	EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_BEV_AT()	0	2	1
A	EPSG_HOTINE_OBLIQUE_MERCATOR_VARIANT_B()	0	2	1
A	EPSG_AMERICAN_POLYCONIC()	0	2	1
A	EPSG_EQUIDISTANT_CYLINDRICAL()	0	2	1
A	EPSG_GENERAL_POLYNOMIAL_OF_DEGREE_6()	0	2	1
A	EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_OSGM15_IRE()	0	2	1
A	EPSG_MARITIME_PROVINCES_POLYNOMIAL_INTERPOLATION()	0	2	1
A	geocentricPositionVectorMolodenskyBadekas()	0	33	1
A	EPSG_NADCON5_3D()	0	2	1
A	EPSG_ORTHOGRAPHIC()	0	2	1
A	EPSG_CARTESIAN_GRID_OFFSETS()	0	2	1
A	EPSG_REVERSIBLE_POLYNOMIAL_OF_DEGREE_3()	0	2	1
A	EPSG_COLOMBIA_URBAN()	0	2	1
A	geographicPositionVectorTransformation()	0	24	1
A	EPSG_MODIFIED_AZIMUTHAL_EQUIDISTANT()	0	2	1
A	EPSG_MERCATOR_VARIANT_C()	0	2	1
A	EPSG_LAMBERT_CYLINDRICAL_EQUAL_AREA_SPHERICAL()	0	2	1
A	EPSG_KROVAK()	0	2	1
A	EPSG_TRANSVERSE_MERCATOR_ZONED_GRID_SYSTEM()	0	2	1

<?php

/**
 * PHPCoord.
 *
 * @author Doug Wright
 */
declare(strict_types=1);

namespace PHPCoord\CoordinateOperation;

use InvalidArgumentException;
use PHPCoord\CoordinateReferenceSystem\CoordinateReferenceSystem;
use PHPCoord\CoordinateReferenceSystem\Geocentric;
use PHPCoord\CoordinateReferenceSystem\Geographic2D;
use PHPCoord\CoordinateReferenceSystem\Geographic3D;
use PHPCoord\Datum\Ellipsoid;
use PHPCoord\GeocentricPoint;
use PHPCoord\GeographicPoint;
use PHPCoord\Point;
use PHPCoord\UnitOfMeasure\Angle\Angle;
use PHPCoord\UnitOfMeasure\Angle\ArcSecond;
use PHPCoord\UnitOfMeasure\Angle\Radian;
use PHPCoord\UnitOfMeasure\Length\Length;
use PHPCoord\UnitOfMeasure\Length\Metre;
use PHPCoord\UnitOfMeasure\Scale\Scale;

class CoordinateOperationMethods
{
    /**
     * Abridged Molodensky
     * This transformation is a truncated Taylor series expansion of a transformation between two geographic coordinate
     * systems, modelled as a set of geocentric translations.
     */
    public const EPSG_ABRIDGED_MOLODENSKY = 9605;

    /**
     * Affine geometric transformation.
     */
    public const EPSG_AFFINE_GEOMETRIC_TRANSFORMATION = 9623;

    /**
     * Affine parametric transformation.
     */
    public const EPSG_AFFINE_PARAMETRIC_TRANSFORMATION = 9624;

    /**
     * Albers Equal Area.
     */
    public const EPSG_ALBERS_EQUAL_AREA = 9822;

    /**
     * American Polyconic
     * See information source for formula and example.
     */
    public const EPSG_AMERICAN_POLYCONIC = 9818;

    /**
     * Axis Order Reversal (2D)
     * This is a parameter-less conversion to reverse the order of the axes of a 2D CRS.
     */
    public const EPSG_AXIS_ORDER_REVERSAL_2D = 9843;

    /**
     * Axis Order Reversal (Geographic3D horizontal)
     * This is a parameter-less conversion to change the order of horizontal coordinates of a geographic 3D CRS.
     */
    public const EPSG_AXIS_ORDER_REVERSAL_GEOGRAPHIC3D_HORIZONTAL = 9844;

    /**
     * Bonne.
     */
    public const EPSG_BONNE = 9827;

    /**
     * Bonne (South Orientated).
     */
    public const EPSG_BONNE_SOUTH_ORIENTATED = 9828;

    /**
     * Cartesian Grid Offsets
     * This transformation allows calculation of coordinates in the target system by adding the parameter value to the
     * coordinate values of the point in the source system.
     */
    public const EPSG_CARTESIAN_GRID_OFFSETS = 9656;

    /**
     * Cartesian Grid Offsets from Form Function
     * Used in German state of Schleswig-Holstein.
     */
    public const EPSG_CARTESIAN_GRID_OFFSETS_FROM_FORM_FUNCTION = 1036;

    /**
     * Cassini-Soldner.
     */
    public const EPSG_CASSINI_SOLDNER = 9806;

    /**
     * Change of Vertical Unit.
     */
    public const EPSG_CHANGE_OF_VERTICAL_UNIT = 1069;

    /**
     * Colombia Urban.
     */
    public const EPSG_COLOMBIA_URBAN = 1052;

    /**
     * Complex polynomial of degree 3
     * Coordinate pairs treated as complex numbers.  This exploits the correlation between the polynomial coefficients
     * and leads to a smaller number of coefficients than the general polynomial of degree 3.
     */
    public const EPSG_COMPLEX_POLYNOMIAL_OF_DEGREE_3 = 9652;

    /**
     * Complex polynomial of degree 4
     * Coordinate pairs treated as complex numbers.  This exploits the correlation between the polynomial coefficients
     * and leads to a smaller number of coefficients than the general polynomial of degree 4.
     */
    public const EPSG_COMPLEX_POLYNOMIAL_OF_DEGREE_4 = 9653;

    /**
     * Coordinate Frame rotation (geocentric domain)
     * This method is a specific case of the Molodensky-Badekas (CF) method (code 1034) in which the evaluation point
     * is at the geocentre with coordinate values of zero. Note the analogy with the Position Vector method (code 1033)
     * but beware of the differences!
     */
    public const EPSG_COORDINATE_FRAME_ROTATION_GEOCENTRIC_DOMAIN = 1032;

    /**
     * Coordinate Frame rotation (geog2D domain)
     * Note the analogy with the Position Vector tfm (code 9606) but beware of the differences!  The Position Vector
     * convention is used by IAG and recommended by ISO 19111. See methods 1032 and 1038 for similar tfms operating
     * between other CRS types.
     */
    public const EPSG_COORDINATE_FRAME_ROTATION_GEOG2D_DOMAIN = 9607;

    /**
     * Coordinate Frame rotation (geog3D domain)
     * Note the analogy with the Position Vector tfm (code 1037) but beware of the differences!  The Position Vector
     * convention is used by IAG and recommended by ISO 19111. See methods 1032 and 9607 for similar tfms operating
     * between other CRS types.
     */
    public const EPSG_COORDINATE_FRAME_ROTATION_GEOG3D_DOMAIN = 1038;

    /**
     * Degree representation conversion: DM to DMSH
     * Applies to 2D and the horizontal component of 3D ellipsoidal  systems.
     * @deprecated
     */
    public const EPSG_DEGREE_REPRESENTATION_CONVERSION_DM_TO_DMSH = 9640;

    /**
     * Degree representation conversion: DMH to DMSH
     * Applies to 2D and the horizontal component of 3D ellipsoidal  systems.
     * @deprecated
     */
    public const EPSG_DEGREE_REPRESENTATION_CONVERSION_DMH_TO_DMSH = 9641;

    /**
     * Degree representation conversion: DMS to DMSH
     * Applies to 2D and the horizontal component of 3D ellipsoidal  systems.
     * @deprecated
     */
    public const EPSG_DEGREE_REPRESENTATION_CONVERSION_DMS_TO_DMSH = 9643;

    /**
     * Degree representation conversion: HDM to DMSH
     * Applies to 2D and the horizontal component of 3D ellipsoidal  systems.
     * @deprecated
     */
    public const EPSG_DEGREE_REPRESENTATION_CONVERSION_HDM_TO_DMSH = 9642;

    /**
     * Degree representation conversion: HDMS to DMSH
     * Applies to 2D and the horizontal component of 3D ellipsoidal  systems.
     * @deprecated
     */
    public const EPSG_DEGREE_REPRESENTATION_CONVERSION_HDMS_TO_DMSH = 9644;

    /**
     * Degree representation conversion: Hdeg to DMSH
     * Applies to 2D and the horizontal component of 3D ellipsoidal  systems.
     * @deprecated
     */
    public const EPSG_DEGREE_REPRESENTATION_CONVERSION_HDEG_TO_DMSH = 9639;

    /**
     * Degree representation conversion: deg to DMSH
     * Applies to 2D and the horizontal component of 3D ellipsoidal  systems.
     * @deprecated
     */
    public const EPSG_DEGREE_REPRESENTATION_CONVERSION_DEG_TO_DMSH = 9637;

    /**
     * Degree representation conversion: degH to DMSH
     * Applies to 2D and the horizontal component of 3D ellipsoidal  systems.
     * @deprecated
     */
    public const EPSG_DEGREE_REPRESENTATION_CONVERSION_DEGH_TO_DMSH = 9638;

    /**
     * Equal Earth.
     */
    public const EPSG_EQUAL_EARTH = 1078;

    /**
     * Equidistant Cylindrical
     * See method code 1029 for spherical development. See also Pseudo Plate Carree, method code 9825.
     */
    public const EPSG_EQUIDISTANT_CYLINDRICAL = 1028;

    /**
     * Equidistant Cylindrical (Spherical)
     * See method code 1028 for ellipsoidal development. If the latitude of natural origin is at the equator, also
     * known as Plate Carrée. See also Pseudo Plate Carree, method code 9825.
     */
    public const EPSG_EQUIDISTANT_CYLINDRICAL_SPHERICAL = 1029;

    /**
     * GNTRANS.
     */
    public const EPSG_GNTRANS = 1040;

    /**
     * General polynomial of degree 2.
     */
    public const EPSG_GENERAL_POLYNOMIAL_OF_DEGREE_2 = 9645;

    /**
     * General polynomial of degree 3.
     */
    public const EPSG_GENERAL_POLYNOMIAL_OF_DEGREE_3 = 9646;

    /**
     * General polynomial of degree 4.
     */
    public const EPSG_GENERAL_POLYNOMIAL_OF_DEGREE_4 = 9647;

    /**
     * General polynomial of degree 6.
     */
    public const EPSG_GENERAL_POLYNOMIAL_OF_DEGREE_6 = 9648;

    /**
     * Geocentric translation by Grid Interpolation (IGN).
     */
    public const EPSG_GEOCENTRIC_TRANSLATION_BY_GRID_INTERPOLATION_IGN = 1087;

    /**
     * Geocentric translations  (geog3D domain)
     * See methods 1031 and 9603 for similar tfms operating between other CRSs types.
     */
    public const EPSG_GEOCENTRIC_TRANSLATIONS_GEOG3D_DOMAIN = 1035;

    /**
     * Geocentric translations (geocentric domain)
     * This method allows calculation of geocentric coords in the target system by adding the parameter values to the
     * corresponding coordinates of the point in the source system. See methods 1035 and 9603 for similar tfms
     * operating between other CRSs types.
     */
    public const EPSG_GEOCENTRIC_TRANSLATIONS_GEOCENTRIC_DOMAIN = 1031;

    /**
     * Geocentric translations (geog2D domain)
     * See methods 1031 and 1035 for similar tfms operating between other CRSs types.
     */
    public const EPSG_GEOCENTRIC_TRANSLATIONS_GEOG2D_DOMAIN = 9603;

    /**
     * Geocentric/topocentric conversions.
     */
    public const EPSG_GEOCENTRIC_TOPOCENTRIC_CONVERSIONS = 9836;

    /**
     * Geog3D to Geog2D+GravityRelatedHeight (US .gtx)
     * Transformation from a Geographic 3D CRS to a Compound CRS consisting of a Geographic 2D CRS and a Vertical CRS,
     * or vice versa. The Geographic 3D and the Geographic 2D CRS must be based on the same Geodetic Datum.
     */
    public const EPSG_GEOG3D_TO_GEOG2D_PLUS_GRAVITYRELATEDHEIGHT_US_GTX = 9635;

    /**
     * Geog3D to Geog2D+Vertical (AUSGeoid v2)
     * Transformation from a Geographic 3D CRS to a Compound CRS consisting of a Geographic 2D CRS and a Vertical CRS,
     * or vice versa. The Geographic 3D CRS and the Geographic 2D CRS must be based on the same Geodetic Datum.
     */
    public const EPSG_GEOG3D_TO_GEOG2D_PLUS_VERTICAL_AUSGEOID_V2 = 1083;

    /**
     * Geographic/geocentric conversions
     * This is a parameter-less conversion. In applications it is often concatenated with the 3- 7- or 10-parameter
     * transformations 9603, 9606, 9607 or 9636 to form a geographic to geographic transformation.
     */
    public const EPSG_GEOGRAPHIC_GEOCENTRIC_CONVERSIONS = 9602;

    /**
     * Geographic/topocentric conversions.
     */
    public const EPSG_GEOGRAPHIC_TOPOCENTRIC_CONVERSIONS = 9837;

    /**
     * Geographic2D offsets
     * This transformation allows calculation of coordinates in the target system by adding the parameter value to the
     * coordinate values of the point in the source system.
     */
    public const EPSG_GEOGRAPHIC2D_OFFSETS = 9619;

    /**
     * Geographic2D with Height Offsets
     * This transformation allows calculation of coordinates in the target system by adding the parameter value to the
     * coordinate values of the point in the source system.
     */
    public const EPSG_GEOGRAPHIC2D_WITH_HEIGHT_OFFSETS = 9618;

    /**
     * Geographic3D offsets
     * This transformation allows calculation of coordinates in the target system by adding the parameter value to the
     * coordinate values of the point in the source system.
     */
    public const EPSG_GEOGRAPHIC3D_OFFSETS = 9660;

    /**
     * Geographic3D to 2D conversion
     * This is a parameter-less conversion.
     */
    public const EPSG_GEOGRAPHIC3D_TO_2D_CONVERSION = 9659;

    /**
     * Geographic3D to GravityRelatedHeight (AUSGeoid v2)
     * The Information Source references software which offers both bi-cubic and bi-linear interpolation methods.
     * Unlike earlier AUSGeoid98 method which used bi-linear interpolation, Ausgeoid v2 uses bi-cubic. See Info Source
     * for file format documentation.
     */
    public const EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_AUSGEOID_V2 = 1048;

    /**
     * Geographic3D to GravityRelatedHeight (AUSGeoid98).
     */
    public const EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_AUSGEOID98 = 9662;

    /**
     * Geographic3D to GravityRelatedHeight (BEV AT).
     */
    public const EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_BEV_AT = 1081;

    /**
     * Geographic3D to GravityRelatedHeight (CGG2013)
     * For consistency with earlier models the Information Source references software which uses bi-quadratic
     * interpolation. Because of denser node spacing in CGG2013 bi-linear interpolation will be sufficient for most
     * uses. See Info Source for file format doc.
     */
    public const EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_CGG2013 = 1060;

    /**
     * Geographic3D to GravityRelatedHeight (CI).
     */
    public const EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_CI = 1050;

    /**
     * Geographic3D to GravityRelatedHeight (EGM)
     * Applies to EGM84 and EGM96 models. For later model see Geographic3D to GravityRelatedHeight (EGM2008), method
     * code 1025.
     */
    public const EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_EGM = 9661;

    /**
     * Geographic3D to GravityRelatedHeight (EGM2008)
     * For earlier EGM84 and EGM96 models see Geographic3D to GravityRelatedHeight (EGM), method code 9661.
     */
    public const EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_EGM2008 = 1025;

    /**
     * Geographic3D to GravityRelatedHeight (Gravsoft).
     */
    public const EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_GRAVSOFT = 1047;

    /**
     * Geographic3D to GravityRelatedHeight (IGN1997)
     * Transformation of the vertical component of a Geographic 3D CRS to a Vertical CRS.
     */
    public const EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_IGN1997 = 9664;

    /**
     * Geographic3D to GravityRelatedHeight (IGN2009)
     * Transformation of the vertical component of a Geographic 3D CRS to a Vertical CRS. File header of method code
     * 9664  (4 lines) has changed in this method (1 line);  recommended interpolation method now in a separate XML
     * file with same name as the grid.
     */
    public const EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_IGN2009 = 1073;

    /**
     * Geographic3D to GravityRelatedHeight (NZgeoid)
     * EPSG initially gave this method the name "Geographic3D to GravityRelatedHeight (NZgeoid2009)". As the same file
     * format was retained for the 2016 geoid, date removed from the method name.
     */
    public const EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_NZGEOID = 1030;

    /**
     * Geographic3D to GravityRelatedHeight (OSGM-GB)
     * Transformation of the vertical component of a Geographic 3D CRS to a Vertical CRS.
     */
    public const EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_OSGM_GB = 9663;

    /**
     * Geographic3D to GravityRelatedHeight (OSGM02-Ire)
     * Transformation of the vertical component of a Geographic 3D CRS to a Vertical CRS.
     */
    public const EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_OSGM02_IRE = 1045;

    /**
     * Geographic3D to GravityRelatedHeight (OSGM15-Ire)
     * Transformation of the vertical component of a Geographic 3D CRS to a Vertical CRS.
     */
    public const EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_OSGM15_IRE = 1072;

    /**
     * Geographic3D to GravityRelatedHeight (PNG)
     * Transformation of the vertical component of a Geographic 3D CRS to a Vertical CRS.
     */
    public const EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_PNG = 1059;

    /**
     * Geographic3D to GravityRelatedHeight (SA 2010)
     * File format: space-separated ascii file, no header, columns: latitude, longitude, separation.
     */
    public const EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_SA_2010 = 1082;

    /**
     * Geographic3D to GravityRelatedHeight (gtx)
     * Transformation of the vertical component of a Geographic 3D CRS to a Vertical CRS. Grid file format: US NGS .gtx
     * (in US sometimes also referred to as 'vdatum format').
     */
    public const EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_GTX = 9665;

    /**
     * Guam Projection
     * Simplified form of Oblique Azimuthal Equidistant projection method.
     */
    public const EPSG_GUAM_PROJECTION = 9831;

    /**
     * Height Depth Reversal
     * This is a parameter-less conversion.
     */
    public const EPSG_HEIGHT_DEPTH_REVERSAL = 1068;

    /**
     * Hotine Oblique Mercator (variant A).
     */
    public const EPSG_HOTINE_OBLIQUE_MERCATOR_VARIANT_A = 9812;

    /**
     * Hotine Oblique Mercator (variant B).
     */
    public const EPSG_HOTINE_OBLIQUE_MERCATOR_VARIANT_B = 9815;

    /**
     * Hyperbolic Cassini-Soldner.
     */
    public const EPSG_HYPERBOLIC_CASSINI_SOLDNER = 9833;

    /**
     * Krovak.
     */
    public const EPSG_KROVAK = 9819;

    /**
     * Krovak (North Orientated).
     */
    public const EPSG_KROVAK_NORTH_ORIENTATED = 1041;

    /**
     * Krovak Modified
     * Incorporates a polynomial transformation which is defined to be exact and for practical purposes is considered
     * to be a map projection.
     */
    public const EPSG_KROVAK_MODIFIED = 1042;

    /**
     * Krovak Modified (North Orientated)
     * Incorporates a polynomial transformation which is defined to be exact and for practical purposes is considered
     * to be a map projection.
     */
    public const EPSG_KROVAK_MODIFIED_NORTH_ORIENTATED = 1043;

    /**
     * Laborde Oblique Mercator.
     */
    public const EPSG_LABORDE_OBLIQUE_MERCATOR = 9813;

    /**
     * Lambert Azimuthal Equal Area
     * This is the ellipsoidal form of the projection.
     */
    public const EPSG_LAMBERT_AZIMUTHAL_EQUAL_AREA = 9820;

    /**
     * Lambert Azimuthal Equal Area (Spherical)
     * This is the spherical form of the projection.  See coordinate operation method Lambert Azimuthal Equal Area
     * (code 9820) for ellipsoidal form.  Differences of several tens of metres result from comparison of the two
     * methods.
     */
    public const EPSG_LAMBERT_AZIMUTHAL_EQUAL_AREA_SPHERICAL = 1027;

    /**
     * Lambert Conic Conformal (1SP).
     */
    public const EPSG_LAMBERT_CONIC_CONFORMAL_1SP = 9801;

    /**
     * Lambert Conic Conformal (2SP Belgium)
     * In 2000 this modification was replaced through use of the regular Lambert Conic Conformal (2SP) method [9802]
     * with appropriately modified parameter values.
     */
    public const EPSG_LAMBERT_CONIC_CONFORMAL_2SP_BELGIUM = 9803;

    /**
     * Lambert Conic Conformal (2SP Michigan).
     */
    public const EPSG_LAMBERT_CONIC_CONFORMAL_2SP_MICHIGAN = 1051;

    /**
     * Lambert Conic Conformal (2SP).
     */
    public const EPSG_LAMBERT_CONIC_CONFORMAL_2SP = 9802;

    /**
     * Lambert Conic Conformal (West Orientated).
     */
    public const EPSG_LAMBERT_CONIC_CONFORMAL_WEST_ORIENTATED = 9826;

    /**
     * Lambert Conic Near-Conformal
     * The Lambert Near-Conformal projection is derived from the Lambert Conformal Conic projection by truncating the
     * series expansion of the projection formulae.
     */
    public const EPSG_LAMBERT_CONIC_NEAR_CONFORMAL = 9817;

    /**
     * Lambert Cylindrical Equal Area
     * This is the ellipsoidal form of the projection.
     */
    public const EPSG_LAMBERT_CYLINDRICAL_EQUAL_AREA = 9835;

    /**
     * Lambert Cylindrical Equal Area (Spherical)
     * This is the spherical form of the projection.  See coordinate operation method Lambert Cylindrical Equal Area
     * (code 9835) for ellipsoidal form.  Differences of several tens of metres result from comparison of the two
     * methods.
     */
    public const EPSG_LAMBERT_CYLINDRICAL_EQUAL_AREA_SPHERICAL = 9834;

    /**
     * Longitude rotation
     * This transformation allows calculation of the longitude of a point in the target system by adding the parameter
     * value to the longitude value of the point in the source system.
     */
    public const EPSG_LONGITUDE_ROTATION = 9601;

    /**
     * Madrid to ED50 polynomial.
     */
    public const EPSG_MADRID_TO_ED50_POLYNOMIAL = 9617;

    /**
     * Maritime Provinces polynomial interpolation
     * This transformation is an executable module within the application NBGeocalc.  It is an adaptation of the ESTPM
     * program developed by Geodetic Survey of Canada.
     */
    public const EPSG_MARITIME_PROVINCES_POLYNOMIAL_INTERPOLATION = 9634;

    /**
     * Mercator (Spherical).
     */
    public const EPSG_MERCATOR_SPHERICAL = 1026;

    /**
     * Mercator (variant A)
     * Note that in these formulas the parameter latitude of natural origin (latO) is not used. However for this
     * Mercator (variant A) method the EPSG dataset includes this parameter, which must have a value of zero, for
     * completeness in CRS labelling.
     */
    public const EPSG_MERCATOR_VARIANT_A = 9804;

    /**
     * Mercator (variant B)
     * Used for most nautical charts.
     */
    public const EPSG_MERCATOR_VARIANT_B = 9805;

    /**
     * Mercator (variant C).
     */
    public const EPSG_MERCATOR_VARIANT_C = 1044;

    /**
     * Modified Azimuthal Equidistant
     * Modified form of Oblique Azimuthal Equidistant projection method developed for Polynesian islands. For the
     * distances over which these projections are used (under 800km) this modification introduces no significant error.
     */
    public const EPSG_MODIFIED_AZIMUTHAL_EQUIDISTANT = 9832;

    /**
     * Molodensky
     * See Abridged Molodensky.
     */
    public const EPSG_MOLODENSKY = 9604;

    /**
     * Molodensky-Badekas (CF geocentric domain)
     * See method codes 1039 and 9636 for this operation in other coordinate domains and method code 1061 for opposite
     * rotation convention in geocentric domain.
     */
    public const EPSG_MOLODENSKY_BADEKAS_CF_GEOCENTRIC_DOMAIN = 1034;

    /**
     * Molodensky-Badekas (CF geog2D domain)
     * See method codes 1034 and 1039 for this operation in other coordinate domains and method code 1063 for the
     * opposite rotation convention in geographic 2D domain.
     */
    public const EPSG_MOLODENSKY_BADEKAS_CF_GEOG2D_DOMAIN = 9636;

    /**
     * Molodensky-Badekas (CF geog3D domain)
     * See method codes 1034 and 9636 for this operation in other coordinate domains and method code 1062 for opposite
     * rotation convention in geographic 3D domain.
     */
    public const EPSG_MOLODENSKY_BADEKAS_CF_GEOG3D_DOMAIN = 1039;

    /**
     * Molodensky-Badekas (PV geocentric domain)
     * See method codes 1062 and 1063 for this operation in other coordinate domains and method code 1034 for opposite
     * rotation convention in geocentric domain.
     */
    public const EPSG_MOLODENSKY_BADEKAS_PV_GEOCENTRIC_DOMAIN = 1061;

    /**
     * Molodensky-Badekas (PV geog2D domain)
     * See method codes 1061 and 1062 for this operation in other coordinate domains and method code 9636 for opposite
     * rotation in geographic 2D domain.
     */
    public const EPSG_MOLODENSKY_BADEKAS_PV_GEOG2D_DOMAIN = 1063;

    /**
     * Molodensky-Badekas (PV geog3D domain)
     * See method codes 1061 and 1063 for this operation in other coordinate domains and method code 1039 for opposite
     * rotation convention in geographic 3D domain.
     */
    public const EPSG_MOLODENSKY_BADEKAS_PV_GEOG3D_DOMAIN = 1062;

    /**
     * NADCON
     * Geodetic transformation operating on geographic coordinate differences by bi-linear interpolation.  Input
     * expects longitudes to be positive west.
     */
    public const EPSG_NADCON = 9613;

    /**
     * NADCON5 (2D)
     * Geodetic transformation operating on geographic coordinate differences by bi-quadratic interpolation.  Input
     * expects longitudes to be positive east in range 0-360° (0° = Greenwich).
     */
    public const EPSG_NADCON5_2D = 1074;

    /**
     * NADCON5 (3D)
     * Geodetic transformation operating on geographic coordinate differences by bi-quadratic interpolation.  Input
     * expects longitudes to be positive east in range 0-360° (0° = Greenwich).
     */
    public const EPSG_NADCON5_3D = 1075;

    /**
     * NTv1
     * Geodetic transformation operating on geographic coordinate differences by bi-linear interpolation. Superseded in
     * 1997 by NTv2 (transformation method code 9615).  Input expects longitudes to be positive west.
     */
    public const EPSG_NTV1 = 9614;

    /**
     * NTv2
     * Geodetic transformation operating on geographic coordinate differences by bi-linear interpolation.  Supersedes
     * NTv1 (transformation method code 9614).  Input expects longitudes to be positive west.
     */
    public const EPSG_NTV2 = 9615;

    /**
     * New Zealand Deformation Model
     * Zip final contains a set of CSV files documentation for a secular tectonic plate velocity model and patches
     * which account for movements due to specific deformation events (earthquakes).
     */
    public const EPSG_NEW_ZEALAND_DEFORMATION_MODEL = 1079;

    /**
     * New Zealand Map Grid.
     */
    public const EPSG_NEW_ZEALAND_MAP_GRID = 9811;

    /**
     * Norway Offshore Interpolation
     * Although in principle this method is not reversible, in practice reversibility is better than 10 cm. For the
     * applications for which it was designed it may be considered reversible.
     */
    public const EPSG_NORWAY_OFFSHORE_INTERPOLATION = 9620;

    /**
     * Oblique Stereographic
     * This is not the same as the projection method of the same name in USGS Professional Paper no. 1395, "Map
     * Projections - A Working Manual" by John P. Snyder.
     */
    public const EPSG_OBLIQUE_STEREOGRAPHIC = 9809;

    /**
     * Ordnance Survey National Transformation
     * Geodetic transformation between ETRS89 (or WGS 84) and OSGB36 / National Grid.  Uses ETRS89 / National Grid as
     * an intermediate coordinate system for bi-linear interpolation of gridded grid coordinate differences.
     */
    public const EPSG_ORDNANCE_SURVEY_NATIONAL_TRANSFORMATION = 9633;

    /**
     * Orthographic
     * If the natural origin of the projection is at the topocentric origin, this is a special case of the Vertical
     * Perspective (orthographic case) (method code 9839) in which the ellipsoid height of all mapped points is zero (h
     * = 0).
     */
    public const EPSG_ORTHOGRAPHIC = 9840;

    /**
     * Point motion (ellipsoidal).
     */
    public const EPSG_POINT_MOTION_ELLIPSOIDAL = 1067;

    /**
     * Point motion (geocen) by grid (INADEFORM).
     */
    public const EPSG_POINT_MOTION_GEOCEN_BY_GRID_INADEFORM = 1086;

    /**
     * Point motion (geocentric Cartesian).
     */
    public const EPSG_POINT_MOTION_GEOCENTRIC_CARTESIAN = 1064;

    /**
     * Point motion by grid (Canada NTv2_Vel)
     * Interpolation within the grid is in the horizontal component of the source CRS.
     */
    public const EPSG_POINT_MOTION_BY_GRID_CANADA_NTV2_VEL = 1070;

    /**
     * Polar Stereographic (variant A)
     * Latitude of natural origin must be either 90 degrees or -90 degrees (or equivalent in alternative angle unit).
     */
    public const EPSG_POLAR_STEREOGRAPHIC_VARIANT_A = 9810;

    /**
     * Polar Stereographic (variant B).
     */
    public const EPSG_POLAR_STEREOGRAPHIC_VARIANT_B = 9829;

    /**
     * Polar Stereographic (variant C).
     */
    public const EPSG_POLAR_STEREOGRAPHIC_VARIANT_C = 9830;

    /**
     * Popular Visualisation Pseudo Mercator
     * Applies spherical formulas to the ellipsoid. As such does not have the properties of a true Mercator projection.
     */
    public const EPSG_POPULAR_VISUALISATION_PSEUDO_MERCATOR = 1024;

    /**
     * Position Vector transformation (geocentric domain)
     * This method is a specific case of the Molodensky-Badekas (PV) method (code 1061) in which the evaluation point
     * is the geocentre with coordinate values of zero. Note the analogy with the Coordinate Frame method (code 1032)
     * but beware of the differences!
     */
    public const EPSG_POSITION_VECTOR_TRANSFORMATION_GEOCENTRIC_DOMAIN = 1033;

    /**
     * Position Vector transformation (geog2D domain)
     * Note the analogy with the Coordinate Frame rotation (code 9607) but beware of the differences!  The Position
     * Vector convention is used by IAG and recommended by ISO 19111. See methods 1033 and 1037 for similar tfms
     * operating between other CRS types.
     */
    public const EPSG_POSITION_VECTOR_TRANSFORMATION_GEOG2D_DOMAIN = 9606;

    /**
     * Position Vector transformation (geog3D domain)
     * Note the analogy with the Coordinate Frame rotation (code 1038) but beware of the differences!  The Position
     * Vector convention is used by IAG and recommended by ISO 19111. See methods 1033 and 9606 for similar tfms
     * operating between other CRS types.
     */
    public const EPSG_POSITION_VECTOR_TRANSFORMATION_GEOG3D_DOMAIN = 1037;

    /**
     * Pseudo Plate Carree
     * Used only for depiction of graticule (latitude/longitude) coordinates on a computer display. The axes units are
     * decimal degrees and of variable scale. The origin is at Lat = 0, Long = 0. See Equidistant Cylindrical, code
     * 1029, for proper Plate Carrée.
     */
    public const EPSG_PSEUDO_PLATE_CARREE = 9825;

    /**
     * Reversible polynomial of degree 13.
     */
    public const EPSG_REVERSIBLE_POLYNOMIAL_OF_DEGREE_13 = 9654;

    /**
     * Reversible polynomial of degree 2
     * Reversibility is subject to constraints.  See Guidance Note 7 for conditions and clarification.
     */
    public const EPSG_REVERSIBLE_POLYNOMIAL_OF_DEGREE_2 = 9649;

    /**
     * Reversible polynomial of degree 3
     * Reversibility is subject to constraints.  See Guidance Note 7 for conditions and clarification.
     */
    public const EPSG_REVERSIBLE_POLYNOMIAL_OF_DEGREE_3 = 9650;

    /**
     * Reversible polynomial of degree 4
     * Reversibility is subject to constraints.  See Guidance Note 7 for conditions and clarification.
     */
    public const EPSG_REVERSIBLE_POLYNOMIAL_OF_DEGREE_4 = 9651;

    /**
     * Similarity transformation
     * Defined for two-dimensional coordinate systems.
     */
    public const EPSG_SIMILARITY_TRANSFORMATION = 9621;

    /**
     * Swiss Oblique Cylindrical
     * Can be accommodated by Oblique Mercator method (code 9815), for which this method is an approximation (see BfL
     * document swissprojectionen.pdf at www.swisstopo.com).
     */
    public const EPSG_SWISS_OBLIQUE_CYLINDRICAL = 9814;

    /**
     * Time-dependent Coordinate Frame rotation (geocen)
     * Note the analogy with the Time-dependent Position Vector transformation (code 1053) but beware of the
     * differences!  The Position Vector convention is used by IAG. See method codes 1057 and 1058 for similar methods
     * operating between other CRS types.
     */
    public const EPSG_TIME_DEPENDENT_COORDINATE_FRAME_ROTATION_GEOCEN = 1056;

    /**
     * Time-dependent Coordinate Frame rotation (geog2D)
     * Note the analogy with the Time-dependent Position Vector transformation (code 1054) but beware of the
     * differences!  The Position Vector convention is used by IAG. See methods 1056 and 1058 for similar tfms
     * operating between other CRS types.
     */
    public const EPSG_TIME_DEPENDENT_COORDINATE_FRAME_ROTATION_GEOG2D = 1057;

    /**
     * Time-dependent Coordinate Frame rotation (geog3D)
     * Note the analogy with the Time-dependent Position Vector transformation (code 1055) but beware of the
     * differences!  The Position Vector convention is used by IAG. See method codes 1056 and 1057 for similar methods
     * operating between other CRS types.
     */
    public const EPSG_TIME_DEPENDENT_COORDINATE_FRAME_ROTATION_GEOG3D = 1058;

    /**
     * Time-dependent Position Vector tfm (geocentric)
     * Note the analogy with the Time-dependent Coordinate Frame rotation (code 1056) but beware of the differences!
     * The Position Vector convention is used by IAG. See method codes 1054 and 1055 for similar methods operating
     * between other CRS types.
     */
    public const EPSG_TIME_DEPENDENT_POSITION_VECTOR_TFM_GEOCENTRIC = 1053;

    /**
     * Time-dependent Position Vector tfm (geog2D)
     * Note the analogy with the Time-dependent Coordinate Frame rotation (code 1057) but beware of the differences!
     * The Position Vector convention is used by IAG. See method codes 1053 and 1055 for similar methods operating
     * between other CRS types.
     */
    public const EPSG_TIME_DEPENDENT_POSITION_VECTOR_TFM_GEOG2D = 1054;

    /**
     * Time-dependent Position Vector tfm (geog3D)
     * Note the analogy with the Coordinate Frame rotation (code 1058) but beware of the differences!  The Position
     * Vector convention is used by IAG. See method codes 1053 and 1054 for similar methods operating between other CRS
     * types.
     */
    public const EPSG_TIME_DEPENDENT_POSITION_VECTOR_TFM_GEOG3D = 1055;

    /**
     * Time-specific Coordinate Frame rotation (geocen)
     * Note the analogy with the Time-specific Position Vector method (code 1065) but beware of the differences!
     */
    public const EPSG_TIME_SPECIFIC_COORDINATE_FRAME_ROTATION_GEOCEN = 1066;

    /**
     * Time-specific Position Vector transform (geocen)
     * Note the analogy with the Time-specifc Coordinate Frame method (code 1066) but beware of the differences!
     */
    public const EPSG_TIME_SPECIFIC_POSITION_VECTOR_TRANSFORM_GEOCEN = 1065;

    /**
     * Transverse Mercator.
     */
    public const EPSG_TRANSVERSE_MERCATOR = 9807;

    /**
     * Transverse Mercator (South Orientated).
     */
    public const EPSG_TRANSVERSE_MERCATOR_SOUTH_ORIENTATED = 9808;

    /**
     * Transverse Mercator Zoned Grid System
     * If locations fall outwith the fixed zones the general Transverse Mercator method (code 9807) must be used for
     * each zone.
     */
    public const EPSG_TRANSVERSE_MERCATOR_ZONED_GRID_SYSTEM = 9824;

    /**
     * Vertical Offset
     * This transformation allows calculation of height (or depth) in the target system by adding the parameter value
     * to the height (or depth)-value of the point in the source system.
     */
    public const EPSG_VERTICAL_OFFSET = 9616;

    /**
     * Vertical Offset and Slope
     * This transformation allows calculation of height in the target system by applying the parameter values to the
     * height value of the point in the source system.
     */
    public const EPSG_VERTICAL_OFFSET_AND_SLOPE = 1046;

    /**
     * Vertical Offset by Grid Interpolation (BEV AT).
     */
    public const EPSG_VERTICAL_OFFSET_BY_GRID_INTERPOLATION_BEV_AT = 1080;

    /**
     * Vertical Offset by Grid Interpolation (NZLVD).
     */
    public const EPSG_VERTICAL_OFFSET_BY_GRID_INTERPOLATION_NZLVD = 1071;

    /**
     * Vertical Offset by Grid Interpolation (VERTCON)
     * Any NAD realization may be used as the Interpolation CRS; bi-linear interpolation is used. Input expects
     * longitudes to be positive west.
     */
    public const EPSG_VERTICAL_OFFSET_BY_GRID_INTERPOLATION_VERTCON = 9658;

    /**
     * Vertical Offset by Grid Interpolation (asc).
     */
    public const EPSG_VERTICAL_OFFSET_BY_GRID_INTERPOLATION_ASC = 1085;

    /**
     * Vertical Offset by Grid Interpolation (gtx).
     */
    public const EPSG_VERTICAL_OFFSET_BY_GRID_INTERPOLATION_GTX = 1084;

    /**
     * Vertical Perspective
     * For a viewing point height approaching or at infinity, see the Vertical Perspective (orthographic case) (method
     * code 9839).
     */
    public const EPSG_VERTICAL_PERSPECTIVE = 9838;

    /**
     * Vertical Perspective (Orthographic case)
     * This is a special case of the general Vertical Perspective (method code 9838) in which the viewing point at
     * infinity.
     */
    public const EPSG_VERTICAL_PERSPECTIVE_ORTHOGRAPHIC_CASE = 9839;

    public function operationByEPSGCode(int $epsgCode)

The parameter $epsgCode is not used and could be removed.

    {
        return [
            self::EPSG_ABRIDGED_MOLODENSKY => 'geographicAbridgedMolodensky',
            self::EPSG_MOLODENSKY => 'geographicMolodensky',
            self::EPSG_GEOGRAPHIC_GEOCENTRIC_CONVERSIONS => 'geographicGeocentricConversion',
            self::EPSG_COORDINATE_FRAME_ROTATION_GEOCENTRIC_DOMAIN => 'geocentricCoordinateFrameRotation',
            self::EPSG_COORDINATE_FRAME_ROTATION_GEOG2D_DOMAIN => 'geographicCoordinateFrameRotation',
            self::EPSG_COORDINATE_FRAME_ROTATION_GEOG3D_DOMAIN => 'geographicCoordinateFrameRotation',
            self::EPSG_POSITION_VECTOR_TRANSFORMATION_GEOCENTRIC_DOMAIN => 'geocentricPositionVectorTransformation',
            self::EPSG_POSITION_VECTOR_TRANSFORMATION_GEOG2D_DOMAIN => 'geographicPositionVectorTransformation',
            self::EPSG_POSITION_VECTOR_TRANSFORMATION_GEOG3D_DOMAIN => 'geographicPositionVectorTransformation',
            self::EPSG_MOLODENSKY_BADEKAS_CF_GEOCENTRIC_DOMAIN => 'geocentricCoordinateFrameMolodenskyBadekas',
            self::EPSG_MOLODENSKY_BADEKAS_CF_GEOG2D_DOMAIN => 'geographicCoordinateFrameMolodenskyBadekas',
            self::EPSG_MOLODENSKY_BADEKAS_CF_GEOG3D_DOMAIN => 'geographicCoordinateFrameMolodenskyBadekas',
            self::EPSG_MOLODENSKY_BADEKAS_PV_GEOCENTRIC_DOMAIN => 'geocentricPositionVectorMolodenskyBadekas',
            self::EPSG_MOLODENSKY_BADEKAS_PV_GEOG2D_DOMAIN => 'geographicPositionVectorMolodenskyBadekas',
            self::EPSG_MOLODENSKY_BADEKAS_PV_GEOG3D_DOMAIN => 'geographicPositionVectorMolodenskyBadekas',
        ];
    }

    /**
     * Abridged Molodensky
     * This transformation is a truncated Taylor series expansion of a transformation between two geographic coordinate
     * systems, modelled as a set of geocentric translations.
     */
    public function geographicAbridgedMolodensky(
        GeographicPoint $from,
        CoordinateReferenceSystem $to,
        Length $translationX,
        Length $translationY,
        Length $translationZ,
        Length $differenceInSemiMajorAxis,
        Scale $differenceInFlattening
    ): GeographicPoint {
        $fromLatitude = $from->getLatitude()->asRadians()->getValue();
        $fromLongitude = $from->getLongitude()->asRadians()->getValue() - $from->getCRS()->getDatum()->getPrimeMeridian()->getGreenwichLongitude()->asRadians()->getValue();
        $fromHeight = $from->getHeight() ? $from->getHeight()->asMetres()->getValue() : 0;
        $tx = $translationX->asMetres()->getValue();
        $ty = $translationY->asMetres()->getValue();
        $tz = $translationZ->asMetres()->getValue();
        $da = $differenceInSemiMajorAxis->asMetres()->getValue();
        $df = $differenceInFlattening->asUnity()->getValue();

        $a = $from->getCRS()->getDatum()->getEllipsoid()->getSemiMajorAxis()->getValue();
        $e2 = $from->getCRS()->getDatum()->getEllipsoid()->getEccentricitySquared();

        $rho = $a * (1 - $e2) / (1 - $e2 * sin($fromLatitude) ** 2) ** (3 / 2);
        $nu = $a / ((1 - $e2 * (sin($fromLatitude) ** 2)) ** 0.5);

        $f = $from->getCRS()->getDatum()->getEllipsoid()->getInverseFlattening();

        $dLatitude = ((-$tx * sin($fromLatitude) * cos($fromLongitude)) - ($ty * sin($fromLatitude) * sin($fromLongitude)) + ($tz * cos($fromLatitude)) + ((($a * $df) + ($from->getCRS()->getDatum()->getEllipsoid()->getInverseFlattening() * $da)) * sin(2 * $fromLatitude))) / ($rho * sin((new ArcSecond(1))->asRadians()->getValue()));
        $dLongitude = (-$tx * sin($fromLongitude) + $ty * cos($fromLongitude)) / (($nu * cos($fromLatitude)) * sin((new ArcSecond(1))->asRadians()->getValue()));
        $dHeight = ($tx * cos($fromLatitude) * cos($fromLongitude)) + ($ty * cos($fromLatitude) * sin($fromLongitude)) + ($tz * sin($fromLatitude)) + (($a * $df + $f * $da) * (sin($fromLatitude) ** 2)) - $da;

        $toLatitude = $fromLatitude + (new ArcSecond($dLatitude))->asRadians()->getValue();
        $toLongitude = $fromLongitude + (new ArcSecond($dLongitude))->asRadians()->getValue();
        $toHeight = $fromHeight + $dHeight;

        return new GeographicPoint(new Radian($toLatitude), new Radian($toLongitude), $to instanceof Geographic3D ? new Metre($toHeight) : null, $to);
    }

    /**
     * Affine geometric transformation.
     */
    public function EPSG_AFFINE_GEOMETRIC_TRANSFORMATION(): Point
    {
    }

    /**
     * Affine parametric transformation.
     */
    public function EPSG_AFFINE_PARAMETRIC_TRANSFORMATION(): Point
    {
    }

    /**
     * Albers Equal Area.
     */
    public function EPSG_ALBERS_EQUAL_AREA(): Point
    {
    }

    /**
     * American Polyconic
     * See information source for formula and example.
     */
    public function EPSG_AMERICAN_POLYCONIC(): Point
    {
    }

    /**
     * Axis Order Reversal (2D)
     * This is a parameter-less conversion to reverse the order of the axes of a 2D CRS.
     */
    public function EPSG_AXIS_ORDER_REVERSAL_2D(): Point
    {
    }

    /**
     * Axis Order Reversal (Geographic3D horizontal)
     * This is a parameter-less conversion to change the order of horizontal coordinates of a geographic 3D CRS.
     */
    public function EPSG_AXIS_ORDER_REVERSAL_GEOGRAPHIC3D_HORIZONTAL(): Point
    {
    }

    /**
     * Bonne.
     */
    public function EPSG_BONNE(): Point
    {
    }

    /**
     * Bonne (South Orientated).
     */
    public function EPSG_BONNE_SOUTH_ORIENTATED(): Point
    {
    }

    /**
     * Cartesian Grid Offsets
     * This transformation allows calculation of coordinates in the target system by adding the parameter value to the
     * coordinate values of the point in the source system.
     */
    public function EPSG_CARTESIAN_GRID_OFFSETS(): Point
    {
    }

    /**
     * Cartesian Grid Offsets from Form Function
     * Used in German state of Schleswig-Holstein.
     */
    public function EPSG_CARTESIAN_GRID_OFFSETS_FROM_FORM_FUNCTION(): Point
    {
    }

    /**
     * Cassini-Soldner.
     */
    public function EPSG_CASSINI_SOLDNER(): Point
    {
    }

    /**
     * Change of Vertical Unit.
     */
    public function EPSG_CHANGE_OF_VERTICAL_UNIT(): Point
    {
    }

    /**
     * Colombia Urban.
     */
    public function EPSG_COLOMBIA_URBAN(): Point
    {
    }

    /**
     * Complex polynomial of degree 3
     * Coordinate pairs treated as complex numbers.  This exploits the correlation between the polynomial coefficients
     * and leads to a smaller number of coefficients than the general polynomial of degree 3.
     */
    public function EPSG_COMPLEX_POLYNOMIAL_OF_DEGREE_3(): Point
    {
    }

    /**
     * Complex polynomial of degree 4
     * Coordinate pairs treated as complex numbers.  This exploits the correlation between the polynomial coefficients
     * and leads to a smaller number of coefficients than the general polynomial of degree 4.
     */
    public function EPSG_COMPLEX_POLYNOMIAL_OF_DEGREE_4(): Point
    {
    }

    /**
     * Coordinate Frame rotation (geocentric domain)
     * This method is a specific case of the Molodensky-Badekas (CF) method (code 1034) in which the evaluation point
     * is at the geocentre with coordinate values of zero. Note the analogy with the Position Vector method (code 1033)
     * but beware of the differences!
     */
    public function geocentricCoordinateFrameRotation(
        GeocentricPoint $from,
        Geocentric $to,
        Length $translationX,
        Length $translationY,
        Length $translationZ,
        Angle $rotationX,
        Angle $rotationY,
        Angle $rotationZ,
        Scale $scaleFactor
    ): GeocentricPoint {
        return $this->geocentricCoordinateFrameMolodenskyBadekas(
            $from,
            $to,
            $translationX,
            $translationY,
            $translationZ,
            $rotationX,
            $rotationY,
            $rotationZ,
            $scaleFactor,
            new Metre(0),
            new Metre(0),
            new Metre(0)
        );
    }

    /**
     * Coordinate Frame rotation (geog2D/geog3D domain)
     * Note the analogy with the Position Vector tfm (codes 9606/1037) but beware of the differences!  The Position Vector
     * convention is used by IAG and recommended by ISO 19111. See methods 1032/1038/9607 for similar tfms operating
     * between other CRS types.
     * @param Geographic2D|Geographic3D $to
     */
    public function geographicCoordinateFrameRotation(
        GeographicPoint $from,
        CoordinateReferenceSystem $to,
        Length $translationX,
        Length $translationY,
        Length $translationZ,
        Angle $rotationX,
        Angle $rotationY,
        Angle $rotationZ,
        Scale $scaleFactor
    ): GeographicPoint {
        return $this->geographicCoordinateFrameMolodenskyBadekas(
            $from,
            $to,
            $translationX,
            $translationY,
            $translationZ,
            $rotationX,
            $rotationY,
            $rotationZ,
            $scaleFactor,
            new Metre(0),
            new Metre(0),
            new Metre(0)
        );
    }

    /**
     * Degree representation conversion: DM to DMSH
     * Applies to 2D and the horizontal component of 3D ellipsoidal  systems.
     * @deprecated
     */
    public function EPSG_DEGREE_REPRESENTATION_CONVERSION_DM_TO_DMSH(): Point
    {
    }

    /**
     * Degree representation conversion: DMH to DMSH
     * Applies to 2D and the horizontal component of 3D ellipsoidal  systems.
     * @deprecated
     */
    public function EPSG_DEGREE_REPRESENTATION_CONVERSION_DMH_TO_DMSH(): Point
    {
    }

    /**
     * Degree representation conversion: DMS to DMSH
     * Applies to 2D and the horizontal component of 3D ellipsoidal  systems.
     * @deprecated
     */
    public function EPSG_DEGREE_REPRESENTATION_CONVERSION_DMS_TO_DMSH(): Point
    {
    }

    /**
     * Degree representation conversion: HDM to DMSH
     * Applies to 2D and the horizontal component of 3D ellipsoidal  systems.
     * @deprecated
     */
    public function EPSG_DEGREE_REPRESENTATION_CONVERSION_HDM_TO_DMSH(): Point
    {
    }

    /**
     * Degree representation conversion: HDMS to DMSH
     * Applies to 2D and the horizontal component of 3D ellipsoidal  systems.
     * @deprecated
     */
    public function EPSG_DEGREE_REPRESENTATION_CONVERSION_HDMS_TO_DMSH(): Point
    {
    }

    /**
     * Degree representation conversion: Hdeg to DMSH
     * Applies to 2D and the horizontal component of 3D ellipsoidal  systems.
     * @deprecated
     */
    public function EPSG_DEGREE_REPRESENTATION_CONVERSION_HDEG_TO_DMSH(): Point
    {
    }

    /**
     * Degree representation conversion: deg to DMSH
     * Applies to 2D and the horizontal component of 3D ellipsoidal  systems.
     * @deprecated
     */
    public function EPSG_DEGREE_REPRESENTATION_CONVERSION_DEG_TO_DMSH(): Point
    {
    }

    /**
     * Degree representation conversion: degH to DMSH
     * Applies to 2D and the horizontal component of 3D ellipsoidal  systems.
     * @deprecated
     */
    public function EPSG_DEGREE_REPRESENTATION_CONVERSION_DEGH_TO_DMSH(): Point
    {
    }

    /**
     * Equal Earth.
     */
    public function EPSG_EQUAL_EARTH(): Point
    {
    }

    /**
     * Equidistant Cylindrical
     * See method code 1029 for spherical development. See also Pseudo Plate Carree, method code 9825.
     */
    public function EPSG_EQUIDISTANT_CYLINDRICAL(): Point
    {
    }

    /**
     * Equidistant Cylindrical (Spherical)
     * See method code 1028 for ellipsoidal development. If the latitude of natural origin is at the equator, also
     * known as Plate Carrée. See also Pseudo Plate Carree, method code 9825.
     */
    public function EPSG_EQUIDISTANT_CYLINDRICAL_SPHERICAL(): Point
    {
    }

    /**
     * GNTRANS.
     */
    public function EPSG_GNTRANS(): Point
    {
    }

    /**
     * General polynomial of degree 2.
     */
    public function EPSG_GENERAL_POLYNOMIAL_OF_DEGREE_2(): Point
    {
    }

    /**
     * General polynomial of degree 3.
     */
    public function EPSG_GENERAL_POLYNOMIAL_OF_DEGREE_3(): Point
    {
    }

    /**
     * General polynomial of degree 4.
     */
    public function EPSG_GENERAL_POLYNOMIAL_OF_DEGREE_4(): Point
    {
    }

    /**
     * General polynomial of degree 6.
     */
    public function EPSG_GENERAL_POLYNOMIAL_OF_DEGREE_6(): Point
    {
    }

    /**
     * Geocentric translation by Grid Interpolation (IGN).
     */
    public function EPSG_GEOCENTRIC_TRANSLATION_BY_GRID_INTERPOLATION_IGN(): Point
    {
    }

    /**
     * Geocentric translations  (geog3D domain)
     * See methods 1031 and 9603 for similar tfms operating between other CRSs types.
     */
    public function EPSG_GEOCENTRIC_TRANSLATIONS_GEOG3D_DOMAIN(): Point
    {
    }

    /**
     * Geocentric translations (geocentric domain)
     * This method allows calculation of geocentric coords in the target system by adding the parameter values to the
     * corresponding coordinates of the point in the source system. See methods 1035 and 9603 for similar tfms
     * operating between other CRSs types.
     */
    public function EPSG_GEOCENTRIC_TRANSLATIONS_GEOCENTRIC_DOMAIN(): Point
    {
    }

    /**
     * Geocentric translations (geog2D domain)
     * See methods 1031 and 1035 for similar tfms operating between other CRSs types.
     */
    public function EPSG_GEOCENTRIC_TRANSLATIONS_GEOG2D_DOMAIN(): Point
    {
    }

    /**
     * Geocentric/topocentric conversions.
     */
    public function EPSG_GEOCENTRIC_TOPOCENTRIC_CONVERSIONS(): Point
    {
    }

    /**
     * Geog3D to Geog2D+GravityRelatedHeight (US .gtx)
     * Transformation from a Geographic 3D CRS to a Compound CRS consisting of a Geographic 2D CRS and a Vertical CRS,
     * or vice versa. The Geographic 3D and the Geographic 2D CRS must be based on the same Geodetic Datum.
     */
    public function EPSG_GEOG3D_TO_GEOG2D_PLUS_GRAVITYRELATEDHEIGHT_US_GTX(): Point
    {
    }

    /**
     * Geog3D to Geog2D+Vertical (AUSGeoid v2)
     * Transformation from a Geographic 3D CRS to a Compound CRS consisting of a Geographic 2D CRS and a Vertical CRS,
     * or vice versa. The Geographic 3D CRS and the Geographic 2D CRS must be based on the same Geodetic Datum.
     */
    public function EPSG_GEOG3D_TO_GEOG2D_PLUS_VERTICAL_AUSGEOID_V2(): Point
    {
    }

    /**
     * Geographic/geocentric conversions
     * This is a parameter-less conversion. In applications it is often concatenated with the 3- 7- or 10-parameter
     * transformations 9603, 9606, 9607 or 9636 to form a geographic to geographic transformation.
     * @param  GeocentricValue|GeographicValue $point
     * @return GeocentricValue|GeographicValue
     */
    public function geographicGeocentricConversion($point, CoordinateReferenceSystem $crs)
    {
        $a = $crs->getDatum()->getEllipsoid()->getSemiMajorAxis()->getValue();
        $e2 = $crs->getDatum()->getEllipsoid()->getEccentricitySquared();

        if ($point instanceof GeocentricValue) {
            $x = $point->getX()->asMetres()->getValue();
            $y = $point->getY()->asMetres()->getValue();
            $z = $point->getZ()->asMetres()->getValue();

            $longitude = atan2($y, $x);
            $p = sqrt(($x ** 2) + ($y ** 2));

            $latitude = atan($z / ($p * (1 - $e2)));

            do {
                $phi1 = $latitude;
                $v = $a / sqrt(1 - $e2 * (sin($latitude) ** 2));
                $latitude = atan(($z + ($e2 * $v * sin($latitude))) / $p);
            } while (abs($latitude - $phi1) >= 0.0000000001);

            $latitude += $crs->getDatum()->getPrimeMeridian()->getGreenwichLongitude()->asRadians()->getValue();
            $h = $p / cos($latitude) - $v;

            return new GeographicValue(new Radian($latitude), new Radian($longitude), new Metre($h));
        }

        if ($point instanceof GeographicValue) {

$point is always a sub-type of PHPCoord\CoordinateOperation\GeographicValue.

            $latitude = $point->getLatitude()->asRadians()->getValue();
            $longitude = $point->getLongitude()->asRadians()->getValue() - $crs->getDatum()->getPrimeMeridian()->getGreenwichLongitude()->asRadians()->getValue();
            $h = $point->getHeight()->asMetres()->getValue();

            $nu = $a / sqrt(1 - $e2 * (sin($latitude) ** 2));
            $x = ($nu + $h) * cos($latitude) * cos($longitude);
            $y = ($nu + $h) * cos($latitude) * sin($longitude);
            $z = ((1 - $e2) * $nu + $h) * sin($latitude);

            return new GeocentricValue(new Metre($x), new Metre($y), new Metre($z));
        }

        throw new InvalidArgumentException('Wrong point type');
    }

    /**
     * Geographic/topocentric conversions.
     */
    public function EPSG_GEOGRAPHIC_TOPOCENTRIC_CONVERSIONS(): Point
    {
    }

    /**
     * Geographic2D offsets
     * This transformation allows calculation of coordinates in the target system by adding the parameter value to the
     * coordinate values of the point in the source system.
     */
    public function EPSG_GEOGRAPHIC2D_OFFSETS(): Point
    {
    }

    /**
     * Geographic2D with Height Offsets
     * This transformation allows calculation of coordinates in the target system by adding the parameter value to the
     * coordinate values of the point in the source system.
     */
    public function EPSG_GEOGRAPHIC2D_WITH_HEIGHT_OFFSETS(): Point
    {
    }

    /**
     * Geographic3D offsets
     * This transformation allows calculation of coordinates in the target system by adding the parameter value to the
     * coordinate values of the point in the source system.
     */
    public function EPSG_GEOGRAPHIC3D_OFFSETS(): Point
    {
    }

    /**
     * Geographic3D to 2D conversion
     * This is a parameter-less conversion.
     */
    public function EPSG_GEOGRAPHIC3D_TO_2D_CONVERSION(): Point
    {
    }

    /**
     * Geographic3D to GravityRelatedHeight (Ausgeoid v2)
     * The Information Source references software which offers both bi-cubic and bi-linear interpolation methods.
     * Unlike earlier Ausgeoid98 method which used bi-linear interpolation, Ausgeoid v2 uses bi-cubic. See Info Source
     * for file format documentation.
     */
    public function EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_AUSGEOID_V2(): Point
    {
    }

    /**
     * Geographic3D to GravityRelatedHeight (Ausgeoid98).
     */
    public function EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_AUSGEOID98(): Point
    {
    }

    /**
     * Geographic3D to GravityRelatedHeight (BEV AT).
     */
    public function EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_BEV_AT(): Point
    {
    }

    /**
     * Geographic3D to GravityRelatedHeight (CGG2013)
     * For consistency with earlier models the Information Source references software which uses bi-quadratic
     * interpolation. Because of denser node spacing in CGG2013 bi-linear interpolation will be sufficient for most
     * uses. See Info Source for file format doc.
     */
    public function EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_CGG2013(): Point
    {
    }

    /**
     * Geographic3D to GravityRelatedHeight (CI).
     */
    public function EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_CI(): Point
    {
    }

    /**
     * Geographic3D to GravityRelatedHeight (EGM)
     * Applies to EGM84 and EGM96 models. For later model see Geographic3D to GravityRelatedHeight (EGM2008), method
     * code 1025.
     */
    public function EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_EGM(): Point
    {
    }

    /**
     * Geographic3D to GravityRelatedHeight (EGM2008)
     * For earlier EGM84 and EGM96 models see Geographic3D to GravityRelatedHeight (EGM), method code 9661.
     */
    public function EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_EGM2008(): Point
    {
    }

    /**
     * Geographic3D to GravityRelatedHeight (Gravsoft).
     */
    public function EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_GRAVSOFT(): Point
    {
    }

    /**
     * Geographic3D to GravityRelatedHeight (IGN1997)
     * Transformation of the vertical component of a Geographic 3D CRS to a Vertical CRS.
     */
    public function EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_IGN1997(): Point
    {
    }

    /**
     * Geographic3D to GravityRelatedHeight (IGN2009)
     * Transformation of the vertical component of a Geographic 3D CRS to a Vertical CRS. File header of method code
     * 9664  (4 lines) has changed in this method (1 line)(){}  recommended interpolation method now in a separate XML
     * file with same name as the grid.
     */
    public function EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_IGN2009(): Point
    {
    }

    /**
     * Geographic3D to GravityRelatedHeight (NZgeoid)
     * EPSG initially gave this method the name "Geographic3D to GravityRelatedHeight (NZgeoid2009)". As the same file
     * format was retained for the 2016 geoid, date removed from the method name.
     */
    public function EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_NZGEOID(): Point
    {
    }

    /**
     * Geographic3D to GravityRelatedHeight (OSGM-GB)
     * Transformation of the vertical component of a Geographic 3D CRS to a Vertical CRS.
     */
    public function EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_OSGM_GB(): Point
    {
    }

    /**
     * Geographic3D to GravityRelatedHeight (OSGM02-Ire)
     * Transformation of the vertical component of a Geographic 3D CRS to a Vertical CRS.
     */
    public function EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_OSGM02_IRE(): Point
    {
    }

    /**
     * Geographic3D to GravityRelatedHeight (OSGM15-Ire)
     * Transformation of the vertical component of a Geographic 3D CRS to a Vertical CRS.
     */
    public function EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_OSGM15_IRE(): Point
    {
    }

    /**
     * Geographic3D to GravityRelatedHeight (PNG)
     * Transformation of the vertical component of a Geographic 3D CRS to a Vertical CRS.
     */
    public function EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_PNG(): Point
    {
    }

    /**
     * Geographic3D to GravityRelatedHeight (SA 2010)
     * File format: space-separated ascii file, no header, columns: latitude, longitude, separation.
     */
    public function EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_SA_2010(): Point
    {
    }

    /**
     * Geographic3D to GravityRelatedHeight (gtx)
     * Transformation of the vertical component of a Geographic 3D CRS to a Vertical CRS. Grid file format: US NGS .gtx
     * (in US sometimes also referred to as 'vdatum format').
     */
    public function EPSG_GEOGRAPHIC3D_TO_GRAVITYRELATEDHEIGHT_GTX(): Point
    {
    }

    /**
     * Guam Projection
     * Simplified form of Oblique Azimuthal Equidistant projection method.
     */
    public function EPSG_GUAM_PROJECTION(): Point
    {
    }

    /**
     * Height Depth Reversal
     * This is a parameter-less conversion.
     */
    public function EPSG_HEIGHT_DEPTH_REVERSAL(): Point
    {
    }

    /**
     * Hotine Oblique Mercator (variant A).
     */
    public function EPSG_HOTINE_OBLIQUE_MERCATOR_VARIANT_A(): Point
    {
    }

    /**
     * Hotine Oblique Mercator (variant B).
     */
    public function EPSG_HOTINE_OBLIQUE_MERCATOR_VARIANT_B(): Point
    {
    }

    /**
     * Hyperbolic Cassini-Soldner.
     */
    public function EPSG_HYPERBOLIC_CASSINI_SOLDNER(): Point
    {
    }

    /**
     * Krovak.
     */
    public function EPSG_KROVAK(): Point
    {
    }

    /**
     * Krovak (North Orientated).
     */
    public function EPSG_KROVAK_NORTH_ORIENTATED(): Point
    {
    }

    /**
     * Krovak Modified
     * Incorporates a polynomial transformation which is defined to be exact and for practical purposes is considered
     * to be a map projection.
     */
    public function EPSG_KROVAK_MODIFIED(): Point
    {
    }

    /**
     * Krovak Modified (North Orientated)
     * Incorporates a polynomial transformation which is defined to be exact and for practical purposes is considered
     * to be a map projection.
     */
    public function EPSG_KROVAK_MODIFIED_NORTH_ORIENTATED(): Point
    {
    }

    /**
     * Laborde Oblique Mercator.
     */
    public function EPSG_LABORDE_OBLIQUE_MERCATOR(): Point
    {
    }

    /**
     * Lambert Azimuthal Equal Area
     * This is the ellipsoidal form of the projection.
     */
    public function EPSG_LAMBERT_AZIMUTHAL_EQUAL_AREA(): Point
    {
    }

    /**
     * Lambert Azimuthal Equal Area (Spherical)
     * This is the spherical form of the projection.  See coordinate operation method Lambert Azimuthal Equal Area
     * (code 9820) for ellipsoidal form.  Differences of several tens of metres result from comparison of the two
     * methods.
     */
    public function EPSG_LAMBERT_AZIMUTHAL_EQUAL_AREA_SPHERICAL(): Point
    {
    }

    /**
     * Lambert Conic Conformal (1SP).
     */
    public function EPSG_LAMBERT_CONIC_CONFORMAL_1SP(): Point
    {
    }

    /**
     * Lambert Conic Conformal (2SP Belgium)
     * In 2000 this modification was replaced through use of the regular Lambert Conic Conformal (2SP) method [9802]
     * with appropriately modified parameter values.
     */
    public function EPSG_LAMBERT_CONIC_CONFORMAL_2SP_BELGIUM(): Point
    {
    }

    /**
     * Lambert Conic Conformal (2SP Michigan).
     */
    public function EPSG_LAMBERT_CONIC_CONFORMAL_2SP_MICHIGAN(): Point
    {
    }

    /**
     * Lambert Conic Conformal (2SP).
     */
    public function EPSG_LAMBERT_CONIC_CONFORMAL_2SP(): Point
    {
    }

    /**
     * Lambert Conic Conformal (West Orientated).
     */
    public function EPSG_LAMBERT_CONIC_CONFORMAL_WEST_ORIENTATED(): Point
    {
    }

    /**
     * Lambert Conic Near-Conformal
     * The Lambert Near-Conformal projection is derived from the Lambert Conformal Conic projection by truncating the
     * series expansion of the projection formulae.
     */
    public function EPSG_LAMBERT_CONIC_NEAR_CONFORMAL(): Point
    {
    }

    /**
     * Lambert Cylindrical Equal Area
     * This is the ellipsoidal form of the projection.
     */
    public function EPSG_LAMBERT_CYLINDRICAL_EQUAL_AREA(): Point
    {
    }

    /**
     * Lambert Cylindrical Equal Area (Spherical)
     * This is the spherical form of the projection.  See coordinate operation method Lambert Cylindrical Equal Area
     * (code 9835) for ellipsoidal form.  Differences of several tens of metres result from comparison of the two
     * methods.
     */
    public function EPSG_LAMBERT_CYLINDRICAL_EQUAL_AREA_SPHERICAL(): Point
    {
    }

    /**
     * Longitude rotation
     * This transformation allows calculation of the longitude of a point in the target system by adding the parameter
     * value to the longitude value of the point in the source system.
     */
    public function EPSG_LONGITUDE_ROTATION(): Point
    {
    }

    /**
     * Madrid to ED50 polynomial.
     */
    public function EPSG_MADRID_TO_ED50_POLYNOMIAL(): Point
    {
    }

    /**
     * Maritime Provinces polynomial interpolation
     * This transformation is an executable module within the application NBGeocalc.  It is an adaptation of the ESTPM
     * program developed by Geodetic Survey of Canada.
     */
    public function EPSG_MARITIME_PROVINCES_POLYNOMIAL_INTERPOLATION(): Point
    {
    }

    /**
     * Mercator (Spherical).
     */
    public function EPSG_MERCATOR_SPHERICAL(): Point
    {
    }

    /**
     * Mercator (variant A)
     * Note that in these formulas the parameter latitude of natural origin (latO) is not used. However for this
     * Mercator (variant A) method the EPSG dataset includes this parameter, which must have a value of zero, for
     * completeness in CRS labelling.
     */
    public function EPSG_MERCATOR_VARIANT_A(): Point
    {
    }

    /**
     * Mercator (variant B)
     * Used for most nautical charts.
     */
    public function EPSG_MERCATOR_VARIANT_B(): Point
    {
    }

    /**
     * Mercator (variant C).
     */
    public function EPSG_MERCATOR_VARIANT_C(): Point
    {
    }

    /**
     * Modified Azimuthal Equidistant
     * Modified form of Oblique Azimuthal Equidistant projection method developed for Polynesian islands. For the
     * distances over which these projections are used (under 800km) this modification introduces no significant error.
     */
    public function EPSG_MODIFIED_AZIMUTHAL_EQUIDISTANT(): Point
    {
    }

    /**
     * Molodensky
     * See Abridged Molodensky.
     */
    public function geographicMolodensky(
        GeographicPoint $from,
        CoordinateReferenceSystem $to,
        Length $translationX,
        Length $translationY,
        Length $translationZ,
        Length $differenceInSemiMajorAxis,
        Scale $differenceInFlattening
    ): GeographicPoint {
        $fromLatitude = $from->getLatitude()->asRadians()->getValue();
        $fromLongitude = $from->getLongitude()->asRadians()->getValue() - $from->getCRS()->getDatum()->getPrimeMeridian()->getGreenwichLongitude()->asRadians()->getValue();
        $fromHeight = $from->getHeight() ? $from->getHeight()->asMetres()->getValue() : 0;
        $tx = $translationX->asMetres()->getValue();
        $ty = $translationY->asMetres()->getValue();
        $tz = $translationZ->asMetres()->getValue();
        $da = $differenceInSemiMajorAxis->asMetres()->getValue();
        $df = $differenceInFlattening->asUnity()->getValue();

        $a = $from->getCRS()->getDatum()->getEllipsoid()->getSemiMajorAxis()->getValue();
        $b = $from->getCRS()->getDatum()->getEllipsoid()->getSemiMinorAxis()->getValue();
        $e2 = $from->getCRS()->getDatum()->getEllipsoid()->getEccentricitySquared();

        $rho = $a * (1 - $e2) / (1 - $e2 * sin($fromLatitude) ** 2) ** (3 / 2);
        $nu = $a / ((1 - $e2 * (sin($fromLatitude) ** 2)) ** 0.5);

        $f = $from->getCRS()->getDatum()->getEllipsoid()->getInverseFlattening();

The assignment to $f is dead and can be removed.

        $dLatitude = ((-$tx * sin($fromLatitude) * cos($fromLongitude)) - ($ty * sin($fromLatitude) * sin($fromLongitude)) + ($tz * cos($fromLatitude)) + ($da * ($nu * $e2 * sin($fromLatitude) * cos($fromLatitude)) / $a + $df * ($rho * ($a / $b) + $nu * ($b / $a)) * sin($fromLatitude) * cos($fromLatitude))) / (($rho + $fromHeight) * sin((new ArcSecond(1))->asRadians()->getValue()));
        $dLongitude = (-$tx * sin($fromLongitude) + $ty * cos($fromLongitude)) / ((($nu + $fromHeight) * cos($fromLatitude)) * sin((new ArcSecond(1))->asRadians()->getValue()));
        $dHeight = ($tx * cos($fromLatitude) * cos($fromLongitude)) + ($ty * cos($fromLatitude) * sin($fromLongitude)) + ($tz * sin($fromLatitude)) - $da * $a / $nu + $df * $b / $a * $nu * sin($fromLatitude) ** 2;

        $toLatitude = $fromLatitude + (new ArcSecond($dLatitude))->asRadians()->getValue();
        $toLongitude = $fromLongitude + (new ArcSecond($dLongitude))->asRadians()->getValue();
        $toHeight = $fromHeight + $dHeight;

        return new GeographicPoint(new Radian($toLatitude), new Radian($toLongitude), $to instanceof Geographic3D ? new Metre($toHeight) : null, $to);
    }

    /**
     * Molodensky-Badekas (CF geocentric domain)
     * See method codes 1039 and 9636 for this operation in other coordinate domains and method code 1061 for opposite
     * rotation convention in geocentric domain.
     */
    public function geocentricCoordinateFrameMolodenskyBadekas(
        GeocentricPoint $from,
        Geocentric $to,
        Length $translationX,
        Length $translationY,
        Length $translationZ,
        Angle $rotationX,
        Angle $rotationY,
        Angle $rotationZ,
        Scale $scaleFactor,
        Length $pivotX,
        Length $pivotY,
        Length $pivotZ
    ): GeocentricPoint {
        $xs = $from->getX()->asMetres()->getValue();
        $ys = $from->getY()->asMetres()->getValue();
        $zs = $from->getZ()->asMetres()->getValue();
        $tx = $translationX->asMetres()->getValue();
        $ty = $translationY->asMetres()->getValue();
        $tz = $translationZ->asMetres()->getValue();
        $rx = $rotationX->asRadians()->getValue();
        $ry = $rotationY->asRadians()->getValue();
        $rz = $rotationZ->asRadians()->getValue();
        $M = 1 + $scaleFactor->asUnity()->getValue();
        $xp = $pivotX->asMetres()->getValue();
        $yp = $pivotY->asMetres()->getValue();
        $zp = $pivotZ->asMetres()->getValue();

        $xt = $M * ((($xs - $xp) * 1) + (($ys - $yp) * $rz) + (($zs - $zp) * -$ry)) + $tx + $xp;
        $yt = $M * ((($xs - $xp) * -$rz) + (($ys - $yp) * 1) + (($zs - $zp) * $rx)) + $ty + $yp;
        $zt = $M * ((($xs - $xp) * $ry) + (($ys - $yp) * -$rx) + (($zs - $zp) * 1)) + $tz + $zp;

        return new GeocentricPoint(new Metre($xt), new Metre($yt), new Metre($zt), $to);
    }

    /**
     * Molodensky-Badekas (CF geog2D/geog3D domain)
     * See method codes 1034 and 1039/9636 for this operation in other coordinate domains and method code 1062/1063 for the
     * opposite rotation convention in geographic 2D domain.
     */
    public function geographicCoordinateFrameMolodenskyBadekas(
        GeographicPoint $from,
        CoordinateReferenceSystem $to,
        Length $translationX,
        Length $translationY,
        Length $translationZ,
        Angle $rotationX,
        Angle $rotationY,
        Angle $rotationZ,
        Scale $scaleFactor,
        Length $pivotX,
        Length $pivotY,
        Length $pivotZ
    ): GeographicPoint {
        $geographicValue = new GeographicValue($from->getLatitude(), $from->getLongitude(), $from->getHeight());
        $asGeocentric = $this->geographicGeocentricConversion($geographicValue, $from->getCRS());

        $xs = $asGeocentric->getX()->asMetres()->getValue();
        $ys = $asGeocentric->getY()->asMetres()->getValue();
        $zs = $asGeocentric->getZ()->asMetres()->getValue();
        $tx = $translationX->asMetres()->getValue();
        $ty = $translationY->asMetres()->getValue();
        $tz = $translationZ->asMetres()->getValue();
        $rx = $rotationX->asRadians()->getValue();
        $ry = $rotationY->asRadians()->getValue();
        $rz = $rotationZ->asRadians()->getValue();
        $M = 1 + $scaleFactor->asUnity()->getValue();
        $xp = $pivotX->asMetres()->getValue();
        $yp = $pivotY->asMetres()->getValue();
        $zp = $pivotZ->asMetres()->getValue();

        $xt = $M * ((($xs - $xp) * 1) + (($ys - $yp) * $rz) + (($zs - $zp) * -$ry)) + $tx + $xp;
        $yt = $M * ((($xs - $xp) * -$rz) + (($ys - $yp) * 1) + (($zs - $zp) * $rx)) + $ty + $yp;
        $zt = $M * ((($xs - $xp) * $ry) + (($ys - $yp) * -$rx) + (($zs - $zp) * 1)) + $tz + $zp;
        $newGeocentric = new GeocentricValue(new Metre($xt), new Metre($yt), new Metre($zt));
        $newGeographic = $this->geographicGeocentricConversion($newGeocentric, $to);

        return new GeographicPoint($newGeographic->getLatitude(), $newGeographic->getLongitude(), $to instanceof Geographic3D ? $newGeographic->getHeight() : null, $to);
    }

    /**
     * Molodensky-Badekas (PV geocentric domain)
     * See method codes 1062 and 1063 for this operation in other coordinate domains and method code 1034 for opposite
     * rotation convention in geocentric domain.
     */
    public function geocentricPositionVectorMolodenskyBadekas(
        GeocentricPoint $from,
        Geocentric $to,
        Length $translationX,
        Length $translationY,
        Length $translationZ,
        Angle $rotationX,
        Angle $rotationY,
        Angle $rotationZ,
        Scale $scaleFactor,
        Length $pivotX,
        Length $pivotY,
        Length $pivotZ
    ): GeocentricPoint {
        $xs = $from->getX()->asMetres()->getValue();
        $ys = $from->getY()->asMetres()->getValue();
        $zs = $from->getZ()->asMetres()->getValue();
        $tx = $translationX->asMetres()->getValue();
        $ty = $translationY->asMetres()->getValue();
        $tz = $translationZ->asMetres()->getValue();
        $rx = $rotationX->asRadians()->getValue();
        $ry = $rotationY->asRadians()->getValue();
        $rz = $rotationZ->asRadians()->getValue();
        $M = 1 + $scaleFactor->asUnity()->getValue();
        $xp = $pivotX->asMetres()->getValue();
        $yp = $pivotY->asMetres()->getValue();
        $zp = $pivotZ->asMetres()->getValue();

        $xt = $M * ((($xs - $xp) * 1) + (($ys - $yp) * -$rz) + (($zs - $zp) * $ry)) + $tx + $xp;
        $yt = $M * ((($xs - $xp) * $rz) + (($ys - $yp) * 1) + (($zs - $zp) * -$rx)) + $ty + $yp;
        $zt = $M * ((($xs - $xp) * -$ry) + (($ys - $yp) * $rx) + (($zs - $zp) * 1)) + $tz + $zp;

        return new GeocentricPoint(new Metre($xt), new Metre($yt), new Metre($zt), $to);
    }

    /**
     * Molodensky-Badekas (PV geog2D/geog3D domain)
     * See method codes 1061 and 1062/1063 for this operation in other coordinate domains and method code 1039/9636 for opposite
     * rotation in geographic 2D/3D domain.
     */
    public function geographicPositionVectorMolodenskyBadekas(
        GeographicPoint $from,
        CoordinateReferenceSystem $to,
        Length $translationX,
        Length $translationY,
        Length $translationZ,
        Angle $rotationX,
        Angle $rotationY,
        Angle $rotationZ,
        Scale $scaleFactor,
        Length $pivotX,
        Length $pivotY,
        Length $pivotZ
    ): GeographicPoint {
        $geographicValue = new GeographicValue($from->getLatitude(), $from->getLongitude(), $from->getHeight());
        $asGeocentric = $this->geographicGeocentricConversion($geographicValue, $from->getCRS());

        $xs = $asGeocentric->getX()->asMetres()->getValue();
        $ys = $asGeocentric->getY()->asMetres()->getValue();
        $zs = $asGeocentric->getZ()->asMetres()->getValue();
        $tx = $translationX->asMetres()->getValue();
        $ty = $translationY->asMetres()->getValue();
        $tz = $translationZ->asMetres()->getValue();
        $rx = $rotationX->asRadians()->getValue();
        $ry = $rotationY->asRadians()->getValue();
        $rz = $rotationZ->asRadians()->getValue();
        $M = 1 + $scaleFactor->asUnity()->getValue();
        $xp = $pivotX->asMetres()->getValue();
        $yp = $pivotY->asMetres()->getValue();
        $zp = $pivotZ->asMetres()->getValue();

        $xt = $M * ((($xs - $xp) * 1) + (($ys - $yp) * -$rz) + (($zs - $zp) * $ry)) + $tx + $xp;
        $yt = $M * ((($xs - $xp) * $rz) + (($ys - $yp) * 1) + (($zs - $zp) * -$rx)) + $ty + $yp;
        $zt = $M * ((($xs - $xp) * -$ry) + (($ys - $yp) * $rx) + (($zs - $zp) * 1)) + $tz + $zp;
        $newGeocentric = new GeocentricValue(new Metre($xt), new Metre($yt), new Metre($zt));
        $newGeographic = $this->geographicGeocentricConversion($newGeocentric, $to);

        return new GeographicPoint($newGeographic->getLatitude(), $newGeographic->getLongitude(), $to instanceof Geographic3D ? $newGeographic->getHeight() : null, $to);
    }

    /**
     * NADCON
     * Geodetic transformation operating on geographic coordinate differences by bi-linear interpolation.  Input
     * expects longitudes to be positive west.
     */
    public function EPSG_NADCON(): Point
    {
    }

    /**
     * NADCON5 (2D)
     * Geodetic transformation operating on geographic coordinate differences by bi-quadratic interpolation.  Input
     * expects longitudes to be positive east in range 0-360° (0° = Greenwich).
     */
    public function EPSG_NADCON5_2D(): Point
    {
    }

    /**
     * NADCON5 (3D)
     * Geodetic transformation operating on geographic coordinate differences by bi-quadratic interpolation.  Input
     * expects longitudes to be positive east in range 0-360° (0° = Greenwich).
     */
    public function EPSG_NADCON5_3D(): Point
    {
    }

    /**
     * NTv1
     * Geodetic transformation operating on geographic coordinate differences by bi-linear interpolation. Superseded in
     * 1997 by NTv2 (transformation method code 9615).  Input expects longitudes to be positive west.
     */
    public function EPSG_NTV1(): Point
    {
    }

    /**
     * NTv2
     * Geodetic transformation operating on geographic coordinate differences by bi-linear interpolation.  Supersedes
     * NTv1 (transformation method code 9614).  Input expects longitudes to be positive west.
     */
    public function EPSG_NTV2(): Point
    {
    }

    /**
     * New Zealand Deformation Model
     * Zip final contains a set of CSV files documentation for a secular tectonic plate velocity model and patches
     * which account for movements due to specific deformation events (earthquakes).
     */
    public function EPSG_NEW_ZEALAND_DEFORMATION_MODEL(): Point
    {
    }

    /**
     * New Zealand Map Grid.
     */
    public function EPSG_NEW_ZEALAND_MAP_GRID(): Point
    {
    }

    /**
     * Norway Offshore Interpolation
     * Although in principle this method is not reversible, in practice reversibility is better than 10 cm. For the
     * applications for which it was designed it may be considered reversible.
     */
    public function EPSG_NORWAY_OFFSHORE_INTERPOLATION(): Point
    {
    }

    /**
     * Oblique Stereographic
     * This is not the same as the projection method of the same name in USGS Professional Paper no. 1395, "Map
     * Projections - A Working Manual" by John P. Snyder.
     */
    public function EPSG_OBLIQUE_STEREOGRAPHIC(): Point
    {
    }

    /**
     * Ordnance Survey National Transformation
     * Geodetic transformation between ETRS89 (or WGS 84) and OSGB36 / National Grid.  Uses ETRS89 / National Grid as
     * an intermediate coordinate system for bi-linear interpolation of gridded grid coordinate differences.
     */
    public function EPSG_ORDNANCE_SURVEY_NATIONAL_TRANSFORMATION(): Point
    {
    }

    /**
     * Orthographic
     * If the natural origin of the projection is at the topocentric origin, this is a special case of the Vertical
     * Perspective (orthographic case) (method code 9839) in which the ellipsoid height of all mapped points is zero (h
     * = 0).
     */
    public function EPSG_ORTHOGRAPHIC(): Point
    {
    }

    /**
     * Point motion (ellipsoidal).
     */
    public function EPSG_POINT_MOTION_ELLIPSOIDAL(): Point
    {
    }

    /**
     * Point motion (geocen) by grid (INADEFORM).
     */
    public function EPSG_POINT_MOTION_GEOCEN_BY_GRID_INADEFORM(): Point
    {
    }

    /**
     * Point motion (geocentric Cartesian).
     */
    public function EPSG_POINT_MOTION_GEOCENTRIC_CARTESIAN(): Point
    {
    }

    /**
     * Point motion by grid (Canada NTv2_Vel)
     * Interpolation within the grid is in the horizontal component of the source CRS.
     */
    public function EPSG_POINT_MOTION_BY_GRID_CANADA_NTV2_VEL(): Point
    {
    }

    /**
     * Polar Stereographic (variant A)
     * Latitude of natural origin must be either 90 degrees or -90 degrees (or equivalent in alternative angle unit).
     */
    public function EPSG_POLAR_STEREOGRAPHIC_VARIANT_A(): Point
    {
    }

    /**
     * Polar Stereographic (variant B).
     */
    public function EPSG_POLAR_STEREOGRAPHIC_VARIANT_B(): Point
    {
    }

    /**
     * Polar Stereographic (variant C).
     */
    public function EPSG_POLAR_STEREOGRAPHIC_VARIANT_C(): Point
    {
    }

    /**
     * Popular Visualisation Pseudo Mercator
     * Applies spherical formulas to the ellipsoid. As such does not have the properties of a true Mercator projection.
     */
    public function EPSG_POPULAR_VISUALISATION_PSEUDO_MERCATOR(): Point
    {
    }

    /**
     * Position Vector transformation (geocentric domain)
     * This method is a specific case of the Molodensky-Badekas (PV) method (code 1061) in which the evaluation point
     * is the geocentre with coordinate values of zero. Note the analogy with the Coordinate Frame method (code 1032)
     * but beware of the differences!
     */
    public function geocentricPositionVectorTransformation(
        GeocentricPoint $from,
        Geocentric $to,
        Length $translationX,
        Length $translationY,
        Length $translationZ,
        Angle $rotationX,
        Angle $rotationY,
        Angle $rotationZ,
        Scale $scaleFactor
    ): GeocentricPoint {
        return $this->geocentricPositionVectorMolodenskyBadekas(
            $from,
            $to,
            $translationX,
            $translationY,
            $translationZ,
            $rotationX,
            $rotationY,
            $rotationZ,
            $scaleFactor,
            new Metre(0),
            new Metre(0),
            new Metre(0)
        );
    }

    /**
     * Position Vector transformation (geog2D/geog3D domain)
     * Note the analogy with the Coordinate Frame rotation (code 9607/1038) but beware of the differences!  The Position
     * Vector convention is used by IAG and recommended by ISO 19111. See methods 1033/1037/9606 for similar tfms
     * operating between other CRS types.
     * @param Geographic2D|Geographic3D $to
     */
    public function geographicPositionVectorTransformation(
        GeographicPoint $from,
        CoordinateReferenceSystem $to,
        Length $translationX,
        Length $translationY,
        Length $translationZ,
        Angle $rotationX,
        Angle $rotationY,
        Angle $rotationZ,
        Scale $scaleFactor
    ): GeographicPoint {
        return $this->geographicPositionVectorMolodenskyBadekas(
            $from,
            $to,
            $translationX,
            $translationY,
            $translationZ,
            $rotationX,
            $rotationY,
            $rotationZ,
            $scaleFactor,
            new Metre(0),
            new Metre(0),
            new Metre(0)
        );
    }

    /**
     * Pseudo Plate Carree
     * Used only for depiction of graticule (latitude/longitude) coordinates on a computer display. The axes units are
     * decimal degrees and of variable scale. The origin is at Lat = 0, Long = 0. See Equidistant Cylindrical, code
     * 1029, for proper Plate Carrée.
     */
    public function EPSG_PSEUDO_PLATE_CARREE(): Point
    {
    }

    /**
     * Reversible polynomial of degree 13.
     */
    public function EPSG_REVERSIBLE_POLYNOMIAL_OF_DEGREE_13(): Point
    {
    }

    /**
     * Reversible polynomial of degree 2
     * Reversibility is subject to constraints.  See Guidance Note 7 for conditions and clarification.
     */
    public function EPSG_REVERSIBLE_POLYNOMIAL_OF_DEGREE_2(): Point
    {
    }

    /**
     * Reversible polynomial of degree 3
     * Reversibility is subject to constraints.  See Guidance Note 7 for conditions and clarification.
     */
    public function EPSG_REVERSIBLE_POLYNOMIAL_OF_DEGREE_3(): Point
    {
    }

    /**
     * Reversible polynomial of degree 4
     * Reversibility is subject to constraints.  See Guidance Note 7 for conditions and clarification.
     */
    public function EPSG_REVERSIBLE_POLYNOMIAL_OF_DEGREE_4(): Point
    {
    }

    /**
     * Similarity transformation
     * Defined for two-dimensional coordinate systems.
     */
    public function EPSG_SIMILARITY_TRANSFORMATION(): Point
    {
    }

    /**
     * Swiss Oblique Cylindrical
     * Can be accommodated by Oblique Mercator method (code 9815), for which this method is an approximation (see BfL
     * document swissprojectionen.pdf at www.swisstopo.com).
     */
    public function EPSG_SWISS_OBLIQUE_CYLINDRICAL(): Point
    {
    }

    /**
     * Time-dependent Coordinate Frame rotation (geocen)
     * Note the analogy with the Time-dependent Position Vector transformation (code 1053) but beware of the
     * differences!  The Position Vector convention is used by IAG. See method codes 1057 and 1058 for similar methods
     * operating between other CRS types.
     */
    public function EPSG_TIME_DEPENDENT_COORDINATE_FRAME_ROTATION_GEOCEN(): Point
    {
    }

    /**
     * Time-dependent Coordinate Frame rotation (geog2D)
     * Note the analogy with the Time-dependent Position Vector transformation (code 1054) but beware of the
     * differences!  The Position Vector convention is used by IAG. See methods 1056 and 1058 for similar tfms
     * operating between other CRS types.
     */
    public function EPSG_TIME_DEPENDENT_COORDINATE_FRAME_ROTATION_GEOG2D(): Point
    {
    }

    /**
     * Time-dependent Coordinate Frame rotation (geog3D)
     * Note the analogy with the Time-dependent Position Vector transformation (code 1055) but beware of the
     * differences!  The Position Vector convention is used by IAG. See method codes 1056 and 1057 for similar methods
     * operating between other CRS types.
     */
    public function EPSG_TIME_DEPENDENT_COORDINATE_FRAME_ROTATION_GEOG3D(): Point
    {
    }

    /**
     * Time-dependent Position Vector tfm (geocentric)
     * Note the analogy with the Time-dependent Coordinate Frame rotation (code 1056) but beware of the differences!
     * The Position Vector convention is used by IAG. See method codes 1054 and 1055 for similar methods operating
     * between other CRS types.
     */
    public function EPSG_TIME_DEPENDENT_POSITION_VECTOR_TFM_GEOCENTRIC(): Point
    {
    }

    /**
     * Time-dependent Position Vector tfm (geog2D)
     * Note the analogy with the Time-dependent Coordinate Frame rotation (code 1057) but beware of the differences!
     * The Position Vector convention is used by IAG. See method codes 1053 and 1055 for similar methods operating
     * between other CRS types.
     */
    public function EPSG_TIME_DEPENDENT_POSITION_VECTOR_TFM_GEOG2D(): Point
    {
    }

    /**
     * Time-dependent Position Vector tfm (geog3D)
     * Note the analogy with the Coordinate Frame rotation (code 1058) but beware of the differences!  The Position
     * Vector convention is used by IAG. See method codes 1053 and 1054 for similar methods operating between other CRS
     * types.
     */
    public function EPSG_TIME_DEPENDENT_POSITION_VECTOR_TFM_GEOG3D(): Point
    {
    }

    /**
     * Time-specific Coordinate Frame rotation (geocen)
     * Note the analogy with the Time-specific Position Vector method (code 1065) but beware of the differences!
     */
    public function EPSG_TIME_SPECIFIC_COORDINATE_FRAME_ROTATION_GEOCEN(): Point
    {
    }

    /**
     * Time-specific Position Vector transform (geocen)
     * Note the analogy with the Time-specifc Coordinate Frame method (code 1066) but beware of the differences!
     */
    public function EPSG_TIME_SPECIFIC_POSITION_VECTOR_TRANSFORM_GEOCEN(): Point
    {
    }

    /**
     * Transverse Mercator.
     */
    public function EPSG_TRANSVERSE_MERCATOR(): Point
    {
    }

    /**
     * Transverse Mercator (South Orientated).
     */
    public function EPSG_TRANSVERSE_MERCATOR_SOUTH_ORIENTATED(): Point
    {
    }

    /**
     * Transverse Mercator Zoned Grid System
     * If locations fall outwith the fixed zones the general Transverse Mercator method (code 9807) must be used for
     * each zone.
     */
    public function EPSG_TRANSVERSE_MERCATOR_ZONED_GRID_SYSTEM(): Point
    {
    }

    /**
     * Vertical Offset
     * This transformation allows calculation of height (or depth) in the target system by adding the parameter value
     * to the height (or depth)-value of the point in the source system.
     */
    public function EPSG_VERTICAL_OFFSET(): Point
    {
    }

    /**
     * Vertical Offset and Slope
     * This transformation allows calculation of height in the target system by applying the parameter values to the
     * height value of the point in the source system.
     */
    public function EPSG_VERTICAL_OFFSET_AND_SLOPE(): Point
    {
    }

    /**
     * Vertical Offset by Grid Interpolation (BEV AT).
     */
    public function EPSG_VERTICAL_OFFSET_BY_GRID_INTERPOLATION_BEV_AT(): Point
    {
    }

    /**
     * Vertical Offset by Grid Interpolation (NZLVD).
     */
    public function EPSG_VERTICAL_OFFSET_BY_GRID_INTERPOLATION_NZLVD(): Point
    {
    }

    /**
     * Vertical Offset by Grid Interpolation (VERTCON)
     * Any NAD realization may be used as the Interpolation CRS(){} bi-linear interpolation is used. Input expects
     * longitudes to be positive west.
     */
    public function EPSG_VERTICAL_OFFSET_BY_GRID_INTERPOLATION_VERTCON(): Point
    {
    }

    /**
     * Vertical Offset by Grid Interpolation (asc).
     */
    public function EPSG_VERTICAL_OFFSET_BY_GRID_INTERPOLATION_ASC(): Point
    {
    }

    /**
     * Vertical Offset by Grid Interpolation (gtx).
     */
    public function EPSG_VERTICAL_OFFSET_BY_GRID_INTERPOLATION_GTX(): Point
    {
    }

    /**
     * Vertical Perspective
     * For a viewing point height approaching or at infinity, see the Vertical Perspective (orthographic case) (method
     * code 9839).
     */
    public function EPSG_VERTICAL_PERSPECTIVE(): Point
    {
    }

    /**
     * Vertical Perspective (Orthographic case)
     * This is a special case of the general Vertical Perspective (method code 9838) in which the viewing point at
     * infinity.
     */
    public function EPSG_VERTICAL_PERSPECTIVE_ORTHOGRAPHIC_CASE(): Point
    {
    }
}