def test_sinusoidal_operation__defaults(): aeop = SinusoidalConversion() assert aeop.name == "unknown" assert aeop.method_name == "Sinusoidal" assert _to_dict(aeop) == { "Longitude of natural origin": 0.0, "False easting": 0.0, "False northing": 0.0, }
def _sinusoidal(cf_params): """ http://cfconventions.org/cf-conventions/cf-conventions.html#_sinusoidal """ return SinusoidalConversion( longitude_natural_origin=cf_params.get("longitude_of_projection_origin", 0.0), false_easting=cf_params.get("false_easting", 0.0), false_northing=cf_params.get("false_northing", 0.0), )
def test_sinusoidal_operation(): aeop = SinusoidalConversion(longitude_natural_origin=2, false_easting=3, false_northing=4) assert aeop.name == "unknown" assert aeop.method_name == "Sinusoidal" assert _to_dict(aeop) == { "Longitude of natural origin": 2.0, "False easting": 3.0, "False northing": 4.0, }
def test_sinusoidal(): crs = ProjectedCRS(conversion=SinusoidalConversion(0, 1, 2)) expected_cf = { "semi_major_axis": 6378137.0, "semi_minor_axis": crs.ellipsoid.semi_minor_metre, "inverse_flattening": crs.ellipsoid.inverse_flattening, "reference_ellipsoid_name": "WGS 84", "longitude_of_prime_meridian": 0.0, "prime_meridian_name": "Greenwich", "horizontal_datum_name": "World Geodetic System 1984", "grid_mapping_name": "sinusoidal", "longitude_of_projection_origin": 0.0, "false_easting": 1.0, "false_northing": 2.0, } cf_dict = crs.to_cf() assert cf_dict.pop("crs_wkt").startswith("PROJCRS[") assert cf_dict == expected_cf # test roundtrip _test_roundtrip(expected_cf, "PROJCRS[")