Application: gvSIG desktop

/*

 * Created on 12-dic-2004

 */

package org.cresques.io;

import org.cresques.px.Extent;

/**

 * Consulta el servidor de mapas de grafcan.

 * Implementa una llamada del tipo:

 * http://www.dgtcan.com/MapServer/Mapa.asp?

 *   XMin=219900&YMax=3151900&XMax=221000&YMin=3151300&

 *   Ancho=1000&Alto=1000&NumeroCapas=6&Capas=2&Servicio=NavGRAFCAN2

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

 */

public class GrafCanMapServerClient extends MapServerClient {

	public final static int CAPA_ORTOFOTO    = 2;

	public final static int CAPA_CARTOGRAFIA = 1;

	private int numeroCapas	= 6;

	private int capas		= CAPA_CARTOGRAFIA;

	private String servicio = "NavGRAFCAN2";

	public GrafCanMapServerClient(String serverName) {

		super(serverName);

		urlBase = "http://www.dgtcan.com/MapServer/Mapa.asp";

/**/	setExtent(new Extent(180020,3270000,660000,3035000));

	}

	public void setCapas(int c) { capas = c; }

	public void setNumeroCapas(int c) { numeroCapas = c; }

	public void setServicio(String s) { servicio = s; }

	public String getUrl() {

		String url = urlBase+"?"+

			"XMin="+xMin+"&"+ "XMax="+xMax+"&"+

			"YMin="+yMin+"&"+ "YMax="+yMax+"&"+

			"Ancho="+ancho+"&"+"Alto="+alto+"&"+

			"NumeroCapas="+numeroCapas+"&"+"Capas="+capas+"&"+

			"Servicio="+servicio;

		return url;

	}

}

/*

 * Created on 25-sep-2004

 *

 * TODO To change the template for this generated file go to

 * Window - Preferences - Java - Code Style - Code Templates

 */

package org.cresques.io;

import java.awt.Color;

import java.awt.Image;

import java.awt.geom.Point2D;

import java.awt.image.BufferedImage;

import java.io.IOException;

import java.util.Vector;

import org.cresques.cts.ICoordTrans;

import org.cresques.cts.IProjection;

import org.cresques.io.GeoFile;

import org.cresques.io.GeoRasterFile;

import org.cresques.px.Extent;

import es.gva.cit.jbsb.BsbException;

/**

 * Soporte 'nativo' para Bsb.

 * Este conjunto de funcionalidades est? tomado de manera casi literal

 * del soporte para ECW de ermapper.<br>

 * 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

 * desde Java.<br><br>   

 * @author Luis W. Sevilla.

 */

class BsbNative extends es.gva.cit.jbsb.Bsb {

	// Polilinea con extent

	class Contour extends Vector {

		public double minX = Double.MAX_VALUE, minY = Double.MAX_VALUE;

		public double maxX = -Double.MAX_VALUE, maxY = -Double.MAX_VALUE;

		public Contour() {

			super();

		}

		public void add(Point2D pt) {

			super.add(pt);

			if (pt.getX() > maxX) maxX = pt.getX();

			if (pt.getX() < minX) minX = pt.getX();

			if (pt.getY() > maxY) maxY = pt.getY();

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

		}

	}

	private int [] paleta = null;

	/**

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

	 */

	public Contour esq = new Contour();

	public int width = 0, height = 0;

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

	public String version = "";

	public BsbNative(String fName) throws BsbException {

		super();

		try {

			init(fName);

		} catch (IOException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

	}

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

		bsbOpen(fName);

		width = g_psInfo.nXSize;

		height = g_psInfo.nYSize;

		Vector pal = new Vector();

		String buf = null;

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

  			buf = g_psInfo.papszHeader[i];

	  		if (buf.length()<1) continue;

	  		if (buf.startsWith("RGB/")) {

	  			String [] dat = buf.split(",");

	  			pal.add(new Color(Integer.parseInt(dat[1]),

	  				Integer.parseInt(dat[2]),

					Integer.parseInt(dat[3])));

	  		} else if (buf.startsWith("VER/")) {

	  			version = buf.substring(4);

	  		} /*else if (buf.startsWith("PLY/")) { // Marco de hoja

	  			String [] dat = buf.split(",");

	  			esq.add(new Point2D.Double(Double.parseDouble(dat[2]),

	  					Double.parseDouble(dat[1])));

	  		}*/

  		}

  		if (true) { //version.startsWith("1")) {

  			esq.add(new Point2D.Double(0D,height));

  			esq.add(new Point2D.Double(0D,0D));

  			esq.add(new Point2D.Double(width,0D));

  			esq.add(new Point2D.Double(width,height));

  		}

  		paleta = new int[pal.size()];

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

  			paleta[i] = ((Color) pal.get(i)).getRGB();

  		originX = esq.minX;

