Application: gvSIG desktop » gvSIG projections

/**

 * Clase para manejar la libreria proj4

 * @author	David Hern?ndez L?pez (david.hernandez@uclm.es)

 * @author  Jos? Luis G?mez Mart?nez (JoseLuis.Gomez@uclm.es)

 *

 */

        public static final Logger logger = LoggerFactory.getLogger(Proj4.class);

	private ArrayList projectionNameList= new ArrayList();

	private ArrayList projectionParameterUnitList= new ArrayList();

	private ArrayList projectionAcronymList= new ArrayList();

	private ArrayList projectionParameterAcronymList= new ArrayList();

	int divider=10000;

	private static double angularTolerance=1.0/3600.0;

	/**

     * Small tolerance factor for rounding errors.

     */

    private static final double EPS = 1E-8;

    /**

     * EPSG codes of the Swiss Oblique Mercator CRSs

     */

    private static final Integer[] somerc_codes = { 2056, 21780, 21781, 21782, 23700 };

	public Proj4() throws CrsException

	{

		//defineProjectionNames();

		defineProjectionParameterList();

		defineProjections();

	}

	private void defineUnitNameList() throws CrsException

	{

		int count=0;

		{

			String[] unitName={"Angular"};

			unitNameList.add(count,unitName);

		}

		count++;

		{

			String[] unitName={"Linear"};

			unitNameList.add(count,unitName);

		}

		count++;

		{

			String[] unitName={"Unitless"};

			unitNameList.add(count,unitName);

			addUnitName(count,"Adimensional");

		}

	}

	private void defineProjectionParameterList() throws CrsException

	{

		int count=0;

		{ // azimuth

			String[] parameterName={"azimuth"};

			projectionParameterNameList.add(count,parameterName);

			String[] parameterUnit={"Angular"};

			projectionParameterUnitList.add(count,parameterUnit);

		{ // central_meridian

			String[] parameterName={"central_meridian"};

			projectionParameterNameList.add(count,parameterName);

			addProjectionParameterName(count,"Longitude of projection center");

			addProjectionParameterName(count,"Longitude_of_center");

			addProjectionParameterName(count,"ProjCenterLong");

			String[] parameterUnit={"Angular"};

			projectionParameterUnitList.add(count,parameterUnit);

		}

		count++;

		{ // false_easting

			String[] parameterName={"false_easting"};

			projectionParameterNameList.add(count,parameterName);

			String[] parameterUnit={"Linear"};

			projectionParameterUnitList.add(count,parameterUnit);

		}

		count++;

		{ // false_northing

			String[] parameterName={"false_northing"};

			projectionParameterNameList.add(count,parameterName);

			String[] parameterUnit={"Linear"};

			projectionParameterUnitList.add(count,parameterUnit);

		{ // latitude_of_center

			String[] parameterName={"latitude_of_center"};

			projectionParameterNameList.add(count,parameterName);

			addProjectionParameterName(count,"FalseOriginLat");

			addProjectionParameterName(count,"Latitude of natural origin");

			addProjectionParameterName(count,"Latitude of projection center");

			addProjectionParameterName(count,"Latitude of projection centre");

			addProjectionParameterName(count,"ProjCenterLat");

			addProjectionParameterName(count,"Spherical_latitude_of_origin");

			addProjectionParameterName(count,"Central_Parallel");

		{ // Latitude_Of_1st_Point

			String[] parameterName={"Latitude_Of_1st_Point"};

			projectionParameterNameList.add(count,parameterName);

			String[] parameterUnit={"Angular"};

			projectionParameterUnitList.add(count,parameterUnit);

		}

		{ // Latitude_Of_2nd_Point

			String[] parameterName={"Latitude_Of_2nd_Point"};

			projectionParameterNameList.add(count,parameterName);

			String[] parameterUnit={"Angular"};

			projectionParameterUnitList.add(count,parameterUnit);

		}

		{ // latitude_of_origin

			String[] parameterName={"latitude_of_origin"};

			projectionParameterNameList.add(count,parameterName);

			addProjectionParameterName(count,"FalseOriginLat");

			addProjectionParameterName(count,"Latitude of natural origin");

			addProjectionParameterName(count,"Latitude of projection center");

			addProjectionParameterName(count,"Latitude of projection centre");

			addProjectionParameterName(count,"ProjCenterLat");

			String[] parameterUnit={"Angular"};

			projectionParameterUnitList.add(count,parameterUnit);

		}

		{ // latitude_of_origin

			String[] parameterName={"latitude_of_standard_parallel"};

			projectionParameterNameList.add(count,parameterName);

			addProjectionParameterName(count,"FalseOriginLat");

			addProjectionParameterName(count,"Latitude of natural origin");

			addProjectionParameterName(count,"Latitude of projection center");

			addProjectionParameterName(count,"Latitude of projection centre");

			addProjectionParameterName(count,"ProjCenterLat");

			String[] parameterUnit={"Angular"};

			projectionParameterUnitList.add(count,parameterUnit);

		}

		{ // longitude_of_center

			String[] parameterName={"longitude_of_center"};

			projectionParameterNameList.add(count,parameterName);

			addProjectionParameterName(count,"Spherical_latitude_of_origin");

			String[] parameterUnit={"Angular"};

		{ // Longitude_Of_1st_Point

			String[] parameterName={"Longitude_Of_1st_Point"};

			projectionParameterNameList.add(count,parameterName);

			String[] parameterUnit={"Angular"};

			projectionParameterUnitList.add(count,parameterUnit);

		}

		{ // Longitude_Of_2nd_Point

			String[] parameterName={"Longitude_Of_2nd_Point"};

			projectionParameterNameList.add(count,parameterName);

			String[] parameterUnit={"Angular"};

			projectionParameterUnitList.add(count,parameterUnit);

		}

		{ // pseudo_standard_parallel_1

			String[] parameterName={"pseudo_standard_parallel_1"};

			projectionParameterNameList.add(count,parameterName);

			String[] parameterUnit={"Angular"};

			projectionParameterUnitList.add(count,parameterUnit);

		}

		count++;

		{ // satellite_height

			String[] parameterName={"rectified_grid_angle"};

			projectionParameterNameList.add(count,parameterName);

			String[] parameterUnit={"Linear"};

			projectionParameterUnitList.add(count,parameterUnit);

		}

		count++;

		{ // satellite_height

			String[] parameterName={"satellite_height"};

			projectionParameterNameList.add(count,parameterName);

			String[] parameterUnit={"Linear"};

			projectionParameterUnitList.add(count,parameterUnit);

		}

		count++;

		{ // scale_factor

			String[] parameterName={"scale_factor"};

			projectionParameterNameList.add(count,parameterName);

			String[] parameterUnit={"Unitless"};

			projectionParameterUnitList.add(count,parameterUnit);

		{ // standard_parallel_1

			String[] parameterName={"standard_parallel_1"};

			projectionParameterNameList.add(count,parameterName);

			String[] parameterUnit={"Angular"};

			projectionParameterUnitList.add(count,parameterUnit);

		{ // standard_parallel_2

			String[] parameterName={"standard_parallel_2"};

			projectionParameterNameList.add(count,parameterName);

			String[] parameterUnit={"Angular"};

			projectionParameterUnitList.add(count,parameterUnit);

		{ // semi_major

			String[] parameterName={"semi_major"};

			projectionParameterNameList.add(count,parameterName);

			String[] parameterUnit={"Linear"};

			projectionParameterUnitList.add(count,parameterUnit);

		{ // semi_minor

			String[] parameterName={"semi_minor"};

			projectionParameterNameList.add(count,parameterName);

			String[] parameterUnit={"Linear"};

			projectionParameterUnitList.add(count,parameterUnit);

		{ // height

			String[] parameterName={"Height"};

			projectionParameterNameList.add(count,parameterName);

			String[] parameterUnit={"Linear"};

			projectionParameterUnitList.add(count,parameterUnit);

		}

	}

	private void defineProjections() throws CrsException

	{

		int count=0;

		{// Aitoff

			String[] projectionName={"Aitoff"};

			projectionNameList.add(count,projectionName);

			String[] projectionAcronym={"aitoff"};

			projectionAcronymList.add(count,projectionAcronym);

		}

		count++;

		{// Albers Equal-Area Conic

			String[] projectionName={"Albers_Conic_Equal_Area"};

			addProjectionName(count,"Albers Equal-Area Conic");

			addProjectionName(count,"Albers Equal Area");

			addProjectionName(count,"9822");

			String[] parameterName={"standard_parallel_1"};

			projectionParameterList.add(count,parameterName);

			addProjectionParameter(count,"standard_parallel_2");

			String[] projectionAcronym={"aea"};

			projectionAcronymList.add(count,projectionAcronym);

		}

		count++;

		{// Azimuthal_Equidistantt

			String[] projectionName={"Azimuthal_Equidistant"};

			addProjectionName(count,"Zenithal Equidistant");

			addProjectionName(count,"Zenithal-Equidistant");

			addProjectionName(count,"Zenithal_Equidistant");

			String[] parameterName={"latitude_of_center"};

			projectionParameterList.add(count,parameterName);

			addProjectionParameter(count,"longitude_of_center");

			String[] projectionAcronym={"cass"};

			projectionAcronymList.add(count,projectionAcronym);

		}

		count++;

		{// Craster Parabolic

			String[] projectionName={"Craster_Parabolic"};

			String[] projectionAcronym={"eck6"};

			projectionAcronymList.add(count,projectionAcronym);

		}

		count++;

		{// Equidistant_Conic

			String[] projectionName={"Equidistant_Conic"};

			projectionNameList.add(count,projectionName);

			addProjectionName(count,"Equidistant Conic");

			String[] parameterName={"latitude_of_center"};

			projectionParameterList.add(count,parameterName);

			addProjectionParameter(count,"longitude_of_center");

			String[] projectionAcronym={"eqc"};

			projectionAcronymList.add(count,projectionAcronym);

		}

		count++;

		{// McBryde-Thomas Flat-Polar Quartic

			String[] projectionName={"McBryde_Thomas_Flat_Polar_Quartic"};

			String[] projectionAcronym={"goode"};

			projectionAcronymList.add(count,projectionAcronym);

		}

		count++;

		{// Hammer-Aitoff

			String[] projectionName={"Hammer_Aitoff"};

			String[] projectionAcronym={"hammer"};

			projectionAcronymList.add(count,projectionAcronym);

		}

		count++;

		{// Krovak

			String[] projectionName={"Krovak"};

			String[] projectionAcronym={"krovak"};

			projectionAcronymList.add(count,projectionAcronym);

		}

		count++;

		{// Laborde Madagascar

			String[] projectionName={"Laborde_Madagascar"};

			String[] projectionAcronym={"labrd"};

			projectionAcronymList.add(count,projectionAcronym);

			/*String[] parameterName={"azimuth"};

			projectionParameterList.add(count,parameterName);

			String[] projectionAcronym={"lcca"};

			projectionAcronymList.add(count,projectionAcronym);

		}

		count++;

		{// Lambert_Conformal_Conic_1SP

			String[] projectionName={"Lambert_Conformal_Conic_1SP"};

			String[] projectionAcronym={"lcc"};

			projectionAcronymList.add(count,projectionAcronym);

		}

		count++;

		{// Lambert_Conformal_Conic_2SP

			String[] projectionName={"Lambert_Conformal_Conic_2SP"};

			String[] projectionAcronym={"lcc"};

			projectionAcronymList.add(count,projectionAcronym);

		}

		count++;

		{// Lambert_Conformal_Conic

			String[] projectionName={"Lambert_Conformal_Conic"};

			String[] projectionAcronym={"lcc"};

			projectionAcronymList.add(count,projectionAcronym);

		}

		count++;

		{// Lambert_Conformal_Conic_2SP_Belgium

			String[] projectionName={"Lambert_Conformal_Conic_2SP_Belgium"};

			String[] projectionAcronym={"loxim"};

			projectionAcronymList.add(count,projectionAcronym);

		}

		count++;

		{// Mercator_1SP

			String[] projectionName={"Mercator_1SP"};

			String[] projectionAcronym={"merc"};

			projectionAcronymList.add(count,projectionAcronym);

		}

		count++;

		{// Mercator_2SP

			String[] projectionName={"Mercator_2SP"};

			String[] projectionAcronym={"mill"};

			projectionAcronymList.add(count,projectionAcronym);

		}

		count++;

		{// Mollweide

			String[] projectionName={"Mollweide"};

			String[] projectionAcronym={"moll"};

			projectionAcronymList.add(count,projectionAcronym);

		}

		count++;

		{// Near-Sided Perspective

			String[] projectionName={"Near_Sided_Perspective"};

			addProjectionName(count,"Vertical Near Side Perspective");

			addProjectionName(count,"Vertical-Near-Side-Perspective");

			addProjectionName(count,"Vertical_Near_Side_Perspective");

			String[] parameterName={"latitude_of_center"};

			projectionParameterList.add(count,parameterName);

			addProjectionParameter(count,"longitude_of_center");

			String[] projectionAcronym={"wintri"};

			projectionAcronymList.add(count,projectionAcronym);

		}

		count++;

		{// Popular_Visualisation_Pseudo_Mercator (based on Mercator_1SP)

			String[] projectionName={"Popular_Visualisation_Pseudo_Mercator"};

			projectionAcronymList.add(count,projectionAcronym);

		}

        }

	public void addProjectionName(int pos,String projectionName) throws CrsException

	{

		if(pos<0||pos>(projectionNameList.size()-1))

		projectionNameList.remove(pos);

		projectionNameList.add(pos,newProjectionNames);

	}

	public void addUnitName(int pos,String unitName) throws CrsException

	{

		if(pos<0||pos>(unitNameList.size()-1))

		unitNameList.remove(pos);

		unitNameList.add(pos,newUnitNames);

	}

	public void addProjectionParameterName(int pos,String projectionParameterName) throws CrsException

	{

		if(pos<0||pos>(projectionParameterNameList.size()-1))

		projectionParameterNameList.remove(pos);

		projectionParameterNameList.add(pos,newProjectionParameterNames);

	}

	public void addProjectionParameter(int pos,String projectionParameter) throws CrsException

	{

		if(pos<0||pos>(projectionParameterList.size()-1))

		projectionParameterList.remove(pos);

		projectionParameterList.add(pos,newProjectionParameters);

	}

	public void addProjectionParameterDefaultValue(int pos,String projectionParameterDefaultValue) throws CrsException

	{

		if(pos<0||pos>(projectionParameterDefaultValueList.size()-1))

		projectionParameterDefaultValueList.remove(pos);

		projectionParameterDefaultValueList.add(pos,newProjectionParameterDefaultValues);

	}

	public void addProjectionParameterMaxValue(int pos,String projectionParameterMaxValue) throws CrsException

	{

		if(pos<0||pos>(projectionParameterMaxValueList.size()-1))

		projectionParameterMaxValueList.remove(pos);

		projectionParameterMaxValueList.add(pos,newProjectionParameterMaxValues);

	}

	public void addProjectionParameterMinValue(int pos,String projectionParameterMinValue) throws CrsException

	{

		if(pos<0||pos>(projectionParameterMinValueList.size()-1))

		projectionParameterMinValueList.remove(pos);

		projectionParameterMinValueList.add(pos,newProjectionParameterMinValues);

	}

	public void addProjectionParameterAcronymList(int pos,String projectionParameterAcronym) throws CrsException

	{

		if(pos<0||pos>(projectionParameterAcronymList.size()-1))

		projectionParameterAcronymList.remove(pos);

		projectionParameterAcronymList.add(pos,newProjectionParameterAcronyms);

	}

	public int findProjection(String projectionName)

	{

		for(int i=0;i<projectionNameList.size();i++)

					return i;

			}

		}

		return -1;

	}

	public int findProjectionParameter(String parameterName)

	{

		for(int i=0;i<projectionParameterNameList.size();i++)

					return i;

			}

		}

		return -1;

	}

	public int findProjectionParameters(String parameterName1, String parameterName2)

	{

		for(int i=0;i<projectionParameterNameList.size();i++)

		{

			boolean existsParameter1=false;

			boolean existsParameter2=false;

			String[] parameterNames=(String[]) projectionParameterNameList.get(i);

			if(parameterNames[0].toLowerCase().replaceAll(" ","").equals(parameterName2.toLowerCase().replaceAll(" ","")))

			{

			if(existsParameter1&&existsParameter2)

				return i;

		}

		return -1;

	}

	public int findUnit(String unitName)

	{

		for(int i=0;i<unitNameList.size();i++)

					return i;

			}

		}

		return -1;

	}

	{

		if(pos<0||pos>(unitNameList.size()-1))

		{

		}

		return ((String[]) unitNameList.get(pos))[0];

	}

	{

		if(pos<0||pos>(projectionNameList.size()-1))

		{

		}

		return ((String[]) projectionNameList.get(pos))[0];

	}

	{

		// Ojo decir a Jose Luis que quite los parametros del elipsoide

		if(pos<0||pos>(projectionParameterNameList.size()-1))

		}

		return ((String[]) projectionParameterNameList.get(pos))[0];

	}

	{

		if(pos<0||pos>(projectionParameterList.size()-1))

				throw new CrsException(new Exception());

		String[] parameterList=(String[])projectionParameterList.get(pos);

		}

		return parameters;

	}

	{

		if(pos<0||pos>(projectionParameterDefaultValueList.size()-1))

		{

			throw new CrsException(new Exception("Error retrieving default value to projection parameter. Position out of valid limits."));

		}

		return parameterDefaultValues;

	}

	{

		if(pos<0||pos>(projectionParameterMaxValueList.size()-1))

		{

			throw new CrsException(new Exception("Error retrieving max value value to projection parameter. Position out of valid limits."));

		}

		return parameterMaxValues;

	}

	{

		if(pos<0||pos>(projectionParameterMinValueList.size()-1))

		{

			throw new CrsException(new Exception("Error retrieving min value value to projection parameter. Position out of valid limits."));

		}

		return parameterMinValues;

	}

	{

		if(pos<0||pos>(projectionParameterAcronymList.size()-1))

		{

			throw new CrsException(new Exception("Error retrieving projection acronym. Position out of valid limits."));

		}

		return parameterAcronyms;

	}

	public String getProjectionParameterUnitList(int pos) throws CrsException {

		if(pos<0||pos>(projectionParameterUnitList.size()-1))

		{

		strDatumName=crs.getCrsWkt().getDatumName();

        int intCodDatum=0;//Integer.parseInt(codDatum);

        strProj4Datum=datumToProj4(strDatumName,intCodDatum);

		String strProj=crs.getCrsWkt().getProjcs();

		if(strProj.equals(""))

		{

	        else

	        	strProj4ToMeter="+units=m ";

	        */

			//System.out.println("Projection      = "+strProj);

			String strProjName=crs.getCrsWkt().getProjection();

			//System.out.println("Projection Name = "+strProjName);

						throw new CrsException(new Exception("In proj4 the 'Oblique_Mercator' projection doesn't allow different azimut and spin axis"));

					}

				}

				// Necesito buscar un caso de ejemplo

			}

			if(projectionName.equals("Hotine_Oblique_Mercator_Two_Point_Center"))

			{

				isLaborde=true;

			}

			strProj4=strProj4+strProjectionAcronym+" ";

			for(int i=0;i<parameterNames.size();i++)

							}

						}

					}

					if(gtParameterAcronymn.equalsIgnoreCase("lat_ts"))

					{

						if(exists_sf)

					if(parameterAcronym.equals("lonc"))

						parameterAcronym="lon_0";

					if(parameterAcronym.equals("k"))

				}

				if(isOmerc)

				{

			strProj4+=strExtraProj4;

			//getProj4ProjectionName();

		}

		// Que pasa si no hay elipsoide? -> OGR por defecto pone WGS84

		double a=0;

		double inv_f=0;

			inv_f = Double.POSITIVE_INFINITY;

                }

		String strEllipseAcronym=ellipseToProj4(a,inv_f);

		if(strEllipseAcronym.equals(""))

		{

			if (!Double.isInfinite(inv_f))

		String strWkt=crs.getWKT();

		if(!strProj4ToMeter.equals(""))

			strProj4+=strProj4ToMeter;

		strProj4 = strProj4.replace(".0 ", " ");

		//System.out.println("- Cadena proj4: "+strProj4);

	}

	public String exportToProj4(CoordinateReferenceSystem crs) throws CrsException

	{

		String strProj4="+proj=";

				strDatumName=val[1];

	        for (Iterator iter =((DefaultProjectedCRS)crs).getDatum().getIdentifiers().iterator();iter.hasNext();) {

	            Identifier element = (Identifier) iter.next();

	        }

	        int intCodDatum=Integer.parseInt(codDatum);

	        strProj4Datum=datumToProj4(strDatumName,intCodDatum);

				strDatumName=val[1];

	        for (Iterator iter =((DefaultGeographicCRS)crs).getDatum().getIdentifiers().iterator();iter.hasNext();) {

	            Identifier element = (Identifier) iter.next();

	        }

	        int intCodDatum=Integer.parseInt(codDatum);

	        strProj4Datum=datumToProj4(strDatumName,intCodDatum);

		else {

			throw(new CrsException(new Exception("the crs is not a geographic nor projected CRS")));

		}

		String primeMeridianName=primeMeridian[0];

		double primeMeridianValue = -1;

		if (primeMeridian[1] == null) {

			throw(new CrsException(new Exception("Null value not allowed for prime meridiam parameters")));

		} else {

			primeMeridianValue=Double.parseDouble(primeMeridian[1]); // -> hacer el cambio de unidades

		String[] strPrimeMeridianProj4=primeMeridianToProj4(primeMeridianName,primeMeridianValue);

		primeMeridianValue=Double.parseDouble(strPrimeMeridianProj4[1]);

		primeMeridianName=strPrimeMeridianProj4[0];

				if(existsParameter)

				{

					parameterValues.set(i,parameterValue);

			}

			boolean isSomerc=false;

			boolean isOmerc=false;

						throw new CrsException(new Exception("The proj4 projection 'Oblique Mercator' not allow diferent azimut and spin axis."));

					}

				}

				// Necesito buscar un caso de ejemplo

			}

			if(projectionName.equals("Hotine_Oblique_Mercator_Two_Point_Center"))

			strProj4+=strExtraProj4;

			//getProj4ProjectionName();

		}

		// Que pasa si no hay elipsoide? -> OGR por defecto pone WGS84

		String strEllipseAcronym=ellipseToProj4(a,inv_f);

		String strEllipse="";

		if(!strProj4ToMeter.equals(""))

			strProj4+=strProj4ToMeter;

		String strWkt=crs.toWKT();

		strProj4 = strProj4.replace(".0 ", " ");

		//System.out.println("- Cadena proj4: "+strProj4);

	}

	private String[] primeMeridianToProj4(String pmName,double pmValue) throws CrsException

            pmValue=dfFromGreenwich;

            acronym="+pm="+pszPM+" ";

        }

        if(pmName.equalsIgnoreCase("paris")||(Math.abs(dfFromGreenwich-pmValue)<tolerance))

        {

            pszPM="paris";

            pmValue=dfFromGreenwich;

            acronym="+pm="+pszPM+" ";

        }

        dfFromGreenwich=-(74+4/60.0+51.3/3600.0);

        if(pmName.equalsIgnoreCase("bogota")||(Math.abs(dfFromGreenwich-pmValue)<tolerance))

        {

            pmValue=dfFromGreenwich;

            acronym="+pm="+pszPM+" ";

        }

        dfFromGreenwich=-(3+41/60.0+16.58/3600.0);  // mal en ogr los segundos pone 16.48

        if(pmName.equalsIgnoreCase("madrid")||(Math.abs(dfFromGreenwich-pmValue)<tolerance))

        {

            pmValue=dfFromGreenwich;

            acronym="+pm="+pszPM+" ";

        }

        dfFromGreenwich=(12+27/60.0+8.4/3600.0);

        if(pmName.equalsIgnoreCase("rome")||(Math.abs(dfFromGreenwich-pmValue)<tolerance))

        {

            pmValue=dfFromGreenwich;

            acronym="+pm="+pszPM+" ";

        }

        dfFromGreenwich=(7+26/60.0+22.5/3600.0);

        if(pmName.equalsIgnoreCase("bern")||(Math.abs(dfFromGreenwich-pmValue)<tolerance))

        {

            pmValue=dfFromGreenwich;

            acronym="+pm="+pszPM+" ";

        }

        dfFromGreenwich=(106+48/60.0+27.79/3600.0);

        if(pmName.equalsIgnoreCase("jakarta")||(Math.abs(dfFromGreenwich-pmValue)<tolerance))

        {

            pmValue=dfFromGreenwich;

            acronym="+pm="+pszPM+" ";

        }

        dfFromGreenwich=-(17+40/60.0+0.0/3600.0);

        if(pmName.equalsIgnoreCase("ferro")||(Math.abs(dfFromGreenwich-pmValue)<tolerance))

        {

            pmValue=dfFromGreenwich;

            acronym="+pm="+pszPM+" ";

        }

        dfFromGreenwich=(4+22/60.0+4.71/3600.0);

        if(pmName.equalsIgnoreCase("brussels")||(Math.abs(dfFromGreenwich-pmValue)<tolerance))

        {

            pmValue=dfFromGreenwich;

            acronym="+pm="+pszPM+" ";

        }

        dfFromGreenwich=(18+3/60.0+29.8/3600.0);

        if(pmName.equalsIgnoreCase("stockholm")||(Math.abs(dfFromGreenwich-pmValue)<tolerance))

        {

            pmValue=dfFromGreenwich;

            acronym="+pm="+pszPM+" ";

        }

        dfFromGreenwich=(23+42/60.0+58.815/3600.0);

        if(pmName.equalsIgnoreCase("athens")||(Math.abs(dfFromGreenwich-pmValue)<tolerance))

        {

            pmValue=dfFromGreenwich;

            acronym="+pm="+pszPM+" ";

        }

        dfFromGreenwich=(10+43/60.0+22.5/3600.0);

        if(pmName.equalsIgnoreCase("oslo")||(Math.abs(dfFromGreenwich-pmValue)<tolerance))

        {

            pmValue=dfFromGreenwich;

            acronym="+pm="+pszPM+" ";

        }

        dfFromGreenwich=(0.0);

        if(pmName.equalsIgnoreCase("Greenwich")||(Math.abs(dfFromGreenwich-pmValue)<tolerance))

        {

	    else if( Math.abs(dfSemiMajor-6378165.0) < 0.01

	             && Math.abs(dfInvFlattening-298.3) < 0.0001 )

	    {

	    }

	    else if( Math.abs(dfSemiMajor-6378145.0) < 0.01

	             && Math.abs(dfInvFlattening-298.25) < 0.0001 )

	    {

	    }

	    else if( Math.abs(dfSemiMajor-6378135.0) < 0.01

	             && Math.abs(dfInvFlattening-298.26) < 0.0001 )

	    {

	    }

	    else if( Math.abs(dfSemiMajor-6378137.0) < 0.01

	             && Math.abs(dfInvFlattening-298.257223563) < 0.000001 )

	    /*

	    else if( EQUAL(pszDatum,"North_American_Datum_1927") )

	    {

	        pszPROJ4Ellipse = "clrk66";

	    }

	    else if( EQUAL(pszDatum,"North_American_Datum_1983") )

	    {

	        pszPROJ4Ellipse = "GRS80";

	    }

	    */

	private String datumToProj4(String datumName,int epsgCode)

	{

	    String datumProj4="";

	    if(datumName.equals(""))

	    	datumProj4="";

		else if(datumName.equalsIgnoreCase(SRS_DN_NAD27) || epsgCode == 6267 )

	    else if(datumName.equalsIgnoreCase(SRS_DN_NAD83) || epsgCode == 6269 )

	    	datumProj4 = "+datum=NAD83 ";

	    	datumProj4 = "+datum=WGS84 ";

	    else if( epsgCode == 6314 )

	// - Transverse Mercator

	// - Cuando el meridiano origen no es 0 se pone pm=madrid y lon_0=

	// +datum

	// Funciones privadas necesarias para el proceso

	private String getName(Identifier name) {

		String[] correctName = name.toString().split(":");

			return correctName[0];

		else

			return correctName[1];

	}

	private String[] Spheroid (Ellipsoid ellips) {

		String[] spheroid = new String[3];

		Unit u = ellips.getAxisUnit();

		spheroid[2] = String.valueOf(inv_f);

		return spheroid;

	}

	private String[] Primem (PrimeMeridian prim) {

		String[] primem = new String[2];

		DefaultPrimeMeridian pm = (DefaultPrimeMeridian) prim;

		primem[1] = String.valueOf(value);

		return primem;

	}

	public double convert(double value, String measure) throws ConversionException {

			value *= this.divider;

	        int deg,min;

	        deg = (int) (value/10000); value -= 10000*deg;

	        value = ((value/60) + min)/60 + deg;

	        return value;

		}

			return value;

			return (value*0.3048/1.0);

		throw new ConversionException("Conversion no contemplada: "+measure);

    }

	public ArrayList getProjectionNameList(){

		return projectionNameList;

	}