Application: gvSIG desktop » gvSIG raster

<assembly>

</assembly>

<projectDescription>

  <name>org.gvsig.raster.gdal.io</name>

  <comment>Gdal data provider for gvSIG</comment>

  <projects>

    <project>org.gvsig.raster.cache.lib.api</project>

    <project>org.gvsig.raster.lib.api</project>

    <project>org.gvsig.raster.lib.impl</project>

  </projects>

  <buildSpec>

    <buildCommand>

      <name>org.eclipse.jdt.core.javabuilder</name>

    </buildCommand>

  </buildSpec>

  <natures>

    <nature>org.eclipse.jdt.core.javanature</nature>

  </natures>

</projectDescription>

/* gvSIG. Geographic Information System of the Valencian Government

 *

 * Copyright (C) 2007-2008 Infrastructures and Transports Department

 * of the Valencian Government (CIT)

 *

 * 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., 51 Franklin Street, Fifth Floor, Boston,

 * MA  02110-1301, USA.

 *

 */

package org.gvsig.raster.gdal.io;

import org.gvsig.fmap.dal.DALFileLibrary;

import org.gvsig.fmap.dal.DALLibrary;

import org.gvsig.raster.impl.store.AbstractRasterFileDataParameters;

import org.gvsig.tools.ToolsLibrary;

import org.gvsig.tools.library.AbstractLibrary;

import org.gvsig.tools.library.Library;

import org.gvsig.tools.library.LibraryException;

/**

 *

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

 */

public class DefaultGdalIOLibrary extends AbstractLibrary {	

	public DefaultGdalIOLibrary() {

		super(DefaultGdalIOLibrary.class,Library.TYPE.IMPL);

		require(ToolsLibrary.class);

		require(DALLibrary.class);

		require(DALFileLibrary.class);

	}

	@Override

	protected void doInitialize() throws LibraryException {

		//RasterLibrary.wakeUp();

	}

	@Override

	protected void doPostInitialize() throws LibraryException {

		AbstractRasterFileDataParameters.registerDynClass();

		// Registro de los drivers de lectura

		GdalProvider.register();

		// Registro de los drivers de escritura

		GdalWriter.register();

		JpegWriter.register();

		PngWriter.register();

	}

}

	/* gvSIG. Geographic Information System of the Valencian Government

 *

 * Copyright (C) 2007-2008 Infrastructures and Transports Department

 * of the Valencian Government (CIT)

 *

 * 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., 51 Franklin Street, Fifth Floor, Boston,

 * MA  02110-1301, USA.

 *

 */

package org.gvsig.raster.gdal.io;

import java.awt.geom.AffineTransform;

import java.awt.geom.NoninvertibleTransformException;

import java.awt.geom.Point2D;

import java.io.IOException;

import org.gvsig.fmap.dal.coverage.RasterLibrary;

import org.gvsig.fmap.dal.coverage.RasterLocator;

import org.gvsig.fmap.dal.coverage.dataset.Buffer;

import org.gvsig.fmap.dal.coverage.datastruct.BandList;

import org.gvsig.fmap.dal.coverage.datastruct.Extent;

import org.gvsig.fmap.dal.coverage.exception.ProcessInterruptedException;

import org.gvsig.raster.impl.datastruct.ExtentImpl;

import org.gvsig.raster.impl.process.RasterTask;

import org.gvsig.raster.impl.process.RasterTaskQueue;

import org.gvsig.raster.impl.store.properties.DataStoreColorInterpretation;

import org.gvsig.raster.impl.store.properties.DataStoreColorTable;

import org.gvsig.raster.impl.store.properties.DataStoreMetadata;

import org.gvsig.raster.impl.store.properties.DataStoreTransparency;

import es.gva.cit.jgdal.Gdal;

import es.gva.cit.jgdal.GdalBuffer;

import es.gva.cit.jgdal.GdalException;

import es.gva.cit.jgdal.GdalRasterBand;

import es.gva.cit.jgdal.GeoTransform;

/**

 * Soporte 'nativo' para ficheros desde GDAL.

 * 

 * @author Luis W. Sevilla (sevilla_lui@gva.es)

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

 */

public class GdalNative extends Gdal {

	private String                       fileName                = null;

	private String                       shortName               = "";

	public 	GeoTransform                 trans                   = null;

	public int                           width                   = 0, height = 0;

	public double                        originX                 = 0D, originY = 0D;

	public String                        version                 = "";

	protected int                        rBandNr                 = 1, gBandNr = 2, bBandNr = 3, aBandNr = 4;

	private int[]                        dataType                = null;

	DataStoreMetadata                    metadata                = null;

	protected boolean                    georeferenced           = true;

	/**

	 * Vectores que contiene los desplazamientos de un pixel cuando hay supersampling.

	 * , es decir el n�mero de pixels de pantalla que tiene un pixel de imagen. Como todos

	 * los pixeles no tienen el mismo ancho y alto ha de meterse en un array y no puede ser

	 * una variable. Adem�s hay que tener en cuenta que el primer y �ltimo pixel son de 

	 * distinto tama�o que el resto.   

	 */

	public int[]                          stepArrayX             = null;

	public int[]                          stepArrayY             = null;

	protected GdalRasterBand[]            gdalBands              = null;

	private double                        lastReadLine           = -1;

	private int                           currentFullWidth       = -1;

	private int                           currentFullHeight      = -1;

	private int                           currentViewWidth       = -1;

	private int                           currentViewHeight      = -1;

	private double                        currentViewX           = 0D;

	private double                        viewportScaleX         = 0D;

	private double                        viewportScaleY         = 0D;

	private double                        stepX                  = 0D;

	private double                        stepY                  = 0D;

	public boolean                        isSupersampling        = false;

	private boolean                       open                   = false;

	/**

	 * Estado de transparencia del raster.

	 */

	DataStoreTransparency                 fileTransparency       = null;

	DataStoreColorTable                   palette                = null;

	DataStoreColorInterpretation          colorInterpr           = null;

	AffineTransform                       ownTransformation      = null;

	AffineTransform                       externalTransformation = new AffineTransform();

	/**

	 * Overview usada en el ?ltimo setView

	 */

	int currentOverview = -1;

	public GdalNative(String fName) throws GdalException, IOException {

		super();

		init(fName);

	}