  		originY = esq.maxY;

	}

	double lastReadLine = -1;

	int currentViewWidth = -1;

	int currentViewHeight = -1;

	double currentViewX = 0D;

	double viewportScale = 0D;

	double step = 0D;

	public Point2D worldToRaster(Point2D pt) {

		double x = (((double) width)/(esq.maxX-esq.minX))*(pt.getX()-esq.minX);

		double y = (((double) height)/(esq.maxY-esq.minY))*(esq.maxY-pt.getY());

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

		return ptRes;

	}

	public int setView(double dWorldTLX, double dWorldTLY,

            double dWorldBRX, double dWorldBRY,

            int nWidth, int nHeight) {

		int err = 0;

		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;

		currentViewX = tl.getX();

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

		step = 1D/viewportScale;

		lastReadLine = tl.getY();

		System.out.println("BsbFile: TL=("+dWorldTLX+","+dWorldTLY+

			"); BR=("+dWorldBRX+","+dWorldBRY+")\n"+

			"BsbFile: escala="+viewportScale+"; lastReadLine="+lastReadLine);

		return err;

	}

	public int readLineRGBA(int [] line) throws BsbException {

		int err = 0;

  		bsbReadLine((int) lastReadLine);

  		lastReadLine += step;

  		byte [] l = g_buffer.pabyScanLineBuf;

  		int white = Color.BLUE.getRGB();

  		float j =(float) currentViewX;

  		for (int i=0; i<currentViewWidth && j<l.length; i++, j+=step) {

  			line[i] = paleta[l[(int) j]-1];

  		}

		//for (int i=0; i<currentViewWidth; i++) line[i] = 128+128*256+128*256*256;

		return err;

	}

	void pintaInfo() {

  		System.out.println("Fichero BSB: Info ("+g_psInfo.papszHeader.length+" l?neas):");

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

	  		if (g_psInfo.papszHeader[i].length()<1) continue;

	  		System.out.println("info["+i+"] = "+g_psInfo.papszHeader[i]);

  		}

	}

	void pintaPaleta() {

  		System.out.println("Fichero BSB: Paleta ("+paleta.length+" colores):");

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

			System.out.println("color("+i+")="+paleta[i]);

	}

}

/**

 * @author Luis W. Sevilla

 */

public class BsbFile extends GeoRasterFile {

	private BsbNative file = null;

	private Extent v = null;

	static {

		GeoRasterFile.registerExtension("nos", BsbFile.class);

		GeoRasterFile.registerExtension("kap", BsbFile.class);

	}

	public BsbFile(IProjection proj, String fName) {

		super(proj, fName);

		extent = new Extent();

		try {

			file = new BsbNative(fName);

			showOnOpen();

			load();

		} catch(BsbException e){

		  	System.out.println("Error en BSBOpen");

		  	e.printStackTrace();

		  	file = null;

		}

	}

	public GeoFile load() {

		/*double minX, minY, maxX, maxY;

		minX = minY = 0D;

		maxX = (double) width;

		maxY = (double) height;*/

		extent = new Extent(file.esq.minX, file.esq.minY, file.esq.maxX, file.esq.maxY);

		return this;

	}

	public void setView(Extent e) { v = new Extent(e); }

	public Extent getView() { return v; }

	public int getWidth() {	return file.width; }

	public int getHeight() { return file.height;}

	/* (non-Javadoc)

	 * @see org.cresques.io.GeoRasterFile#reProject(org.cresques.cts.ICoordTrans)

	 */

	public void reProject(ICoordTrans rp) {

		// TODO Auto-generated method stub

	}

	/* (non-Javadoc)

	 * @see org.cresques.io.GeoRasterFile#setTransparency(boolean)

	 */

	public void setTransparency(boolean t) {

		// TODO Auto-generated method stub

	}

	public void setTransparency(int t) {

		// TODO Auto-generated method stub

	}

	/* (non-Javadoc)

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

	 */

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

		double dFileAspect, dWindowAspect;

		int line, pRGBArray[] = null;

		Image image = null;

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

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

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

		if (dFileAspect > dWindowAspect) {

		  height =(int)((double)width/dFileAspect);

		} else {

		  width = (int)((double)height*dFileAspect);

		}

		// Set the view

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

			width, height);

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

		pRGBArray = new int[width];

		try {

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

				file.readLineRGBA(pRGBArray);

				((BufferedImage)image).setRGB(0, line, width, 1, pRGBArray, 0, width);

			}

		} catch (Exception e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

		return image;

	}

	private void showOnOpen() {

  		// Report en la apertura (quitar)

  		System.out.println("Fichero BSB '"+getName()+"' abierto.");

  		System.out.println("Version = "+file.version);

  		System.out.println("   Size = ("+file.width+","+file.height+")");

  		file.pintaInfo();

  		file.pintaPaleta();

	}

	/* (non-Javadoc)

	 * @see org.cresques.io.GeoRasterFile#close()

	 */

	public void close() {

	}

	/* (non-Javadoc)

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

	 */

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

		// TODO Auto-generated method stub

		return null;

	}

	/* (non-Javadoc)

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

	 */

	public Object getData(int x, int y, int band) {

		// TODO Auto-generated method stub

		return null;

	}

}

/*

 * Created on 25-sep-2004 (20-oct-2004)

 */

package org.cresques.io;

import java.awt.Color;

import java.awt.Image;

import java.awt.Point;

import java.awt.geom.Point2D;

import java.awt.image.BandedSampleModel;

import java.awt.image.BufferedImage;

import java.awt.image.DataBuffer;

import java.awt.image.Raster;

import java.awt.image.WritableRaster;

import java.io.IOException;

import java.util.Vector;

import org.cresques.cts.ICoordTrans;

