[geotk] Transform from EPSG 4979 to 3D CRS

Aaron Braeckel braeckel at ucar.edu
Tue May 3 00:12:01 EDT 2011


Hello,

I am trying to transform between a 3D CRS (EPSG:4979) and a 3D Compound
CRS I am creating from scratch:

    Hints.putSystemDefault( Hints.LENIENT_DATUM_SHIFT, Boolean.TRUE );
    CoordinateReferenceSystem crs4979 = CRS.decode( "EPSG:4979" ); //3D
CRS - lat,lon,ellipsoidal height
    ProjectedCRS lambert = createLambertConformal( 25.0, 25.0, 25.0,
-95.0 );  //latitudes = 25.0, longitude = -95.0
    CRS.findMathTransform( CRS.getHorizontalCRS( crs4979 ), lambert ); 
//this works
    CoordinateReferenceSystem vertCRS = CRS.decode( "EPSG:5714"); 
//mean sea level height
    CoordinateReferenceSystem compoundCRS = new DefaultCompoundCRS(
"compound", lambert, vertCRS );
    CRS.findMathTransform( crs4979, compoundCRS );  //this fails with
"No transformation available from system "GeographicCRS[WGS 84]" to
"CompoundCRS[ CompoundCRS ]"."

Perhaps I am making a wrong assumption here.  Any thoughts on why I am
not able to find a transform between these two CRSs?

Thanks,
Aaron


More information about the Geotoolkit mailing list