	private void init(String fName) throws GdalException, IOException {

		fileName = fName;

		open(fName, GA_ReadOnly);

		open = true;

		if (getPtro() == -1)

			throw new GdalException("Error en la apertura del fichero. El fichero no tiene un formato v?lido.");

//		ext = RasterUtilities.getExtensionFromFileName(fName);

		width = getRasterXSize();

		height = getRasterYSize();

		int[] dt = new int[getRasterCount()];

		for (int i = 0; i < getRasterCount(); i++)

			dt[i] = this.getRasterBand(i + 1).getRasterDataType();

		setDataType(dt);

		shortName = getDriverShortName();

		fileTransparency = new DataStoreTransparency();

		colorInterpr = new DataStoreColorInterpretation();

		metadata = new DataStoreMetadata(getMetadata(), colorInterpr);

		// Asignamos la interpretaci?n de color leida por gdal a cada banda. Esto

		// nos sirve para saber que banda de la imagen va asignada a cada banda de

		// visualizaci?n (ARGB)

		colorInterpr.initColorInterpretation(getRasterCount());

		metadata.initNoDataByBand(getRasterCount());

		for (int i = 0; i < getRasterCount(); i++) {

			GdalRasterBand rb = getRasterBand(i + 1);

			String colorInt = getColorInterpretationName(rb.getRasterColorInterpretation());

			metadata.setNoDataValue(i, (rb.getRasterNoDataValue() == -9999.0) ? RasterLibrary.defaultNoDataValue : rb.getRasterNoDataValue());

			metadata.setNoDataEnabled(rb.existsNoDataValue());

			colorInterpr.setColorInterpValue(i, colorInt);

			if (colorInt.equals("Alpha"))

				fileTransparency.setTransparencyBand(i);

			if (rb.getRasterColorTable() != null && palette == null) {

				palette = new DataStoreColorTable();

				palette.createPaletteFromGdalColorTable(rb.getRasterColorTable());

//				fileTransparency.setTransparencyRangeList(palette.getTransparencyRange());

			}

		}

		fileTransparency.setTransparencyByPixelFromMetadata(metadata);

		try {

			trans = getGeoTransform();

			boolean isCorrect = false;

			for (int i = 0; i < trans.adfgeotransform.length; i++)

				if (trans.adfgeotransform[i] != 0)

					isCorrect = true;

			if (!isCorrect)

				throw new GdalException("");

			ownTransformation = new AffineTransform(trans.adfgeotransform[1], trans.adfgeotransform[4], trans.adfgeotransform[2], trans.adfgeotransform[5], trans.adfgeotransform[0], trans.adfgeotransform[3]);

			externalTransformation = (AffineTransform) ownTransformation.clone();

			currentFullWidth = width;

			currentFullHeight = height;

			this.georeferenced = true;

		} catch (GdalException exc) {

			// Transformaci�n para ficheros sin georreferenciaci�n. Se invierte la Y

			// ya que las WC decrecen de

			// arriba a abajo y los pixeles crecen de arriba a abajo

			ownTransformation = new AffineTransform(1, 0, 0, -1, 0, height);

			externalTransformation = (AffineTransform) ownTransformation.clone();

			currentFullWidth = width;

			currentFullHeight = height;

			this.georeferenced = false;

		}

	}

	/**

	 * Returns true if this provider is open and false if don't

	 * @return

	 */

	public boolean isOpen() {

		return open;

	}

	/*

	 * (non-Javadoc)

	 * @see es.gva.cit.jgdal.Gdal#close()

	 */

	public void close() throws GdalException {

		open = false;

		super.close();

	}

	/**

	 * Obtiene el flag que informa de si el raster tiene valor no data o no.

	 * Consultar� todas las bandas del mismo y si alguna tiene valor no data

	 * devuelve true sino devolver� false.

	 * @return true si tiene valor no data y false si no lo tiene

	 * @throws GdalException

	 */

	public boolean existsNoDataValue() throws GdalException {

		for (int i = 0; i < getRasterCount(); i++) {

			GdalRasterBand rb = getRasterBand(i + 1);

			if (rb.existsNoDataValue())

				return true;

		}

		return false;

	}

	/**

	 * Obtiene el flag que informa de si el raster tiene valor no data o no

	 * en una banda concreta.

	 * @return true si tiene valor no data en esa banda y false si no lo tiene

	 * @param band Posici�n de la banda a consultar (0..n)

	 * @throws GdalException

	 */

	public boolean existsNoDataValue(int band) throws GdalException {

		GdalRasterBand rb = getRasterBand(band + 1);

		return rb.existsNoDataValue();

	}

	/**

	 * Devuelve el valor NoData en caso de existir, sino existe devuelve null.

	 * @return

	 */

	public double getNoDataValue() {

		if (metadata == null)

			return RasterLibrary.defaultNoDataValue;

		if (metadata.getNoDataValue().length == 0)

			return RasterLibrary.defaultNoDataValue;

		return metadata.getNoDataValue()[0];

	}

	/**

	 * Asigna el tipo de dato

	 * @param dt entero que representa el tipo de dato

	 */

	public void setDataType(int[] dt) { 

		dataType = dt; 

	}

	/**

	 * Obtiene el tipo de dato

	 * @return entero que representa el tipo de dato

	 */

	public int[] getDataType() { 

		return dataType; 

	}

	/**

	 * Obtiene un punto 2D con las coordenadas del raster a partir de uno en coordenadas

	 * del punto real.

	 * Supone rasters no girados

	 * @param pt	punto en coordenadas del punto real

	 * @return	punto en coordenadas del raster

	 */

	public Point2D worldToRasterWithoutRot(Point2D pt) {

		Point2D p = new Point2D.Double();

		AffineTransform at = new AffineTransform(	externalTransformation.getScaleX(), 0, 

													0, externalTransformation.getScaleY(), 

													externalTransformation.getTranslateX(), externalTransformation.getTranslateY());

		try {

			at.inverseTransform(pt, p);

		} catch (NoninvertibleTransformException e) {

			return pt;

		}

		return p;

	}

	/**

	 * Obtiene un punto 2D con las coordenadas del raster a partir de uno en coordenadas

	 * del punto real.

	 * Supone rasters no girados

	 * @param pt	punto en coordenadas del punto real

	 * @return	punto en coordenadas del raster

	 */

	public Point2D worldToRaster(Point2D pt) {

		Point2D p = new Point2D.Double();

		try {

			externalTransformation.inverseTransform(pt, p);

		} catch (NoninvertibleTransformException e) {

			return pt;

		}

		return p;

	}

	/**

	 * Obtiene un punto del raster en coordenadas pixel a partir de un punto en coordenadas

	 * reales. 

	 * @param pt Punto en coordenadas reales

	 * @return Punto en coordenadas pixel.

	 */

	public Point2D rasterToWorld(Point2D pt) {

		Point2D p = new Point2D.Double();

		externalTransformation.transform(pt, p);

		return p;

	}

	/**

	 * Calcula el overview a usar. Hay que tener en cuenta que tenemos que tener calculadas las variables

	 * viewPortScale, currentFullWidth y currentFulHeight

	 * @param coordenada pixel expresada en double que indica la posici�n superior izquierda

	 * @throws GdalException

	 */

	private void calcOverview(Point2D tl, Point2D br) throws GdalException {

		gdalBands[0] = getRasterBand(1);

		currentOverview = -1;

		if (gdalBands[0].getOverviewCount() > 0) {

			GdalRasterBand ovb = null;

			for (int i = gdalBands[0].getOverviewCount() - 1; i > 0; i--) {

				ovb = gdalBands[0].getOverview(i);

				if (ovb.getRasterBandXSize() > getRasterXSize() * viewportScaleX) {

					currentOverview = i;

					viewportScaleX *= ((double) width / (double) ovb.getRasterBandXSize());

					viewportScaleY *= ((double) height / (double) ovb.getRasterBandYSize());

					stepX = 1D / viewportScaleX;

					stepY = 1D / viewportScaleY;

					currentFullWidth = ovb.getRasterBandXSize();

					currentFullHeight = ovb.getRasterBandYSize();

					currentViewX = Math.min(tl.getX(), br.getX());

					lastReadLine = Math.min(tl.getY(), br.getY());

					break;

				}

			}

		}

	}