import org.cresques.cts.IProjection;

import org.cresques.px.Extent;

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.

 * Este conjunto de funcionalidades est? tomado de manera casi literal

 * del soporte para ECW de ermapper.<br>

 * 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

 * desde Java.<br><br>   

 * @author Luis W. Sevilla.

 */

class GdalNative extends Gdal {

	static boolean WITH_OVERVIEWS = true;

	// Polilinea con extent

	class Contour extends Vector {

		public double minX = Double.MAX_VALUE, minY = Double.MAX_VALUE;

		public double maxX = -Double.MAX_VALUE, maxY = -Double.MAX_VALUE;

		public Contour() {

			super();

		}

		public void add(Point2D pt) {

			super.add(pt);

			if (pt.getX() > maxX) maxX = pt.getX();

			if (pt.getX() < minX) minX = pt.getX();

			if (pt.getY() > maxY) maxY = pt.getY();

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

		}

	}

	/**

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

	 */

	public Contour esq = new Contour();

	public int width = 0, height = 0;

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

	public String version = "";

	private int alpha = 0;

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

	private int dataType = GDT_Byte;

	public GdalNative(String fName) throws GdalException, IOException {

		super();

		init(fName);

	}

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

		open(fName,GA_ReadOnly);

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

		if (ext.compareTo("tif") == 0)

			WITH_OVERVIEWS = false;

		width = getRasterXSize();

		height = getRasterYSize();

