GeotoolsTransformation.java

Go to the documentation of this file.

/* *********************************************************************** *
 * project: org.matsim.*
 * GeotoolsTransformation.java
 *                                                                         *
 * *********************************************************************** *
 *                                                                         *
 * copyright       : (C) 2007 by the members listed in the COPYING,        *
 *                   LICENSE and WARRANTY file.                            *
 * email           : info at matsim dot org                                *
 *                                                                         *
 * *********************************************************************** *
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *   See also COPYING, LICENSE and WARRANTY file                           *
 *                                                                         *
 * *********************************************************************** */

package org.matsim.core.utils.geometry.transformations;

import org.geotools.api.referencing.FactoryException;
import org.geotools.api.referencing.crs.CoordinateReferenceSystem;
import org.geotools.api.referencing.operation.MathTransform;
import org.geotools.api.referencing.operation.TransformException;
import org.geotools.geometry.jts.JTS;
import org.geotools.referencing.CRS;
import org.locationtech.jts.geom.Point;
import org.matsim.api.core.v01.Coord;
import org.matsim.core.utils.geometry.CoordUtils;
import org.matsim.core.utils.geometry.CoordinateTransformation;
import org.matsim.core.utils.geometry.geotools.MGC;

public class GeotoolsTransformation implements CoordinateTransformation {

    private MathTransform transform;

    public GeotoolsTransformation(final String from, final String to) {
        CoordinateReferenceSystem sourceCRS = MGC.getCRS(from);
        CoordinateReferenceSystem targetCRS = MGC.getCRS(to);

        try {
            this.transform = CRS.findMathTransform(sourceCRS, targetCRS,true);
        } catch (FactoryException e) {
            throw new RuntimeException(e);
        }
    }

    @Override
    public Coord transform(final Coord coord) {
        Point p = null;
        try {
            p = (Point) JTS.transform(MGC.coord2Point(coord), this.transform);
        } catch (TransformException e) {
            throw new RuntimeException(e);
        }
        if(coord.hasZ()){
            return CoordUtils.createCoord(p.getX(), p.getY(), coord.getZ());
        } else{
            return MGC.point2Coord(p);
        }
    }



}