	public void setView(double dWorldTLX, double dWorldTLY,

						double dWorldBRX, double dWorldBRY,

						int nWidth, int nHeight) throws GdalException {

		currentFullWidth = width;

		currentFullHeight = height;

		Point2D tl = worldToRaster(new Point2D.Double(dWorldTLX, dWorldTLY));

		Point2D br = worldToRaster(new Point2D.Double(dWorldBRX, dWorldBRY));

		// Calcula cual es la primera l�nea a leer;

		currentViewWidth = nWidth;

		currentViewHeight = nHeight;

//		wcWidth = Math.abs(br.getX() - tl.getX());

		currentViewX = Math.min(tl.getX(), br.getX());

		viewportScaleX = (double) currentViewWidth / (br.getX() - tl.getX());

		viewportScaleY = (double) currentViewHeight / (br.getY() - tl.getY());

		stepX = 1D / viewportScaleX;

		stepY = 1D / viewportScaleY;

		lastReadLine = Math.min(tl.getY(), br.getY());

		//Para lectura del renderizado (ARGB). readWindow selecciona las bandas que necesita.

		// calcula el overview a usar

		gdalBands = new GdalRasterBand[4];

		calcOverview(tl, br);

		// Selecciona las bandas y los overviews necesarios

		/*gdalBands[0] = getRasterBand(rBandNr);

		gdalBands[1] = gdalBands[0]; 

		gdalBands[2] = gdalBands[1]; 

		if(getRasterCount() >= 2) {

			gdalBands[1] = getRasterBand(gBandNr);

			gdalBands[2] = gdalBands[1]; 

		}

		if(this.getRasterCount() >= 3) 

			gdalBands[2] = getRasterBand(bBandNr);

		if(colorInterpr.isAlphaBand())

			gdalBands[3] = getRasterBand(aBandNr);			

		assignDataTypeFromGdalRasterBands(gdalBands);

		if (currentOverview > 0) {

			gdalBands[0] = gdalBands[0].getOverview(currentOverview);

			if(getRasterCount() >= 2) {

				gdalBands[1] = gdalBands[1].getOverview(currentOverview);

			}

			if(this.getRasterCount() >= 3) 

				gdalBands[2] = gdalBands[2].getOverview(currentOverview);

			if(colorInterpr.isAlphaBand())

				gdalBands[3] = gdalBands[3].getOverview(currentOverview);			

		}*/

	}

	/**

	 * Selecciona bandas y overview en el objeto GdalRasterBand[] para el n�mero de bandas solicitado.

	 * @param nbands N�mero de bandas solicitado.

	 * @throws GdalException

	 */

	public void selectGdalBands(int nbands) throws GdalException {

		gdalBands = new GdalRasterBand[nbands];

		// Selecciona las bandas y los overviews necesarios

		gdalBands[0] = getRasterBand(1);

		for (int i = 0; i < nbands; i++)

			gdalBands[i] = gdalBands[0];

		assignDataTypeFromGdalRasterBands(gdalBands);

//		setDataType(gdalBands[0].getRasterDataType());

		for (int i = 2; i <= nbands; i++) {

			if (getRasterCount() >= i) {

				gdalBands[i - 1] = getRasterBand(i);

				for (int j = i; j < nbands; j++)

					gdalBands[j] = gdalBands[i - 1];

			}

		}

		if (currentOverview > 0) {

			gdalBands[0] = gdalBands[0].getOverview(currentOverview);

			for (int i = 2; i <= nbands; i++) {

				if (getRasterCount() >= i)

					gdalBands[i - 1] = gdalBands[i - 1].getOverview(currentOverview);

			}

		}

	}

	int lastY = -1;

	/**

	 * Lee una l�nea de bytes

	 * @param line Buffer donde se cargan los datos

	 * @param initOffset Desplazamiento inicial desde el margen inzquierdo. Esto es necesario para cuando

	 * se supersamplea ya que cada pixel de imagen ocupa muchos pixeles de pantalla y puede empezar a dibujarse

	 * por la izquierda a mitad de pixel

	 * @param gdalBuffer Buffer con la l�nea de datos original

	 */

	private void readLine(byte[][] line, double initOffset, GdalBuffer[] gdalBuffer) {

		double j = 0D;

		int i = 0;

		for (int iBand = 0; iBand < gdalBuffer.length; iBand++) {

			for (i = 0, j = initOffset; i < currentViewWidth && j < gdalBuffer[0].getSize(); i++, j += stepX) {

				line[iBand][i] = gdalBuffer[iBand].buffByte[(int) j];

			}

		}

	}

	/**

	 * Lee una l�nea de shorts

	 * @param line Buffer donde se cargan los datos

	 * @param initOffset Desplazamiento inicial desde el margen inzquierdo. Esto es necesario para cuando

	 * se supersamplea ya que cada pixel de imagen ocupa muchos pixeles de pantalla y puede empezar a dibujarse

	 * por la izquierda a mitad de pixel

	 * @param gdalBuffer Buffer con la l�nea de datos original

	 */

	private void readLine(short[][] line, double initOffset, GdalBuffer[] gdalBuffer) {

		double j = 0D;

		int i = 0;

		for (int iBand = 0; iBand < gdalBuffer.length; iBand++) {

			for (i = 0, j = initOffset; i < currentViewWidth && j < gdalBuffer[0].getSize(); i++, j += stepX) {

				line[iBand][i] = (short) (gdalBuffer[iBand].buffShort[(int) j] & 0xffff);

			}

		}

	}

	/**

	 * Lee una l�nea de ints

	 * @param line Buffer donde se cargan los datos

	 * @param initOffset Desplazamiento inicial desde el margen inzquierdo. Esto es necesario para cuando

	 * se supersamplea ya que cada pixel de imagen ocupa muchos pixeles de pantalla y puede empezar a dibujarse

	 * por la izquierda a mitad de pixel

	 * @param gdalBuffer Buffer con la l�nea de datos original

	 */

	private void readLine(int[][] line, double initOffset, GdalBuffer[] gdalBuffer) {

		double j = 0D;

		int i = 0;

		for (int iBand = 0; iBand < gdalBuffer.length; iBand++) {

			for (i = 0, j = initOffset; i < currentViewWidth && j < gdalBuffer[0].getSize(); i++, j += stepX) {

				line[iBand][i] = (gdalBuffer[iBand].buffInt[(int) j] & 0xffffffff);

			}

		}

	}

	/**

	 * Lee una l�nea de float

	 * @param line Buffer donde se cargan los datos

	 * @param initOffset Desplazamiento inicial desde el margen izquierdo. Esto es necesario para cuando

	 * se supersamplea ya que cada pixel de imagen ocupa muchos pixeles de pantalla y puede empezar a dibujarse

	 * por la izquierda a mitad de pixel

	 * @param gdalBuffer Buffer con la l�nea de datos original

	 */

	private void readLine(float[][] line, double initOffset, GdalBuffer[] gdalBuffer) {

		double j = 0D;

		int i = 0;

		for (int iBand = 0; iBand < gdalBuffer.length; iBand++) {

			for (i = 0, j = initOffset; i < currentViewWidth && j < gdalBuffer[0].getSize(); i++, j += stepX) {

				line[iBand][i] = gdalBuffer[iBand].buffFloat[(int) j];

			}

		}

	}

	/**

	 * Lee una l�nea de doubles

	 * @param line Buffer donde se cargan los datos

	 * @param initOffset Desplazamiento inicial desde el margen inzquierdo. Esto es necesario para cuando

	 * se supersamplea ya que cada pixel de imagen ocupa muchos pixeles de pantalla y puede empezar a dibujarse

	 * por la izquierda a mitad de pixel

	 * @param gdalBuffer Buffer con la l�nea de datos original

	 */

	private void readLine(double[][] line, double initOffset, GdalBuffer[] gdalBuffer) {

		double j = 0D;

		int i = 0;

		for (int iBand = 0; iBand < gdalBuffer.length; iBand++) {

			for (i = 0, j = initOffset; i < currentViewWidth && j < gdalBuffer[0].getSize(); i++, j += stepX) {

				line[iBand][i] = gdalBuffer[iBand].buffDouble[(int) j];

			}

		}

	}

	/**

	 * Lee una l�nea completa del raster y devuelve un array del tipo correcto. Esta funci�n es util

	 * para una lectura rapida de todo el fichero sin necesidad de asignar vista. 

	 * @param nLine N�mero de l�nea a leer

	 * @param band Banda requerida

	 * @return Object que es un array unidimendional del tipo de datos del raster

	 * @throws GdalException

	 */

