Trying to remember the details; dropping an axis (to take things down to 2d) is
a bit of a trick as you need to determine which axis is "vertical". There was
some fancy word for it in the referencing module. There was also the idea of
using a Feature 2d hint when you perform the query. It was specifically done in
order to map 3d datasets.
Andrea will remember the details; or you can search the email archive.
Sorry I cannot be more help.
--
Jody Garnett
On Tuesday, 15 February 2011 at 1:14 AM, Farrukh Najmi wrote:
>
> I am using a PostGres/PostGIS 8.3 database with gt-epsg-postgresql
> plugin. The EPSG Geodetic Parameter Dataset has been loaded into
> postgres from EPSG version 7.6.
>
> When I try to transform a geometry from EPSG:4326 (2D CRS) to EPSG:4979
> (3D CRS) I get the following:
>
> Caused by: org.opengis.referencing.operation.OperationNotFoundException:
> No transformation available from system "EllipsoidalCS[Ellipsoidal 2D
> CS. Axes: latitude, longitude. Orientations: north, east. UoM: degree]"
> to "EllipsoidalCS[Ellipsoidal 3D CS. Axes: latitude, longitude,
> ellipsoidal height. Orientations: north, east, up. UoM: degree, degree,
> metre.]".
> at
> org.geotools.referencing.operation.AbstractCoordinateOperationFactory.swapAndScaleAxis(AbstractCoordinateOperationFactory.java:278)
> at
> org.geotools.referencing.operation.DefaultCoordinateOperationFactory.swapAndScaleAxis(DefaultCoordinateOperationFactory.java:591)
> at
> org.geotools.referencing.operation.DefaultCoordinateOperationFactory.createOperationStep(DefaultCoordinateOperationFactory.java:780)
> at
> org.geotools.referencing.operation.DefaultCoordinateOperationFactory.createOperation(DefaultCoordinateOperationFactory.java:239)
> at
> org.geotools.referencing.operation.BufferedCoordinateOperationFactory.createOperation(BufferedCoordinateOperationFactory.java:254)
> at org.geotools.referencing.CRS.findMathTransform(CRS.java:1097)
> ... 50 more
> Caused by: java.lang.IllegalArgumentException: No source axis match UP.
> at
> org.geotools.referencing.operation.matrix.GeneralMatrix.<init>(GeneralMatrix.java:314)
> at
> org.geotools.referencing.operation.matrix.GeneralMatrix.<init>(GeneralMatrix.java:221)
> at
> org.geotools.referencing.cs.AbstractCS.swapAndScaleAxis(AbstractCS.java:331)
> at
> org.geotools.referencing.operation.AbstractCoordinateOperationFactory.swapAndScaleAxis(AbstractCoordinateOperationFactory.java:276)
>
> It seems to me that 2D to 3D transform should be simple (just add a 0
> value for Z axis). Is there a way to do this? If not, shall I file an
> RFE? Thanks for your help.
>
> --
> Regards,
> Farrukh Najmi
>
> Web: http://www.wellfleetsoftware.com
>
>
>
> ------------------------------------------------------------------------------
> The ultimate all-in-one performance toolkit: Intel(R) Parallel Studio XE:
> Pinpoint memory and threading errors before they happen.
> Find and fix more than 250 security defects in the development cycle.
> Locate bottlenecks in serial and parallel code that limit performance.
> http://p.sf.net/sfu/intel-dev2devfeb
> _______________________________________________
> Geotools-gt2-users mailing list
> Geotools-gt2-users@lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/geotools-gt2-users
>
------------------------------------------------------------------------------
The ultimate all-in-one performance toolkit: Intel(R) Parallel Studio XE:
Pinpoint memory and threading errors before they happen.
Find and fix more than 250 security defects in the development cycle.
Locate bottlenecks in serial and parallel code that limit performance.
http://p.sf.net/sfu/intel-dev2devfeb
_______________________________________________
Geotools-gt2-users mailing list
Geotools-gt2-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/geotools-gt2-users
