org.cts.crs.Geographic3DCRS Maven / Gradle / Ivy
/*
* Coordinate Transformations Suite (abridged CTS) is a library developped to
* perform Coordinate Transformations using well known geodetic algorithms
* and parameter sets.
* Its main focus are simplicity, flexibility, interoperability, in this order.
*
* This library has been originally developed by Michaël Michaud under the JGeod
* name. It has been renamed CTS in 2009 and shared to the community from
* the Atelier SIG code repository.
*
* Since them, CTS is supported by the Atelier SIG team in collaboration with Michaël
* Michaud.
* The new CTS has been funded by the French Agence Nationale de la Recherche
* (ANR) under contract ANR-08-VILL-0005-01 and the regional council
* "Région Pays de La Loire" under the projet SOGVILLE (Système d'Orbservation
* Géographique de la Ville).
*
* CTS 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 3 of the License, or (at your option) any later
* version.
*
* CTS is distributed in the hope that it will be useful, but WITHOUT ANY
* WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
* A PARTICULAR PURPOSE. See the GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License along with
* CTS. If not, see A Geographic CoordinateReferenceSystem is a reference system based on a
* GeodeticDatum and a 2D or 3D Ellipsoidal Coordinate System.
*
* @author Michaël Michaud
*/
public class Geographic3DCRS extends GeodeticCRS {
/**
* A 3D {@link CoordinateSystem} whose first {@link Axis} contains latitude,
* second {@link Axis} contains longitude and third axis contains
* ellipsoidal height. The units used by these axes are radian and meter.
*/
public static CoordinateSystem LATLONH_RRM_CS = new CoordinateSystem(
new Axis[]{LATITUDE, LONGITUDE, HEIGHT}, new Unit[]{RADIAN,
RADIAN, METER});
/**
* A 3D {@link CoordinateSystem} whose first {@link Axis} contains
* longitude, second {@link Axis} contains latitude and third axis contains
* ellipsoidal height. The units used by these axes are radian and meter.
*/
public static CoordinateSystem LONLATH_RRM_CS = new CoordinateSystem(
new Axis[]{LONGITUDE, LATITUDE, HEIGHT}, new Unit[]{RADIAN,
RADIAN, METER});
/**
* A 3D {@link CoordinateSystem} whose first {@link Axis} contains latitude,
* second {@link Axis} contains longitude and third axis contains
* ellipsoidal height. The units used by these axes are decimal degree and
* meter.
*/
public static CoordinateSystem LATLONH_DDM_CS = new CoordinateSystem(
new Axis[]{LATITUDE, LONGITUDE, HEIGHT}, new Unit[]{DEGREE,
DEGREE, METER});
/**
* A 3D {@link CoordinateSystem} whose first {@link Axis} contains
* longitude, second {@link Axis} contains latitude and third axis contains
* ellipsoidal height. The units used by these axes are decimal degree and
* meter.
*/
public static CoordinateSystem LONLATH_DDM_CS = new CoordinateSystem(
new Axis[]{LONGITUDE, LATITUDE, HEIGHT}, new Unit[]{DEGREE,
DEGREE, METER});
/**
* A 3D {@link CoordinateSystem} whose first {@link Axis} contains latitude,
* second {@link Axis} contains longitude and third axis contains
* ellipsoidal height. The units used by these axes are grad and meter.
*/
public static CoordinateSystem LATLONH_GGM_CS = new CoordinateSystem(
new Axis[]{LATITUDE, LONGITUDE, HEIGHT}, new Unit[]{GRAD,
GRAD, METER});
/**
* A 3D {@link CoordinateSystem} whose first {@link Axis} contains
* longitude, second {@link Axis} contains latitude and third axis contains
* ellipsoidal height. The units used by these axes are grad and meter.
*/
public static CoordinateSystem LONLATH_GGM_CS = new CoordinateSystem(
new Axis[]{LONGITUDE, LATITUDE, HEIGHT}, new Unit[]{GRAD,
GRAD, METER});
/**
* Create a new Geographic3DCRS.
*
* @param identifier the identifier of the Geographic3DCRS
* @param datum the datum associated with the Geographic3DCRS
* @param coordSys the coordinate system associated with the Geographic3DCRS
*/
public Geographic3DCRS(Identifier identifier, GeodeticDatum datum,
CoordinateSystem coordSys) {
super(identifier, datum, coordSys);
}
/**
* Create a new Geographic2DCRS. The first {@link Axis} of the associated
* {@link CoordinateSystem} contains latitude, the second {@link Axis}
* contains longitude and the third contains the ellipsoidal height in
* meters.
*
* @param identifier the identifier of the Geographic3DCRS
* @param datum the datum associated with the Geographic3DCRS
* @param unit the angular unit to use for the two first axis of the
* coordinate system associated with the Geographic3DCRS
*/
public Geographic3DCRS(Identifier identifier, GeodeticDatum datum, Unit unit) {
super(identifier, datum, LATLONH_DDM_CS);
if (unit == RADIAN) {
this.coordinateSystem = LATLONH_RRM_CS;
} else if (unit == DEGREE) {
this.coordinateSystem = LATLONH_DDM_CS;
} else if (unit == GRAD) {
this.coordinateSystem = LATLONH_GGM_CS;
} else;
}
/**
* Create a new Geographic3DCRS. The first {@link Axis} of the associated
* {@link CoordinateSystem} contains latitude, the second {@link Axis}
* contains longitude and the third contains the ellipsoidal height in
* meters.
*
* @param identifier the identifier of the Geographic2DCRS
* @param datum the datum associated with the Geographic2DCRS
*/
public Geographic3DCRS(Identifier identifier, GeodeticDatum datum) {
super(identifier, datum, LATLONH_DDM_CS);
}
/**
* @see GeodeticCRS#getType()
*/
@Override
public Type getType() {
return Type.GEOGRAPHIC3D;
}
/**
* @see GeodeticCRS#getProjection()
*/
@Override
public Projection getProjection() {
return null;
}
/**
* @see GeodeticCRS#toGeographicCoordinateConverter()
*/
@Override
public CoordinateOperation toGeographicCoordinateConverter() {
List ops = new ArrayList();
for (int i = 0; i < 3; i++) {
if (getCoordinateSystem().getAxis(i).getDirection() == SOUTH
|| getCoordinateSystem().getAxis(i).getDirection() == WEST
|| getCoordinateSystem().getAxis(i).getDirection() == DOWN) {
ops.add(new OppositeCoordinate(i));
}
}
// Convert from source unit to radians
ops.add(UnitConversion.createUnitConverter(getCoordinateSystem().getUnit(0), RADIAN,
getCoordinateSystem().getUnit(2), METER));
// switch from LON/LAT to LAT/LON coordinate if necessary
if (getCoordinateSystem().getAxis(0).getDirection() == EAST
|| getCoordinateSystem().getAxis(0).getDirection() == WEST) {
ops.add(CoordinateSwitch.SWITCH_LAT_LON);
}
return new CoordinateOperationSequence(new Identifier(
CoordinateOperationSequence.class), ops);
}
/**
* @see GeodeticCRS#fromGeographicCoordinateConverter()
*/
@Override
public CoordinateOperation fromGeographicCoordinateConverter() {
List ops = new ArrayList();
// switch from LON/LAT to LAT/LON coordinate if necessary
if (getCoordinateSystem().getAxis(0).getDirection() == EAST
|| getCoordinateSystem().getAxis(0).getDirection() == WEST) {
ops.add(CoordinateSwitch.SWITCH_LAT_LON);
}
// Convert from radian to this coordinate system's units
ops.add(UnitConversion.createUnitConverter(RADIAN, getCoordinateSystem().getUnit(0),
METER, getCoordinateSystem().getUnit(2)));
for (int i = 0; i < 3; i++) {
if (getCoordinateSystem().getAxis(i).getDirection() == SOUTH
|| getCoordinateSystem().getAxis(i).getDirection() == WEST
|| getCoordinateSystem().getAxis(i).getDirection() == DOWN) {
ops.add(new OppositeCoordinate(i));
}
}
return new CoordinateOperationSequence(new Identifier(
CoordinateOperationSequence.class), ops);
}
/**
* Returns a WKT representation of the geographic 3D CRS.
*
*/
public String toWKT() {
StringBuilder w = new StringBuilder();
w.append("GEOGCS[\"");
w.append(this.getName());
w.append("\",");
w.append(this.getDatum().toWKT());
w.append(',');
w.append(this.getDatum().getPrimeMeridian().toWKT());
w.append(',');
w.append(this.getCoordinateSystem().getUnit(0).toWKT());
for (int i = 0; i < this.getCoordinateSystem().getDimension(); i++) {
w.append(',');
w.append(this.getCoordinateSystem().getAxis(i).toWKT());
}
if (!this.getAuthorityName().startsWith(Identifiable.LOCAL)) {
w.append(',');
w.append(this.getIdentifier().toWKT());
}
w.append(']');
return w.toString();
}
}