	public Object readCompleteLine(int nLine, int band) throws GdalException {

		GdalRasterBand gdalBand = super.getRasterBand(band + 1);

		GdalBuffer gdalBuf = null;

		gdalBuf = gdalBand.readRaster(0, nLine, getRasterXSize(), 1, getRasterXSize(), 1, dataType[band]);

		if (dataType[band] == GDT_Byte)

			return gdalBuf.buffByte;

		if (dataType[band] == GDT_Int16 || dataType[band] == GDT_UInt16)

			return gdalBuf.buffShort;

		if (dataType[band] == GDT_Int32 || dataType[band] == GDT_UInt32)

			return gdalBuf.buffInt;

		if (dataType[band] == GDT_Float32)

			return gdalBuf.buffFloat;

		if (dataType[band] == GDT_Float64)

			return gdalBuf.buffDouble;

		if (dataType[band] == GDT_CInt16 || dataType[band] == GDT_CInt32 ||

				dataType[band] == GDT_CFloat32 || dataType[band] == GDT_CFloat64)

			return null;

		return null;

	}

	/**

	 * Lee una bloque completo del raster y devuelve un array tridimensional del tipo correcto. Esta funci�n es util

	 * para una lectura rapida de todo el fichero sin necesidad de asignar vista. 

	 * @param nLine N�mero de l�nea a leer

	 * @param band Banda requerida

	 * @return Object que es un array unidimendional del tipo de datos del raster

	 * @throws GdalException

	 */

	public Object readBlock(int pos, int blockHeight) throws GdalException, ProcessInterruptedException {

		bBandNr = super.getRasterCount();

		int nX = getRasterXSize();

		RasterTask task = RasterTaskQueue.get(Thread.currentThread().toString());

		GdalRasterBand[] gdalBand = new GdalRasterBand[bBandNr];

		for (int iBand = 0; iBand < gdalBand.length; iBand++) 

			gdalBand[iBand] = super.getRasterBand(iBand + 1);

		GdalBuffer[] gdalBuf = new GdalBuffer[bBandNr];

		if (dataType[0] == GDT_Byte) {

			byte[][][] buf = new byte[bBandNr][blockHeight][getRasterXSize()];

			for (int iBand = 0; iBand < gdalBuf.length; iBand++) {

				gdalBuf[iBand] = gdalBand[iBand].readRaster(0, pos, nX, blockHeight, nX, blockHeight, dataType[0]);

				for (int iRow = 0; iRow < blockHeight; iRow++) {

					for (int iCol = 0; iCol < nX; iCol++) 

						buf[iBand][iRow][iCol] = gdalBuf[iBand].buffByte[iRow * nX + iCol];

					if(task.getEvent() != null)

						task.manageEvent(task.getEvent());

				}

			}	

			return buf;

		} else if (dataType[0] == GDT_CInt16 || dataType[0] == GDT_Int16  || dataType[0] == GDT_UInt16) {

			short[][][] buf = new short[bBandNr][blockHeight][getRasterXSize()];

			for (int iBand = 0; iBand < gdalBuf.length; iBand++) {

				gdalBuf[iBand] = gdalBand[iBand].readRaster(0, pos, nX, blockHeight, nX, blockHeight, dataType[0]);

				for (int iRow = 0; iRow < blockHeight; iRow++) {

					for (int iCol = 0; iCol < nX; iCol++) 

						buf[iBand][iRow][iCol] = gdalBuf[iBand].buffShort[iRow * nX + iCol];

					if(task.getEvent() != null)

						task.manageEvent(task.getEvent());

				}

			}	

			return buf;

		} else if (dataType[0] == GDT_CInt32 || dataType[0] == GDT_Int32  || dataType[0] == GDT_UInt32) {

			int[][][] buf = new int[bBandNr][blockHeight][getRasterXSize()];

			for (int iBand = 0; iBand < gdalBuf.length; iBand++) {

				gdalBuf[iBand] = gdalBand[iBand].readRaster(0, pos, nX, blockHeight, nX, blockHeight, dataType[0]);

				for (int iRow = 0; iRow < blockHeight; iRow++) { 

					for (int iCol = 0; iCol < nX; iCol++)

						buf[iBand][iRow][iCol] = gdalBuf[iBand].buffInt[iRow * nX + iCol];

					if(task.getEvent() != null)

						task.manageEvent(task.getEvent());

				}

			}	

			return buf;

		} else if(dataType[0] == GDT_Float32 || dataType[0] == GDT_CFloat32) {

			float[][][] buf = new float[bBandNr][blockHeight][getRasterXSize()];

			for (int iBand = 0; iBand < gdalBuf.length; iBand++) {

				gdalBuf[iBand] = gdalBand[iBand].readRaster(0, pos, nX, blockHeight, nX, blockHeight, dataType[0]);

				for (int iRow = 0; iRow < blockHeight; iRow++) {

					for (int iCol = 0; iCol < nX; iCol++)

						buf[iBand][iRow][iCol] = gdalBuf[iBand].buffFloat[iRow * nX + iCol];

					if(task.getEvent() != null)

						task.manageEvent(task.getEvent());

				}

			}	

			return buf;

		} else if(dataType[0] == GDT_Float64 || dataType[0] == GDT_CFloat64) {

			double[][][] buf = new double[bBandNr][blockHeight][getRasterXSize()];

			for (int iBand = 0; iBand < gdalBuf.length; iBand++) {

				gdalBuf[iBand] = gdalBand[iBand].readRaster(0, pos, nX, blockHeight, nX, blockHeight, dataType[0]);

				for (int iRow = 0; iRow < blockHeight; iRow++) {

					for (int iCol = 0; iCol < nX; iCol++) 

						buf[iBand][iRow][iCol] = gdalBuf[iBand].buffDouble[iRow * nX + iCol];

					if(task.getEvent() != null)

						task.manageEvent(task.getEvent());

				}

			}		

			return buf;

		}

			return null;

	}

	/**

	 * Lectura de una l�nea de datos.

	 * @param line

	 * @throws GdalException

	 */

	public void readLine(Object line) throws GdalException {

		int w = (int) (Math.ceil(((double)currentViewWidth)*stepX) + 1);

		int x = (int) (currentViewX);

		int y = (int) (lastReadLine);

		GdalBuffer r = null, g = null, b = null;

		GdalBuffer a = new GdalBuffer();

		while(y >= gdalBands[0].getRasterBandYSize())

			y--;

		if (x+w > gdalBands[0].getRasterBandXSize()) 

			w = gdalBands[0].getRasterBandXSize()-x;

		if(gdalBands[0].getRasterColorTable() != null) {

			palette = new DataStoreColorTable();

			palette.createPaletteFromGdalColorTable(gdalBands[0].getRasterColorTable());

			r = gdalBands[0].readRaster(x, y, w, 1, w, 1, dataType[0]);

		} else {

			a.buffByte = new byte[w];

			r = gdalBands[0].readRaster(x, y, w, 1, w, 1, dataType[0]);

			g = b = r;

			if (getRasterCount() > 1 && gdalBands[1] != null)

				g = gdalBands[1].readRaster(x, y, w, 1, w, 1, dataType[0]);

			if (getRasterCount() > 2 && gdalBands[2] != null)

				b = gdalBands[2].readRaster(x, y, w, 1, w, 1, dataType[0]);

		}

		lastReadLine += stepY;

		double initOffset =  Math.abs(currentViewX - ((int)currentViewX));

		GdalBuffer[] bands = {r, g, b};

		if (dataType[0] == GDT_Byte)

			readLine((byte[][])line, initOffset, bands);

		else if (dataType[0] == GDT_CInt16 || dataType[0] == GDT_Int16  || dataType[0] == GDT_UInt16)

			readLine((short[][])line, initOffset, bands);

		else if (dataType[0] == GDT_CInt32 || dataType[0] == GDT_Int32  || dataType[0] == GDT_UInt32)

			readLine((int[][])line, initOffset, bands);

		else if(dataType[0] == GDT_Float32 || dataType[0] == GDT_CFloat32)

			readLine((float[][])line, initOffset, bands);

		else if(dataType[0] == GDT_Float64 || dataType[0] == GDT_CFloat64)

			readLine((double[][])line, initOffset, bands);

		return;

	}

