Trouble using CRS from ESRI ArcGIS 9.3 (esri:102003)

View: New views
3 Messages — Rating Filter:   Alert me  

Trouble using CRS from ESRI ArcGIS 9.3 (esri:102003)

by Benjamin Wright :: Rate this Message:

Reply to Author | View Threaded | Show Only this Message

I'm having trouble using a CRS I'm getting from an ESRI ArcGIS 9.3 server. (MapServer Service)

The error I get is:

Caused by: org.opengis.referencing.operation.OperationNotFoundException: Bursa wolf parameters required.
        at org.geotools.referencing.operation.DefaultCoordinateOperationFactory.createOperationStep(DefaultCoordinateOperationFactory.java:1072)
        at org.geotools.referencing.operation.DefaultCoordinateOperationFactory.createOperationStep(DefaultCoordinateOperationFactory.java:1136)
        at org.geotools.referencing.operation.DefaultCoordinateOperationFactory.createOperationStep(DefaultCoordinateOperationFactory.java:881)
        at org.geotools.referencing.operation.DefaultCoordinateOperationFactory.createOperationStep(DefaultCoordinateOperationFactory.java:998)
        at org.geotools.referencing.operation.DefaultCoordinateOperationFactory.createOperation(DefaultCoordinateOperationFactory.java:267)
        at org.geotools.referencing.operation.BufferedCoordinateOperationFactory.createOperation(BufferedCoordinateOperationFactory.java:254)
        at org.geotools.referencing.CRS.findMathTransform(CRS.java:1001)
        at org.geotools.referencing.CRS.findMathTransform(CRS.java:969)

The call made into this is:
CRS.findMathTransform(CRS_IN_QUESTION,DefaultGeographicCRS.WGS84);

