GeoTools

Geotools-generated prj file incompatible with ESRI ArcMap

Details

  • Type: Bug Bug
  • Status: Open Open
  • Priority: Minor Minor
  • Resolution: Unresolved
  • Affects Version/s: 2.4.4
  • Fix Version/s: 2.7.6
  • Component/s: referencing, shapefile plugin
  • Labels:
    None
  • Environment:
    Linux OS, JDK 1.6

Description

The prj file provided by the ArcMap distribution for WGS 84 / UTM Zone 35S is this (line breaks added for easier reading):

PROJCS["WGS_1984_UTM_Zone_35S",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137,298.257223563]],
PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["Central_Meridian",27],
PARAMETER["Latitude_Of_Origin",0],PARAMETER["Scale_Factor",0.9996],PARAMETER["False_Easting",500000],PARAMETER["False_Northing",10000000],
UNIT["Meter",1]

Geotools ShapefileDataStore generates this prj file when given a CRS with ShapefileDataStore.forceSchemaCRS():

PROJCS["WGS 84 / UTM zone 35S",
  GEOGCS["WGS 84",
    DATUM["World Geodetic System 1984",
      SPHEROID["WGS 84", 6378137.0, 298.257223563, AUTHORITY["EPSG","7030"]],
      AUTHORITY["EPSG","6326"]],
    PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]],
    UNIT["degree", 0.017453292519943295],
    AXIS["Geodetic latitude", NORTH],
    AXIS["Geodetic longitude", EAST],
    AUTHORITY["EPSG","4326"]],
  PROJECTION["Transverse Mercator", AUTHORITY["EPSG","9807"]],
  PARAMETER["central_meridian", 27.0],
  PARAMETER["latitude_of_origin", 0.0],
  PARAMETER["scale_factor", 0.9996],
  PARAMETER["false_easting", 500000.0],
  PARAMETER["false_northing", 10000000.0],
  UNIT["m", 1.0],
  AXIS["Easting", EAST],
  AXIS["Northing", NORTH],
  AUTHORITY["EPSG","32735"]]

Among other dissimilarities, ArcMap doesn't recognize the Geotools-generated prj definition because units are "m" instead of "Meter", and the datum name is "World Geodetic System 1984" versus "D_WGS_1984". ArcMap 9.2 shows the units as "Unknown", and complains about the projection being different when the shapefile is layered on other shapefiles with the same projection.

See GEOS-830 for a similar complaint.

Issue Links

depends upon

Task - A task that needs to be done. GEOT-2117 Review GeoToolkit fixes

  • Minor - Minor loss of function, or other problem where easy workaround is present.
  • Open - The issue is open and ready for the assignee to start work on it.
is depended upon by

Bug - A problem which impairs or prevents the functions of the product. GEOS-830 shapefile .prj output should not have spaces

  • Major - Major loss of function.
  • Closed - The issue is considered finished, the resolution is correct. Issues which are not closed can be reopened.
is superceded by

Improvement - An improvement or enhancement to an existing feature or task. GEOS-4503 Make SHAPE-ZIP WFS output format capable of sending .prj files in ESRI WKT format

  • Major - Major loss of function.
  • Closed - The issue is considered finished, the resolution is correct. Issues which are not closed can be reopened.
relates to

Improvement - An improvement or enhancement to an existing feature or task. GEOT-3571 Generate prj files with ESRI friendly nomenclature

  • Major - Major loss of function.
  • Open - The issue is open and ready for the assignee to start work on it.

Activity

Ascending order - Click to sort in descending order
Hide
Permalink
Martin Desruisseaux added a comment -
Geotidy now formats "m" as "meter" if the authority is ESRI or ORACLE, and "metre" otherwise. The list of authorities for which to format as "meter" contains for now ESRI based on this report, and ORACLE based on the example provided in GEOT-1164. More authorities may be added to this list as we get feedback (note that "metre" is the international name while "meter" is US-specific).

The Datum should now formats using the ESRI name if the CRS is formatted using the toWKT(Citations.ESRI, ...) method call. I would expect this call to work in GeoTools 2.x, but I have not verified. However I verified that it works in Geotidy and added a test case for ensuring that.
Show
Martin Desruisseaux added a comment - Geotidy now formats "m" as "meter" if the authority is ESRI or ORACLE, and "metre" otherwise. The list of authorities for which to format as "meter" contains for now ESRI based on this report, and ORACLE based on the example provided in GEOT-1164 . More authorities may be added to this list as we get feedback (note that "metre" is the international name while "meter" is US-specific). The Datum should now formats using the ESRI name if the CRS is formatted using the toWKT(Citations.ESRI, ...) method call. I would expect this call to work in GeoTools 2.x, but I have not verified. However I verified that it works in Geotidy and added a test case for ensuring that.
Hide
Permalink
Andrea Aime added a comment -
After the patch in GEOT-3571 GeoTools can generate the following representation:

{panel}
PROJCS["WGS 84 / UTM zone 35S", GEOGCS["WGS 84", DATUM["D_WGS_1984", SPHEROID["D_WGS_1984", 6378137.0, 298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]], UNIT["Degree", 0.017453292519943295], AXIS["Geodetic longitude", EAST], AXIS["Geodetic latitude", NORTH], AUTHORITY["EPSG","4326"]], PROJECTION["Transverse_Mercator", AUTHORITY["EPSG","9807"]], PARAMETER["Longitude_Of_Center", 27.0], PARAMETER["Latitude_Of_Origin", 0.0], PARAMETER["Scale_Factor", 0.9996], PARAMETER["false_easting", 500000.0], PARAMETER["False_Northing", 10000000.0], UNIT["Meter", 1.0], AXIS["Easting", EAST], AXIS["Northing", NORTH], AUTHORITY["EPSG","32735"]]
{panel}

which is as close as I could get it, in particular not only the unit but also the projection name and parameter name appear to match. What is missing is a flag in the store to enable an "ESRI compatibility mode"
Show
Andrea Aime added a comment - After the patch in GEOT-3571 GeoTools can generate the following representation: {panel} PROJCS["WGS 84 / UTM zone 35S", GEOGCS["WGS 84", DATUM["D_WGS_1984", SPHEROID["D_WGS_1984", 6378137.0, 298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]], UNIT["Degree", 0.017453292519943295], AXIS["Geodetic longitude", EAST], AXIS["Geodetic latitude", NORTH], AUTHORITY["EPSG","4326"]], PROJECTION["Transverse_Mercator", AUTHORITY["EPSG","9807"]], PARAMETER["Longitude_Of_Center", 27.0], PARAMETER["Latitude_Of_Origin", 0.0], PARAMETER["Scale_Factor", 0.9996], PARAMETER["false_easting", 500000.0], PARAMETER["False_Northing", 10000000.0], UNIT["Meter", 1.0], AXIS["Easting", EAST], AXIS["Northing", NORTH], AUTHORITY["EPSG","32735"]] {panel} which is as close as I could get it, in particular not only the unit but also the projection name and parameter name appear to match. What is missing is a flag in the store to enable an "ESRI compatibility mode"

People

  • Assignee:
    Andrea Aime
    Reporter:
    Larry Reeder
Vote (0)
Watch (0)

Dates

  • Created:
    Updated: