Application: gvSIG desktop

/* gvSIG. Sistema de Informaci?n Geogr?fica 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.

 */

package org.gvsig.georeferencing.utils;

import java.awt.geom.Point2D;

import java.awt.geom.Rectangle2D;

import java.io.File;

import java.io.FileInputStream;

import java.io.FileNotFoundException;

import java.io.FileOutputStream;

import java.io.IOException;

import java.io.InputStream;

import java.io.OutputStream;

import com.iver.cit.gvsig.fmap.ViewPort;

public class GeoUtils {

        public static double INCHESMTR = 39.73007874D;

        public static double INCHESCM = 0.3973007874D;

        public static double INCHESMM = 0.03973007874D;

        public static double MTRSINCH = 0.0254D;

        public static double MMSINCH = 25.4D;

        public static double CMSINCH = 2.54D;

        /**

     * Recorta el n?mero de decimales a n del n?mero pasado por par?metro

     * @param num Valor a recortar el n?mero de decimales

     * @param n N?mero de decimales a recortar

     * @return N?mero recortado

     */

    public static double tailDecimals(double num, int n){

            long m = (long)Math.pow(10, n);

        num *= m;

        long aux = ((long)num);

        num = (double)((double)aux / (double)m);

        return num;

    }

    /**

     * Recorta el n?mero de decimales a n del n?mero pasado por par?metro

     * @param num Valor a recortar el n?mero de decimales

     * @param n N?mero de decimales a recortar

     * @return N?mero recortado

     */

    public static String tailDecimals(String num, int n){

            int index = num.indexOf(".");

            if(index != -1){

                    try{

                            num = num.substring(0, index + n);

                    }catch(IndexOutOfBoundsException e){

                            //No hacemos nada pq retornar? num

                    }

            }

            return num;

    }

        /**

         * Convierte pixels en milimetros dependiendo del los puntos por pulgada

         * @param pixels N?mero de pixels

         * @param ppp Puntos por pulgada

         * @return N?mero de milimetros

         */

        public static double convertPixelsToMms(int pixels, int ppp){

                return convertPixelsToInches(pixels, ppp) * MMSINCH;

        }

        /**

         * Convierte pixels en centimetros dependiendo del los metros por pixel

         * @param pixels N?mero de pixels

         * @param mtsPixel Metros por pixel

         * @return N?mero de centimetros

         */

        public static double convertPixelsToCms(int pixels, int ppp){

                return convertPixelsToInches(pixels, ppp) * CMSINCH;

        }

        /**

         * Convierte pixels en metros dependiendo del los metros por pixel

         * @param pixels N?mero de pixels

         * @param mtsPixel Metros por pixel

         * @return N?mero de metros

         */

        public static double convertPixelsToMts(int pixels, int ppp){

                return convertPixelsToInches(pixels, ppp) * MTRSINCH;

        }

        /**

         * Convierte pixels en pulgadas dependiendo del los metros por pixel

         * @param pixels N?mero de pixels

         * @param mtsPixel Metros por pixel

         * @return N?mero de pulgadas

         */

        public static double convertPixelsToInches(int pixels, int ppp){

                return (double)((double)pixels / (double)ppp);

        }

        /**

         * Convierte metros en pixels dependiendo de los puntos por pulgada

         * @param mts N?mero de metros

         * @param ppp Puntos por pulgada

         * @return N?mero de pixeles

         */

        public static int convertMtsToPixels(double mts, int ppp){

                return (int)(mts * GeoUtils.INCHESMTR * ppp);

        }

        /**

         * Convierte pulgadas en pixels dependiendo de los puntos por pulgada

         * @param mts N?mero de pulgadas

         * @param ppp Puntos por pulgada

         * @return N?mero de pixeles

         */

        public static int convertInchesToPixels(double inches, int ppp){

                return (int)(inches * ppp);

        }

        /**

         * Convierte centimetros en pixels dependiendo de los puntos por pulgada

         * @param mts N?mero de centmetros

         * @param ppp Puntos por pulgada

         * @return N?mero de pixeles

         */

        public static int convertCmsToPixels(double cms, int ppp){

                return (int)(cms * GeoUtils.INCHESCM * ppp);

        }

        /**

         * Convierte milimetros en pixels dependiendo de los puntos por pulgada

         * @param mts N?mero de milimetros

         * @param ppp Puntos por pulgada

         * @return N?mero de pixeles

         */