	/**

	 * Cuando se hace una petici�n de carga de buffer la extensi�n pedida puede

	 * estar ajustada a la extensi�n del raster o no estarlo. En caso de no

	 * estarlo los pixeles del buffer que caen fuera de la extensi�n del raster

	 * tendr�n valor de NoData. Esta funci�n calcula en que pixel del buffer hay

	 * que empezar a escribir en caso de que este sea mayor que los datos a leer.

	 * 

	 * @param dWorldTLX Posici�n X superior izquierda en coord reales

	 * @param dWorldTLY Posici�n Y superior izquierda en coord reales

	 * @param dWorldBRX Posici�n X inferior derecha en coord reales

	 * @param dWorldBRY Posici�n Y inferior derecha en coord reales

	 * @param nWidth Ancho en pixeles del buffer

	 * @param nHeight Alto en pixeles del buffer

	 * @return desplazamiento dentro del buffer en X e Y

	 */ 

	private int[] calcStepBuffer(Extent dataExtent, int nWidth, int nHeight, int[] stpBuffer) {

		Extent imageExtent = getExtentWithoutRot();

		Extent ajustDataExtent = RasterLocator.getManager().getRasterUtils().calculateAdjustedView(dataExtent, imageExtent);

		if(!RasterLocator.getManager().getRasterUtils().compareExtents(dataExtent, ajustDataExtent)){

			Point2D p1 = worldToRasterWithoutRot(new Point2D.Double(ajustDataExtent.minX(), ajustDataExtent.maxY()));

			Point2D p2 = worldToRasterWithoutRot(new Point2D.Double(ajustDataExtent.maxX(), ajustDataExtent.minY()));

			Point2D p3 = worldToRasterWithoutRot(new Point2D.Double(dataExtent.minX(), dataExtent.maxY()));

			//    		Point2D p4 = worldToRasterWithoutRot(new Point2D.Double(dataExtent.maxX(), dataExtent.minY()));

			//Ese es el ancho y alto q tendr�a el buffer en caso de haberse ajustado

			int w = (int)Math.abs(Math.ceil(p2.getX()) - Math.floor(p1.getX())); 

			int h = (int)Math.abs(Math.floor(p1.getY()) - Math.ceil(p2.getY()));

			stpBuffer[0] = (int)(p1.getX() + (-p3.getX()));

			stpBuffer[1] = (int)(p1.getY() + (-p3.getY()));

			stpBuffer[2] = stpBuffer[0] + w; 

			stpBuffer[3] = stpBuffer[1] + h;

			return new int[]{w, h};

		}

		return new int[]{nWidth, nHeight};

	}

	/**

	 * Lee una ventana de datos sin resampleo a partir de coordenadas reales.

	 * @param buf Buffer donde se almacenan los datos

	 * @param bandList Lista de bandas que queremos leer y sobre que bandas del buffer de destino queremos escribirlas

	 * @param dWorldTLX Posici�n X superior izquierda en coord reales

	 * @param dWorldTLY Posici�n Y superior izquierda en coord reales

	 * @param dWorldBRX Posici�n X inferior derecha en coord reales

	 * @param dWorldBRY Posici�n Y inferior derecha en coord reales

	 * @param nWidth Ancho en pixeles del buffer

	 * @param nHeight Alto en pixeles del buffer

	 * @throws GdalException

	 */

	public void readWindow(Buffer buf, BandList bandList, double ulx, double uly,double lrx, double lry,

			int nWidth, int nHeight, boolean adjustToExtent) throws GdalException, ProcessInterruptedException {

		Extent petExtent = new ExtentImpl(ulx, uly, lrx, lry);

		setView(ulx, uly, lrx, lry, nWidth, nHeight);

		Point2D tl = worldToRaster(new Point2D.Double(ulx, uly));

		Point2D br = worldToRaster(new Point2D.Double(lrx, lry));

		if(tl.getX() > br.getX())

			tl.setLocation(tl.getX() - 1, tl.getY());

		else

			br.setLocation(br.getX() - 1, br.getY());

		if(tl.getY() > br.getY())

			tl.setLocation(tl.getX(), tl.getY() - 1);

		else

			br.setLocation(br.getX(), br.getY() - 1);

		if(gdalBands.length == 0)

			return;

		selectGdalBands(/*buf.getBandCount()*/getRasterCount());

		int x = (int) Math.round(Math.min(tl.getX(), br.getX()));

		int y = (int) Math.round(Math.min(tl.getY(), br.getY()));

		int[] stpBuffer = new int[]{0, 0 , buf.getWidth(), buf.getHeight()};

		//Si el buffer no se ajusta al extent entonces calculamos en que posici�n comienza a escribirse dentro del buffer

		//ya que lo que cae fuera ser�n valores NoData

		if(!adjustToExtent){

			int[] wh = calcStepBuffer(petExtent, nWidth, nHeight, stpBuffer);

			if(x < 0)

				x  = 0;

			if(y < 0)

				y  = 0;

			readData(buf, bandList, x, y, wh[0], wh[1], wh[0], wh[1], 0, 0, stpBuffer);

			return;

		}

		readData(buf, bandList, x, y, nWidth, nHeight, nWidth, nHeight, 0, 0, stpBuffer);

	}

	/**

	 * Lee una ventana de datos con resampleo a partir de coordenadas reales. Este m�todo lee la

	 * ventana de una vez cargando los datos de un golpe en el buffer. Las coordenadas se solicitan

	 * en coordenadas del mundo real por lo que estas pueden caer en cualquier parte de un pixel.

	 * Esto se hace m�s evidente cuando supersampleamos en la petici�n, es decir el buffer de de 

	 * mayor tama�o que el n�mero de pixels solicitado.

	 * 

	 * Para resolver esto escribiremos con la funci�n readRaster los datos sobre un buffer mayor 

	 * que el solicitado. Despu�s calcularemos el desplazamiento en pixels dentro de este buffer 

	 * de mayor tama�o hasta llegar a la coordenada real donde comienza la petici�n real que ha

	 * hecho el usuario. Esto es as� porque cuando supersampleamos no queremos los pixeles del 

	 * raster de disco completos sino que en los bordes del buffer quedan cortados.  

	 *  

	 * @param buf Buffer donde se almacenan los datos

	 * @param bandList Lista de bandas que queremos leer y sobre que bandas del buffer de destino queremos escribirlas

	 * @param dWorldTLX Posici�n X superior izquierda en coord reales

	 * @param dWorldTLY Posici�n Y superior izquierda en coord reales

	 * @param dWorldBRX Posici�n X inferior derecha en coord reales

	 * @param dWorldBRY Posici�n Y inferior derecha en coord reales

	 * @param nWidth Ancho en pixeles de la petici�n

	 * @param nHeight Alto en pixeles de la petici�n

	 * @param bufWidth Ancho del buffer

	 * @param bufHeight Alto del buffer

	 * @throws GdalException

	 */

