Martin Desruisseaux created SIS-472:
---------------------------------------

             Summary: Unexpected lost of vertical coordinate
                 Key: SIS-472
                 URL: https://issues.apache.org/jira/browse/SIS-472
             Project: Spatial Information Systems
          Issue Type: Bug
          Components: Referencing
    Affects Versions: 1.0
            Reporter: Martin Desruisseaux


Conversion from projected (x,y,h) to geographic (φ,λ,H) coordinates causes the 
lost of vertical coordinates, which is set to NaN. Actually given that 
operation implies a transformation from geoidal to ellipsoidal height and 
SIS-362 is not yet fixed, we would rather expect an exception to be thrown. 
Code to reproduce:

{code:java}
import org.apache.sis.referencing.CRS;
import org.apache.sis.referencing.CommonCRS;
import org.opengis.referencing.crs.CoordinateReferenceSystem;
import org.opengis.referencing.operation.CoordinateOperation;

final class Test {
    public static void main(String[] args) throws Exception {
        CoordinateReferenceSystem sourceCRS = 
CRS.compound(CRS.forCode("EPSG:2154"), CommonCRS.Vertical.MEAN_SEA_LEVEL.crs());
        CoordinateReferenceSystem targetCRS = CommonCRS.WGS84.geographic3D();
        CoordinateOperation op = CRS.findOperation(sourceCRS, targetCRS, null);
        System.out.println(op);
        System.out.println();
        System.out.println(op.getMathTransform());
    }
}
{code}




--
This message was sent by Atlassian Jira
(v8.3.2#803003)

Reply via email to