| gvSIG desktop 1

/*

 * Cresques Mapping Suite. Graphic Library for constructing mapping applications.

 *

 * 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

 * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307,USA.

 *

 * For more information, contact:

 * cresques@gmail.com

 */

package org.cresques.io;

 * Probablemente esto deber?a formar parte del JNI que recubre a la

 * librer?a en C extraida de gdal.<br>

 * Lo pongo aqu? a manera de ejemplo de como atacar un formato binario

 * @author Luis W. Sevilla.

 */

class GdalNative extends Gdal {

	static boolean 				WITH_OVERVIEWS = true;

	private String 				ext = "";

	public 	GeoTransform 		trans = null;

	/**

	 * Contorno en coordenadas geogr?ficas. (y Extent del raster).

	public Contour 				bBoxRot = new Contour();

	/**

	 * Contorno en coordenadas geogr?ficas sin rotaci?n aplicada. Esto es util para poder

	 */

	public Contour				bBoxWithoutRot = new Contour();

	public int 					width = 0, height = 0;

	private Metadata			metadata = null;

	private boolean 			georeferenced = true;

	// Polilinea con extent

	public class Contour extends Vector {

		final private static long serialVersionUID = -3370601314380922368L;

			if (pt.getY() < minY) minY = pt.getY();

		}

	}

	public GdalNative(String fName) throws GdalException, IOException {

		super();

		init(fName);

	}

	/**

	 * <P>

	 * Calcula la bounding box en la que est? metido el raster teniendo en cuenta

	 * leidos por gdal en el vector de 6 elementos adfGeoTransform donde cada elemento

	 * del vector represnta los siguientes valores

	 * </P>

	 * PtoY = originY + shearY * x + pixelSizeY * y;<BR>

	 * Aplicandolo a las cuatro esquinas sustituimos en cada una de ellas por.

	 * </P>

	 * <LI>Esquina superior izquierda: x = 0; y = 0;</LI>

	 * <LI>Esquina superior derecha: x = MaxX; y = 0;</LI>

	 * <LI>Esquina inferior izquierda: x = 0; y = MaxY;</LI>

	 * <LI>Esquina inferior derecha: x = MaxX; y = MaxY;</LI>

	 * <P>

	 * quedandonos en los cuatro casos:

	 * </P>

	 * <LI>Esquina superior izquierda: originX; originY;</LI>

	 * <LI>Esquina superior derecha: PtoX = originX + pixelSizeX * x; PtoY = originY + shearY * x;</LI>

	 * <LI>Esquina inferior izquierda:  PtoX = originX + shearX * y; PtoY = originY + pixelSizeY * y;</LI>

	 * <LI>Esquina inferior derecha: PtoX = originX + pixelSizeX * x + shearX * y; PtoY = originY + shearY * x + pixelSizeY * y;</LI>

	 * </UL>

	 */

	private void boundingBoxFromGeoTransform(){

		double geoX = 0D, geoY = 0D;

		geoX = trans.adfgeotransform[0] +  trans.adfgeotransform[2] * height;

		geoY = trans.adfgeotransform[3] +  trans.adfgeotransform[5] * height;

		bBoxRot.add(new Point2D.Double(geoX, geoY));

		//Upper right corner

		geoX = trans.adfgeotransform[0] + trans.adfgeotransform[1] * width;

		geoY = trans.adfgeotransform[3] + trans.adfgeotransform[4] * width;

		bBoxRot.add(new Point2D.Double(geoX, geoY));

		//Lower right corner

		geoX = trans.adfgeotransform[0] + trans.adfgeotransform[1] * width + trans.adfgeotransform[2] * height;

		geoY = trans.adfgeotransform[3] + trans.adfgeotransform[4] * width + trans.adfgeotransform[5] * height;

		bBoxRot.add(new Point2D.Double(geoX, geoY));

		//TODO: ?OJO! con coordenadas geogr?ficas

	}

	/**

	 * Calcula la bounding box en la que est? metido el raster teniendo en cuenta

	 */

	private void boundingBoxWithoutRotation(){

		double ox = trans.adfgeotransform[0];

		double oy = trans.adfgeotransform[3];

		double resx = trans.adfgeotransform[1];

		double resy = trans.adfgeotransform[5];

		bBoxWithoutRot.add(new Point2D.Double(ox, oy));

		bBoxWithoutRot.add(new Point2D.Double(ox + resx * width, oy));

		bBoxWithoutRot.add(new Point2D.Double(ox, oy + resy * height));

		//TODO: ?OJO! con coordenadas geogr?ficas

	}

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

		open(fName,GA_ReadOnly);

		ext = fName.toLowerCase().substring(fName.lastIndexOf('.')+1);

		width = getRasterXSize();

		height = getRasterYSize();

		setDataType(this.getRasterBand(1).getRasterDataType());

		metadata = new Metadata(getMetadata());

		//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)

		metadata.initColorInterpretation(getRasterCount());

		metadata.initNoDataByBand(getRasterCount());

		try{

			trans = getGeoTransform();

			boundingBoxWithoutRotation();

			boundingBoxFromGeoTransform();

			this.georeferenced = true;

		}catch(GdalException exc){

			bBoxRot.add(new Point2D.Double(0, 0));

			this.georeferenced = false;

		}

	}

	public void setAlpha(int a) { alpha = a; }

	public void setDataType(int dt) { dataType = dt; }

	public int getDataType() { return dataType; }

	double lastReadLine = -1;

	int currentFullWidth = -1;

	int currentFullHeight = -1;

	double stepX = 0D;

	double stepY = 0D;

	int currentOverview = -1;

	protected GdalRasterBand bandR = null, bandG = null, bandB = null, bandA = null;

	private boolean[] orientation;

	/**

	 * Devuelve la banda actualmente en uso para el color especificado.

	 * @param color 0=Rojo, 1=Green, 2=Blue.

	 * @return

	 */

	public GdalRasterBand getCurrentBand(int color) {

			return bandR;

			return bandG;

		return bandB;

	}

	//Supone rasters no girados

	public Point2D worldToRaster(Point2D pt) {

		double x = (((double) currentFullWidth) / (bBoxWithoutRot.maxX - bBoxWithoutRot.minX)) * (pt.getX() - bBoxWithoutRot.minX);

		Point2D ptRes = new Point2D.Double(x, y);

		return ptRes;

	}

	/**

	 * Si el tama?o de pixel en X es menor que 0 entonces la imagen se orienta al contrario en X por lo que en los zooms

	 * siguiente forma:

	 * Nuevo_punto_inicialX = (Ancho_total_raster - punto_inicial_del_zoomX) - Ancho_de_petici?n

	 *

	 * siguiente forma:

	 * Nuevo_punto_inicialY = (Alto_total_raster - punto_inicial_del_zoomY) - Alto_de_petici?n

	 * @param dWorldTLX

	 * @param dWorldTLY

	 * @param dWorldBRX

	 * @return

	 */

	public int setView(double dWorldTLX, double dWorldTLY,

		int err = 0;

		this.orientation = orientation;

		currentFullWidth = width;

		currentViewWidth = nWidth;

		currentViewHeight = nHeight;

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

		if(!orientation[0]) //Invierte la orientaci?n en X

			currentViewX = (width - tl.getX()) - (br.getX()-tl.getX());

		else

			currentViewX = tl.getX();

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

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

		stepX = 1D/viewportScaleX;

			lastReadLine = (height - tl.getY()) - (br.getY()-tl.getY());

		else

			lastReadLine = tl.getY();

		try {

			// calcula el overview a usar

			bandR = getRasterBand(1);

			currentOverview = -1;

			if (WITH_OVERVIEWS && bandR.getOverviewCount() > 0) {

				GdalRasterBand ovb = null;

					ovb = bandR.getOverview(i);

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

						currentOverview = i;

					}

				}

			}

			// Selecciona las bandas y los overviews necesarios

			bandR = getRasterBand(rBandNr);

			setDataType(bandR.getRasterDataType());

			if (this.getRasterCount() > 1) {

				bandG = getRasterBand(gBandNr);

				if(metadata.isAlphaBand())

			}

			if (currentOverview > 0) {

				bandR = bandR.getOverview(currentOverview);

					bandG = bandG.getOverview(currentOverview);

					bandB = bandB.getOverview(currentOverview);

					if(metadata.isAlphaBand())

		}

		return err;

	}

	int lastY = -1;

	public void readLine(int[][] line) throws GdalException {

			r = bandR.readRaster(x, y, w, 1, w, 1, dataType);

			if (bandG != null)

			if (bandB != null)

		lastReadLine += stepY;

			for (int k=0; k<4; k++){

			}

		return;

	}

	//int liney = 0;

	int readLineRGBA(int [] line) throws GdalException {

		int err = 0;

			if (bandG != null)

			if (bandB != null)

		int alpha = (this.alpha & 0xff) << 24;

		if(orientation[0]){ //Pixel size en X positivo

		}else{ //Pixel size en X negativo

			if (dataType == GDT_Byte){

		}

		return err;

	}

	/**

	 * Lee una franja de la imagen.

	 * @param bandH Altura de la franja

	 */

	public int readBandRGBA(int bandH, int bufH, int [] buf) throws GdalException {

		int err = 0;

			if (bandG != null)

			if (bandB != null)

		int alpha = (this.alpha & 0xff) << 24;

		return err;

	}

	/* (non-Javadoc)

			return null;

		}

	}

	void pintaInfo() {

		try {

			//System.out.println("Origin = "+originX+","+originY);

				System.out.println(metadata[i]);

			}

		} catch (GdalException e) {

		}

	}

	void pintaPaleta() {

	}

	public int getBlockSize(){

		return this.getBlockSize();

	}

	protected GdalNative 		file = null;

	private Extent v = null;

	public GdalFile(IProjection proj, String fName){

		super(proj, fName);

		extent = new Extent();

			file = new GdalNative(fName);

			load();

			readGeoInfo(fName);

			if ( bandCount > 2) {

				setBand(RED_BAND,   0);

				setBand(GREEN_BAND, 1);

			} else

				setBand(RED_BAND|GREEN_BAND|BLUE_BAND, 0);

		} catch(Exception e){

		}

		switch(file.getDataType()){

		case 4://GDT_UInt32

		case 5:setDataType(DataBuffer.TYPE_INT);break;//GDT_Int32

		case 6:setDataType(DataBuffer.TYPE_FLOAT);break;//GDT_Float32

		case 10:setDataType(DataBuffer.TYPE_UNDEFINED);break;//GDT_CFloat32

		case 11:setDataType(DataBuffer.TYPE_UNDEFINED);break;//GDT_CFloat64

		}

	}

	/**

	 * Obtenemos o calculamos el extent de la imagen.

	 */

		requestExtent = new Extent(file.bBoxWithoutRot.minX, file.bBoxWithoutRot.minY, file.bBoxWithoutRot.maxX, file.bBoxWithoutRot.maxY);

		return this;

	}

	/**

	 * Cierra el fichero de imagen

	 */

			e.printStackTrace();

		}

	}

	/**

	 * Asigna a cada banda R,G o B una banda de la imagen

	 */

		if ((flag & GeoRasterFile.GREEN_BAND) == GeoRasterFile.GREEN_BAND) file.gBandNr = bandNr+1;

		if ((flag & GeoRasterFile.BLUE_BAND) == GeoRasterFile.BLUE_BAND) file.bBandNr = bandNr+1;

	}

	/**

	 * Asigna el extent de la vista actual. existe un fichero .rmf debemos hacer una transformaci?n

	 * de la vista asignada ya que la petici?n viene en coordenadas del fichero .rmf y la vista (v)

	 * ha de estar en coordenadas del fichero.

	 */

		if(rmfExists){

			Point2D.Double petInit = null, petEnd = null;

			try{

				petInit = new Point2D.Double(e.minX(),  e.maxY());

			}catch(NoninvertibleTransformException ex){}

			double h = file.bBoxWithoutRot.maxY - file.bBoxWithoutRot.minY;

			if(!file.isGeoreferenced())

			else

				v = new Extent(	petInit.getX(), petInit.getY(), petEnd.getX(), petEnd.getY());

		}else

	}

	 /**

	 * Calcula la transformaci?n que se produce sobre la vista cuando la imagen tiene un fichero .rmf

	 * asociado. En Gdal el origen de coordenadas en Y es el valor m?nimo y crece hasta el m?ximo. De la

	 * @param originX Origen de la imagen en la coordenada X

	 * @param originY Origen de la imagen en la coordenada Y

	 */

		transformRMF.setToTranslation(originX, originY);

		transformRMF.scale(psX, psY);

			transformTFW.setToTranslation(file.trans.adfgeotransform[0], file.trans.adfgeotransform[3]);

			transformTFW.scale(file.trans.adfgeotransform[1], file.trans.adfgeotransform[5]);

		}

	}

	/**

	 * Obtiene extent de la vista actual

	 */

	}

	/**

	 * Obtiene la anchura del fichero

	 */

	}

	/**

	 * Obtiene la altura del fichero

	 */

		return file.height;

	}

	public void reProject(ICoordTrans rp) {

		// TODO Auto-generated method stub

	}

	/**

	 * Obtiene la orientaci?n de la imagen a partir del signo del tama?o de pixel para poder

	 * en este caso las X de menor a mayor y las Y de mayor a menor. Los casos posibles son:

	 * <UL>

	 * <LI><B>X > 0; Y < 0;</B> {true, false}</LI>

	 * <LI><B>X < 0; Y > 0;</B> {false, true}</LI>

	 * <LI><B>X < 0; Y < 0;</B> {false, false}</LI>

	 * </UL>

	 * @return

	 */

	private boolean[] getOrientation(){

		}

		return orientation;

	}

	/* (non-Javadoc)

	 * @see org.cresques.io.GeoRasterFile#updateImage(int, int, org.cresques.cts.ICoordTrans)

	 */

	public Image updateImage(int width, int height, ICoordTrans rp) {

		int line, pRGBArray[] = null;

		Image image = null;

		if (mustVerifySize()) {

			// Work out the correct aspect for the setView call.

			double dFileAspect = (double)v.width()/(double)v.height();

			double dWindowAspect = (double)width /(double)height;

			if (dFileAspect > dWindowAspect) {

			} else {

			}

		}

		file.setView(v.minX(), v.maxY(), v.maxX(), v.minY(),

			width, height, getOrientation());

		if(width<=0)width=1;

		if(height<=0)height=1;

		image = new BufferedImage(width, height, BufferedImage.TYPE_INT_ARGB);

		//image = new BufferedImage(width, height, BufferedImage.TYPE_INT_RGB);

		pRGBArray = new int[width/**BAND_HEIGHT*/];

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

		return image;

	}

	public RasterBuf getRaster(int width, int height, ICoordTrans rp) {

		int line;

		RasterBuf raster = null;

		if(mustVerifySize()){

			// Work out the correct aspect for the setView call.

			double dFileAspect = (double)v.width()/(double)v.height();

			double dWindowAspect = (double)width /(double)height;

			if (dFileAspect > dWindowAspect) {

			} else {

			}

		}

		// Set the view

		boolean[] orientation = getOrientation();

		file.setView(v.minX(), v.maxY(), v.maxX(), v.minY(),

			width, height, orientation);

		raster = new RasterBuf(DataBuffer.TYPE_INT, width, height, 4, new Point(0,0));

		try {

		} catch (Exception e) {

			e.printStackTrace();

		}

		return raster;

	}

	/**

	 * vector de enteros.

	 * @param image	imagen con los datos actuales

	 * @param startX	inicio de la posici?n en X dentro de la imagen

	 * @param offset	desplazamiento

	 * @param scansize	tama?o de imagen recorrida por cada p

	 */

			 int offset, int scansize) {

		image.setRGB(startX, startY, w, h, rgbArray, offset, scansize);

	}

	/**

	 * con los dato de la imagen contenidos en el vector de enteros. De los valores RGB

	 * que ya contiene se mantienen las bandas que no coinciden con el valor de flags. La

	 * banda correspondiente a flags es sustituida por los datos del vector.

	 * @param scansize	tama?o de imagen recorrida por cada paso

	 * @param flags	banda que se va a sustituir (Ctes de GeoRasterFile)

	 */

			 int offset, int scansize, int flags) {

		image.getRGB(startX, startY, w, h, line, offset, scansize);

		if (flags == GeoRasterFile.RED_BAND)

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

				line[i] = (line[i] & 0x00ffff00) | (rgbArray[i] & 0xff0000ff);

		image.setRGB(startX, startY, w, h, line, offset, scansize);

	}

	/**

	 * con los dato de la imagen contenidos en el vector de enteros. De los valores RGB

	 * que ya contiene se mantienen las bandas que no coinciden con el valor de flags. La

	 * banda correspondiente a flags es sustituida por los datos del vector.

	 * @param origBand	Banda origen del GeoRasterFile

	 * @param destBandFlag	banda que se va a sustituir (Ctes de GeoRasterFile)

	 */

			 int offset, int scansize, int origBand, int destBandFlag) {

		image.getRGB(startX, startY, w, h, line, offset, scansize);

		if (origBand == 0 && destBandFlag == GeoRasterFile.RED_BAND)

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

		else if (origBand == 2 && destBandFlag == GeoRasterFile.BLUE_BAND)

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

				line[i] = (line[i] & 0x00ffff00) | (rgbArray[i] & 0xff0000ff);

		else if (origBand == 0 && destBandFlag == GeoRasterFile.GREEN_BAND)

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

				line[i] = (line[i] & 0xffff00ff) | ((rgbArray[i] & 0x00ff0000) >> 8) ;

		else if (origBand == 1 && destBandFlag == GeoRasterFile.RED_BAND)

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

				line[i] = (line[i] & 0xff00ffff) | ((rgbArray[i] & 0x0000ff00) << 8);

		else if (origBand == 1 && destBandFlag == GeoRasterFile.BLUE_BAND)

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

				line[i] = (line[i] & 0xffffff00) | ((rgbArray[i] & 0x0000ff00) >> 8);

				line[i] = (line[i] & 0xffff00ff) | ((rgbArray[i] & 0x000000ff) << 8);

		image.setRGB(startX, startY, w, h, line, offset, scansize);

	}

	private void showOnOpen() {

			System.out.println("   NumBands = ("+file.getRasterCount()+")");

		} catch (GdalException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

	}

	/* (non-Javadoc)

	 * @see org.cresques.io.GeoRasterFile#updateImage(int, int, org.cresques.cts.ICoordTrans, java.awt.Image, int, int)

	 */

	public Image updateImage(int width, int height, ICoordTrans rp, Image img, int origBand, int destBandFlag)throws SupersamplingNotSupportedException{

		int line, pRGBArray[] = null;

		if(mustVerifySize()){

			// Work out the correct aspect for the setView call.

			double dFileAspect = (double)v.width()/(double)v.height();

			double dWindowAspect = (double)width /(double)height;

			if (dFileAspect > dWindowAspect) {

			} else {

			}

		}

		// Set the view

		boolean[] orientation = getOrientation();

		file.setView(v.minX(), v.maxY(), v.maxX(), v.minY(),

			width, height, orientation);

		if(width<=0)width=1;

		if(height<=0)height=1;

		pRGBArray = new int[width];

		try {

			setBand(RED_BAND,   rBandNr);

					for (line=0; line < height; line++) {

						file.readLineRGBA(pRGBArray);

						setRGBLine((BufferedImage) img, 0, height - 1 - line, width, 1, pRGBArray, 0, width, origBand, destBandFlag);

				}else{

					for (line=0; line < height; line++) {

						file.readLineRGBA(pRGBArray);

					for (line=0; line < height; line++) {

						file.readLineRGBA(pRGBArray);

						setRGBLine((BufferedImage) image, 0, height - 1 - line, width, 1, pRGBArray, 0, width);

				}else{

					for (line=0; line < height; line++) {

						file.readLineRGBA(pRGBArray);

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

		return img;

	}

	/* (non-Javadoc)

	 * @see org.cresques.io.GeoRasterFile#getData(int, int, int)

	 */

		}

		return null;

	}

	/**

	 * Devuelve los datos de una ventana solicitada

	 * @param ulX	coordenada X superior izda.

	 * @param band	Banda solicitada.

	 */

	public byte[] getWindow(int ulX, int ulY, int sizeX, int sizeY, int band){

		return null;

	}

	/**

	 * Obtiene la zona (Norte / Sur)

	 * @return true si la zona es norte y false si es sur

	 */

	public boolean getZone(){

		return false;

	}

	/**

	 *Devuelve el n?mero de zona UTM

	 */

	public int getUTM(){

	}

	/**

	 * Obtiene el sistema de coordenadas geograficas

	 * @return Sistema de coordenadas geogr?ficas

	 */

	public String getGeogCS(){

	}

	/**

	 * Devuelve el tama?o de bloque

	 * @return Tama?o de bloque

	 */

	public int getBlockSize(){

	}

	/**

	 * Obtiene el objeto que contiene los metadatos

	 */

	public Metadata getMetadata() {

		if(file != null)

			return file.getMetadataJavaObject();

			return null;

	}

	/**

	 * Obtiene el flag que dice si la imagen est? o no georreferenciada

	 * @return true si est? georreferenciada y false si no lo est?.

	public boolean isGeoreferenced() {

		return file.isGeoreferenced();

	}

	/**

	 * Obtiene los par?metros de la transformaci?n af?n que corresponde con los elementos de

	 * un fichero tfw.

	 * <LI>[1]tama?o de pixel en X</LI>

	 * <LI>[2]rotaci?n en X</LI>

	 * <LI>[4]rotaci?n en Y</LI>

	 * </UL>

	 * Este m?todo debe ser reimplementado por el driver si tiene esta informaci?n. En principio

	 * Gdal es capaz de proporcionarla de esta forma.

	 * En caso de que exista fichero .rmf asociado al raster pasaremos de la informaci?n de georreferenciaci?n

	 * del .tfw y devolveremos la que est? asociada al rmf

	 * @return vector de double con los elementos de la transformaci?n af?n.

			return file.trans.adfgeotransform;

		else{

			if(this.rmfExists){

										rmfTransform.getScaleX(),

										rmfTransform.getTranslateY(),

										rmfTransform.getShearY(),

										rmfTransform.getScaleY()};

			}

			return null;

		}

	}

}