		if (true) { //ext.compareTo("sid") == 0) {

			String [] metadata = getMetadata();

			double ox=0D, oy=0D, resx=0D, resy=0D;

			GeoTransform trans = getGeoTransform();

			ox = trans.adfgeotransform[0];

			oy = trans.adfgeotransform[3];

			resx = trans.adfgeotransform[1];

			resy = trans.adfgeotransform[5];

			System.out.println("Origin = ("+ox+","+oy+")");

			System.out.println("Pixel Size = ("+resx+","+resy+")");

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

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

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

  			esq.add(new Point2D.Double(ox+resx*width, oy+resy*height));

		} else  { //version.startsWith("1")) {

  			//double [] transParam = getGeoTransform();

  			esq.add(new Point2D.Double(0D,height));

  			esq.add(new Point2D.Double(0D,0D));

  			esq.add(new Point2D.Double(width,0D));

  			esq.add(new Point2D.Double(width,height));

  		}

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

	}

	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;

	int currentViewWidth = -1;

	int currentViewHeight = -1;

	double currentViewX = 0D;

	double viewportScale = 0D;

	double step = 0D;

	int currentOverview = -1;

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

	/**

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

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

	 * @return

	 */

	public GdalRasterBand getCurrentBand(int color) {

		if (color == 0) return bandR;

		else if (color == 1) return bandG;

		return bandB;

	}

	// Supone rasters no girados

	public Point2D worldToRaster(Point2D pt) {

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

		double y = (((double) currentFullHeight)/(esq.maxY-esq.minY))*(esq.maxY-pt.getY());

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

		return ptRes;

	}

	public int setView(double dWorldTLX, double dWorldTLY,

            double dWorldBRX, double dWorldBRY,

            int nWidth, int nHeight) {

		int err = 0;

		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;

		currentViewX = tl.getX();

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

		step = 1D/viewportScale;

		lastReadLine = tl.getY();

		try {

			// calcula el overview a usar

			bandR = getRasterBand(1);

			currentOverview = -1;

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

				GdalRasterBand ovb = null;

				for (int i=bandR.getOverviewCount()-1; i>0; i--) {              

					ovb = bandR.getOverview(i);

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

						currentOverview = i;

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

	            		step = 1D/viewportScale;

	            		currentFullWidth = ovb.getRasterBandXSize();

	            		currentFullHeight = ovb.getRasterBandYSize();

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

	            		currentViewX = tl.getX();

	            		lastReadLine = tl.getY();

	            		break;

					}

				}

			}

			// Selecciona las bandas y los overviews necesarios

			bandR = getRasterBand(rBandNr);

			setDataType(bandR.getRasterDataType());

			if (this.getRasterCount() > 1) {

				bandG = getRasterBand(gBandNr);

				bandB = getRasterBand(bBandNr);

			}

			if (currentOverview > 0) {

				bandR = bandR.getOverview(currentOverview);

    			if (this.getRasterCount() > 1) {

					bandG = bandG.getOverview(currentOverview);

					bandB = bandB.getOverview(currentOverview);

				}

			}

			//System.out.println(band.)

		} catch (GdalException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

		System.out.println("GdalFile: TL=("+dWorldTLX+","+dWorldTLY+

			"); BR=("+dWorldBRX+","+dWorldBRY+")\n"+

			"GdalFile: escala="+viewportScale+"; lastReadLine="+lastReadLine+"\n"+

		    "Actual Raster Size="+currentFullWidth+"x"+currentFullHeight+

			"\nDataType="+dataType);

		return err;

	}

	int lastY = -1;

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

		int err = 0;

		int nbands = getRasterCount();

		GdalRasterBand band = null;

        int w = (int)(((double)currentViewWidth)*step);

        int x = (int) currentViewX;

        int y = (int) lastReadLine;

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

        //if (alpha > 0) a = alpha << 24;

        if (x+w > bandR.getRasterBandXSize()) 

        	w = bandR.getRasterBandXSize()-x;

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

		if (bandG != null)

    		g = bandG.readRaster(x, y, w, 1, w, 1, dataType);

		if (bandB != null)

    		b = bandB.readRaster(x, y, w, 1, w, 1, dataType);

  		lastReadLine += step;

  		int white = Color.BLUE.getRGB(), i2=0, i=0;

  		float j =0F;

  		if (dataType == GDT_CInt16 || dataType == GDT_Int16  || dataType == GDT_UInt16)

  			if (g == null) // Sibgle Band (Typical DEM)

	      		/*for (i=0, j=0F, i2 = 1; i<currentViewWidth && i2<r.length;

	      			i++, j+=step, i2 = (((int) j)*2)+1) {

	      			line[i] = a + ((r[i2-1]) << 8) + r[i2];

	      		}*/

  		  		for (i=0, j=0F; i<currentViewWidth && j<r.getSize(); i++, j+=step) {

	      			line[0][i] = r.buffShort[(int) j] & 0xffff;

	      			line[1][i] = r.buffShort[(int) j] & 0xffff;

	      			line[2][i] = r.buffShort[(int) j] & 0xffff;

	      		}

  			else { // Multiband

  				short px;

  				//System.err.println("readLine(): Raster 16bits multibanda");

	      		for (i=0, j=0F; i<currentViewWidth && j<r.getSize(); i++, j+=step) {

	      			line[0][i] = r.buffShort[(int) j] & 0xffff;

	      			line[1][i] = g.buffShort[(int) j] & 0xffff;

	      			line[2][i] = b.buffShort[(int) j] & 0xffff;

	      		}

  			}

	    //for (int i=0; i<currentViewWidth; i++) line[i] = 128+128*256+128*256*256;

		return;

	}

	public int readLineRGBA(int [] line) throws GdalException {

		int err = 0;

		int nbands = getRasterCount();

		GdalRasterBand band = null;

        int w = (int)(((double)currentViewWidth)*step);

        int x = (int) currentViewX;

        int y = (int) lastReadLine;

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

        //if (alpha > 0) a = alpha << 24;

        if (x+w > bandR.getRasterBandXSize()) 

        	w = bandR.getRasterBandXSize()-x;

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

		if (bandG != null)

    		g = bandG.readRaster(x, y, w, 1, w, 1, dataType);

		if (bandB != null)

    		b = bandB.readRaster(x, y, w, 1, w, 1, dataType);

  		lastReadLine += step;

  		int white = Color.BLUE.getRGB(), i2=0, i=0;

  		float j =0F;

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

  		if (dataType == GDT_Byte)

	  		if (g != null)

	      		for (i=0; i<currentViewWidth && j<r.getSize(); i++, j+=step) {

	      			line[i] = alpha + ((r.buffByte[(int) j] & 0xff) << 16) + ((g.buffByte[(int) j] & 0xff) << 8) + (b.buffByte[(int) j] & 0xff);

	      		}

	      	else

	      		for (i=0; i<currentViewWidth && j<r.getSize(); i++, j+=step) {

	      			line[i] = alpha + ((r.buffByte[(int) j] & 0xff) << 16) +

						((r.buffByte[(int) j] & 0xff) << 8) + (r.buffByte[(int) j] & 0xff);

	      		}

  		else if (dataType == GDT_CInt16 || dataType == GDT_Int16  || dataType == GDT_UInt16)

  			if (g == null) // Sibgle Band (Typical DEM)

	      		/*for (i=0, j=0F, i2 = 1; i<currentViewWidth && i2<r.length;

	      			i++, j+=step, i2 = (((int) j)*2)+1) {

	      			line[i] = a + ((r[i2-1]) << 8) + r[i2];

	      		}*/

  		  		for (i=0, j=0F; i<currentViewWidth && j<r.getSize(); i++, j+=step) {

	      			line[i] = alpha + r.buffShort[(int) j];

	      		}

  			else { // Multiband

  				System.err.println("Raster 16bits multibanda");

	      		for (i=0, j=0F; i<currentViewWidth && j<r.getSize(); i++, j+=step) {

	      			line[i] = alpha | (((r.buffShort[(int) j] & 0xfff0) << 12) & 0xff0000 ) | 

									  (((g.buffShort[(int) j] & 0xfff0) << 4 ) & 0xff00 ) |

									  (((b.buffShort[(int) j] & 0xfff0) >> 4 ) & 0xff );

	      		}

  			}

	    //for (int i=0; i<currentViewWidth; i++) line[i] = 128+128*256+128*256*256;

		return err;

	}

	/**

	 * Lee una franja de la imagen.

	 * @param bandH Altura de la franja

	 * @param bufH	Altura del buffer

	 * @param buf	Buffer con la franja (retorno)

	 * @return

	 * @throws GdalException

	 */

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

		int err = 0;

		int nbands = getRasterCount();

		GdalRasterBand band = null;

        int w = (int)(((double)currentViewWidth)*step);

        int x = (int)(((double)currentViewX)*step);

        int y = (int) lastReadLine;

        int h = (int) (((double)bandH)*step);

        System.out.println("Leyendo "+y);

        if (x+w > bandR.getRasterBandXSize()) 

        	w = bandR.getRasterBandXSize()-x;

		GdalBuffer r = bandR.readRaster(x, y, w, h, w, h, GDT_Byte);

		GdalBuffer g = bandG.readRaster(x, y, w, h, w, h, GDT_Byte);

		GdalBuffer b = bandB.readRaster(x, y, w, h, w, h, GDT_Byte);

  		lastReadLine += ((double)bandH)*step;

  		// TODO Acabar de implementarlo

  		float k=0F;

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

  		for (int j=0, t=0; j<bandH; j++) {

  			k = j*w; t=j*currentViewWidth;

  			for (int i=0; i<currentViewWidth && k<r.getSize(); i++, k+=step) {

  				buf[t+i] = alpha + ((r.buffByte[(int) k]) << 16) + ((g.buffByte[(int) k]) << 8) + b.buffByte[(int) k];

  			}

  		}

		//for (int i=0; i<currentViewWidth; i++) line[i] = 128+128*256+128*256*256;

		return err;

	}

	void pintaInfo() {

		try {

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

			//System.out.println("Origin = "+this.);

			System.out.println("GeoTransform:");

			GeoTransform trans = getGeoTransform();

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

				System.out.println("  param["+i+"]="+trans.adfgeotransform[i]);

			System.out.println("Metadata:");

			String [] metadata = getMetadata();

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

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

			}

		} catch (GdalException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

	}

	void pintaPaleta() {

	}

}

