[SIS-472] Unexpected lost of vertical coordinate

• Type: Bug
• Status: Open
• Priority: Major
• Resolution: Unresolved
• Affects Version/s: 1.0
• Component/s: Referencing
• Assignee: Unassigned
• Reporter: Martin Desruisseaux
• Created: 11/Sep/19 10:13
• Updated: 11/Sep/19 10:13

Description

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:

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());
    }
}

Issue Links

• depends upon: SIS-362 Is there a way to convert geoidal altitudes to ellipsoidal altitudes?