Application: gvSIG desktop

 /**********************************************************************

 * $Id: crsJniProj.c,v 1.9 2006/09/07 10:38:13

 *

 * Name:     crsJniProj.c

 * Project:  Proj-4.4.9. Cartographic Projections library(Gerald Evenden).

 * Purpose:  Projections Basic Funcions.

 * Author:   Miguel Garcia Jimenez, garciajimenez.miguel@gmail.com

 * Collaborator: Jose Luis Gomez Martinez, jolugomar@gmail.com

 *               David Hernandez Lopez , dhernan@agr-ab.uclm.es

 *

 **********************************************************************/

/* gvSIG. Sistema de Informacion Geografica de la Generalitat Valenciana

*

* Copyright (C) 2004 IVER T.I. and Generalitat Valenciana.

*

* 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.

*

* This program 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 this program; if not, write to the Free Software

* Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307,USA.

*

*/

#include "projects.h"

#include <jni.h>

#include <proj_api.h>

#define arraysize 300

//Realizar una operacion entre CRS, pasando como parametros tres arrays de doubles(puntos)

JNIEXPORT int JNICALL Java_org_gvsig_crs_proj_JNIBaseOperation_operation

  (JNIEnv * env, jobject parent, jdoubleArray firstcoord, jdoubleArray secondcoord, jdoubleArray values, jlong src, jlong dest)

{

 int i;

 projPJ src_pj, dst_pj; // comprobar si se puede quitar

 src_pj = src;

 dst_pj = dest;

 double *xcoord = (* env)-> GetDoubleArrayElements(env, firstcoord, NULL);

 double *ycoord = (* env) -> GetDoubleArrayElements(env, secondcoord, NULL);

 double *zcoord = (* env) -> GetDoubleArrayElements(env, values, NULL);

 jint sizeofdata = (*env)->GetArrayLength(env, firstcoord);

 if(src_pj->is_latlong == 1)

 {

     for(i = 0;i<sizeofdata;i++)

     {

      *xcoord = *xcoord / (180/PI);

      *ycoord = *ycoord / (180/PI);

      pj_transform( src_pj, dst_pj,1,1, xcoord, ycoord, zcoord);

      if(dst_pj->is_latlong == 1)

      {

           *xcoord = *xcoord * (180/PI);

           *ycoord = *ycoord * (180/PI);

       }

       xcoord++;

       ycoord++;

       zcoord++;

      }

 }

 if(src_pj->is_latlong == 0)

  {

    for(i = 0;i<sizeofdata;i++)

     {

      pj_transform( src_pj, dst_pj,1,1, xcoord, ycoord, zcoord);

      if(dst_pj->is_latlong == 1)

      {

         *xcoord = *xcoord * (180/PI);

         *ycoord = *ycoord * (180/PI);

      }

      xcoord++;

      ycoord++;

      zcoord++;

     }

  }

 xcoord = xcoord - sizeofdata;

 ycoord = ycoord - sizeofdata;

 zcoord = zcoord - sizeofdata;

 (* env)->ReleaseDoubleArrayElements(env,firstcoord,(jdouble *) xcoord,JNI_COMMIT);

 (* env)->ReleaseDoubleArrayElements(env,secondcoord,(jdouble *) ycoord,JNI_COMMIT);

 (* env)->ReleaseDoubleArrayElements(env,values,(jdouble *) zcoord,JNI_COMMIT);

 return 1; // pendiente de analizar error

}

//Realizar una operacion entre CRS, pasando como parametros tres doubles(puntos)

JNIEXPORT int JNICALL Java_org_gvsig_crs_proj_JNIBaseOperation_operationSimple

  (JNIEnv * env, jobject parent, jdouble firstcoord, jdouble secondcoord, jdouble values, jlong src, jlong dest)

{

 int i;

 projPJ src_pj, dst_pj; // comprobar si se puede quitar

 src_pj = src;

 dst_pj = dest;

 double x,y,z;

 x=firstcoord;

 y=secondcoord;

 z=values;

 double *xcoord;

 double *ycoord;

 double *zcoord;

 *xcoord = firstcoord;

 *ycoord = secondcoord;

 *zcoord = values;

 if(src_pj->is_latlong == 1)

 {

   *xcoord = *xcoord / (180/PI);

   *ycoord = *ycoord / (180/PI);

 }

 pj_transform( src_pj, dst_pj,1,1, xcoord,ycoord,zcoord);

 if(dst_pj->is_latlong == 1)

 {

   *xcoord = *xcoord * (180/PI);

   *ycoord = *ycoord * (180/PI);

 }

 firstcoord =*xcoord;

 secondcoord =*ycoord;

 values =*zcoord;

 return 1; // pendiente de analizar error

}

//Realizar una operacion entre CRS, pasando como parametro un array de tres dobles(puntos)

JNIEXPORT int JNICALL Java_org_gvsig_crs_proj_JNIBaseOperation_operationArraySimple

  (JNIEnv * env, jobject parent, jdoubleArray coord, jlong src, jlong dest)

{

 projPJ src_pj, dst_pj; // comprobar si se puede quitar

 src_pj = src;

 dst_pj = dest;

 double *zcoord = (* env)-> GetDoubleArrayElements(env, coord, NULL);

 double *xcoord = zcoord++;

 double *ycoord = zcoord++;

 if(src_pj->is_latlong == 1)

 {

   *xcoord = *xcoord / (180/PI);

   *ycoord = *ycoord / (180/PI);

 }

 pj_transform( src_pj, dst_pj,1,1, xcoord, ycoord, zcoord);

 if(dst_pj->is_latlong == 1)

 {

    *xcoord = *xcoord * (180/PI);

    *ycoord = *ycoord * (180/PI);

 }

 (* env)->ReleaseDoubleArrayElements(env,coord,(jdouble *) xcoord,JNI_COMMIT);

 return 1; // pendiente de analizar error

}