	public void readWindow(Buffer buf, BandList bandList, double ulx, double uly, double lrx, double lry,

										double nWidth, double nHeight, int bufWidth, int bufHeight, boolean adjustToExtent) throws GdalException, ProcessInterruptedException {

		Extent petExtent = new ExtentImpl(ulx, uly, lrx, lry);

		setView(ulx, uly, lrx, lry, bufWidth, bufHeight);

		Point2D ul = worldToRaster(new Point2D.Double(ulx, uly));

		Point2D lr = worldToRaster(new Point2D.Double(lrx, lry));

		ul.setLocation(ul.getX() < 0 ? 1 : ul.getX(), ul.getY() < 0 ? 1 : ul.getY());

		lr.setLocation(lr.getX() < 0 ? 1 : lr.getX(), lr.getY() < 0 ? 1 : lr.getY());

		ul.setLocation(ul.getX() - 0.5, ul.getY() - 0.5);

		lr.setLocation(lr.getX() - 0.5, lr.getY() - 0.5);

		/*if(tl.getX() > br.getX())

			tl.setLocation(tl.getX() - 1, tl.getY());

		else

			br.setLocation(br.getX() - 1, br.getY());

		if(tl.getY() > br.getY())

			tl.setLocation(tl.getX(), tl.getY() - 1);

		else

			br.setLocation(br.getX(), br.getY() - 1);*/

		adjustPoints(ul, lr);

		if(gdalBands.length == 0)

			return;

		selectGdalBands(/*buf.getBandCount()*/getRasterCount());

		int x = (int) Math.min(ul.getX(), lr.getX());

		int y = (int) Math.min(ul.getY(), lr.getY());

		//int endX = (int) Math.ceil(Math.max(br.getX(), tl.getX()));

		//int endY = (int) Math.ceil(Math.max(br.getY(), tl.getY()));

		int stpX = 0;

		int stpY = 0;

		/*if(bufWidth > Math.ceil(nWidth)){

			stpX = (int)(((tl.getX() - x) * bufWidth) / nWidth);

			bufWidth = (int)((Math.abs(endX - x) * bufWidth) / nWidth);

		}

		if(bufHeight > Math.ceil(nHeight)){

			stpY = (int)(((tl.getY() - y) * bufHeight) / nHeight);

			bufHeight = (int)((Math.abs(endY - y) * bufHeight) / nHeight);

		}

		nWidth = (int)Math.abs(endX - x);

		nHeight = (int)Math.abs(endY - y);*/

		nWidth = (nWidth * currentFullWidth) / width;

		nHeight = (nHeight * currentFullHeight) / height;

		x = (int)(((long)x * (long)currentFullWidth) / (long)width);

		y = (int) (((long)y * (long)currentFullHeight) / (long)height);

		int[] stpBuffer = new int[]{0, 0 , buf.getWidth(), buf.getHeight()};

		//Si el buffer no se ajusta al extent entonces calculamos en que posici�n comienza a escribirse dentro del buffer

		//ya que lo que cae fuera ser�n valores NoData

		if(!adjustToExtent){

			int[] wh = calcStepBuffer(petExtent, bufWidth, bufHeight, stpBuffer);

			if(x < 0)

				x  = 0;

			if(y < 0)

				y  = 0;

			stpBuffer[0] = (int)((stpBuffer[0] * bufWidth) / nWidth);

			stpBuffer[1] = (int)((stpBuffer[1] * bufHeight) / nHeight);

			stpBuffer[2] = (int)((stpBuffer[2] * bufWidth) / nWidth);

			stpBuffer[3] = (int)((stpBuffer[3] * bufHeight) / nHeight);

			bufWidth = (int)Math.abs(stpBuffer[2] - stpBuffer[0]);

			bufHeight = (int)Math.abs(stpBuffer[3] - stpBuffer[1]);

			readData(buf, bandList, x, y, wh[0], wh[1], bufWidth, bufHeight, 0, 0, stpBuffer);

			return;

		}

		if ((x + nWidth) > gdalBands[0].getRasterBandXSize()) 

			nWidth = gdalBands[0].getRasterBandXSize() - x;

		if ((y + nHeight) > gdalBands[0].getRasterBandYSize()) 

			nHeight = gdalBands[0].getRasterBandYSize() - y;

		readData(buf, bandList, x, y, (int)nWidth, (int)nHeight, bufWidth, bufHeight, stpX, stpY, stpBuffer);

	}

	private void adjustPoints(Point2D ul, Point2D lr) {

		double a = (ul.getX() - (int)ul.getX());

		double b = (ul.getY() - (int)ul.getY());

		ul.setLocation(	(a > 0.95 || a < 0.05) ? Math.round(ul.getX()) : ul.getX(), 

						(b > 0.95 || b < 0.05) ? Math.round(ul.getY()) : ul.getY());

		lr.setLocation(	(a > 0.95 || a < 0.05) ? Math.round(lr.getX()) : lr.getX(), 

						(b > 0.95 || b < 0.05) ? Math.round(lr.getY()) : lr.getY());

	}

	/**

	 * Lee una ventana de datos sin resampleo a partir de coordenadas en pixeles.

	 * @param buf Buffer donde se almacenan los datos

	 * @param bandList Lista de bandas que queremos leer y sobre que bandas del buffer de destino queremos escribirlas

	 * @param x Posici�n X en pixeles

	 * @param y Posici�n Y en pixeles

	 * @param w Ancho en pixeles

	 * @param h Alto en pixeles

	 * @throws GdalException

	 */

	public void readWindow(Buffer buf, BandList bandList, int x, int y, int w, int h) 

		throws GdalException, ProcessInterruptedException {

		gdalBands = new GdalRasterBand[getRasterCount()];

		isSupersampling = false;

		if(gdalBands.length == 0)

			return;

		// Selecciona las bandas

		gdalBands[0] = getRasterBand(1);

		for(int iBand = 1; iBand < gdalBands.length; iBand++)

			gdalBands[iBand] = getRasterBand(iBand + 1);

		assignDataTypeFromGdalRasterBands(gdalBands);

		int[] stepBuffer = new int[]{0, 0, w, h};

		readData(buf, bandList, x, y, x + w, y + h, w, h, 0, 0, stepBuffer);

		//int yMax = y + h;

		//readDataByLine(buf, bandList, x, y, w, yMax);

	}

	/**

	 * Lee una ventana de datos con resampleo a partir de coordenadas en pixeles. Este m�todo lee la

	 * ventana de una vez cargando los datos de un golpe en el buffer.

	 * @param buf Buffer donde se almacenan los datos

	 * @param bandList Lista de bandas que queremos leer y sobre que bandas del buffer de destino queremos escribirlas

	 * @param x Posici�n X en pixeles

	 * @param y Posici�n Y en pixeles

	 * @param w Ancho en pixeles

	 * @param h Alto en pixeles

	 * @param bufWidth Ancho del buffer

	 * @param bufHeight Alto del buffer

	 * @throws GdalException

	 */

	public void readWindow(Buffer buf, BandList bandList, int x, int y, int w, int h, int bufWidth, int bufHeight) throws GdalException, ProcessInterruptedException {

		gdalBands = new GdalRasterBand[getRasterCount()];

		if(gdalBands.length == 0)

			return;

		// Selecciona las bandas

		gdalBands[0] = getRasterBand(1);

		for(int iBand = 1; iBand < gdalBands.length; iBand++)

			gdalBands[iBand] = getRasterBand(iBand + 1);

		assignDataTypeFromGdalRasterBands(gdalBands);

		int[] stpBuffer = new int[]{0, 0 , buf.getWidth(), buf.getHeight()};

		readData(buf, bandList, x, y, w, h, bufWidth, bufHeight, 0, 0, stpBuffer);

	}

	/**

	 * Asigna el tipo de datos de las bandas a partir de una lista de GdalRasterBands

	 * @param gdalBands

	 * @throws GdalException

	 */

	private void assignDataTypeFromGdalRasterBands(GdalRasterBand[] gdalBands) throws GdalException {

		int[] dt = new int[gdalBands.length];

		for (int i = 0; i < gdalBands.length; i++) {

			if(gdalBands[i] != null)

				dt[i] = gdalBands[i].getRasterDataType();

		}

		setDataType(dt);

	}

