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;

	}

}

/*

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

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

import com.iver.cit.gvsig.fmap.layers.FLyrVect;

/**

 * Container of MappedPosition objects, that represents digitized vector error

 * for spatial adjustment process.

 * 

 * @author Alvaro Zabala

 *

 */

public interface MappedPositionContainer {

	public void addMappedPosition(MappedPosition mappedPosition);

	public MappedPosition getMappedPosition(int idx);

	public int getCount();

	public List<MappedPosition> getAsList();

	public boolean existsLinksLyr();

	public FLyrVect getLinkLyr();

}