diff options
Diffstat (limited to 'jniwrap/org/proj4/Projections.java')
| -rw-r--r-- | jniwrap/org/proj4/Projections.java | 268 |
1 files changed, 268 insertions, 0 deletions
diff --git a/jniwrap/org/proj4/Projections.java b/jniwrap/org/proj4/Projections.java new file mode 100644 index 00000000..6245eb68 --- /dev/null +++ b/jniwrap/org/proj4/Projections.java @@ -0,0 +1,268 @@ +/**===================================================================================== + + FILE: Projections.java + + DESCRIPTION: + + NOTES: --- + AUTHOR: Antonello Andrea + EMAIL: andrea.antonello@hydrologis.com + COMPANY: HydroloGIS / Engineering, University of Trento / CUDAM + COPYRIGHT: Copyright (C) 2004 HydroloGIS / University of Trento / CUDAM, ITALY, GPL + VERSION: $version$ + CREATED OR MODIFIED: Oct 15, 2004 + REVISION: --- + =====================================================================================*/ + +/* + This library is free software; you can redistribute it and/or + modify it under the terms of the GNU Library General Public + License as published by the Free Software Foundation; either + version 2 of the License, or (at your option) any later version. + + This library 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 + Library General Public License for more details. + + You should have received a copy of the GNU Library General Public + License along with this library; if not, write to the Free + Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 + USA + + 1. Redistributions of source code must retain the above copyright + notice, this list of conditions and the following disclaimer. + 2. Redistributions in binary form must reproduce the above copyright + notice, this list of conditions and the following disclaimer in the + documentation and/or other materials provided with the distribution. + */ +package org.proj4; + +import java.util.LinkedHashMap; + +/** + * superclass of all the projections the main proj methods are held here + */ +public abstract class Projections +{ + + double degtorad = .0174532925199432958; + double radtodeg = 57.29577951308232; + + String proj = null; + /** + * Comment for <code>projParameters</code> + * it is a hashmap that contains all the parameters that + * define the projection. Usually that can be:<BR> + * +init<BR> + * +proj<BR> + * +zone<BR> + * +ellps<BR> + * +datum<BR> + * +units<BR> + * +no_defs<BR> + * +ellps<BR> + * +towgs84<BR> + *<BR> + * ELLIPSOID INFO:<BR> + * name<BR> + * a<BR> + * e<BR> + * es<BR> + * ra<BR> + * one_es<BR> + * rone_es<BR> + * lam0<BR> + * phi0<BR> + * x0<BR> + * y0<BR> + * k0<BR> + * to_meter<BR> + * fr_meter<BR> + */ + private LinkedHashMap projParameters = null; + + static + { + System.loadLibrary("jproj"); + } + + /** + * get all the projection informations needed from the + * projection code (reproduces pj_get_def() of the proj api) + * + * @param the proj code or options + * @return the info String + */ + protected native String getProjInfo(String proj); + + /** + * get the ellipsoid parameters from the + * projection code + * + * @param the proj code or options + * @return the info String + */ + protected native String getEllipsInfo(String proj); + + /** + * native call to the reprojections routines of proj + * + * @param firstCoord array of x + * @param secondCoord array of y + * @param values array of z + * @param srcCodeString source projection code or option + * @param destCodeString destination projection code or option + * @param pointcount + * @param pointoffset + */ + protected native void transform(double[] firstCoord, double[] secondCoord, + double[] values, String srcCodeString, String destCodeString, + long pointcount, int pointoffset); + + /** + * public method to call the native getProjInfo + * + * @return quoting the proj api:"Returns the PROJ.4 command string that + * would produce this definition expanded as much as possible. For + * instance, +init= calls and +datum= defintions would be expanded" + */ + public String getProjInfo() + { + return getProjInfo(proj); + } + + /** + * public method to call the native getEllispdInfo + * + * @return the list of ellipsoid parameters + */ + public String getEllipseInfo() + { + return getEllipsInfo(proj); + } + + /** + * print to standard output the proj info in a nice format + */ + public void printProjInfo() + { + String projinfo = getProjInfo(); + + System.out.println(); + System.out + .println("******************************************************"); + System.out.println("* PROJECTION INFO:"); + System.out.println("*"); + String[] infos = projinfo.split("\\+"); + for (int i = 0; i < infos.length - 1; i++) + { + System.out.println("* +" + infos[i + 1].trim()); + } + System.out.println("*"); + System.out.println("* ELLIPSOID INFO:"); + System.out.println("*"); + String ellipsinfo = getEllipseInfo(); + String[] ellipsoid = ellipsinfo.split(";"); + for (int i = 0; i < ellipsoid.length; i++) + { + System.out.println("* " + ellipsoid[i].trim()); + } + + System.out + .println("******************************************************"); + System.out.println(); + } + + /** + * return all the proj info into a Linked Hashmap + */ + public LinkedHashMap mapProjInfo() + { + projParameters = new LinkedHashMap(); + String projinfo = getProjInfo(); + + String[] infos = projinfo.split("\\+"); + for (int i = 0; i < infos.length - 1; i++) + { + String[] pairs = infos[i + 1].split("="); + projParameters.put(pairs[0].trim(), pairs[1].trim()); + } + + String ellipsinfo = getEllipseInfo(); + String[] ellipsoid = ellipsinfo.split(";"); + for (int i = 0; i < ellipsoid.length; i++) + { + String[] pairs2 = ellipsoid[i].split(":"); + projParameters.put(pairs2[0].trim(), pairs2[1].trim()); + } + + return projParameters; + } + + /** + * if there is some operation to perform on the input data, this is + * the right moment (ex. latlong from degree to radiant) + * -> i.e. do whatever is needed + * + * @param dataTP + */ + public abstract void prepareData(ProjectionData dataTP); + + /** + * do the transform. The srcProjection is passed to the destination proj and + * the transformation takes place. Then the resulting transformed data + * are passed to the destProj, so that the destProj in case can take care of + * final transformation of data (ex. if the destProj is latlong, the values + * have to be set beck to degrees + * + * @param srcProj object holding the source projection + * @param dataTP the data set + * @param point_count + * @param point_offset + */ + protected void doTheTransform(Projections srcProj, ProjectionData dataTP, + long point_count, int point_offset) + { + transform(dataTP.x, dataTP.y, dataTP.z, srcProj.proj, proj, point_count, + point_offset); + } + + /** + * this takes care that the reprojected data are in the correct format + * (ex. latlong has to be transformed back to radiant) + * + * @param dataTP the data set + */ + public abstract void prepareTransformedData(ProjectionData dataTP); + + /** + * transform latitude and longitude from degree to radiant format + * + * @param la + * @param lo + */ + protected void degreeToRadiant(double[] la, double[] lo) + { + for (int i = 0; i < la.length; i++) + { + la[i] = la[i] * degtorad; + lo[i] = lo[i] * degtorad; + } + } + + /** + * transform latitude and longitude from radiant to degree format + * + * @param la + * @param lo + */ + protected void radiantToDegree(double[] la, double[] lo) + { + for (int i = 0; i < la.length; i++) + { + la[i] = la[i] * radtodeg; + lo[i] = lo[i] * radtodeg; + } + } +} |