def test_make_geographic_3d_crs(): gcrs = GeographicCRS(ellipsoidal_cs=Ellipsoidal3DCS()) assert gcrs.type_name == "Geographic 3D CRS" expected_authority = ("IGNF", "WGS84GEODD") if PROJ_GTE_901: expected_authority = ("OGC", "CRS84h") assert gcrs.to_authority() == expected_authority
def test_ellipsoidal_3d_cs(axis, name_0, direction_0, name_1, direction_1): vcs = Ellipsoidal3DCS(axis=axis) assert len(vcs.axis_list) == 3 assert vcs.axis_list[0].direction == direction_0 assert vcs.axis_list[0].name == name_0 assert vcs.axis_list[0].unit_name == "degree" assert vcs.axis_list[1].direction == direction_1 assert vcs.axis_list[1].name == name_1 assert vcs.axis_list[1].unit_name == "degree" assert vcs.axis_list[2].direction == "up" assert vcs.axis_list[2].name == "Ellipsoidal height" assert vcs.axis_list[2].unit_name == "metre"
def test_ellipsoidal3dcs_to_cf(): ecs = Ellipsoidal3DCS(axis=Ellipsoidal3DCSAxis.LONGITUDE_LATITUDE_HEIGHT) assert ecs.to_cf() == [ { "standard_name": "longitude", "long_name": "longitude coordinate", "units": "degrees_east", "axis": "X", }, { "standard_name": "latitude", "long_name": "latitude coordinate", "units": "degrees_north", "axis": "Y", }, { "standard_name": "height_above_reference_ellipsoid", "long_name": "Ellipsoidal height", "units": "metre", "positive": "up", "axis": "Z", }, ]
def test_make_geographic_3d_crs(): gcrs = GeographicCRS(ellipsoidal_cs=Ellipsoidal3DCS()) assert gcrs.type_name == "Geographic 3D CRS" assert gcrs.to_authority() == ("IGNF", "WGS84GEODD")