/**

 * @author Luis W. Sevilla

 */

public class GdalFile extends GeoRasterFile {

	public final static int BAND_HEIGHT = 64;

	protected GdalNative file = null;

	private Extent v = null;

	static {

		//GeoRasterFile.registerExtension("sid", GdalFile.class);

		GeoRasterFile.registerExtension("img", GdalFile.class);

		GeoRasterFile.registerExtension("tif", GdalFile.class);

	}

	public GdalFile(IProjection proj, String fName) {

		super(proj, fName);

		extent = new Extent();

		try {

			file = new GdalNative(fName);

			showOnOpen();

			load();

			bandCount = file.getRasterCount(); 

			if ( bandCount > 2) {

				setBand(RED_BAND,   0);

				setBand(GREEN_BAND, 1);

				setBand(BLUE_BAND,  2);

			} else

				setBand(RED_BAND|GREEN_BAND|BLUE_BAND, 0);

		} catch(Exception e){

		  	System.out.println("Error en GdalOpen");

		  	e.printStackTrace();

		  	file = null;

		}

		int dataType = file.getDataType();

  		if (dataType == Gdal.GDT_Byte)

  			setDataType(DataBuffer.TYPE_BYTE);

  		else if (dataType == Gdal.GDT_CInt16 || dataType == Gdal.GDT_Int16)

  			setDataType(DataBuffer.TYPE_SHORT);

  		else if(dataType == Gdal.GDT_UInt16)

  			setDataType(DataBuffer.TYPE_USHORT);

	}

	public GeoFile load() {

		/*double minX, minY, maxX, maxY;

		minX = minY = 0D;

		maxX = (double) width;

		maxY = (double) height;*/

		extent = new Extent(file.esq.minX, file.esq.minY, file.esq.maxX, file.esq.maxY);

		return this;

	}

	public void close() {

		try {

			file.close();

		} catch (GdalException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

	}

	public void setBand(int flag, int bandNr) {

		super.setBand(flag, bandNr);

		if ((flag & GeoRasterFile.RED_BAND) == GeoRasterFile.RED_BAND) file.rBandNr = bandNr+1;

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

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

	}

	public void setView(Extent e) { v = new Extent(e); }

	public Extent getView() { return v; }

	public int getWidth() {	return file.width; }

	public int getHeight() { return file.height;}

	/* (non-Javadoc)

	 * @see org.cresques.io.GeoRasterFile#reProject(org.cresques.cts.ICoordTrans)

	 */

	public void reProject(ICoordTrans rp) {

		// TODO Auto-generated method stub

	}

	/* (non-Javadoc)

	 * @see org.cresques.io.GeoRasterFile#setTransparency(boolean)

	 * /

	public void setTransparency(boolean t) {

		// TODO Auto-generated method stub

	}

	public void setTransparency(int t) {

		// TODO Auto-generated method stub

	}*/

	/* (non-Javadoc)

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

	 */

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

		double dFileAspect, dWindowAspect;

		int line, pRGBArray[] = null;

		Image image = null;

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

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

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

		if (dFileAspect > dWindowAspect) {

		  height =(int)((double)width/dFileAspect);

		} else {

		  width = (int)((double)height*dFileAspect);

		}

		// Set the view

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

			width, height);

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

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

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