	/**

	 * Lee una ventana de datos sin resampleo a partir de coordenadas en pixeles. Esta funci�n es usuada por

	 * readWindow para coordenadas reales y readWindow en coordenadas pixel. 

	 * @param buf Buffer donde se almacenan los datos

	 * @param bandList Lista de bandas que queremos leer y sobre que bandas del buffer de destino queremos escribirlas

	 * @param x Posici�n X en pixeles

	 * @param y Posici�n Y en pixeles

	 * @param w Ancho en pixeles

	 * @param h Alto en pixeles

	 * @param bufWidth Ancho del buffer

	 * @param bufHeight Alto del buffer

	 * @param stepX Desplazamiento en pixeles en X a partir de la posici�n x. Este desplazamiento es util cuando hay un 

	 * supersampleo ya que puede ser que de los pixeles que est�n en el borde izquierdo de la petici�n solo queramos una

	 * parte de ellos. 

	 * @param stepY Desplazamiento en pixeles en Y a partir de la posici�n y. Este desplazamiento es util cuando hay un 

	 * supersampleo ya que puede ser que de los pixeles que est�n en el borde superior de la petici�n solo queramos una

	 * parte de ellos.

	 * @param stepBuffer El buffer puede empezar a escribirse a partir de un pixel determinado y acabar de escribir antes 

	 * de fin de buffer. Este par�metro indica el desplazamiento desde el inicio del buffer y la posici�n final.

	 * <UL>

	 * <LI>stepBuffer[0]:Desplazamiento en X desde el inicio</LI>

	 * <LI>stepBuffer[1]:Desplazamiento en Y desde el inicio</LI>

	 * <LI>stepBuffer[2]:Posici�n X final</LI>

	 * <LI>stepBuffer[3]:Posici�n Y final</LI>

	 * </UL>

	 * @throws GdalException

	 */

	private void readData(Buffer buf, BandList bandList, int x, int y, int w, int h, 

			int bufWidth, int bufHeight, int stpX, int stpY, int[] stepBuffer) throws GdalException, ProcessInterruptedException {

		RasterTask task = RasterTaskQueue.get(Thread.currentThread().toString()); 

		GdalBuffer gdalBuf = null;

		for(int iBand = 0; iBand < gdalBands.length; iBand++) {

			int[] drawableBands = bandList.getBufferBandToDraw(fileName, iBand);

			if(drawableBands == null || (drawableBands.length == 1 && drawableBands[0] == -1))

				continue;	

			int init = (int)((bufWidth * stpY) + stpX); //Pos inicial. Desplazamos stpX pixels hacia la derecha y bajamos stpY lineas

			int pos = init;

			gdalBuf = gdalBands[iBand].readRaster(x, y, w, h, bufWidth, bufHeight, dataType[iBand]);

			if(dataType[iBand] == Gdal.GDT_Byte){

				for (int line = stepBuffer[1]; line < stepBuffer[3]/*buf.getHeight()*/; line++) {

					pos = (int)((bufWidth * (line - stepBuffer[0])) + init);

					for (int col = stepBuffer[0]; col < stepBuffer[2]/*buf.getWidth()*/; col ++) {

						for (int i = 0; i < drawableBands.length; i++) 

							buf.setElem(line, col, drawableBands[i], gdalBuf.buffByte[pos]);

						pos ++;

					}

					if(task.getEvent() != null)

						task.manageEvent(task.getEvent());

				}

			}else if((dataType[iBand] == Gdal.GDT_UInt16) || (dataType[iBand] == Gdal.GDT_Int16) || (dataType[iBand] == Gdal.GDT_CInt16)){

				for (int line = stepBuffer[1]; line < stepBuffer[3]; line++) {

					pos = (int)((bufWidth * (line - stepBuffer[0])) + init);

					for (int col = stepBuffer[0]; col < stepBuffer[2]; col ++) {

						for (int i = 0; i < drawableBands.length; i++)

							buf.setElem(line, col, drawableBands[i], gdalBuf.buffShort[pos]);

						pos ++;

					}

					if(task.getEvent() != null)

						task.manageEvent(task.getEvent());

				}

			}else if((dataType[iBand] == Gdal.GDT_UInt32) || (dataType[iBand] == Gdal.GDT_Int32) || (dataType[iBand] == Gdal.GDT_CInt32)){

				for (int line = stepBuffer[1]; line < stepBuffer[3]; line++) {

					pos = (int)((bufWidth * (line - stepBuffer[0])) + init);

					for (int col = stepBuffer[0]; col < stepBuffer[2]; col ++) {

						for (int i = 0; i < drawableBands.length; i++)

							buf.setElem(line, col, drawableBands[i], gdalBuf.buffInt[pos]);

						pos ++;

					}

					if(task.getEvent() != null)

						task.manageEvent(task.getEvent());

				}

			}else if(dataType[iBand] == Gdal.GDT_Float32){

				for (int line = stepBuffer[1]; line < stepBuffer[3]; line++) {

					pos = (int)((bufWidth * (line - stepBuffer[0])) + init);

					for (int col = stepBuffer[0]; col < stepBuffer[2]; col ++) {

						for (int i = 0; i < drawableBands.length; i++)

							buf.setElem(line, col, drawableBands[i], gdalBuf.buffFloat[pos]);

						pos ++;

					}

					if(task.getEvent() != null)

						task.manageEvent(task.getEvent());

				}

			}else if(dataType[iBand] == Gdal.GDT_Float64){

				for (int line = stepBuffer[1]; line < stepBuffer[3]; line++) {

					pos = (int)((bufWidth * (line - stepBuffer[0])) + init);

					for (int col = stepBuffer[0]; col < stepBuffer[2]; col ++) {

						for (int i = 0; i < drawableBands.length; i++)

							buf.setElem(line, col, drawableBands[i], gdalBuf.buffDouble[pos]);

						pos ++;

					}

					if(task.getEvent() != null)

						task.manageEvent(task.getEvent());

				}

			}

		}

	}

	/**

	 * Lee una ventana de datos sin resampleo a partir de coordenadas en pixeles. Esta funci�n es usuada por

	 * readWindow para coordenadas reales y readWindow en coordenadas pixel. 

	 * @param buf Buffer donde se almacenan los datos

	 * @param bandList Lista de bandas que queremos leer y sobre que bandas del buffer de destino queremos escribirlas

	 * @param x Posici�n X en pixeles

	 * @param y Posici�n Y en pixeles

	 * @param w Ancho en pixeles

	 * @param yMax altura m�xima de y

	 * @throws GdalException

	 */

	@SuppressWarnings("unused")

	private void readDataByLine(Buffer buf, BandList bandList, int x, int y, int w, int yMax) throws GdalException, ProcessInterruptedException {

		GdalBuffer gdalBuf = null;

		int rasterBufLine;

		RasterTask task = RasterTaskQueue.get(Thread.currentThread().toString());

		for(int iBand = 0; iBand < gdalBands.length; iBand++) {

			int[] drawableBands = bandList.getBufferBandToDraw(fileName, iBand);

			if(drawableBands == null || (drawableBands.length == 1 && drawableBands[0] == -1))

				continue;	

			if(dataType[iBand] == Gdal.GDT_Byte) {

				for (int line = y; line < yMax; line++) {

					gdalBuf = gdalBands[iBand].readRaster(x, line, w, 1, w, 1, dataType[iBand]);

					rasterBufLine = line - y;

					for (int i = 0; i < drawableBands.length; i++) {

						buf.setLineInBandByte(gdalBuf.buffByte, rasterBufLine, drawableBands[i]);

					}

					if(task.getEvent() != null)

						task.manageEvent(task.getEvent());

				}

			}else if((dataType[iBand] == Gdal.GDT_UInt16) || (dataType[iBand] == Gdal.GDT_Int16) || (dataType[iBand] == Gdal.GDT_CInt16)) {

				for (int line = y; line < yMax; line++) {

					gdalBuf = gdalBands[iBand].readRaster(x, line, w, 1, w, 1, dataType[iBand]);

					rasterBufLine = line - y;

					for (int i = 0; i < drawableBands.length; i++) {

						buf.setLineInBandShort(gdalBuf.buffShort, rasterBufLine, drawableBands[i]);

					}

					if(task.getEvent() != null)

						task.manageEvent(task.getEvent());

				}

			}else if((dataType[iBand] == Gdal.GDT_UInt32) || (dataType[iBand] == Gdal.GDT_Int32) || (dataType[iBand] == Gdal.GDT_CInt32)) {

				for (int line = y; line < yMax; line++) {

					gdalBuf = gdalBands[iBand].readRaster(x, line, w, 1, w, 1, dataType[iBand]);

					rasterBufLine = line - y;

					for (int i = 0; i < drawableBands.length; i++) {

						buf.setLineInBandInt(gdalBuf.buffInt, rasterBufLine, drawableBands[i]);

					}

					if(task.getEvent() != null)

						task.manageEvent(task.getEvent());

				}

			}else if(dataType[iBand] == Gdal.GDT_Float32){

				for (int line = y; line < yMax; line++) {

					gdalBuf = gdalBands[iBand].readRaster(x, line, w, 1, w, 1, dataType[iBand]);

					rasterBufLine = line - y;

					for (int i = 0; i < drawableBands.length; i++) {

						buf.setLineInBandFloat(gdalBuf.buffFloat, rasterBufLine, drawableBands[i]);

					}

					if(task.getEvent() != null)

						task.manageEvent(task.getEvent());

				}

			}else if(dataType[iBand] == Gdal.GDT_Float64){

				for (int line = y; line < yMax; line++) {

					gdalBuf = gdalBands[iBand].readRaster(x, line, w, 1, w, 1, dataType[iBand]);

					rasterBufLine = line - y;

					for (int i = 0; i < drawableBands.length; i++) {

						buf.setLineInBandDouble(gdalBuf.buffDouble, rasterBufLine, drawableBands[i]);

					}

					if(task.getEvent() != null)

						task.manageEvent(task.getEvent());

				}

			}

		}

	}

