Application: gvSIG desktop

/*

 * Created on 10-abr-2006

 *

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

 *

 * For more information, contact:

 *

 *  Generalitat Valenciana

 *   Conselleria d'Infraestructures i Transport

 *   Av. Blasco Ib??ez, 50

 *   46010 VALENCIA

 *   SPAIN

 *

 *      +34 963862235

 *   gvsig@gva.es

 *      www.gvsig.gva.es

 *

 *    or

 *

 *   IVER T.I. S.A

 *   Salamanca 50

 *   46005 Valencia

 *   Spain

 *

 *   +34 963163400

 *   dac@iver.es

 */

/* CVS MESSAGES:

*

* $Id: 

* $Log: 

*/

package org.gvsig.referencing;

import java.util.ArrayList;

import java.util.List;

import java.util.Map;

import javax.vecmath.MismatchedSizeException;

import org.geotools.referencefork.geometry.DirectPosition2D;

import org.geotools.referencefork.referencing.operation.builder.AdvancedAffineBuilder;

import org.geotools.referencefork.referencing.operation.builder.AffineTransformBuilder;

import org.geotools.referencefork.referencing.operation.builder.BursaWolfTransformBuilder;

import org.geotools.referencefork.referencing.operation.builder.MappedPosition;

import org.geotools.referencefork.referencing.operation.builder.MathTransformBuilder;

import org.geotools.referencefork.referencing.operation.builder.ProjectiveTransformBuilder;

import org.geotools.referencefork.referencing.operation.builder.RubberSheetBuilder;

import org.geotools.referencefork.referencing.operation.builder.SimilarTransformBuilder;

import org.gvsig.util.ReferencingUtil;

import org.opengis.referencing.FactoryException;

import org.opengis.referencing.operation.TransformException;

import org.opengis.spatialschema.geometry.DirectPosition;

import org.opengis.spatialschema.geometry.MismatchedDimensionException;

import org.opengis.spatialschema.geometry.MismatchedReferenceSystemException;

public class TransformBuilderFactory {

        //least squares methods

        public static final int PROJECTIVE_TRANSFORM = 0;

        public static final int AFFINE_TRANSFORM = 1;

        public static final int SIMILAR_TRANSFORM = 2;

        public static final int BURSA_WOLF_TRANSFORM = 3;

        //TIN 

        public static final int RUBBER_SHEET_TRANSFORM = 4;

        //ADVANCED TRANSFORMATIONS (GRID BASED)

        public static final int ADVANCED_AFFINE_TRANSFORM = 5;

//        public static final int IDW_WARP_GRID = 6;

//        public static final int NADCON_GRID = 7;

//        public static final int RS_TRANSFORM = 8;

//        public static final int TPS_GRID = 9;

        public static MathTransformBuilder createTransformBuilder(int transformType, 

                                                                                                        MappedPositionContainer mappedPositionsContainer, 

                                                                                                        Map<String, Object> params) 

                                throws MismatchedDimensionException, 

                                           MismatchedReferenceSystemException, 

                                           MismatchedSizeException, 

                                           FactoryException, 

                                           TransformException{

                MathTransformBuilder solution = null;

                List<MappedPosition> mappedPositions = mappedPositionsContainer.getAsList();

                //TODO SUSTITUIR EL SWITCH POR UN HASHMAP, Y POR EL REGISTRO EST?TICO

                //DE M?TODOS DE TRANSFORMACION MEDIANTE PUNTOS DE EXTENSION

                switch(transformType){

                case PROJECTIVE_TRANSFORM:

                        solution = new ProjectiveTransformBuilder(mappedPositions);

                        break;

                case AFFINE_TRANSFORM:

                        solution = new AffineTransformBuilder(mappedPositions);

                        break;

                case SIMILAR_TRANSFORM:

                        solution = new SimilarTransformBuilder(mappedPositions);

                        break;

                case BURSA_WOLF_TRANSFORM:

                        solution = new BursaWolfTransformBuilder(mappedPositions);

                        break;

                case RUBBER_SHEET_TRANSFORM:

                        List<DirectPosition> roi = createRoi(mappedPositions);

                        solution = new RubberSheetBuilder(mappedPositions, roi);

                        break;

                case ADVANCED_AFFINE_TRANSFORM:

                        solution = new AdvancedAffineBuilder(mappedPositions);

                        //aqu? habr?a que a?adir los diferentes constre?imientos con

                        //solution.setConstraint(String, Object)

                        break;

//                case IDW_WARP_GRID:

//                        solution = new IDWGridBuilder(mappedPositions, 0d, 0d, null);

//                        break;

//                case NADCON_GRID:

//                        solution = new NADCONBuilder(mappedPositions, 0d, 0d, null);

//                        break;

//                case RS_TRANSFORM:

//                        solution = new RSGridBuilder(mappedPositions, 0d, 0d, null, null);

//                        break;

//                case TPS_GRID:

//                        solution = new TPSGridBuilder(mappedPositions, 0d, 0d, null, null);

//                        break;

                default:

                        solution = new AffineTransformBuilder(mappedPositions);

                        break;

                }

                return solution;

        }

        public static List<DirectPosition> createRoi(List<MappedPosition> mappedPositions){

                List<DirectPosition> solution = new ArrayList<DirectPosition>();

                double minX = Double.MAX_VALUE, minY = Double.MAX_VALUE, maxX = Double.MIN_VALUE, maxY = Double.MIN_VALUE;

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

                        MappedPosition position = mappedPositions.get(i);

                        DirectPosition source = position.getSource();

                        double[] sourceCoords = source.getCoordinates();

                        if(sourceCoords[0] > maxX)

                                maxX = sourceCoords[0];

                        if(sourceCoords[0] < minX)

                                minX = sourceCoords[0];

                        if(sourceCoords[1] > maxY)

                                maxY = sourceCoords[1];

                        if(sourceCoords[1] < minY)

                                minY = sourceCoords[1];

                        DirectPosition target = position.getTarget();

                        double[] targetCoords = target.getCoordinates();

                        if(targetCoords[0] > maxX)

                                maxX = targetCoords[0];

                        if(targetCoords[0] < minX)

                                minX = targetCoords[0];

                        if(targetCoords[1] > maxY)

                                maxY = targetCoords[1];

                        if(targetCoords[1] < minY)

                                minY = targetCoords[1];

                }

                DirectPosition a = new DirectPosition2D(minX, minY);

                DirectPosition b = new DirectPosition2D(minX, maxY);

                DirectPosition c = new DirectPosition2D(maxX, minY);

                DirectPosition d = new DirectPosition2D(maxX, maxY);

                //azo: this code has been modified for compatibility problem between

                //different versions of geoapi (2.0 used in libFMap, and 2.2 needed for

                //topology

//                DirectPosition a = ReferencingUtil.getInstance().create(new double[]{minX, minY}, null);

//                DirectPosition b = ReferencingUtil.getInstance().create(new double[]{minX, maxY}, null);

//                DirectPosition c = ReferencingUtil.getInstance().create(new double[]{maxX, minY}, null);

//                DirectPosition d = ReferencingUtil.getInstance().create(new double[]{maxX, maxY}, null);

                solution.add(a);

                solution.add(b);

                solution.add(c);

                solution.add(d);

                return solution;

        }

}