Application: gvSIG desktop

/* gvSIG. Sistema de Informaci?n Geogr?fica de la Generalitat Valenciana

         *

         * Copyright (C) 2006 Instituto de Desarrollo Regional 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.

         *

         * For more information, contact:

         *

         *  Generalitat Valenciana

         *   Conselleria d'Infraestructures i Transport

         *   Av. Blasco Iba?ez, 50

         *   46010 VALENCIA

         *   SPAIN

         *

         *      +34 963862235

         *   gvsig@gva.es

         *      www.gvsig.gva.es

         *

         *    or

         *

         *   Instituto de Desarrollo Regional (Universidad de Castilla La-Mancha)

         *   Campus Universitario s/n

         *   02071 Alabacete

         *   Spain

         *

         *   +34 967 599 200

         */

package org.gvsig.georeferencing.process;

import java.awt.geom.AffineTransform;

import java.io.IOException;

import org.gvsig.fmap.raster.layers.FLyrRasterSE;

import org.gvsig.georeferencing.process.geotransform.GeoTransformDataResult;

import org.gvsig.georeferencing.process.geotransform.GeoTransformProcess;

import org.gvsig.raster.IProcessActions;

import org.gvsig.raster.RasterProcess;

import org.gvsig.raster.buffer.BufferFactory;

import org.gvsig.raster.buffer.RasterBuffer;

import org.gvsig.raster.buffer.RasterBufferInvalidException;

import org.gvsig.raster.dataset.GeoRasterWriter;

import org.gvsig.raster.dataset.IBuffer;

import org.gvsig.raster.dataset.IRasterDataSource;

import org.gvsig.raster.dataset.NotSupportedExtensionException;

import org.gvsig.raster.dataset.io.RasterDriverException;

import org.gvsig.raster.datastruct.GeoPoint;

import org.gvsig.raster.grid.Grid;

import org.gvsig.raster.grid.GridException;

import org.gvsig.raster.grid.GridExtent;

import org.gvsig.raster.grid.OutOfGridException;

import org.gvsig.raster.util.RasterToolsUtil;

import com.iver.andami.PluginServices;

/**

 *  Clase que representa una proceso de georreferenciacion de un raster.

 *  

 *  @author Alejandro Mu?oz Sanchez (alejandro.munoz@uclm.es)

 *         @version 10/2/2008

 **/

public class GeoreferencingProcess extends RasterProcess implements IProcessActions{

        //Capa a georreferenciar

        private FLyrRasterSE rasterSE                        = null;

        //Grid resultante de georreferenciacion

        private Grid imageGrid = null;

        //Fichero de salida

        private String filename = null;

        // resutado de la transformacion 

        private GeoTransformDataResult transformData=null;

        // Lista puntos de control

        private GeoPoint gpcs[] =null;

        //Extend de imagen corregida

        GridExtent newExtend=null;

        // Metodo de resampleado utilizado

        private int rMethod= 0;

        //Indicador de progreso

        private int percent=0;

        // Grid resultado

        private Grid gridResult=null;

        WriterBufferServer writerBufferServer=null;

        /** Metodo que recoge los parametros del proceso georreferenciacion de un raster

        * <LI>rasterSE: capa a georeferenciar</LI>

        * <LI>filename: path con el fichero de salida</LI>

        * <LI>method: metodo de resampleo </LI>

        */

        public void init() {

                rasterSE= (FLyrRasterSE)getParam("fLayer");

                filename= (String)getParam("filename");

                rMethod=(int)getIntParam("method");

                gpcs = (GeoPoint[])getParam("gpcs");

                // Inicializacion del grid correspondiente a la imagen a corregir

                IRasterDataSource dsetCopy = null; 

                dsetCopy = rasterSE.getDataSource().newDataset();

                BufferFactory bufferFactory = new BufferFactory(dsetCopy);

                if (!RasterBuffer.loadInMemory(dsetCopy))

                        bufferFactory.setReadOnly(true);        

                try {

                                imageGrid = new Grid(bufferFactory);        

                }catch (RasterBufferInvalidException e) {

                        e.printStackTrace();                        

                }

        }

