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")