		try {

			//int nLin = height % BAND_HEIGHT;

			file.setAlpha(getAlpha());

			setBand(RED_BAND,   rBandNr);

			setBand(GREEN_BAND, gBandNr);

			setBand(BLUE_BAND,  bBandNr);

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

				//int bandH = Math.min(BAND_HEIGHT, height-line);

				//file.readBandRGBA(bandH, BAND_HEIGHT, pRGBArray);

				file.readLineRGBA(pRGBArray);

				setRGBLine((BufferedImage) image, 0, line, width, 1/*bandH*/, pRGBArray, 0, width);

			}

		} catch (Exception e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

		return image;

	}

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

		double dFileAspect, dWindowAspect;

		int line, pRGBArray[][] = null;

		WritableRaster raster = null;

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

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

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

		if (dFileAspect > dWindowAspect) {

		  height =(int)((double)width/dFileAspect);

		} else {

		  width = (int)((double)height*dFileAspect);

		}

		// Set the view

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

			width, height);

		raster = Raster.createBandedRaster(getDataType(), width, height, 3, new Point(0,0));

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

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

		try {

			//int nLin = height % BAND_HEIGHT;

			file.setAlpha(getAlpha());

			setBand(RED_BAND,   rBandNr);

			setBand(GREEN_BAND, gBandNr);

			setBand(BLUE_BAND,  bBandNr);

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

				//int bandH = Math.min(BAND_HEIGHT, height-line);

				//file.readBandRGBA(bandH, BAND_HEIGHT, pRGBArray);

				file.readLine(pRGBArray);

				raster.setSamples(0, line, width, 1, 0, pRGBArray[0]);

				raster.setSamples(0, line, width, 1, 1, pRGBArray[1]);

				raster.setSamples(0, line, width, 1, 2, pRGBArray[2]);

//				setRGBLine((BufferedImage) image, 0, line, width, 1/*bandH*/, pRGBArray, 0, width);

			}

		} catch (Exception e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

		return raster;

	}

	protected void setRGBLine(BufferedImage image, int startX, int startY, int w, int h, int[] rgbArray, 

			 int offset, int scansize) {

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

	}

	protected void setRGBLine(BufferedImage image, int startX, int startY, int w, int h, int[] rgbArray, 

			 int offset, int scansize, int flags) {

		int [] line = new int[rgbArray.length]; 

		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] & 0x0000ffff) | (rgbArray[i] & 0xffff0000);

		else if (flags == GeoRasterFile.GREEN_BAND)

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

				line[i] = (line[i] & 0x00ff00ff) | (rgbArray[i] & 0xff00ff00);

		else if (flags == GeoRasterFile.BLUE_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);

	}

	private void showOnOpen() {

  		// Report en la apertura (quitar)

  		System.out.println("Fichero GDAL '"+getName()+"' abierto.");

  		System.out.println("Version = "+file.version);

  		System.out.println("   Size = ("+file.width+","+file.height+")");

  		try {

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

		} catch (GdalException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

  		file.pintaInfo();

  		file.pintaPaleta();

	}

	/* (non-Javadoc)

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

	 */

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

		double dFileAspect, dWindowAspect;

		int line, pRGBArray[] = null;

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

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

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

		if (dFileAspect > dWindowAspect) {

		  height =(int)((double)width/dFileAspect);

		} else {

		  width = (int)((double)height*dFileAspect);

		}

		// Set the view

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

			width, height);

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

		try {

			//int nLin = height % BAND_HEIGHT;

			setBand(RED_BAND,   rBandNr);

			setBand(GREEN_BAND, gBandNr);

			setBand(BLUE_BAND,  bBandNr);

			file.setAlpha(getAlpha());

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

				//int bandH = Math.min(BAND_HEIGHT, height-line);

				//file.readBandRGBA(bandH, BAND_HEIGHT, pRGBArray);

				file.readLineRGBA(pRGBArray);

				setRGBLine((BufferedImage) img, 0, line, width, 1/*bandH*/, pRGBArray, 0, width, flags);

			}

		} catch (Exception e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

		return img;

	}

	/* (non-Javadoc)

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

	 */

	public Object getData(int x, int y, int band) {

		// TODO Auto-generated method stub

		return null;

	}

}

/*

 * Created on 16-dic-2004

 */

package org.cresques.io;

import org.cresques.px.Extent;

/**

 * Consulta el servidor de mapas de Alia de Sevilla.

 * (hecho por tracasa)

 * ver http://sevilla.tracasa.es/navegar/?lang=cas

 * Implementa una llamada del tipo:

 * http://sevilla.tracasa.es/navegar/asp/mime.aspx?

 *   TIPO=0&Xmin=225916.666666667&Ymin=4133000&

 *   Xmax=254583.333333333&Ymax=4150000&PixAncho=430&PixAlto=255

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

 */

public class AliaMapServerClient extends MapServerClient {

	public final static int CAPA_ORTOFOTO_2000  = 0;

	public final static int CAPA_ORTOFOTO_7500  = 1;

	public final static int CAPA_GUIA_URBANA	= 2;

	public final static String [] nomCapas = {

		"Ortofoto 1:2.000 (Noviembre 2.001)",

		"Ortofoto 1:7.500 (Enero 1.999)",

		"Guía Urbana (Marzo 1.999)"};

	private int capas		= CAPA_ORTOFOTO_2000;

	public AliaMapServerClient(String serverName) {

		super(serverName);

		urlBase = "http://sevilla.tracasa.es/navegar/asp/mime.aspx";

		setExtent(new Extent(225916.67, 4150000.0, 254583.3, 4133000.0));

		setMaxViewSize(650,400);

	}

	public void setCapas(int c) { capas = c; }