        public void process() throws InterruptedException {

                //COMENTAR:         Podiamos ahorrarlo si le pasamos  ya un  geoTransformDataResult

                GeoTransformProcess transform = new GeoTransformProcess();

                transform.setActions(this);

                transform.addParam("gpcs",gpcs);

                transform.addParam("orden",new Integer(1));

                transform.run();                

                // Obtenida la transformacion la aplicamos a los puntos extremos de la imagen

                double p1[]=transformData.getCoordMap(imageGrid.getGridExtent().getLLX(),imageGrid.getGridExtent().getLLY());

                double p2[]=transformData.getCoordMap(imageGrid.getGridExtent().getULX(),imageGrid.getGridExtent().getULY());

                double p3[]=transformData.getCoordMap(imageGrid.getGridExtent().getURX(),imageGrid.getGridExtent().getURY());

                double p4[]=transformData.getCoordMap(imageGrid.getGridExtent().getLRX(),imageGrid.getGridExtent().getLRY());

                double xmin=Math.min( Math.min(p1[0],p2[0]),Math.min(p3[0],p4[0]));

                double ymin=Math.min( Math.min(p1[1],p2[1]),Math.min(p3[1],p4[1]));

                double xmax=Math.max( Math.max(p1[0],p2[0]),Math.max(p3[0],p4[0]));        

                double ymax=Math.max( Math.max(p1[1],p2[1]),Math.max(p3[1],p4[1]));

                // OJO CON EL TAMA?O DE CELDA: CONSULTAR PARA EJEMPLO DE JOSE ES 25

                newExtend= new GridExtent(xmin,ymin,xmax,ymax,25);

                int datatype= rasterSE.getBufferFactory().getRasterBuf().getDataType();

                int bands[]= rasterSE.getBufferFactory().getDrawableBands();

                try {

                        gridResult= new Grid(newExtend,newExtend,datatype,bands);

                } catch (RasterBufferInvalidException e) {

                        RasterToolsUtil.messageBoxError("error_grid", this, e);

                }

                double minPointX=gridResult.getGridExtent().getMin().getX();

                double minPointY=gridResult.getGridExtent().getMin().getY();

                double cellsizeX=gridResult.getCellSize();

                double cellsizeY=gridResult.getCellSize();

                double coord[]=null;

                try {

                        // Caso imagen tipo byte

                        if(datatype==IBuffer.TYPE_BYTE)

                        {

                                for(int fil=0; fil<gridResult.getLayerNX(); fil++){

                                        for(int col=0; col<gridResult.getLayerNY();col++)

                                        {

                                                coord=transform.getCoordPixel(fil*cellsizeX+minPointX, col*cellsizeY+minPointY);        

                                                double value;

                                                value = imageGrid.getCellValueAsByte((int)coord[0],(int)coord[1]);

                                                gridResult.setCellValue(fil,gridResult.getLayerNY()-col-1,(byte)value);

                                        }

                                        percent= (int)fil*100/gridResult.getLayerNX();        

                                }        

                        }

                        // Caso imagen tipo short

                        if(datatype==IBuffer.TYPE_SHORT)

                        {

                                for(int fil=0; fil<gridResult.getLayerNX(); fil++){

                                        for(int col=0; col<gridResult.getLayerNY();col++)

                                        {

                                                coord=transform.getCoordPixel(fil*cellsizeX+minPointX, col*cellsizeY+minPointY);        

                                                short value;

                                                value = imageGrid.getCellValueAsByte((int)coord[0],(int)coord[1]);

                                                gridResult.setCellValue(fil,gridResult.getLayerNY()-col-1,(short)value);

                                        }

                                        percent= (int)fil*100/gridResult.getLayerNX();        

                                }        

                        }

                        //        Caso imagen tipo int

                        if(datatype==IBuffer.TYPE_INT)

                        {

                                for(int fil=0; fil<gridResult.getLayerNX(); fil++){

                                        for(int col=0; col<gridResult.getLayerNY();col++)

                                        {

                                                coord=transform.getCoordPixel(fil*cellsizeX+minPointX, col*cellsizeY+minPointY);        

                                                short value;

                                                value = imageGrid.getCellValueAsByte((int)coord[0],(int)coord[1]);

                                                gridResult.setCellValue(fil,gridResult.getLayerNY()-col-1,(int)value);

                                        }

                                        percent= (int)fil*100/gridResult.getLayerNX();        

                                }        

                        }

                        // Caso imagen tipo float

                        if(datatype==IBuffer.TYPE_FLOAT)

                        {

                                for(int fil=0; fil<gridResult.getLayerNX(); fil++){

                                        for(int col=0; col<gridResult.getLayerNY();col++)

                                        {

                                                coord=transform.getCoordPixel(fil*cellsizeX+minPointX, col*cellsizeY+minPointY);        

                                                double value;

                                                value = imageGrid.getCellValueAsByte((int)coord[0],(int)coord[1]);

                                                gridResult.setCellValue(fil,gridResult.getLayerNY()-col-1,(float)value);

                                        }

                                        percent= (int)fil*100/gridResult.getLayerNX();        

                                }        

                        }

                        // Caso imagen tipo double

                        if(datatype==IBuffer.TYPE_FLOAT)

                        {

                                for(int fil=0; fil<gridResult.getLayerNX(); fil++){

                                        for(int col=0; col<gridResult.getLayerNY();col++)

                                        {

                                                coord=transform.getCoordPixel(fil*cellsizeX+minPointX, col*cellsizeY+minPointY);        

                                                double value;

                                                value = imageGrid.getCellValueAsByte((int)coord[0],(int)coord[1]);

                                                gridResult.setCellValue(fil,gridResult.getLayerNY()-col-1,(double)value);

                                        }

                                        percent= (int)fil*100/gridResult.getLayerNX();        

                                }        

                        }

                } catch (GridException e) {

                                e.printStackTrace();        

                } catch (OutOfGridException e) {

                        e.printStackTrace();

                }

                generateLayer();

                if(externalActions!=null)

                        externalActions.end(filename);

        }