//Crea un CRS a partir de una cadena proj4, devolviendo la direccion de memoria del CRS

JNIEXPORT projPJ JNICALL Java_org_gvsig_crs_proj_JNIBaseCrs_loadCrs

  (JNIEnv * env, jobject parent, jstring strCrs)

{

        projPJ src_pj;

        char * srcproj_def = (char *) (*env)->GetStringUTFChars (env, strCrs, 0);

        if (!(src_pj = pj_init_plus(srcproj_def))) return 0;

        return src_pj;

}

//Libera la direccion de memoria del CRS creado

JNIEXPORT void JNICALL Java_org_gvsig_crs_proj_JNIBaseCrs_freeCrs

  (JNIEnv * env, jobject parent, jlong crs )

{

    projPJ proj_pj;

    proj_pj = crs;

        pj_free( proj_pj );

}

//Devuelve un 1 si el CRS es latlong,  un 0 en caso contrario o un codigo de error (negativo) si el crs no está bien creado.

JNIEXPORT int JNICALL Java_org_gvsig_crs_proj_JNIBaseCrs_isLatlong

  (JNIEnv * env, jobject parent, jlong crs)

  {

          projPJ proj_pj = crs;

    int *errno = pj_get_errno_ref();

    if(*errno>=0)

        return proj_pj->is_latlong;

    else

        return *errno;

  }

//Devuelve el codigo de error correspondiente al ultimo crs creado.

JNIEXPORT int JNICALL Java_org_gvsig_crs_proj_JNIBaseCrs_getErrno

  (JNIEnv * env, jobject parent)

  {

    int *errno = pj_get_errno_ref();

    return *errno;

  }

  //Traduce un codigo de error a su correspondiente mensaje.

  JNIEXPORT jstring JNICALL Java_org_gvsig_crs_proj_JNIBaseCrs_strErrno

  (JNIEnv * env, jobject parent, jint errno)

  {

    int err = errno;

    char * strerr = pj_strerrno(err);

    return (*env)->NewStringUTF(env,strerr);

  }

/*!

 * \brief

 * retrieves projection parameters

 *

 * JNI informations:

 * Class:     org_proj4_Projections

 * Method:    getProjInfo

 * Signature: (Ljava/lang/String;)Ljava/lang/String;

 *

 *

 * \param env - parameter used by jni (see JNI specification)

 * \param parent - parameter used by jni (see JNI specification)

 * \param projdefinition - definition of the projection

*/

/*

JNIEXPORT jstring JNICALL Java_org_proj4_Projections_getProjInfo

  (JNIEnv * env, jobject parent, jstring projdefinition)

{

 PJ *pj;

 char * pjdesc;

 char info[arraysize];

 char * proj_def = (char *) (*env)->GetStringUTFChars (env, projdefinition, 0);

 if (!(pj = pj_init_plus(proj_def)))

  exit(1);

 // put together all the info of the projection and free the pointer to pjdesc

 pjdesc = pj_get_def(pj, 0);

 strcpy(info,pjdesc);

 pj_dalloc(pjdesc);

 return (*env)->NewStringUTF(env,info);

}

*/

/*!

 * \brief

 * retrieves ellipsoid parameters

 *

 * JNI informations:

 * Class:     org_proj4_Projections

 * Method:    getEllipsInfo

 * Signature: (Ljava/lang/String;)Ljava/lang/String;

 *

 *

 * \param env - parameter used by jni (see JNI specification)

 * \param parent - parameter used by jni (see JNI specification)

 * \param projdefinition - definition of the projection

*/

/*

JNIEXPORT jstring JNICALL Java_org_proj4_Projections_getEllipsInfo

  (JNIEnv * env, jobject parent, jstring projdefinition)

{

 PJ *pj;

 char * pjdesc;

 char ellipseinfo[arraysize];

 char temp[50];

 char * proj_def = (char *) (*env)->GetStringUTFChars (env, projdefinition, 0);

 if (!(pj = pj_init_plus(proj_def)))

  exit(1);

 // put together all the info of the ellipsoid

//  sprintf(temp,"name: %s;", pj->descr);

 sprintf(temp,"name: not available;");

 strcpy(ellipseinfo,temp);

 sprintf(temp,"a: %lf;", pj->a);

 strcat(ellipseinfo,temp);

 sprintf(temp,"e: %lf;", pj->e);

 strcat(ellipseinfo,temp);

 sprintf(temp,"es: %lf;", pj->es);

 strcat(ellipseinfo,temp);

 sprintf(temp,"ra: %lf;", pj->ra);

 strcat(ellipseinfo,temp);

 sprintf(temp,"one_es: %lf;", pj->one_es);

 strcat(ellipseinfo,temp);

 sprintf(temp,"rone_es: %lf;", pj->rone_es);

 strcat(ellipseinfo,temp);

 sprintf(temp,"lam0: %lf;", pj->lam0);

 strcat(ellipseinfo,temp);

 sprintf(temp,"phi0: %lf;", pj->phi0);

 strcat(ellipseinfo,temp);

 sprintf(temp,"x0: %lf;", pj->x0);

 strcat(ellipseinfo,temp);

 sprintf(temp,"y0: %lf;", pj->y0);

 strcat(ellipseinfo,temp);

 sprintf(temp,"k0: %lf;", pj->k0);

 strcat(ellipseinfo,temp);

 sprintf(temp,"to_meter: %lf;", pj->to_meter);

 strcat(ellipseinfo,temp);

 sprintf(temp,"fr_meter: %lf;", pj->fr_meter);

 strcat(ellipseinfo,temp);

 return (*env)->NewStringUTF(env,ellipseinfo);

}

*/