	public String getUrl() {

		String url = urlBase+"?"+

			"TIPO="+capas+"&"+

			"Xmin="+xMin+"&"+ "Ymin="+yMax+"&"+

			"Xmax="+xMax+"&"+ "Ymax="+yMin+"&"+

			"PixAncho="+ancho+"&"+"PixAlto="+alto;;

		return url;

	}

}

/*

 * Created on 16-dic-2004

 */

package org.cresques.io;

import org.cresques.px.Extent;

/**

 * Consulta el servidor de mapas de Sitna de Navarra.

 * (hecho por tracasa)

 * ver http://sitna.tracasa.es/navegar/?lang=cas

 * Implementa una llamada del tipo:

 * http://sitna.tracasa.es/navegar/asp/mime.aspx?

 *   TIPO=0&Xmin=450321.428571429&Ymin=628000&

 *   Xmax=780678.571428571&Ymax=813000&PixAncho=500&PixAlto=280

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

 */

public class SitnaMapServerClient extends MapServerClient {

	public final static int CAPA_ORTOCOLOR_2K_03  = 0;

	public final static int CAPA_ORTOFOTO_7500  = 1;

	public final static int CAPA_GUIA_URBANA	= 2;

	public final static String [] nomCapas = {

		"Ortofoto Color 1:5.000 - A?o 2.003",

		"Ortofoto Color 1:5.000 - A?os 98-00",

		"Ortofoto Color 1:2.000",

		"Ortofoto Color 1:1.000",

		"Mapa de Relieve Color (25 m)",

		"Mapa de Relieve B/N (7,5 m)",

		"Cartograf?a Topogr?fica",

		"Catastro",

		"Catastro + Ortofoto 1:5.000",

		"LANDSAT-5 27/05/91 Falso Color",

		"LANDSAT-5 27/05/91 Pseudocolor",

		"LANDSAT-5 12/06/91 Falso Color",

		"LANDSAT-5 12/06/91 Pseudocolor",

		"LANDSAT-5 15/08/91 Falso Color",

		"LANDSAT-5 15/08/91 Pseudocolor",

		"LANDSAT-5 13/05/92 Falso Color",

		"LANDSAT-5 13/05/92 Pseudocolor",

		"LANDSAT-5 17/08/92 Falso Color",

		"LANDSAT-5 17/08/92 Pseudocolor",

		"LANDSAT-5 01/06/93 Falso Color",

		"LANDSAT-5 01/06/93 Pseudocolor",

		"LANDSAT-5 20/08/93 Falso Color",

		"LANDSAT-5 20/08/93 Pseudocolor",

		"LANDSAT-5 03/05/94 Falso Color",

		"LANDSAT-5 03/05/94 Pseudocolor",

		"LANDSAT-5 23/08/94 Falso Color",

		"LANDSAT-5 23/08/94 Pseudocolor",

		"LANDSAT-5 04/04/95 Falso Color",

		"LANDSAT-5 04/04/95 Pseudocolor",

		"LANDSAT-5 10/08/95 Falso Color",

		"LANDSAT-5 10/08/95 Pseudocolor",

		"LANDSAT-5 24/05/96 Falso Color",

		"LANDSAT-5 24/05/96 Pseudocolor",

		"LANDSAT-5 11/07/96 Falso Color",

		"LANDSAT-5 11/07/96 Pseudocolor",

		"LANDSAT-5 08/03/97 Falso Color",

		"LANDSAT-5 08/03/97 Pseudocolor",

		"LANDSAT-5 27/03/98 Falso Color",

		"LANDSAT-5 27/03/98 Pseudocolor",

		"LANDSAT-5 30/03/99 Falso Color",

		"LANDSAT-5 30/03/99 Pseudocolor",

		"LANDSAT-5 21/08/99 Falso Color",

		"LANDSAT-5 21/08/99 Pseudocolor",

		"LANDSAT-7 08/03/00 Falso Color",

		"LANDSAT-7 08/03/00 Pseudocolor",

		"LANDSAT-5 07/08/00 Falso Color",

		"LANDSAT-5 07/08/00 Pseudocolor",

		"LANDSAT-7 28/04/01 Falso Color",

		"LANDSAT-7 28/04/01 Pseudocolor",

		"LANDSAT-5 31/07/03 Falso Color",

		"LANDSAT-5 31/07/03 Pseudocolor",

		"IRS-1C 13/04/97 Falso Color",

		"IRS-1C 13/04/97 Pseudocolor",

		"IRS-1C 04/09/97  Falso Color",

		"IRS-1C 04/09/97 Pseudocolor",

		"IRS-1C 19/06/98 Falso Color",

		"IRS-1C 19/06/98 Pseudocolor",

		"IRS-1C 06/08/98 Falso Color",

		"IRS-1C 06/08/98 Pseudocolor",

		"IRS-1D 15/06/99 Falso Color",

		"IRS-1D 15/06/99 Pseudocolor",

		"IRS-1C 08/06/00 Falso Color",

		"IRS-1C 08/06/00 Pseudocolor",

		"IRS-1D 23/08/01 Falso Color",

		"IRS-1D 23/08/01 Pseudocolor",

		"IRS-1D 05/04/02 Falso Color",

		"IRS-1D 05/04/02 Pseudocolor",

		"IRS-1C 26/09/02 Falso Color",

		"IRS-1C 26/09/02 Pseudocolor",

		"IRS-1D 10/05/03 Falso Color",

		"IRS-1D 10/05/03 Pseudocolor"};