        private void generateLayer(){

                GeoRasterWriter grw = null;

                IBuffer buffer= gridResult.getRasterBuf();

            writerBufferServer = new WriterBufferServer(buffer);

                AffineTransform aTransform = new AffineTransform(newExtend.getCellSize(),0.0,0.0,-newExtend.getCellSize(),newExtend.getMin().getX(),newExtend.getMax().getY());

                int endIndex =filename.lastIndexOf(".");

                if (endIndex < 0)

                        endIndex = filename.length();

                try {

                        grw = GeoRasterWriter.getWriter(writerBufferServer, filename,gridResult.getRasterBuf().getBandCount(),aTransform, gridResult.getRasterBuf().getWidth(),gridResult.getRasterBuf().getHeight(), gridResult.getRasterBuf().getDataType(), GeoRasterWriter.getWriter(filename).getParams(),null);

                        grw.dataWrite();

                        grw.setWkt(rasterSE.getWktProjection());

                        grw.writeClose();

                } catch (NotSupportedExtensionException e1) {

                        RasterToolsUtil.messageBoxError(PluginServices.getText(this, "error_writer_notsupportedextension"), this, e1);

                } catch (RasterDriverException e1) {

                        e1.printStackTrace();

                } catch (IOException e) {

                        RasterToolsUtil.messageBoxError(PluginServices.getText(this, "error_writer"), this, e);

                } catch (InterruptedException e) {

                        Thread.currentThread().interrupt();

                }

        }

        // Metodos para la interpolacion

        /*

         * (non-Javadoc)

         * @see org.gvsig.gui.beans.incrementabletask.IIncrementable#getTitle()

         */

        public String getTitle() {

                return PluginServices.getText(this,"georreferenciacion_process");

        }

        /*

         * (non-Javadoc)

         * @see org.gvsig.gui.beans.incrementabletask.IIncrementable#getLabel()

         */

        public int getPercent() {

                if(writerBufferServer==null)

                        return percent;

                else

                        return writerBufferServer.getPercent();

        }

        public void interrupted() {

                // TODO Auto-generated method stub

        }

        public void end(Object param) {        

                transformData = (GeoTransformDataResult)param;        

        }

}