def osr_usgs_2(): srs = osr.SpatialReference() srs.ImportFromWkt("""PROJCS["unnamed",GEOGCS["NAD27",\ DATUM["North_American_Datum_1927",\ SPHEROID["Clarke 1866",6378206.4,294.9786982139006,\ AUTHORITY["EPSG","7008"]],AUTHORITY["EPSG","6267"]],\ PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],\ AUTHORITY["EPSG","4267"]],PROJECTION["Lambert_Conformal_Conic_2SP"],\ PARAMETER["standard_parallel_1",33.90363402777778],\ PARAMETER["standard_parallel_2",33.62529002777778],\ PARAMETER["latitude_of_origin",33.76446202777777],\ PARAMETER["central_meridian",-117.4745428888889],\ PARAMETER["false_easting",0],PARAMETER["false_northing",0],\ UNIT["metre",1,AUTHORITY["EPSG","9001"]]]""") (proj_code, zone, parms, datum_code) = srs.ExportToUSGS() if proj_code != 4 or datum_code != 0 \ or abs(gdal.PackedDMSToDec(parms[2]) - 33.90363403) > 0.0000005 \ or abs(gdal.PackedDMSToDec(parms[3]) - 33.62529003) > 0.0000005 \ or abs(gdal.PackedDMSToDec(parms[4]) - -117.4745429) > 0.0000005 \ or abs(gdal.PackedDMSToDec(parms[5]) - 33.76446203) > 0.0000005: gdaltest.post_reason( 'Can not import Lambert Conformal Conic projection.') return 'fail' return 'success'
def test_osr_usgs_2(): srs = osr.SpatialReference() srs.ImportFromWkt("""PROJCS["unnamed",GEOGCS["NAD27",\ DATUM["North_American_Datum_1927",\ SPHEROID["Clarke 1866",6378206.4,294.9786982139006,\ AUTHORITY["EPSG","7008"]],AUTHORITY["EPSG","6267"]],\ PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],\ AUTHORITY["EPSG","4267"]],PROJECTION["Lambert_Conformal_Conic_2SP"],\ PARAMETER["standard_parallel_1",33.90363402777778],\ PARAMETER["standard_parallel_2",33.62529002777778],\ PARAMETER["latitude_of_origin",33.76446202777777],\ PARAMETER["central_meridian",-117.4745428888889],\ PARAMETER["false_easting",0],PARAMETER["false_northing",0],\ UNIT["metre",1,AUTHORITY["EPSG","9001"]]]""") (proj_code, _, params, datum_code) = srs.ExportToUSGS() assert proj_code == 4 and datum_code == 0 and gdal.PackedDMSToDec(params[2]) == pytest.approx(33.90363403, abs=0.0000005) and gdal.PackedDMSToDec(params[3]) == pytest.approx(33.62529003, abs=0.0000005) and gdal.PackedDMSToDec(params[4]) == pytest.approx(-117.4745429, abs=0.0000005) and gdal.PackedDMSToDec(params[5]) == pytest.approx(33.76446203, abs=0.0000005), \ 'Can not import Lambert Conformal Conic projection.'