I get the CRS from the WKT:
PROJCS["USA_Contiguous_Albers_Equal_Area_Conic",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Albers"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-96.0],PARAMETER["Standard_Parallel_1",29.5],PARAMETER["Standard_Parallel_2",45.5],PARAMETER["Latitude_Of_Origin",37.5],UNIT["Meter",1.0],AUTHORITY["ESRI","102003"]]

Spatial Reference:
http://spatialreference.org/ref/esri/102003/

I'm using Geotools version 2.5.2

What are 'Bursa wolf parameters'?  Do I need them?  Any ideas?  Should this be fixed if I upgrade to 2.6? (working across a few companies - upgrading isn't as easy as it would sound but it is possible)

Re: Trouble using CRS from ESRI ArcGIS 9.3 (esri:102003)

by Andrea Aime-4 :: Rate this Message:

Reply to Author | View Threaded | Show Only this Message

Benjamin Wright ha scritto:

> I'm having trouble using a CRS I'm getting from an ESRI ArcGIS 9.3 server.
> (MapServer Service)
>
> The error I get is:
>
> Caused by: org.opengis.referencing.operation.OperationNotFoundException:
> Bursa wolf parameters required.
> at
> org.geotools.referencing.operation.DefaultCoordinateOperationFactory.createOperationStep(DefaultCoordinateOperationFactory.java:1072)
> at
> org.geotools.referencing.operation.DefaultCoordinateOperationFactory.createOperationStep(DefaultCoordinateOperationFactory.java:1136)
> at
> org.geotools.referencing.operation.DefaultCoordinateOperationFactory.createOperationStep(DefaultCoordinateOperationFactory.java:881)
> at
> org.geotools.referencing.operation.DefaultCoordinateOperationFactory.createOperationStep(DefaultCoordinateOperationFactory.java:998)
> at
> org.geotools.referencing.operation.DefaultCoordinateOperationFactory.createOperation(DefaultCoordinateOperationFactory.java:267)
> at
> org.geotools.referencing.operation.BufferedCoordinateOperationFactory.createOperation(BufferedCoordinateOperationFactory.java:254)
> at org.geotools.referencing.CRS.findMathTransform(CRS.java:1001)
> at org.geotools.referencing.CRS.findMathTransform(CRS.java:969)
>
> The call made into this is:
> CRS.findMathTransform(CRS_IN_QUESTION,DefaultGeographicCRS.WGS84);

Just use
CRS.findMathTransform(CRS_IN_QUESTION,DefaultGeographicCRS.WGS84, true);

instead

> I get the CRS from the WKT:
> PROJCS["USA_Contiguous_Albers_Equal_Area_Conic",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Albers"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-96.0],PARAMETER["Standard_Parallel_1",29.5],PARAMETER["Standard_Parallel_2",45.5],PARAMETER["Latitude_Of_Origin",37.5],UNIT["Meter",1.0],AUTHORITY["ESRI","102003"]]
>
> Spatial Reference:
> http://spatialreference.org/ref/esri/102003/
>
> I'm using Geotools version 2.5.2
>
> What are 'Bursa wolf parameters'?  Do I need them?  

The Bursa Wolf parameter are required to get decent precision (4m
usually) in datum shifts, without them you might get a transformation
error around 100m easily.
In the case of NAD83 -> WGS84 the error should not be so big without
them, but I've never actually measured it

Cheers
Andrea

--
Andrea Aime
OpenGeo - http://opengeo.org
Expert service straight from the developers.

---------------------------------------------------------------------
To unsubscribe from this list, please visit:

    http://xircles.codehaus.org/manage_email



Re: Trouble using CRS from ESRI ArcGIS 9.3 (esri:102003)

by Benjamin Wright :: Rate this Message:

Reply to Author | View Threaded | Show Only this Message

Thanks Andrea, that fixed the problem.

I see what you mean with the inaccurate transformations though.  With this particular transformation we don't see a significant problem, but with another one that was not working we're seeing latitude shifted by several kilometers (in the Eastern US / New England Area).  Could be a problem with what we're doing or with the transformation - we'll have to look into it further.

-Ben Wright


Andrea Aime-4 wrote:
Benjamin Wright ha scritto:
> I'm having trouble using a CRS I'm getting from an ESRI ArcGIS 9.3 server.
> (MapServer Service)
>
> The error I get is:
>
> Caused by: org.opengis.referencing.operation.OperationNotFoundException:
> Bursa wolf parameters required.
> at
> org.geotools.referencing.operation.DefaultCoordinateOperationFactory.createOperationStep(DefaultCoordinateOperationFactory.java:1072)
> at
> org.geotools.referencing.operation.DefaultCoordinateOperationFactory.createOperationStep(DefaultCoordinateOperationFactory.java:1136)
> at
> org.geotools.referencing.operation.DefaultCoordinateOperationFactory.createOperationStep(DefaultCoordinateOperationFactory.java:881)
> at
> org.geotools.referencing.operation.DefaultCoordinateOperationFactory.createOperationStep(DefaultCoordinateOperationFactory.java:998)
> at
> org.geotools.referencing.operation.DefaultCoordinateOperationFactory.createOperation(DefaultCoordinateOperationFactory.java:267)
> at
> org.geotools.referencing.operation.BufferedCoordinateOperationFactory.createOperation(BufferedCoordinateOperationFactory.java:254)
> at org.geotools.referencing.CRS.findMathTransform(CRS.java:1001)
> at org.geotools.referencing.CRS.findMathTransform(CRS.java:969)
>
> The call made into this is:
> CRS.findMathTransform(CRS_IN_QUESTION,DefaultGeographicCRS.WGS84);

Just use
CRS.findMathTransform(CRS_IN_QUESTION,DefaultGeographicCRS.WGS84, true);

instead

> I get the CRS from the WKT:
> PROJCS["USA_Contiguous_Albers_Equal_Area_Conic",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Albers"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-96.0],PARAMETER["Standard_Parallel_1",29.5],PARAMETER["Standard_Parallel_2",45.5],PARAMETER["Latitude_Of_Origin",37.5],UNIT["Meter",1.0],AUTHORITY["ESRI","102003"]]
>
> Spatial Reference:
> http://spatialreference.org/ref/esri/102003/
>
> I'm using Geotools version 2.5.2
>
> What are 'Bursa wolf parameters'?  Do I need them?  

The Bursa Wolf parameter are required to get decent precision (4m
usually) in datum shifts, without them you might get a transformation
error around 100m easily.
In the case of NAD83 -> WGS84 the error should not be so big without
them, but I've never actually measured it

Cheers
Andrea

--
Andrea Aime
OpenGeo - http://opengeo.org
Expert service straight from the developers.

---------------------------------------------------------------------
To unsubscribe from this list, please visit:

    http://xircles.codehaus.org/manage_email