| gvSIG desktop 1

/*

 * Created on 19-jul-2006

 *

 * To change the template for this generated file go to

 * Window>Preferences>Java>Code Generation>Code and Comments

 */

/* 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.cresques.io.data;

import junit.framework.TestCase;

import org.cresques.io.GeoRasterFile;

/**

 * Este test prueba la gesti?n de multifichero de la clase GeoRasterMultiFile

 * y que incorpora un Grid, as? como la gesti?n de bandas que tiene BandList

 * y que incorpora GeoRasterMultiFile.

 * @author Nacho Brodin (nachobrodin@gmail.com)

 *

 */

public class TestDataSourceWithNoDataAndSizeBuffer extends TestCase{

	private String baseDir = "/home/nacho/images/";

	private String path1 = baseDir + "miniRaster25x24.tif";

	private GeoRasterFile f1 = null;

	private Grid g = null;

	public void setUp() {

		System.out.println("TestDataSourceWithNoDataAndSizeBuffer running...");

		f1 = GeoRasterFile.openFile(null, path1);

	}

	public void testStack(){

		int[] drawableBands = {0, 1, 2};

		g = new Grid(f1);

		g.addDrawableBands(drawableBands);

		//El par?metro false indica que no se ajuste el extent al raster sino que se rellene con NoData los huecos

		//TODO: Hacer las llamadas para probar setAreaOfInterest con tama?o de buffer

		g.setAreaOfInterest(645819.2583, 4923849.7, 25.584, 25.2); //ZOOM

		//g.setAreaOfInterest(645818.2583, 4923850.7, 34.584, 33.2);

		//g.setAreaOfInterest(645808.2583, 4923850.7, 44.584, 33.2); //NODATA IZDA

		//g.setAreaOfInterest(645818.2583, 4923860.7, 34.584, 43.2); //NODATA ARRIBA

		//g.setAreaOfInterest(645818.2583, 4923850.7, 34.584, 43.2); //NODATA ABAJO

		//g.setAreaOfInterest(645818.2583, 4923850.7, 44.584, 33.2); //NODATA DCHA

		//g.setAreaOfInterest(645808.2583, 4923860.7, 54.584, 53.2); //NODATA 4 LADOS

		RasterBuf raster = g.getRasterBuf();

		for(int line = 0; line < raster.getHeight(); line++){

			for(int col = 0; col < raster.getWidth(); col++)

				System.out.print("("+raster.getElemByte(line, col, 0)+" "+raster.getElemByte(line, col, 1)+" "+raster.getElemByte(line, col, 2)+")");

			System.out.println();

		}

	}

}

		Extent fitExtent = Utilities.calculateAdjustedView(dataExtent, grmf.getExtent());

		//Caso en que la extensi?n pedida SI se sale fuera del ?rea del raster

		//El pixel a medias entre NoData y datos se pone NoData

		if(!Utilities.isInside(dataExtent, fitExtent)){

			//Upper Left

			double distWcX = Math.abs(fitExtent.getMin().getX() - dataExtent.getMin().getX());

			double distWcY = Math.abs(fitExtent.getMax().getY() - dataExtent.getMax().getY());

			//Pixel inicial del buffer donde se empieza a dibujar. Redondeamos por arriba pq lo que sobra se pone NoData

			double initPxX = Math.ceil((distWcX * bufWidth) / (Math.abs(maxX - minX))); 

			double initPxY = Math.ceil((distWcY * bufHeight) / (Math.abs(maxY - minY)));

			//Obtenemos la coordenada real para el pixel inicial

			double wcXInit = (initPxX * (Math.abs(maxX - minX))) / bufWidth;

			double wcYInit = (initPxY * (Math.abs(maxY - minY))) / bufHeight;

			//Lower Right

			distWcX = Math.abs(fitExtent.getMax().getX() - dataExtent.getMin().getX());

			distWcY = Math.abs(fitExtent.getMin().getY() - dataExtent.getMax().getY());

			//Pixel final del buffer donde se dibuja. Redondeamos por abajo pq lo que sobra se pone NoData

			double endPxX = Math.floor((distWcX * bufWidth) / (Math.abs(maxX - minX))); 

			double endPxY = Math.floor((distWcY * bufHeight) / (Math.abs(maxY - minY)));

			//Obtenemos la coordenada real para el pixel

			double wcXEnd = (endPxX * (Math.abs(maxX - minX))) / bufWidth;

			double wcYEnd = (endPxY * (Math.abs(maxY - minY))) / bufHeight;

			rasterBuf = grmf.getWindowRaster(wcXInit, wcYInit, wcXEnd, wcYEnd, (int)Math.abs(endPxX - initPxX), (int)Math.abs(endPxY - initPxY));

			RasterBuf buf = new RasterBuf(rasterBuf.getDataType(), bufWidth, bufHeight, rasterBuf.getBandCount(), true);

			buf.setToNoData();

			for (int iBand = 0; iBand < rasterBuf.getBandCount(); iBand++)	

				for (int row = 0; row < rasterBuf.getHeight(); row++) 

					for (int col = 0; col < rasterBuf.getWidth(); col++) 

						buf.setElemByte((int)(row + initPxY), (int)(col + initPxX), 

										iBand, 

										rasterBuf.getElemByte(row, col, iBand));

			rasterBuf = buf;

			return null;

		}

		//Caso en que la extensi?n pedida NO se sale fuera del ?rea del raster

		rasterBuf = grmf.getWindowRaster(fitExtent.minX(), fitExtent.maxY(), fitExtent.width(), fitExtent.height(), bufWidth, bufHeight);