Opened 9 years ago
Closed 5 years ago
#5901 closed defect (wontfix)
Latitude of origin vs center for LAEA, AEA and AEQD projections
| 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
LAEA codes EPSG:42106 and EPSG:102017, AEA codes EPSG:42303 and EPSG:42308, and AEQD codes ESRI:102016, ESRI:102019 and ESRI:102020 have different PROJ.4 and WKT from the GTiff driver.
There is a possible confusion between latitude_of_origin and latitude_of_center in WKT, and that +lat_0 values are not preserved in the GTiff file.
See #5897 for test_prj_roundtrip, and run:
for code in [42106, 102017, 42303, 42308, 102016, 102019, 102020]:
print(test_prj_roundtrip(code, True, True))
with output
Differences in PROJ.4:
- +proj=laea +lat_0=5 +lon_0=20 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs
? ^ -
+ +proj=laea +lat_0=0 +lon_0=0 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs
? ^
Differences in WKT:
PROJCS["WGS84 / Lambert Azim Mozambique",
GEOGCS["WGS 84",
DATUM["WGS_1984",
- SPHEROID["Sphere_radius_6370997_m",6370997,0]],
+ SPHEROID["WGS 84",6370997,447.2543030894074,
+ AUTHORITY["EPSG","7030"]],
+ AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
- UNIT["Decimal_Degree",0.0174532925199433]],
? ^^^^^^^^^ -
+ UNIT["degree",0.0174532925199433],
? ^
+ AUTHORITY["EPSG","4326"]],
- PROJECTION["Lambert_Azimuthal_equal_area"],
? ^ ^
+ PROJECTION["Lambert_Azimuthal_Equal_Area"],
? ^ ^
- PARAMETER["latitude_of_origin",5],
? ^ ---- ^
+ PARAMETER["latitude_of_center",0],
? ^^^^^ ^
+ PARAMETER["longitude_of_center",0],
- PARAMETER["central_meridian",20],
- PARAMETER["standard_parallel_1",5],
- PARAMETER["standard_parallel_2",5],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
- UNIT["Meter",1],
? ^ - -
+ UNIT["metre",1,
? ^ +
- AUTHORITY["EPSG","42106"]]? ^^ --
+ AUTHORITY["EPSG","9001"]]]? ++++ ^^^ +
EPSG:102017
Differences in PROJ.4:
- +proj=laea +lat_0=90 +lon_0=0 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs
? -
+ +proj=laea +lat_0=0 +lon_0=0 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs
Differences in WKT:
PROJCS["North_Pole_Lambert_Azimuthal_Equal_Area",
- GEOGCS["GCS_WGS_1984",
? ---- ^^^
+ GEOGCS["WGS 84",
? ^
DATUM["WGS_1984",
- SPHEROID["WGS_1984",6378137,298.257223563]],
? ^^^ --
+ SPHEROID["WGS 84",6378137,298.257223563,
? ^
+ AUTHORITY["EPSG","7030"]],
+ AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
- UNIT["Degree",0.017453292519943295]],
? ^ ^^^^
+ UNIT["degree",0.0174532925199433],
? ^ ^
+ AUTHORITY["EPSG","4326"]],
PROJECTION["Lambert_Azimuthal_Equal_Area"],
+ PARAMETER["latitude_of_center",0],
+ PARAMETER["longitude_of_center",0],
- PARAMETER["False_Easting",0],
? ^ ^
+ PARAMETER["false_easting",0],
? ^ ^
- PARAMETER["False_Northing",0],
? ^ ^
+ PARAMETER["false_northing",0],
? ^ ^
- PARAMETER["Central_Meridian",0],
- PARAMETER["Latitude_Of_Origin",90],
- UNIT["Meter",1],
? ^ - -
+ UNIT["metre",1,
? ^ +
- AUTHORITY["EPSG","102017"]]? ^ - -
+ AUTHORITY["EPSG","9001"]]]? ++++ ^ +
EPSG:42303
Differences in PROJ.4:
- +proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=23 +lon_0=-96 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs
? ^^ ^^^
+ +proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=0 +lon_0=0 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs
? ^ ^
Differences in WKT:
PROJCS["NAD83 / Albers NorthAm",
GEOGCS["NAD83",
DATUM["North_American_Datum_1983",
- SPHEROID["GRS_1980",6378137,298.257222101]],
? ^ ^^
+ SPHEROID["GRS 1980",6378137,298.2572221010002,
? ^ ^^^^
+ AUTHORITY["EPSG","7019"]],
+ AUTHORITY["EPSG","6269"]],
PRIMEM["Greenwich",0],
- UNIT["Decimal_Degree",0.0174532925199433]],
? ^^^^^^^^^ -
+ UNIT["degree",0.0174532925199433],
? ^
+ AUTHORITY["EPSG","4269"]],
- PROJECTION["Albers_conic_equal_area"],
? ^ ^ ^
+ PROJECTION["Albers_Conic_Equal_Area"],
? ^ ^ ^
- PARAMETER["central_meridian",-96.0],
- PARAMETER["latitude_of_origin",23],
PARAMETER["standard_parallel_1",29.5],
PARAMETER["standard_parallel_2",45.5],
+ PARAMETER["latitude_of_center",0],
+ PARAMETER["longitude_of_center",0],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
- UNIT["Meter",1],
? ^ - -
+ UNIT["metre",1,
? ^ +
- AUTHORITY["EPSG","42303"]]? ^^^ ^
+ AUTHORITY["EPSG","9001"]]]? ++++ ^ ^^ +
EPSG:42308
Differences in PROJ.4:
- +proj=aea +lat_1=34 +lat_2=40.5 +lat_0=0 +lon_0=-120 +x_0=0 +y_0=-4000000 +datum=NAD27 +units=m +no_defs
? ---
+ +proj=aea +lat_1=34 +lat_2=40.5 +lat_0=0 +lon_0=0 +x_0=0 +y_0=-4000000 +datum=NAD27 +units=m +no_defs
Differences in WKT:
PROJCS["NAD27 / California Albers",
GEOGCS["NAD27",
DATUM["North_American_Datum_1927",
- SPHEROID["Clarke 1866",6378206.4,294.978698213901]],
? ^^^
+ SPHEROID["Clarke 1866",6378206.4,294.9786982139006,
? ^^
+ AUTHORITY["EPSG","7008"]],
+ AUTHORITY["EPSG","6267"]],
PRIMEM["Greenwich",0],
- UNIT["degree",0.0174532925199433]],
? -
+ UNIT["degree",0.0174532925199433],
+ AUTHORITY["EPSG","4267"]],
- PROJECTION["Albers_conic_equal_area"],
? ^ ^ ^
+ PROJECTION["Albers_Conic_Equal_Area"],
? ^ ^ ^
- PARAMETER["central_meridian",-120.0],
- PARAMETER["latitude_of_origin",0],
PARAMETER["standard_parallel_1",34],
PARAMETER["standard_parallel_2",40.5],
+ PARAMETER["latitude_of_center",0],
+ PARAMETER["longitude_of_center",0],
PARAMETER["false_easting",0],
PARAMETER["false_northing",-4000000],
- UNIT["Meter",1],
? ^ - -
+ UNIT["metre",1,
? ^ +
- AUTHORITY["EPSG","42308"]]? ^^^ ^
+ AUTHORITY["EPSG","9001"]]]? ++++ ^ ^^ +
EPSG:102016
Differences in PROJ.4:
- +proj=aeqd +lat_0=90 +lon_0=0 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs
? -
+ +proj=aeqd +lat_0=0 +lon_0=0 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs
Differences in WKT:
PROJCS["North_Pole_Azimuthal_Equidistant",
- GEOGCS["GCS_WGS_1984",
? ---- ^^^
+ GEOGCS["WGS 84",
? ^
DATUM["WGS_1984",
- SPHEROID["WGS_1984",6378137,298.257223563]],
? ^^^ --
+ SPHEROID["WGS 84",6378137,298.257223563,
? ^
+ AUTHORITY["EPSG","7030"]],
+ AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
- UNIT["Degree",0.017453292519943295]],
? ^ ^^^^
+ UNIT["degree",0.0174532925199433],
? ^ ^
+ AUTHORITY["EPSG","4326"]],
PROJECTION["Azimuthal_Equidistant"],
+ PARAMETER["latitude_of_center",0],
+ PARAMETER["longitude_of_center",0],
- PARAMETER["False_Easting",0],
? ^ ^
+ PARAMETER["false_easting",0],
? ^ ^
- PARAMETER["False_Northing",0],
? ^ ^
+ PARAMETER["false_northing",0],
? ^ ^
- PARAMETER["Central_Meridian",0],
- PARAMETER["Latitude_Of_Origin",90],
- UNIT["Meter",1],
? ^ - -
+ UNIT["metre",1,
? ^ +
- AUTHORITY["EPSG","102016"]]? ^ - -
+ AUTHORITY["EPSG","9001"]]]? ++++ ^ +
EPSG:102019
Differences in PROJ.4:
- +proj=aeqd +lat_0=-90 +lon_0=0 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs
? --
+ +proj=aeqd +lat_0=0 +lon_0=0 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs
Differences in WKT:
PROJCS["South_Pole_Azimuthal_Equidistant",
- GEOGCS["GCS_WGS_1984",
? ---- ^^^
+ GEOGCS["WGS 84",
? ^
DATUM["WGS_1984",
- SPHEROID["WGS_1984",6378137,298.257223563]],
? ^^^ --
+ SPHEROID["WGS 84",6378137,298.257223563,
? ^
+ AUTHORITY["EPSG","7030"]],
+ AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
- UNIT["Degree",0.017453292519943295]],
? ^ ^^^^
+ UNIT["degree",0.0174532925199433],
? ^ ^
+ AUTHORITY["EPSG","4326"]],
PROJECTION["Azimuthal_Equidistant"],
+ PARAMETER["latitude_of_center",0],
+ PARAMETER["longitude_of_center",0],
- PARAMETER["False_Easting",0],
? ^ ^
+ PARAMETER["false_easting",0],
? ^ ^
- PARAMETER["False_Northing",0],
? ^ ^
+ PARAMETER["false_northing",0],
? ^ ^
- PARAMETER["Central_Meridian",0],
- PARAMETER["Latitude_Of_Origin",-90],
- UNIT["Meter",1],
? ^ - -
+ UNIT["metre",1,
? ^ +
- AUTHORITY["EPSG","102019"]]? ^ - -
+ AUTHORITY["EPSG","9001"]]]? ++++ ^ +
EPSG:102020
Differences in PROJ.4:
- +proj=laea +lat_0=-90 +lon_0=0 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs
? --
+ +proj=laea +lat_0=0 +lon_0=0 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs
Differences in WKT:
PROJCS["South_Pole_Lambert_Azimuthal_Equal_Area",
- GEOGCS["GCS_WGS_1984",
? ---- ^^^
+ GEOGCS["WGS 84",
? ^
DATUM["WGS_1984",
- SPHEROID["WGS_1984",6378137,298.257223563]],
? ^^^ --
+ SPHEROID["WGS 84",6378137,298.257223563,
? ^
+ AUTHORITY["EPSG","7030"]],
+ AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
- UNIT["Degree",0.017453292519943295]],
? ^ ^^^^
+ UNIT["degree",0.0174532925199433],
? ^ ^
+ AUTHORITY["EPSG","4326"]],
PROJECTION["Lambert_Azimuthal_Equal_Area"],
+ PARAMETER["latitude_of_center",0],
+ PARAMETER["longitude_of_center",0],
- PARAMETER["False_Easting",0],
? ^ ^
+ PARAMETER["false_easting",0],
? ^ ^
- PARAMETER["False_Northing",0],
? ^ ^
+ PARAMETER["false_northing",0],
? ^ ^
- PARAMETER["Central_Meridian",0],
- PARAMETER["Latitude_Of_Origin",-90],
- UNIT["Meter",1],
? ^ - -
+ UNIT["metre",1,
? ^ +
- AUTHORITY["EPSG","102020"]]? -----
+ AUTHORITY["EPSG","9001"]]]? ++++ +++ +
Change History (2)
comment:1 by , 9 years ago
comment:2 by , 5 years ago
| Milestone: | → closed_because_of_github_migration |
|---|---|
| Resolution: | → wontfix |
| Status: | new → closed |
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.
The issue is not really in the GeoTIFF driver but more in the WKT of those "EPSG" codes that improperly use Latitude_of_origin instead of the 'official' Latitude_of_center (http://www.remotesensing.org/geotiff/proj_list/lambert_azimuthal_equal_area.html)