	/**

	 * Obtiene el valor de un pixel determinado por las coordenadas x e y que se pasan

	 * por par�metro

	 * @param x Coordenada X del pixel

	 * @param y Coordenada Y del pixel

	 * @return Array de Object donde cada posici�n representa una banda y el valor ser� Integer

	 * en caso de ser byte, shot o int, Float en caso de ser float y Double en caso de ser double.

	 */

	public Object[] getData(int x, int y) {

		try {

			Object[] data = new Object[getRasterCount()];

			for(int i = 0; i < getRasterCount(); i++){

				GdalRasterBand rb = getRasterBand(i + 1);

				GdalBuffer r = rb.readRaster(x, y, 1, 1, 1, 1, dataType[i]);

				switch(dataType[i]){

				case 0:	break;									//Sin tipo

				case 1:	data[i] = new Integer(r.buffByte[0]); 	//Buffer byte (8)

						break;

				case 2:											//Buffer short (16)

				case 3:	data[i] = new Integer(r.buffShort[0]);	//Buffer short (16)

						break;

				case 4:											//Buffer int (32)

				case 5: data[i] = new Integer(r.buffInt[0]);	//Buffer int (32)

						break;

				case 6:	data[i] = new Float(r.buffFloat[0]);	//Buffer float (32)

						break;

				case 7:	data[i] = new Double(r.buffDouble[0]);	//Buffer double (64)

						break;

				}

			}

			return data;

		} catch (GdalException e) {

			return null;

		}

	}

	public int getBlockSize(){

		return this.getBlockSize();

	}

	/**

	 * Devuelve la transformaci�n del fichero de georreferenciaci�n

	 * @return AffineTransform

	 */

	public AffineTransform getOwnTransformation() {

		return ownTransformation;

	}

	/**

	 * Calcula el extent en coordenadas del mundo real sin rotaci�n. Solo coordenadas y tama�o de pixel

	 * @return Extent

	 */

	public Extent getExtentWithoutRot() {

		AffineTransform at = new AffineTransform(	externalTransformation.getScaleX(), 0, 

													0, externalTransformation.getScaleY(), 

													externalTransformation.getTranslateX(), externalTransformation.getTranslateY());

		Point2D p1 = new Point2D.Double(0, 0);

		Point2D p2 = new Point2D.Double(width, height);

		at.transform(p1, p1);

		at.transform(p2, p2);

		return new ExtentImpl(p1, p2);

	}

	/**

	 * Asigna una transformaci�n que es aplicada sobre la que ya tiene el propio fichero

	 * @param t

	 */

	public void setExternalTransform(AffineTransform t){

		externalTransformation = t;

	}

	/**

	 * Obtiene el nombre del driver de Gdal

	 * @return Cadena que representa el nombre del driver de gdal

	 */

	public String getGdalShortName() {

		return shortName;

	}

}

/* gvSIG. Geographic Information System of the Valencian Government

 *

 * Copyright (C) 2007-2008 Infrastructures and Transports Department

 * of the Valencian Government (CIT)

 *

 * 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., 51 Franklin Street, Fifth Floor, Boston,

 * MA  02110-1301, USA.

 *

 */

package org.gvsig.raster.gdal.io.features;

import org.gvsig.raster.gdal.io.GdalProvider;

import org.gvsig.raster.gdal.io.GdalWriter;

import org.gvsig.raster.impl.store.WriteFileFormatFeatures;

/**

 * Caracteristicas del formato ILWIS para escritura.

 * Soporta tipos de datos enteros en 8, 16 y 32 bits y coma flotante en 32 y 64 bits

 * Genera un fichero .mpl con la informaci?n de contenido (ficheros), un fichero .grf con la

 * georreferenciaci?n, y dos ficheros por banda. Un .mpr con la informaci?n referente a los

 * datos de la banda y un .mp# con los datos de la banda

 * 

 * @version 04/06/2007

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

 */

public class ILWIS_MprFeatures extends WriteFileFormatFeatures {

	public ILWIS_MprFeatures() {

		super(GdalProvider.FORMAT_ILWIS, "mpl", new int[] { -1 }, new int[] { 0, 1, 2, 3, 4, 5 }, GdalWriter.class);

	}

	/**

	 * Carga los par?metros de este driver.

	 */

	public void loadParams() {

		super.loadParams();

	}

}

/* gvSIG. Geographic Information System of the Valencian Government

 *

 * Copyright (C) 2007-2008 Infrastructures and Transports Department

 * of the Valencian Government (CIT)

 *

 * 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., 51 Franklin Street, Fifth Floor, Boston,

 * MA  02110-1301, USA.

 *

 */

package org.gvsig.raster.gdal.io.features;

import org.gvsig.fmap.dal.coverage.datastruct.Params;

import org.gvsig.raster.gdal.io.JpegWriter;

import org.gvsig.raster.impl.store.WriteFileFormatFeatures;

/**

 * Caracteristicas del formato Jpeg para escritura.

 * Soporta escritura de imagenes 1 banda (monocromo) o 3 (RGB) en 8 bits. 

 * La georreferenciaci?n puede hacerse mediante un fichero .wld.

 * 

 * @version 04/06/2007

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

 *

 */

public class JpegFeatures extends WriteFileFormatFeatures {

	public JpegFeatures() {

		super("Jpeg", "jpg", new int[]{3}, null, JpegWriter.class);

	}

	/**

	 * Carga los par?metros de este driver.

	 */

	public void loadParams() {

		super.loadParams();

		driverParams.setParam("quality",

				new Integer(16),

				Params.SLIDER,

				new String[]{ "10", "100", "16", "5", "30"}); //min, max, valor defecto, intervalo peque?o, intervalo grande;

		driverParams.setParam("progressive", 

				new Boolean("false"), 

				Params.CHECK, 

				null);

	}

}

/* gvSIG. Geographic Information System of the Valencian Government

 *

 * Copyright (C) 2007-2008 Infrastructures and Transports Department

 * of the Valencian Government (CIT)

 *

 * 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., 51 Franklin Street, Fifth Floor, Boston,

 * MA  02110-1301, USA.

 *

 */

package org.gvsig.raster.gdal.io.features;