From b52885bf5359b58c4414bf1c110e90dd705d0b88 Mon Sep 17 00:00:00 2001 From: Elias Carvalho Date: Tue, 29 Oct 2024 11:03:17 -0300 Subject: [PATCH] Add PROJ commands used to generate tests --- test/crs/conversions.jl | 85 ++++++++++++++++++++++++++++++++++++++++- 1 file changed, 84 insertions(+), 1 deletion(-) diff --git a/test/crs/conversions.jl b/test/crs/conversions.jl index d4add46..7849fec 100644 --- a/test/crs/conversions.jl +++ b/test/crs/conversions.jl @@ -974,6 +974,7 @@ @testset "Projected" begin @testset "LatLon <> Mercator" begin + # tested against Proj.Transformation("EPSG:4326", "EPSG:3395") c1 = LatLon(T(45), T(90)) c2 = convert(Mercator, c1) @test allapprox(c2, Mercator(T(10018754.171394622), T(5591295.9185533915))) @@ -1006,6 +1007,7 @@ end @testset "LatLon <> WebMercator" begin + # tested against Proj.Transformation("EPSG:4326", "EPSG:3857") c1 = LatLon(T(45), T(90)) c2 = convert(WebMercator, c1) @test allapprox(c2, WebMercator(T(10018754.171394622), T(5621521.486192066))) @@ -1038,6 +1040,7 @@ end @testset "LatLon <> PlateCarree" begin + # tested against Proj.Transformation("EPSG:4326", "EPSG:32662") c1 = LatLon(T(45), T(90)) c2 = convert(PlateCarree, c1) @test allapprox(c2, PlateCarree(T(10018754.171394622), T(5009377.085697311))) @@ -1070,6 +1073,7 @@ end @testset "LatLon <> Lambert" begin + # tested against Proj.Transformation("EPSG:4326", "ESRI:54034") c1 = LatLon(T(45), T(90)) c2 = convert(Lambert, c1) @test allapprox(c2, Lambert(T(10018754.171394622), T(4489858.8869480025))) @@ -1102,6 +1106,7 @@ end @testset "LatLon <> Behrmann" begin + # tested against Proj.Transformation("EPSG:4326", "ESRI:54017") c1 = LatLon(T(45), T(90)) c2 = convert(Behrmann, c1) @test allapprox(c2, Behrmann(T(8683765.222580686), T(5180102.328839251))) @@ -1134,6 +1139,18 @@ end @testset "LatLon <> GallPeters" begin + # forward tested against Proj.Transformation(""" + # proj=pipeline + # step proj=axisswap order=2,1 + # step proj=unitconvert xy_in=deg xy_out=rad + # step proj=cea lat_ts=45 ellps=WGS84 + # """) + # inverse tested against Proj.Transformation(""" + # proj=pipeline + # step proj=cea inv lat_ts=45 ellps=WGS84 + # step proj=unitconvert xy_in=rad xy_out=deg + # step proj=axisswap order=2,1 + # """) c1 = LatLon(T(45), T(90)) c2 = convert(GallPeters, c1) @test allapprox(c2, GallPeters(T(7096215.158458031), T(6338983.732612475))) @@ -1166,6 +1183,7 @@ end @testset "LatLon <> WinkelTripel" begin + # tested against Proj.Transformation("EPSG:4326", "ESRI:54042") c1 = LatLon(T(45), T(90)) c2 = convert(WinkelTripel, c1) @test allapprox(c2, WinkelTripel(T(7044801.6979576545), T(5231448.051548355))) @@ -1204,6 +1222,7 @@ end @testset "LatLon <> Robinson" begin + # tested against Proj.Transformation("EPSG:4326", "ESRI:54030") c1 = LatLon(T(45), T(90)) c2 = convert(Robinson, c1) @test allapprox(c2, Robinson(T(7620313.925950073), T(4805073.646653474))) @@ -1236,6 +1255,18 @@ end @testset "LatLon <> OrthoNorth" begin + # forward tested against Proj.Transformation(""" + # proj=pipeline + # step proj=axisswap order=2,1 + # step proj=unitconvert xy_in=deg xy_out=rad + # step proj=ortho lat_0=90 ellps=WGS84 + # """) + # inverse tested against Proj.Transformation(""" + # proj=pipeline + # step proj=ortho inv lat_0=90 ellps=WGS84 + # step proj=unitconvert xy_in=rad xy_out=deg + # step proj=axisswap order=2,1 + # """) c1 = LatLon(T(30), T(60)) c2 = convert(OrthoNorth, c1) @test allapprox(c2, OrthoNorth(T(4787610.688267582), T(-2764128.319646418))) @@ -1256,6 +1287,18 @@ end @testset "LatLon <> OrthoSouth" begin + # forward tested against Proj.Transformation(""" + # proj=pipeline + # step proj=axisswap order=2,1 + # step proj=unitconvert xy_in=deg xy_out=rad + # step proj=ortho lat_0=-90 ellps=WGS84 + # """) + # inverse tested against Proj.Transformation(""" + # proj=pipeline + # step proj=ortho inv lat_0=-90 ellps=WGS84 + # step proj=unitconvert xy_in=rad xy_out=deg + # step proj=axisswap order=2,1 + # """) c1 = LatLon(-T(30), T(60)) c2 = convert(OrthoSouth, c1) @test allapprox(c2, OrthoSouth(T(4787610.688267582), T(2764128.319646418))) @@ -1279,6 +1322,7 @@ OrthoNorthSpherical = CoordRefSystems.get(ESRI{102035}) OrthoSouthSpherical = CoordRefSystems.get(ESRI{102037}) + # tested against Proj.Transformation("EPSG:4326", "ESRI:102035") c1 = LatLon(T(30), T(60)) c2 = convert(OrthoNorthSpherical, c1) @test allapprox(c2, OrthoNorthSpherical(T(4783602.75), T(-2761814.335408735))) @@ -1291,6 +1335,7 @@ c3 = convert(LatLon, c2) @test allapprox(c3, c1) + # tested against Proj.Transformation("EPSG:4326", "ESRI:102037") c1 = LatLon(-T(30), T(60)) c2 = convert(OrthoSouthSpherical, c1) @test allapprox(c2, OrthoSouthSpherical(T(4783602.75), T(2761814.335408735))) @@ -1319,6 +1364,18 @@ # link: https://sourceforge.net/projects/geographiclib/files/testdata/TMcoords.dat.gz TM = TransverseMercator{0.9996,0.0°,WGS84Latest} + # forward tested against Proj.Transformation(""" + # proj=pipeline + # step proj=axisswap order=2,1 + # step proj=unitconvert xy_in=deg xy_out=rad + # step proj=tmerc lat_0=0 k_0=0.9996 ellps=WGS84 + # """) + # inverse tested against Proj.Transformation(""" + # proj=pipeline + # step proj=tmerc inv lat_0=0 k_0=0.9996 ellps=WGS84 + # step proj=unitconvert xy_in=rad xy_out=deg + # step proj=axisswap order=2,1 + # """) c1 = LatLon(T(70.579277094557), T(45.599419731762)) c2 = convert(TM, c1) @test allapprox(c2, TM(T(1548706.7916191491794), T(8451449.1987722350778))) @@ -1346,6 +1403,18 @@ # latₒ, lonₒ ≠ 0 TM = CoordRefSystems.shift(TransverseMercator{0.9996,15.0°,WGS84Latest}, lonₒ=25.0°) + # forward tested against Proj.Transformation(""" + # proj=pipeline + # step proj=axisswap order=2,1 + # step proj=unitconvert xy_in=deg xy_out=rad + # step proj=tmerc lat_0=15 lon_0=25 k_0=0.9996 ellps=WGS84 + # """) + # inverse tested against Proj.Transformation(""" + # proj=pipeline + # step proj=tmerc inv lat_0=15 lon_0=25 k_0=0.9996 ellps=WGS84 + # step proj=unitconvert xy_in=rad xy_out=deg + # step proj=axisswap order=2,1 + # """) c1 = LatLon(T(30), T(60)) c2 = convert(TM, c1) @test allapprox(c2, TM(T(3478021.0568801453), T(2237622.8976712096))) @@ -1366,7 +1435,7 @@ end @testset "LatLon <> Albers" begin - # tested aganist Proj.Transformation("EPSG:4269", "EPSG:5070") + # tested against Proj.Transformation("EPSG:4269", "EPSG:5070") AlbersUS = CoordRefSystems.shift(Albers{23.0°,29.5°,45.5°,NAD83}, lonₒ=-96°) c1 = LatLon{NAD83}(T(45), T(90)) c2 = convert(AlbersUS, c1) @@ -1403,12 +1472,14 @@ UTMNorth32 = utmnorth(32) UTMSouth59 = utmsouth(59) + # tested against Proj.Transformation("EPSG:4326", "EPSG:32632") c1 = LatLon(T(56), T(12)) c2 = convert(UTMNorth32, c1) @test allapprox(c2, UTMNorth32(T(687071.439107327), T(6210141.326872105))) c3 = convert(LatLon, c2) @test allapprox(c3, c1) + # tested against Proj.Transformation("EPSG:4326", "EPSG:32759") c1 = LatLon(-T(44), T(174)) c2 = convert(UTMSouth59, c1) @test allapprox(c2, UTMSouth59(T(740526.3210524899), T(5123750.873037999))) @@ -1427,6 +1498,18 @@ end @testset "LatLon <> ShiftedCRS" begin + # forward tested against Proj.Transformation(""" + # proj=pipeline + # step proj=axisswap order=2,1 + # step proj=unitconvert xy_in=deg xy_out=rad + # step proj=merc lon_0=15 x_0=200 y_0=200 ellps=WGS84 + # """) + # inverse tested against Proj.Transformation(""" + # proj=pipeline + # step proj=merc inv lon_0=15 x_0=200 y_0=200 ellps=WGS84 + # step proj=unitconvert xy_in=rad xy_out=deg + # step proj=axisswap order=2,1 + # """) ShiftedMercator = CoordRefSystems.shift(Mercator{WGS84Latest}, lonₒ=15.0°, xₒ=200.0m, yₒ=200.0m) c1 = LatLon(T(45), T(90)) c2 = convert(ShiftedMercator, c1)