	private int capas		= CAPA_ORTOCOLOR_2K_03;

	public SitnaMapServerClient(String serverName) {

		super(serverName);

		urlBase = "http://sitna.tracasa.es/navegar/asp/mime.aspx";

		setExtent(new Extent(450321.428571429, 813000, 780678.571428571, 628000));

		setMaxViewSize(700,420);

	}

	public void setCapas(int c) { capas = c; }

	public String getUrl() {

		String url = urlBase+"?"+

			"TIPO="+capas+"&"+

			"Xmin="+xMin+"&"+ "Ymin="+yMax+"&"+

			"Xmax="+xMax+"&"+ "Ymax="+yMin+"&"+

			"PixAncho="+ancho+"&"+"PixAlto="+alto;;

		return url;

	}

}

	protected byte [] bandR, bandG, bandB;

  				if (rBandNr == 1) bandR = buffer.buf1;

  				else if (rBandNr == 2) bandR = buffer.buf2;

  				else if (rBandNr == 3) bandR = buffer.buf3;

  				if (gBandNr == 1) bandG = buffer.buf1;

  				else if (gBandNr == 2) bandG = buffer.buf2;

  				else if (gBandNr == 3) bandG = buffer.buf3;

  				if (bBandNr == 1) bandB = buffer.buf1;

  				else if (bBandNr == 2) bandB = buffer.buf2;

  				else if (bBandNr == 3) bandB = buffer.buf3;

	  			if (eColorSpace==LTIColorSpace.LTI_COLORSPACE_RGB)

	  							line[kd] = alpha + ((0xff & bandR[k])<<16) + (( 0xff & bandG[k])<<8) + (0xff & bandB[k]);

			if ( bandCount > 2) {

				setBand(RED_BAND,   0);

				setBand(GREEN_BAND, 1);

				setBand(BLUE_BAND,  2);

			} else

				setBand(RED_BAND|GREEN_BAND|BLUE_BAND, 0);

		if ((flag & GeoRasterFile.RED_BAND) == GeoRasterFile.RED_BAND) file.rBandNr = bandNr+1;

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

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

			file.setAlpha(getAlpha());

			setBand(RED_BAND,   rBandNr);

			setBand(GREEN_BAND, gBandNr);

			setBand(BLUE_BAND,  bBandNr);

import java.awt.image.DataBuffer;

	private int dataType = DataBuffer.TYPE_BYTE;

	public int getAlpha() {

		if (tFilter == null)

			return 255;

		return tFilter.getAlpha();

	}

	/**

	 * @return Returns the dataType.

	 */

	public int getDataType() {

		return dataType;

	}

	/**

	 * @param dataType The dataType to set.

	 */

	public void setDataType(int dataType) {

		this.dataType = dataType;

	}

import java.awt.geom.GeneralPath;

import java.awt.geom.Point2D;

import java.awt.image.BufferedImage;

import java.awt.image.DataBuffer;

import java.awt.image.Raster;

import org.cresques.io.GdalFile;

		if (geoFile.getBandCount()>2) {

			setBand(GeoRasterFile.RED_BAND, 0);

			setBand(GeoRasterFile.GREEN_BAND, 1);

			setBand(GeoRasterFile.BLUE_BAND, 2);

		} else {

			setBand(GeoRasterFile.RED_BAND|GeoRasterFile.GREEN_BAND|GeoRasterFile.BLUE_BAND, 0);

		}

		if (fnames.length >= 3 || geoFile.getBandCount()>2) {

			setBand(GeoRasterFile.RED_BAND, 0);

			setBand(GeoRasterFile.GREEN_BAND, 1);

			setBand(GeoRasterFile.BLUE_BAND, 2);

			setBand(GeoRasterFile.RED_BAND|GeoRasterFile.GREEN_BAND|GeoRasterFile.BLUE_BAND, 0);

		if (geoFile.getBandCount()>2) {

			setBand(GeoRasterFile.RED_BAND, 0);

			setBand(GeoRasterFile.GREEN_BAND, 1);

			setBand(GeoRasterFile.BLUE_BAND, 2);

		} else {

			setBand(GeoRasterFile.RED_BAND|GeoRasterFile.GREEN_BAND|GeoRasterFile.BLUE_BAND, 0);

		}

	public GeoRasterFile [] getFiles() {

		if (colorBand != null)

			return colorBand;

		GeoRasterFile [] ret = {geoFile};

		return ret;

	}

		if (colorBand != null) {

			if ((flag & GeoRasterFile.RED_BAND) == GeoRasterFile.RED_BAND)

				rBand = nBand;

			else if ((flag & GeoRasterFile.GREEN_BAND) == GeoRasterFile.GREEN_BAND)

				gBand = nBand;

			else if ((flag & GeoRasterFile.BLUE_BAND) == GeoRasterFile.BLUE_BAND)

				bBand = nBand;

		} else

			geoFile.setBand(flag, nBand);

		if (colorBand != null) {

			if (flag == GeoRasterFile.RED_BAND)

				return rBand;

			else if (flag == GeoRasterFile.GREEN_BAND)

				return gBand;

			else if (flag == GeoRasterFile.BLUE_BAND)

				return bBand;

			else