        public static int convertMmsToPixels(double mms, int ppp){

                return (int)(mms * GeoUtils.INCHESMM * ppp);

        }

        /**

         * Devuelve la traspuesta de una matriz

         * @param M Matriz

         * @param m N?mero de filas

         * @param n N?mero de columnas

         * @return Matriz traspuesta

         */

        public static double[][] transpose(double[][] M, int m, int n){

          double[][]        a = new double[n][m];

          int                        i, j;

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

            for (j = 0; j < m; j++)

              a[i][j] = M[j][i];

          return a;

        }

        /**

         * Producto de matrices (AxB)

         * @param a        Matriz A

         * @param b Matriz B

         * @param M 

         * @param N        

         * @param O 

         * @return

         */

        public static double[][] multmatrix(double[][] a, double[][] b, int M, int N, int O){

          double[][]        c = new double[M][O]; 

          double                sum;

          int                        i, j, k;

          for (i = 0; i < M; i++) {

            for (j = 0; j < O; j++) {

              sum = 0;

              for (k = 0; k < N; k++)

                sum += b[k][j] * a[i][k];

              c[i][j] = sum;

            }

          }

          return (c);

        }

    /**

     * Calcula el zoom dependiendo del factor de escala pasado por

     * par?metro.

     */

    public static ViewPort calcZoom(ViewPort vp, double factor){

            ViewPort viewPortBlock = new ViewPort(vp.getProjection());

                viewPortBlock.setExtent(vp.getExtent());

                viewPortBlock.setImageSize(vp.getImageSize());        

                viewPortBlock.setScale();

            Rectangle2D.Double r = new Rectangle2D.Double();

                double nuevoX = (vp.getExtent().getCenterX()) - ((vp.getExtent().getWidth() * factor) / 2.0);

                double nuevoY = (vp.getExtent().getCenterY()) - ((vp.getExtent().getHeight() * factor) / 2.0);

                r.x = nuevoX;

                r.y = nuevoY;

                r.width = viewPortBlock.getExtent().getWidth() * factor;

                r.height = viewPortBlock.getExtent().getHeight() * factor;

                viewPortBlock.setExtent(r);

                viewPortBlock.setScale();

                return viewPortBlock;

        }

    /**

     * Desplaza el extent al nuevo centro 

     * @param pto

     */

    public static ViewPort shiftExtent(Point2D oldCenter, Point2D newCenter, ViewPort vp){

            double difX = oldCenter.getX() - newCenter.getX();

            double difY = oldCenter.getY() - newCenter.getY();

            ViewPort viewPortBlock = new ViewPort(vp.getProjection());

            double minX = vp.getExtent().getMinX() - difX;

            double minY = vp.getExtent().getMinY() - difY;

            Rectangle2D ex = new Rectangle2D.Double(minX, minY, vp.getExtent().getWidth(), vp.getExtent().getHeight());

                viewPortBlock.setExtent(ex);

                viewPortBlock.setImageSize(vp.getImageSize());        

                viewPortBlock.setScale();

                return viewPortBlock;

    }

    /**

     * Hace una copia de seguridad sobre el fichero que le pasamos por par?metro.

     * El fichero de copia se llama como el original pero acabado en ~. Si ya existiese

     * un fichero nombrado as? se nombraria como filename.rmf~1, filename.rmf~2, y as?

     * sucesivamente.

     * @param fName Nombre del fichero que se quiere hacer un backup.

     */

    public static void fileBackup(String fName){ 

            File backupFile = new File(fName+"~");

            int nBackupFile = 1;

            while(backupFile.exists()){

                    backupFile = new File(fName+"~"+nBackupFile);

                    nBackupFile ++;

            }

            File inFile = new File(fName);

            try{

                    InputStream in = new FileInputStream(inFile);

                    OutputStream out = new FileOutputStream(backupFile);

                    byte[] buf = new byte[1024];

                    int len;

                    while ((len = in.read(buf)) > 0) 

                            out.write(buf, 0, len);

                    in.close();

                    out.close();

            }catch(FileNotFoundException exc){

                    System.err.print("RMF Backup isn't possible: File not found, "+fName);

            }catch(IOException exc){

                    System.err.print("RMF Backup isn't possible: IOException, "+fName);

            }

    }

}