Opened 9 years ago

Closed 5 years ago

#5899 closed defect (wontfix)

SRS with UNIT = GRAD not preserved by GTiff

Reported by: Mike Taves Owned by: warmerdam
Priority: low Milestone: closed_because_of_github_migration
Component: default Version: 1.11.1
Severity: normal Keywords:
Cc:

Description

Lambert Conformal Conic projection codes 27500, 27561, 27562, 27563, 27564, 27571, 27572, 27573, 27574, 27581, 27582, 27583, 27584, 27591, 27592, 27593 and 27594 are not preserved by the GTiff driver.

There are seemingly random differences with various parameters in PROJ.4 and WKT outputs. See #5897 for test_prj_roundtrip, and run:

for code in [27500, 27561, 27562, 27563, 27564, 27571, 27572, 27573, 27574, 27581, 27582, 27583, 27584, 27591, 27592, 27593, 27594]:
    print(test_prj_roundtrip(code, True, True))

with output

EPSG:27500 
Differences in PROJ.4:
- +proj=lcc +lat_1=49.50000000000001 +lat_0=49.50000000000001 +lon_0=5.4 +k_0=0.9995090800000001 +x_0=500000 +y_0=300000 +a=6376523 +b=6355862.933255573 +pm=2.337208333333333 +units=m +no_defs 
?                   ^               ------------------------------------                                                                                       ^^^^^ ^^^^^^^^^
+ +proj=lcc +lat_1=44.55 +lat_0=44.55 +lon_0=4.860000000000001 +k_0=0.9995090800000001 +x_0=500000 +y_0=300000 +a=6376523 +b=6355862.933255573 +pm=2.596898148148148 +units=m +no_defs 
?                   ^  ++++++++++++++++++++++++++                                                                                                    ^^^ ^^^^^^^^^^^
Differences in WKT:
  PROJCS["ATF (Paris) / Nord de Guerre",
      GEOGCS["ATF (Paris)",
          DATUM["Ancienne_Triangulation_Francaise_Paris",
-             SPHEROID["Plessis 1817",6376523,308.64,
?                                                  ^
+             SPHEROID["Plessis 1817",6376523,308.6399999999991,
?                                                  ^^^^^^^^^^^^
                  AUTHORITY["EPSG","7027"]],
              AUTHORITY["EPSG","6901"]],
+         PRIMEM["Paris RGS",2.596898148148148],
-         PRIMEM["Paris RGS",2.337208333333333,
-             AUTHORITY["EPSG","8914"]],
-         UNIT["grad",0.01570796326794897,
+         UNIT["grad",0.01570796326794897],
?                                        +
-             AUTHORITY["EPSG","9105"]],
          AUTHORITY["EPSG","4901"]],
      PROJECTION["Lambert_Conformal_Conic_1SP"],
-     PARAMETER["latitude_of_origin",55],
?                                    ^
+     PARAMETER["latitude_of_origin",49.5],
?                                    ^^^
-     PARAMETER["central_meridian",5.999999999999998],
?                                    ^^^^^^^^^^^^^^^
+     PARAMETER["central_meridian",5.4],
?                                    ^
      PARAMETER["scale_factor",0.99950908],
      PARAMETER["false_easting",500000],
      PARAMETER["false_northing",300000],
      UNIT["metre",1,
          AUTHORITY["EPSG","9001"]],
-     AXIS["X",EAST],
-     AXIS["Y",NORTH],
      AUTHORITY["EPSG","27500"]]
EPSG:27561 
Differences in PROJ.4:
- +proj=lcc +lat_1=49.50000000000001 +lat_0=49.50000000000001 +lon_0=0 +k_0=0.999877341 +x_0=600000 +y_0=200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=paris +units=m +no_defs 
?                   ^  ^^^^^^^^^^^^^         ^  ^^^^^^^^^^^^^                                                                                                             ^^^^^
+ +proj=lcc +lat_1=44.55 +lat_0=44.55 +lon_0=0 +k_0=0.999877341 +x_0=600000 +y_0=200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=2.5969213 +units=m +no_defs 
?                   ^  ^         ^  ^                                                                                                             ^^^^^^^^^
Differences in WKT:
  PROJCS["NTF (Paris) / Lambert Nord France",
      GEOGCS["NTF (Paris)",
          DATUM["Nouvelle_Triangulation_Francaise_Paris",
              SPHEROID["Clarke 1880 (IGN)",6378249.2,293.4660212936265,
                  AUTHORITY["EPSG","7011"]],
              TOWGS84[-168,-60,320,0,0,0,0],
              AUTHORITY["EPSG","6807"]],
-         PRIMEM["Paris",2.33722917,
?                           ^^^^^^^
+         PRIMEM["Paris",2.5969213],
?                          ++++++ ^
-             AUTHORITY["EPSG","8903"]],
-         UNIT["grad",0.01570796326794897,
+         UNIT["grad",0.01570796326794897],
?                                        +
-             AUTHORITY["EPSG","9105"]],
          AUTHORITY["EPSG","4807"]],
      PROJECTION["Lambert_Conformal_Conic_1SP"],
-     PARAMETER["latitude_of_origin",55],
?                                    ^
+     PARAMETER["latitude_of_origin",49.5],
?                                    ^^^
      PARAMETER["central_meridian",0],
      PARAMETER["scale_factor",0.999877341],
      PARAMETER["false_easting",600000],
      PARAMETER["false_northing",200000],
      UNIT["metre",1,
          AUTHORITY["EPSG","9001"]],
-     AXIS["X",EAST],
-     AXIS["Y",NORTH],
      AUTHORITY["EPSG","27561"]]
EPSG:27562 
Differences in PROJ.4:
- +proj=lcc +lat_1=46.8 +lat_0=46.8 +lon_0=0 +k_0=0.99987742 +x_0=600000 +y_0=200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=paris +units=m +no_defs 
?                   ^ ^         ^ ^                                                                                                            ^^^^^
+ +proj=lcc +lat_1=42.12000000000001 +lat_0=42.12000000000001 +lon_0=0 +k_0=0.99987742 +x_0=600000 +y_0=200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=2.5969213 +units=m +no_defs 
?                   ^ ^^^^^^^^^^^^^^         ^ ^^^^^^^^^^^^^^                                                                                                            ^^^^^^^^^
Differences in WKT:
  PROJCS["NTF (Paris) / Lambert Centre France",
      GEOGCS["NTF (Paris)",
          DATUM["Nouvelle_Triangulation_Francaise_Paris",
              SPHEROID["Clarke 1880 (IGN)",6378249.2,293.4660212936265,
                  AUTHORITY["EPSG","7011"]],
              TOWGS84[-168,-60,320,0,0,0,0],
              AUTHORITY["EPSG","6807"]],
-         PRIMEM["Paris",2.33722917,
?                           ^^^^^^^
+         PRIMEM["Paris",2.5969213],
?                          ++++++ ^
-             AUTHORITY["EPSG","8903"]],
-         UNIT["grad",0.01570796326794897,
+         UNIT["grad",0.01570796326794897],
?                                        +
-             AUTHORITY["EPSG","9105"]],
          AUTHORITY["EPSG","4807"]],
      PROJECTION["Lambert_Conformal_Conic_1SP"],
-     PARAMETER["latitude_of_origin",52],
?                                    ^^
+     PARAMETER["latitude_of_origin",46.8],
?                                    ^^^^
      PARAMETER["central_meridian",0],
      PARAMETER["scale_factor",0.99987742],
      PARAMETER["false_easting",600000],
      PARAMETER["false_northing",200000],
      UNIT["metre",1,
          AUTHORITY["EPSG","9001"]],
-     AXIS["X",EAST],
-     AXIS["Y",NORTH],
      AUTHORITY["EPSG","27562"]]
EPSG:27563 
Differences in PROJ.4:
- +proj=lcc +lat_1=44.10000000000001 +lat_0=44.10000000000001 +lon_0=0 +k_0=0.999877499 +x_0=600000 +y_0=200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=paris +units=m +no_defs 
?                  ^^ ^^                    ^^ ^^                                                                                                                         ^^^^^
+ +proj=lcc +lat_1=39.69000000000001 +lat_0=39.69000000000001 +lon_0=0 +k_0=0.999877499 +x_0=600000 +y_0=200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=2.5969213 +units=m +no_defs 
?                  ^^ ^^                    ^^ ^^                                                                                                                         ^^^^^^^^^
Differences in WKT:
  PROJCS["NTF (Paris) / Lambert Sud France",
      GEOGCS["NTF (Paris)",
          DATUM["Nouvelle_Triangulation_Francaise_Paris",
              SPHEROID["Clarke 1880 (IGN)",6378249.2,293.4660212936265,
                  AUTHORITY["EPSG","7011"]],
              TOWGS84[-168,-60,320,0,0,0,0],
              AUTHORITY["EPSG","6807"]],
-         PRIMEM["Paris",2.33722917,
?                           ^^^^^^^
+         PRIMEM["Paris",2.5969213],
?                          ++++++ ^
-             AUTHORITY["EPSG","8903"]],
-         UNIT["grad",0.01570796326794897,
+         UNIT["grad",0.01570796326794897],
?                                        +
-             AUTHORITY["EPSG","9105"]],
          AUTHORITY["EPSG","4807"]],
      PROJECTION["Lambert_Conformal_Conic_1SP"],
-     PARAMETER["latitude_of_origin",49],
?                                     ^
+     PARAMETER["latitude_of_origin",44.1],
?                                     ^^^
      PARAMETER["central_meridian",0],
      PARAMETER["scale_factor",0.999877499],
      PARAMETER["false_easting",600000],
      PARAMETER["false_northing",200000],
      UNIT["metre",1,
          AUTHORITY["EPSG","9001"]],
-     AXIS["X",EAST],
-     AXIS["Y",NORTH],
      AUTHORITY["EPSG","27563"]]
EPSG:27564 
Differences in PROJ.4:
- +proj=lcc +lat_1=42.16500000000001 +lat_0=42.16500000000001 +lon_0=0 +k_0=0.99994471 +x_0=234.358 +y_0=185861.369 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=paris +units=m +no_defs 
?                   ^^^^ -----------         ^^^^ -----------                                                                                                                 ^^^^^
+ +proj=lcc +lat_1=37.9485 +lat_0=37.9485 +lon_0=0 +k_0=0.99994471 +x_0=234.358 +y_0=185861.369 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=2.5969213 +units=m +no_defs 
?                  ++++ ^         ++++ ^                                                                                                                  ^^^^^^^^^
Differences in WKT:
  PROJCS["NTF (Paris) / Lambert Corse",
      GEOGCS["NTF (Paris)",
          DATUM["Nouvelle_Triangulation_Francaise_Paris",
              SPHEROID["Clarke 1880 (IGN)",6378249.2,293.4660212936265,
                  AUTHORITY["EPSG","7011"]],
              TOWGS84[-168,-60,320,0,0,0,0],
              AUTHORITY["EPSG","6807"]],
-         PRIMEM["Paris",2.33722917,
?                           ^^^^^^^
+         PRIMEM["Paris",2.5969213],
?                          ++++++ ^
-             AUTHORITY["EPSG","8903"]],
-         UNIT["grad",0.01570796326794897,
+         UNIT["grad",0.01570796326794897],
?                                        +
-             AUTHORITY["EPSG","9105"]],
          AUTHORITY["EPSG","4807"]],
      PROJECTION["Lambert_Conformal_Conic_1SP"],
-     PARAMETER["latitude_of_origin",46.85],
?                                      --
+     PARAMETER["latitude_of_origin",42.165],
?                                     +++
      PARAMETER["central_meridian",0],
      PARAMETER["scale_factor",0.99994471],
      PARAMETER["false_easting",234.358],
      PARAMETER["false_northing",185861.369],
      UNIT["metre",1,
          AUTHORITY["EPSG","9001"]],
-     AXIS["X",EAST],
-     AXIS["Y",NORTH],
      AUTHORITY["EPSG","27564"]]
EPSG:27571 
Differences in PROJ.4:
- +proj=lcc +lat_1=49.50000000000001 +lat_0=49.50000000000001 +lon_0=0 +k_0=0.999877341 +x_0=600000 +y_0=1200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=paris +units=m +no_defs 
?                   ^  ^^^^^^^^^^^^^         ^  ^^^^^^^^^^^^^                                                                                                              ^^^^^
+ +proj=lcc +lat_1=44.55 +lat_0=44.55 +lon_0=0 +k_0=0.999877341 +x_0=600000 +y_0=1200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=2.5969213 +units=m +no_defs 
?                   ^  ^         ^  ^                                                                                                              ^^^^^^^^^
Differences in WKT:
  PROJCS["NTF (Paris) / Lambert zone I",
      GEOGCS["NTF (Paris)",
          DATUM["Nouvelle_Triangulation_Francaise_Paris",
              SPHEROID["Clarke 1880 (IGN)",6378249.2,293.4660212936265,
                  AUTHORITY["EPSG","7011"]],
              TOWGS84[-168,-60,320,0,0,0,0],
              AUTHORITY["EPSG","6807"]],
-         PRIMEM["Paris",2.33722917,
?                           ^^^^^^^
+         PRIMEM["Paris",2.5969213],
?                          ++++++ ^
-             AUTHORITY["EPSG","8903"]],
-         UNIT["grad",0.01570796326794897,
+         UNIT["grad",0.01570796326794897],
?                                        +
-             AUTHORITY["EPSG","9105"]],
          AUTHORITY["EPSG","4807"]],
      PROJECTION["Lambert_Conformal_Conic_1SP"],
-     PARAMETER["latitude_of_origin",55],
?                                    ^
+     PARAMETER["latitude_of_origin",49.5],
?                                    ^^^
      PARAMETER["central_meridian",0],
      PARAMETER["scale_factor",0.999877341],
      PARAMETER["false_easting",600000],
      PARAMETER["false_northing",1200000],
      UNIT["metre",1,
          AUTHORITY["EPSG","9001"]],
-     AXIS["X",EAST],
-     AXIS["Y",NORTH],
      AUTHORITY["EPSG","27571"]]
EPSG:27572 
Differences in PROJ.4:
- +proj=lcc +lat_1=46.8 +lat_0=46.8 +lon_0=0 +k_0=0.99987742 +x_0=600000 +y_0=2200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=paris +units=m +no_defs 
?                   ^ ^         ^ ^                                                                                                             ^^^^^
+ +proj=lcc +lat_1=42.12000000000001 +lat_0=42.12000000000001 +lon_0=0 +k_0=0.99987742 +x_0=600000 +y_0=2200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=2.5969213 +units=m +no_defs 
?                   ^ ^^^^^^^^^^^^^^         ^ ^^^^^^^^^^^^^^                                                                                                             ^^^^^^^^^
Differences in WKT:
  PROJCS["NTF (Paris) / Lambert zone II",
      GEOGCS["NTF (Paris)",
          DATUM["Nouvelle_Triangulation_Francaise_Paris",
              SPHEROID["Clarke 1880 (IGN)",6378249.2,293.4660212936265,
                  AUTHORITY["EPSG","7011"]],
              TOWGS84[-168,-60,320,0,0,0,0],
              AUTHORITY["EPSG","6807"]],
-         PRIMEM["Paris",2.33722917,
?                           ^^^^^^^
+         PRIMEM["Paris",2.5969213],
?                          ++++++ ^
-             AUTHORITY["EPSG","8903"]],
-         UNIT["grad",0.01570796326794897,
+         UNIT["grad",0.01570796326794897],
?                                        +
-             AUTHORITY["EPSG","9105"]],
          AUTHORITY["EPSG","4807"]],
      PROJECTION["Lambert_Conformal_Conic_1SP"],
-     PARAMETER["latitude_of_origin",52],
?                                    ^^
+     PARAMETER["latitude_of_origin",46.8],
?                                    ^^^^
      PARAMETER["central_meridian",0],
      PARAMETER["scale_factor",0.99987742],
      PARAMETER["false_easting",600000],
      PARAMETER["false_northing",2200000],
      UNIT["metre",1,
          AUTHORITY["EPSG","9001"]],
-     AXIS["X",EAST],
-     AXIS["Y",NORTH],
      AUTHORITY["EPSG","27572"]]
EPSG:27573 
Differences in PROJ.4:
- +proj=lcc +lat_1=44.10000000000001 +lat_0=44.10000000000001 +lon_0=0 +k_0=0.999877499 +x_0=600000 +y_0=3200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=paris +units=m +no_defs 
?                  ^^ ^^                    ^^ ^^                                                                                                                          ^^^^^
+ +proj=lcc +lat_1=39.69000000000001 +lat_0=39.69000000000001 +lon_0=0 +k_0=0.999877499 +x_0=600000 +y_0=3200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=2.5969213 +units=m +no_defs 
?                  ^^ ^^                    ^^ ^^                                                                                                                          ^^^^^^^^^
Differences in WKT:
  PROJCS["NTF (Paris) / Lambert zone III",
      GEOGCS["NTF (Paris)",
          DATUM["Nouvelle_Triangulation_Francaise_Paris",
              SPHEROID["Clarke 1880 (IGN)",6378249.2,293.4660212936265,
                  AUTHORITY["EPSG","7011"]],
              TOWGS84[-168,-60,320,0,0,0,0],
              AUTHORITY["EPSG","6807"]],
-         PRIMEM["Paris",2.33722917,
?                           ^^^^^^^
+         PRIMEM["Paris",2.5969213],
?                          ++++++ ^
-             AUTHORITY["EPSG","8903"]],
-         UNIT["grad",0.01570796326794897,
+         UNIT["grad",0.01570796326794897],
?                                        +
-             AUTHORITY["EPSG","9105"]],
          AUTHORITY["EPSG","4807"]],
      PROJECTION["Lambert_Conformal_Conic_1SP"],
-     PARAMETER["latitude_of_origin",49],
?                                     ^
+     PARAMETER["latitude_of_origin",44.1],
?                                     ^^^
      PARAMETER["central_meridian",0],
      PARAMETER["scale_factor",0.999877499],
      PARAMETER["false_easting",600000],
      PARAMETER["false_northing",3200000],
      UNIT["metre",1,
          AUTHORITY["EPSG","9001"]],
-     AXIS["X",EAST],
-     AXIS["Y",NORTH],
      AUTHORITY["EPSG","27573"]]
EPSG:27574 
Differences in PROJ.4:
- +proj=lcc +lat_1=42.16500000000001 +lat_0=42.16500000000001 +lon_0=0 +k_0=0.99994471 +x_0=234.358 +y_0=4185861.369 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=paris +units=m +no_defs 
?                   ^^^^ -----------         ^^^^ -----------                                                                                                                  ^^^^^
+ +proj=lcc +lat_1=37.9485 +lat_0=37.9485 +lon_0=0 +k_0=0.99994471 +x_0=234.358 +y_0=4185861.369 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=2.5969213 +units=m +no_defs 
?                  ++++ ^         ++++ ^                                                                                                                   ^^^^^^^^^
Differences in WKT:
  PROJCS["NTF (Paris) / Lambert zone IV",
      GEOGCS["NTF (Paris)",
          DATUM["Nouvelle_Triangulation_Francaise_Paris",
              SPHEROID["Clarke 1880 (IGN)",6378249.2,293.4660212936265,
                  AUTHORITY["EPSG","7011"]],
              TOWGS84[-168,-60,320,0,0,0,0],
              AUTHORITY["EPSG","6807"]],
-         PRIMEM["Paris",2.33722917,
?                           ^^^^^^^
+         PRIMEM["Paris",2.5969213],
?                          ++++++ ^
-             AUTHORITY["EPSG","8903"]],
-         UNIT["grad",0.01570796326794897,
+         UNIT["grad",0.01570796326794897],
?                                        +
-             AUTHORITY["EPSG","9105"]],
          AUTHORITY["EPSG","4807"]],
      PROJECTION["Lambert_Conformal_Conic_1SP"],
-     PARAMETER["latitude_of_origin",46.85],
?                                      --
+     PARAMETER["latitude_of_origin",42.165],
?                                     +++
      PARAMETER["central_meridian",0],
      PARAMETER["scale_factor",0.99994471],
      PARAMETER["false_easting",234.358],
      PARAMETER["false_northing",4185861.369],
      UNIT["metre",1,
          AUTHORITY["EPSG","9001"]],
-     AXIS["X",EAST],
-     AXIS["Y",NORTH],
      AUTHORITY["EPSG","27574"]]
EPSG:27581 
Differences in PROJ.4:
- +proj=lcc +lat_1=49.50000000000001 +lat_0=49.50000000000001 +lon_0=0 +k_0=0.999877341 +x_0=600000 +y_0=1200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=paris +units=m +no_defs 
?                   ^  ^^^^^^^^^^^^^         ^  ^^^^^^^^^^^^^                                                                                                              ^^^^^
+ +proj=lcc +lat_1=44.55 +lat_0=44.55 +lon_0=0 +k_0=0.999877341 +x_0=600000 +y_0=1200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=2.5969213 +units=m +no_defs 
?                   ^  ^         ^  ^                                                                                                              ^^^^^^^^^
Differences in WKT:
- PROJCS["NTF (Paris) / France I (deprecated)",
?                               -------------
+ PROJCS["NTF (Paris) / France I",
      GEOGCS["NTF (Paris)",
          DATUM["Nouvelle_Triangulation_Francaise_Paris",
              SPHEROID["Clarke 1880 (IGN)",6378249.2,293.4660212936265,
                  AUTHORITY["EPSG","7011"]],
              TOWGS84[-168,-60,320,0,0,0,0],
              AUTHORITY["EPSG","6807"]],
-         PRIMEM["Paris",2.33722917,
?                           ^^^^^^^
+         PRIMEM["Paris",2.5969213],
?                          ++++++ ^
-             AUTHORITY["EPSG","8903"]],
-         UNIT["grad",0.01570796326794897,
+         UNIT["grad",0.01570796326794897],
?                                        +
-             AUTHORITY["EPSG","9105"]],
          AUTHORITY["EPSG","4807"]],
      PROJECTION["Lambert_Conformal_Conic_1SP"],
-     PARAMETER["latitude_of_origin",55],
?                                    ^
+     PARAMETER["latitude_of_origin",49.5],
?                                    ^^^
      PARAMETER["central_meridian",0],
      PARAMETER["scale_factor",0.999877341],
      PARAMETER["false_easting",600000],
      PARAMETER["false_northing",1200000],
      UNIT["metre",1,
          AUTHORITY["EPSG","9001"]],
-     AXIS["X",EAST],
-     AXIS["Y",NORTH],
      AUTHORITY["EPSG","27581"]]
EPSG:27582 
Differences in PROJ.4:
- +proj=lcc +lat_1=46.8 +lat_0=46.8 +lon_0=0 +k_0=0.99987742 +x_0=600000 +y_0=2200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=paris +units=m +no_defs 
?                   ^ ^         ^ ^                                                                                                             ^^^^^
+ +proj=lcc +lat_1=42.12000000000001 +lat_0=42.12000000000001 +lon_0=0 +k_0=0.99987742 +x_0=600000 +y_0=2200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=2.5969213 +units=m +no_defs 
?                   ^ ^^^^^^^^^^^^^^         ^ ^^^^^^^^^^^^^^                                                                                                             ^^^^^^^^^
Differences in WKT:
- PROJCS["NTF (Paris) / France II (deprecated)",
?                                -------------
+ PROJCS["NTF (Paris) / France II",
      GEOGCS["NTF (Paris)",
          DATUM["Nouvelle_Triangulation_Francaise_Paris",
              SPHEROID["Clarke 1880 (IGN)",6378249.2,293.4660212936265,
                  AUTHORITY["EPSG","7011"]],
              TOWGS84[-168,-60,320,0,0,0,0],
              AUTHORITY["EPSG","6807"]],
-         PRIMEM["Paris",2.33722917,
?                           ^^^^^^^
+         PRIMEM["Paris",2.5969213],
?                          ++++++ ^
-             AUTHORITY["EPSG","8903"]],
-         UNIT["grad",0.01570796326794897,
+         UNIT["grad",0.01570796326794897],
?                                        +
-             AUTHORITY["EPSG","9105"]],
          AUTHORITY["EPSG","4807"]],
      PROJECTION["Lambert_Conformal_Conic_1SP"],
-     PARAMETER["latitude_of_origin",52],
?                                    ^^
+     PARAMETER["latitude_of_origin",46.8],
?                                    ^^^^
      PARAMETER["central_meridian",0],
      PARAMETER["scale_factor",0.99987742],
      PARAMETER["false_easting",600000],
      PARAMETER["false_northing",2200000],
      UNIT["metre",1,
          AUTHORITY["EPSG","9001"]],
-     AXIS["X",EAST],
-     AXIS["Y",NORTH],
      AUTHORITY["EPSG","27582"]]
EPSG:27583 
Differences in PROJ.4:
- +proj=lcc +lat_1=44.10000000000001 +lat_0=44.10000000000001 +lon_0=0 +k_0=0.999877499 +x_0=600000 +y_0=3200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=paris +units=m +no_defs 
?                  ^^ ^^                    ^^ ^^                                                                                                                          ^^^^^
+ +proj=lcc +lat_1=39.69000000000001 +lat_0=39.69000000000001 +lon_0=0 +k_0=0.999877499 +x_0=600000 +y_0=3200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=2.5969213 +units=m +no_defs 
?                  ^^ ^^                    ^^ ^^                                                                                                                          ^^^^^^^^^
Differences in WKT:
- PROJCS["NTF (Paris) / France III (deprecated)",
?                                 -------------
+ PROJCS["NTF (Paris) / France III",
      GEOGCS["NTF (Paris)",
          DATUM["Nouvelle_Triangulation_Francaise_Paris",
              SPHEROID["Clarke 1880 (IGN)",6378249.2,293.4660212936265,
                  AUTHORITY["EPSG","7011"]],
              TOWGS84[-168,-60,320,0,0,0,0],
              AUTHORITY["EPSG","6807"]],
-         PRIMEM["Paris",2.33722917,
?                           ^^^^^^^
+         PRIMEM["Paris",2.5969213],
?                          ++++++ ^
-             AUTHORITY["EPSG","8903"]],
-         UNIT["grad",0.01570796326794897,
+         UNIT["grad",0.01570796326794897],
?                                        +
-             AUTHORITY["EPSG","9105"]],
          AUTHORITY["EPSG","4807"]],
      PROJECTION["Lambert_Conformal_Conic_1SP"],
-     PARAMETER["latitude_of_origin",49],
?                                     ^
+     PARAMETER["latitude_of_origin",44.1],
?                                     ^^^
      PARAMETER["central_meridian",0],
      PARAMETER["scale_factor",0.999877499],
      PARAMETER["false_easting",600000],
      PARAMETER["false_northing",3200000],
      UNIT["metre",1,
          AUTHORITY["EPSG","9001"]],
-     AXIS["X",EAST],
-     AXIS["Y",NORTH],
      AUTHORITY["EPSG","27583"]]
EPSG:27584 
Differences in PROJ.4:
- +proj=lcc +lat_1=42.16500000000001 +lat_0=42.16500000000001 +lon_0=0 +k_0=0.99994471 +x_0=234.358 +y_0=4185861.369 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=paris +units=m +no_defs 
?                   ^^^^ -----------         ^^^^ -----------                                                                                                                  ^^^^^
+ +proj=lcc +lat_1=37.9485 +lat_0=37.9485 +lon_0=0 +k_0=0.99994471 +x_0=234.358 +y_0=4185861.369 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=2.5969213 +units=m +no_defs 
?                  ++++ ^         ++++ ^                                                                                                                   ^^^^^^^^^
Differences in WKT:
- PROJCS["NTF (Paris) / France IV (deprecated)",
?                                -------------
+ PROJCS["NTF (Paris) / France IV",
      GEOGCS["NTF (Paris)",
          DATUM["Nouvelle_Triangulation_Francaise_Paris",
              SPHEROID["Clarke 1880 (IGN)",6378249.2,293.4660212936265,
                  AUTHORITY["EPSG","7011"]],
              TOWGS84[-168,-60,320,0,0,0,0],
              AUTHORITY["EPSG","6807"]],
-         PRIMEM["Paris",2.33722917,
?                           ^^^^^^^
+         PRIMEM["Paris",2.5969213],
?                          ++++++ ^
-             AUTHORITY["EPSG","8903"]],
-         UNIT["grad",0.01570796326794897,
+         UNIT["grad",0.01570796326794897],
?                                        +
-             AUTHORITY["EPSG","9105"]],
          AUTHORITY["EPSG","4807"]],
      PROJECTION["Lambert_Conformal_Conic_1SP"],
-     PARAMETER["latitude_of_origin",46.85],
?                                      --
+     PARAMETER["latitude_of_origin",42.165],
?                                     +++
      PARAMETER["central_meridian",0],
      PARAMETER["scale_factor",0.99994471],
      PARAMETER["false_easting",234.358],
      PARAMETER["false_northing",4185861.369],
      UNIT["metre",1,
          AUTHORITY["EPSG","9001"]],
-     AXIS["X",EAST],
-     AXIS["Y",NORTH],
      AUTHORITY["EPSG","27584"]]
EPSG:27591 
Differences in PROJ.4:
- +proj=lcc +lat_1=49.50000000000001 +lat_0=49.50000000000001 +lon_0=0 +k_0=0.999877341 +x_0=600000 +y_0=200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=paris +units=m +no_defs 
?                   ^  ^^^^^^^^^^^^^         ^  ^^^^^^^^^^^^^                                                                                                             ^^^^^
+ +proj=lcc +lat_1=44.55 +lat_0=44.55 +lon_0=0 +k_0=0.999877341 +x_0=600000 +y_0=200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=2.5969213 +units=m +no_defs 
?                   ^  ^         ^  ^                                                                                                             ^^^^^^^^^
Differences in WKT:
- PROJCS["NTF (Paris) / Nord France (deprecated)",
?                                  -------------
+ PROJCS["NTF (Paris) / Nord France",
      GEOGCS["NTF (Paris)",
          DATUM["Nouvelle_Triangulation_Francaise_Paris",
              SPHEROID["Clarke 1880 (IGN)",6378249.2,293.4660212936265,
                  AUTHORITY["EPSG","7011"]],
              TOWGS84[-168,-60,320,0,0,0,0],
              AUTHORITY["EPSG","6807"]],
-         PRIMEM["Paris",2.33722917,
?                           ^^^^^^^
+         PRIMEM["Paris",2.5969213],
?                          ++++++ ^
-             AUTHORITY["EPSG","8903"]],
-         UNIT["grad",0.01570796326794897,
+         UNIT["grad",0.01570796326794897],
?                                        +
-             AUTHORITY["EPSG","9105"]],
          AUTHORITY["EPSG","4807"]],
      PROJECTION["Lambert_Conformal_Conic_1SP"],
-     PARAMETER["latitude_of_origin",55],
?                                    ^
+     PARAMETER["latitude_of_origin",49.5],
?                                    ^^^
      PARAMETER["central_meridian",0],
      PARAMETER["scale_factor",0.999877341],
      PARAMETER["false_easting",600000],
      PARAMETER["false_northing",200000],
      UNIT["metre",1,
          AUTHORITY["EPSG","9001"]],
-     AXIS["X",EAST],
-     AXIS["Y",NORTH],
      AUTHORITY["EPSG","27591"]]
EPSG:27592 
Differences in PROJ.4:
- +proj=lcc +lat_1=46.8 +lat_0=46.8 +lon_0=0 +k_0=0.99987742 +x_0=600000 +y_0=200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=paris +units=m +no_defs 
?                   ^ ^         ^ ^                                                                                                            ^^^^^
+ +proj=lcc +lat_1=42.12000000000001 +lat_0=42.12000000000001 +lon_0=0 +k_0=0.99987742 +x_0=600000 +y_0=200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=2.5969213 +units=m +no_defs 
?                   ^ ^^^^^^^^^^^^^^         ^ ^^^^^^^^^^^^^^                                                                                                            ^^^^^^^^^
Differences in WKT:
- PROJCS["NTF (Paris) / Centre France (deprecated)",
?                                    -------------
+ PROJCS["NTF (Paris) / Centre France",
      GEOGCS["NTF (Paris)",
          DATUM["Nouvelle_Triangulation_Francaise_Paris",
              SPHEROID["Clarke 1880 (IGN)",6378249.2,293.4660212936265,
                  AUTHORITY["EPSG","7011"]],
              TOWGS84[-168,-60,320,0,0,0,0],
              AUTHORITY["EPSG","6807"]],
-         PRIMEM["Paris",2.33722917,
?                           ^^^^^^^
+         PRIMEM["Paris",2.5969213],
?                          ++++++ ^
-             AUTHORITY["EPSG","8903"]],
-         UNIT["grad",0.01570796326794897,
+         UNIT["grad",0.01570796326794897],
?                                        +
-             AUTHORITY["EPSG","9105"]],
          AUTHORITY["EPSG","4807"]],
      PROJECTION["Lambert_Conformal_Conic_1SP"],
-     PARAMETER["latitude_of_origin",52],
?                                    ^^
+     PARAMETER["latitude_of_origin",46.8],
?                                    ^^^^
      PARAMETER["central_meridian",0],
      PARAMETER["scale_factor",0.99987742],
      PARAMETER["false_easting",600000],
      PARAMETER["false_northing",200000],
      UNIT["metre",1,
          AUTHORITY["EPSG","9001"]],
-     AXIS["X",EAST],
-     AXIS["Y",NORTH],
      AUTHORITY["EPSG","27592"]]
EPSG:27593 
Differences in PROJ.4:
- +proj=lcc +lat_1=44.10000000000001 +lat_0=44.10000000000001 +lon_0=0 +k_0=0.999877499 +x_0=600000 +y_0=200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=paris +units=m +no_defs 
?                  ^^ ^^                    ^^ ^^                                                                                                                         ^^^^^
+ +proj=lcc +lat_1=39.69000000000001 +lat_0=39.69000000000001 +lon_0=0 +k_0=0.999877499 +x_0=600000 +y_0=200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=2.5969213 +units=m +no_defs 
?                  ^^ ^^                    ^^ ^^                                                                                                                         ^^^^^^^^^
Differences in WKT:
- PROJCS["NTF (Paris) / Sud France (deprecated)",
?                                 -------------
+ PROJCS["NTF (Paris) / Sud France",
      GEOGCS["NTF (Paris)",
          DATUM["Nouvelle_Triangulation_Francaise_Paris",
              SPHEROID["Clarke 1880 (IGN)",6378249.2,293.4660212936265,
                  AUTHORITY["EPSG","7011"]],
              TOWGS84[-168,-60,320,0,0,0,0],
              AUTHORITY["EPSG","6807"]],
-         PRIMEM["Paris",2.33722917,
?                           ^^^^^^^
+         PRIMEM["Paris",2.5969213],
?                          ++++++ ^
-             AUTHORITY["EPSG","8903"]],
-         UNIT["grad",0.01570796326794897,
+         UNIT["grad",0.01570796326794897],
?                                        +
-             AUTHORITY["EPSG","9105"]],
          AUTHORITY["EPSG","4807"]],
      PROJECTION["Lambert_Conformal_Conic_1SP"],
-     PARAMETER["latitude_of_origin",49],
?                                     ^
+     PARAMETER["latitude_of_origin",44.1],
?                                     ^^^
      PARAMETER["central_meridian",0],
      PARAMETER["scale_factor",0.999877499],
      PARAMETER["false_easting",600000],
      PARAMETER["false_northing",200000],
      UNIT["metre",1,
          AUTHORITY["EPSG","9001"]],
-     AXIS["X",EAST],
-     AXIS["Y",NORTH],
      AUTHORITY["EPSG","27593"]]
EPSG:27594 
Differences in PROJ.4:
- +proj=lcc +lat_1=42.16500000000001 +lat_0=42.16500000000001 +lon_0=0 +k_0=0.99994471 +x_0=234.358 +y_0=185861.369 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=paris +units=m +no_defs 
?                   ^^^^ -----------         ^^^^ -----------                                                                                                                 ^^^^^
+ +proj=lcc +lat_1=37.9485 +lat_0=37.9485 +lon_0=0 +k_0=0.99994471 +x_0=234.358 +y_0=185861.369 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=2.5969213 +units=m +no_defs 
?                  ++++ ^         ++++ ^                                                                                                                  ^^^^^^^^^
Differences in WKT:
- PROJCS["NTF (Paris) / Corse (deprecated)",
?                            -------------
+ PROJCS["NTF (Paris) / Corse",
      GEOGCS["NTF (Paris)",
          DATUM["Nouvelle_Triangulation_Francaise_Paris",
              SPHEROID["Clarke 1880 (IGN)",6378249.2,293.4660212936265,
                  AUTHORITY["EPSG","7011"]],
              TOWGS84[-168,-60,320,0,0,0,0],
              AUTHORITY["EPSG","6807"]],
-         PRIMEM["Paris",2.33722917,
?                           ^^^^^^^
+         PRIMEM["Paris",2.5969213],
?                          ++++++ ^
-             AUTHORITY["EPSG","8903"]],
-         UNIT["grad",0.01570796326794897,
+         UNIT["grad",0.01570796326794897],
?                                        +
-             AUTHORITY["EPSG","9105"]],
          AUTHORITY["EPSG","4807"]],
      PROJECTION["Lambert_Conformal_Conic_1SP"],
-     PARAMETER["latitude_of_origin",46.85],
?                                      --
+     PARAMETER["latitude_of_origin",42.165],
?                                     +++
      PARAMETER["central_meridian",0],
      PARAMETER["scale_factor",0.99994471],
      PARAMETER["false_easting",234.358],
      PARAMETER["false_northing",185861.369],
      UNIT["metre",1,
          AUTHORITY["EPSG","9001"]],
-     AXIS["X",EAST],
-     AXIS["Y",NORTH],
      AUTHORITY["EPSG","27594"]]

Attachments (1)

5899.patch (24.1 KB ) - added by Even Rouault 9 years ago.

Download all attachments as: .zip

Change History (6)

comment:1 by Even Rouault, 9 years ago

The issue here is related to a UNIT = GRAD. Not specificaly with LCC

comment:2 by Even Rouault, 9 years ago

Summary: Lambert Conformal Conic projections not preserved by GTiffSRS with UNIT = GRAD not preserved by GTiff

Similar ticket #5902

comment:3 by Even Rouault, 9 years ago

And #5898

comment:4 by Even Rouault, 9 years ago

Attached a patch that attempts to resolve the issue. But there are some unknowns relative how to interprete the PRIMEM in GDAL WKT. I believe that currently it is in degrees, whatever the angular unit.

by Even Rouault, 9 years ago

Attachment: 5899.patch added

comment:5 by Even Rouault, 5 years ago

Milestone: closed_because_of_github_migration
Resolution: wontfix
Status: newclosed

This ticket has been automatically closed because Trac is no longer used for GDAL bug tracking, since the project has migrated to GitHub. If you believe this ticket is still valid, you may file it to https://github.com/OSGeo/gdal/issues if it is not already reported there.

Note: See TracTickets for help on using tickets.