| gvSIG desktop 1

/*

 * Created on 18-dic-2004

 */

package org.cresques.io;

import java.awt.Image;

import java.awt.image.BufferedImage;

import java.awt.image.DataBuffer;

import java.awt.image.DataBufferByte;

import java.awt.image.DataBufferShort;

import java.awt.image.DataBufferInt;

import java.awt.image.DataBufferFloat;

import java.awt.image.DataBufferDouble;

import java.io.BufferedReader;

import java.io.File;

import java.io.FileReader;

import java.io.IOException;

import java.io.InputStreamReader;

import java.io.Reader;

import java.util.StringTokenizer;

import java.util.zip.ZipEntry;

import org.cresques.cts.ICoordTrans;

import org.cresques.cts.IProjection;

import org.cresques.px.Extent;

import org.cresques.px.IObjList;

/**

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

 */

public class AscArcGridFile extends GeoRasterFile implements QueryableRaster {

	public final static byte	NODATA_BYTE   = Byte.MIN_VALUE; 

	public final static short	NODATA_SHORT  = Short.MIN_VALUE; 

	public final static int		NODATA_INT    = Integer.MIN_VALUE; 

	public final static float	NODATA_FLOAT  = Float.MIN_VALUE; 

	public final static double	NODATA_DOUBLE = Double.MIN_VALUE; 

	private String formatName = "AscArcGrid";

	private BufferedReader fi;

	private String buf = null;

	private int width = 0, height = 0;

	private int xCnt = 0, yCnt = 0;

	private BufferedImage im = null;

	private double posX, posY, cellSize;

	private Number noData = null;

	private int dataType = TYPE_UNDEFINED;

	private DataBuffer data;

	private boolean toCorner = true;

	/**

	 * Lector de ArcGrid ascii.

	 * 

	 * @param proj

	 * @param name

	 */

	public AscArcGridFile(IProjection proj, String name) {

		super(proj, name);

	}

	/* (non-Javadoc)

	 * @see org.cresques.io.GeoFile#load()

	 */

	public GeoFile load() {

		System.out.println(formatName+": Intento de carga de '" + getName() + "'.");

		long t1, t2;

		try {

			if (FileFolder.isUrl(getName())) {

				ZipFileFolder zFolder = new ZipFileFolder(getName());

				ZipEntry ze = zFolder.getZipEntry(getName());

				setFileSize(ze.getSize());

				// El primero calcula el tama?o, el segundo la imagen

				t1 = getTime();

				load(new InputStreamReader(zFolder.getInputStream(ze)));

				t2 = getTime();

				System.out.println(formatName+": Primer load: " + (t2-t1)/1000 + " seg.");

				/*t1 = getTime();

				load(new InputStreamReader(zFolder.getInputStream(ze)));

				t2 = getTime();

				System.out.println(formatName+": Segundo load: " + (t2-t1)/1000 + " seg.");*/

			} else {

				File f = new File(getName());

				setFileSize(f.length());

				// El primero calcula el tama?o, el segundo la imagen

				t1 = getTime();

				load(new FileReader(getName()));

				t2 = getTime();

				System.out.println(formatName+": Primer load: " + (t2-t1)/1000 + " seg.");

				/*t1 = getTime();

				load(new FileReader(getName()));

				t2 = getTime();

				System.out.println(formatName+": Segundo load: " + (t2-t1)/1000 + " seg.");*/

			}

			//v = new Extent(getExtent());

			return this;

		} catch (Exception e) {

			System.err.println(formatName+": ERROR: "+lineCnt+" lineas leidas ("+bytesReaded+"/"+getFileSize()+") bytes.");

			//if (e instanceof java.lang.ArrayIndexOutOfBoundsException)

			//	System.err.println("     row = "+rowAct+", col ="+colAct); 

			e.printStackTrace();

		}

		return this;

	}

	public GeoFile load(Reader fr) throws IOException {

		lineCnt = 0;

		if (im == null){

			System.out.println(formatName+": Cargando "+getName()+" ...");

		} else

			System.out.println(formatName+": Generando imagen ...");

		fi = new BufferedReader(fr);

		while ((buf = fi.readLine()) != null) {

			bytesReaded += buf.length()+1;

			parseLine(buf);

			lineCnt++;

		}

		fi.close();

		System.out.println(formatName+": Fichero cargado; "+lineCnt+" l?neas procesadas.\n"+

				"  extent = "+getExtent());

		return this;

	}

	public void parseLine(String buf) {

		//System.out.println("lin "+lineCnt+" Length = "+buf.length());

		if (lineCnt < 6) {

			String value = buf.substring(14);

			if (buf.toLowerCase().startsWith("ncols"))

				width = Integer.parseInt(value);

			else if (buf.toLowerCase().startsWith("nrows"))

				height = Integer.parseInt(value);

			/* xllcorner e yllcorner se refiere a la esquina

			 * inferior izquierda de la celda inferior izquierda

			 * del raster.

			 * xllcenter e yllcenter se refiere al centro de la

			 * misma celda.

			 */

			else if (buf.toLowerCase().startsWith("xllcorner")) {

				posX = Double.parseDouble(value);

				toCorner = true;

			} else if (buf.toLowerCase().startsWith("yllcorner")) {

				posY = Double.parseDouble(value);

				toCorner = true;

			} else if (buf.toLowerCase().startsWith("xllcenter")) {

				posX = Double.parseDouble(value);

				toCorner = false;

			} else if (buf.toLowerCase().startsWith("yllcenter")) {

				posY = Double.parseDouble(value);

				toCorner = false;

			} else if (buf.toLowerCase().startsWith("cellsize")) {

				cellSize = Double.parseDouble(value);

				/* TODO WARNING: Este c?lculo solo estar? bien en el hemisferio

				 * norte.

				 */

				double x1 = posX, y1 = posY;

				if (!toCorner) {

					x1 -= cellSize/2.0; y1 -= cellSize/2.0; }

				double x2 = x1 + ((double)width)*cellSize,

					y2 = y1 + ((double) height)*cellSize;

				extent = new Extent (x1, y2, x2, y1);

			} else if (buf.startsWith("NODATA_value")) {

				if (dataType == TYPE_UNDEFINED)

					dataType = TYPE_FLOAT;

				switch (dataType) {

				case TYPE_FLOAT:

					noData = Float.valueOf(value);

					data = new DataBufferFloat(width, height);

					break;

				case TYPE_DOUBLE:

					noData = Double.valueOf(value);

					data = new DataBufferDouble(width, height);

					break;

				case TYPE_INT:

					noData = Integer.valueOf(value);

					data = new DataBufferInt(width, height);

					break;

				case TYPE_BYTE:

					noData = Byte.valueOf(value);

					data = new DataBufferByte(width, height);

					break;

				case TYPE_SHORT:

					noData = Short.valueOf(value);

					data = new DataBufferShort(width, height);

					break;

				}

			}

			return;

		}

		StringTokenizer st = new StringTokenizer(buf, " ");

		xCnt = 0;

		while (st.hasMoreTokens()) {

			setData(xCnt, yCnt, st.nextToken());

			//System.out.print(".");

			xCnt++;

		}

		//System.out.println(yCnt+" ");

/*		int nextIndex=0, lastIndex = 0;

		String value;

		xCnt = 0;

		while (lastIndex < buf.length()) {

			nextIndex = buf.indexOf(' ', lastIndex);

			if (nextIndex == -1) nextIndex = buf.length();

			value = buf.substring(lastIndex, nextIndex);

			System.out.print(".");

			setData(xCnt, yCnt, value);

			lastIndex = nextIndex+1;

			xCnt++;

			System.out.println(xCnt+" ");

		}*/

		yCnt++;

	}

	private void setData(int x, int y, String value) {

		switch (dataType) {

		case TYPE_BYTE:

		case TYPE_SHORT:

		case TYPE_INT:

			//data.setElem(y, x, Integer.parseInt(value));

			data.setElem(y, x, (int) Double.parseDouble(value));

			break;

		case TYPE_FLOAT:

			data.setElemFloat(y, x, Float.parseFloat(value));

			break;

		case TYPE_DOUBLE:

			data.setElemDouble(y, x, Double.parseDouble(value));

			break;

		}

	}

	/* (non-Javadoc)

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

	 */

	public void close() {

		// TODO Auto-generated method stub

	}

	/* (non-Javadoc)

	 * @see org.cresques.io.GeoFile#getObjects()

	 */

	public IObjList getObjects() {

		// TODO Auto-generated method stub

		return null;

	}

	public int getWidth() { return width;}

	public int getHeight() { return height;}

	/* (non-Javadoc)

	 * @see org.cresques.io.GeoRasterFile#setView(org.cresques.px.Extent)

	 */

	public void setView(Extent e) {

		// TODO Auto-generated method stub

	}

	/* (non-Javadoc)

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

	 */

	public Extent getView() {

		// TODO Auto-generated method stub

		return null;

	}

	/* (non-Javadoc)

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

	 */

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

		// TODO Auto-generated method stub

		return null;

	}

	public Object getNoData() { return noData; }

	/* (non-Javadoc)

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

	 */

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

		return getElem(y, x, band);

	}

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

		switch (dataType) {

		case TYPE_FLOAT:

			return new Double(data.getElemFloat(y, x));

		case TYPE_DOUBLE:

			return new Double(data.getElemDouble(y, x));

		}

		return new Integer(data.getElem(y, x));

	}

	public int getElemInt(int x, int y, int band) {

		return data.getElem(y, x);

	}

	public double getElemDouble(int x, int y, int band) {

		return data.getElem(y, x);

	}

	public float getElemFloat(int x, int y, int band) {

		return data.getElem(y, x);

	}

	/* (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.geo.Projected#reProject(org.cresques.cts.ICoordTrans)

	 */

	public void reProject(ICoordTrans rp) {

		// TODO Auto-generated method stub

	}

	/* (non-Javadoc)

	 * @see org.cresques.io.QueryableRaster#getDataType()

	 */

	public int getDataType() {

		return dataType;

	}

	/* (non-Javadoc)

	 * @see org.cresques.io.QueryableRaster#setDataType(int)

	 */

	public void setDataType(int dt) {

		dataType = dt;

	}

	/* (non-Javadoc)

	 * @see org.cresques.io.QueryableRaster#startQuerying()

	 */

	public void startQuerying() {

		// TODO Auto-generated method stub

	}

}

/*

 * Created on 18-dic-2004

 */

package org.cresques.io;

import java.awt.image.DataBuffer;

/**

 * Interface que implementan los rasters que soporten an?lisis.

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

 */

public interface QueryableRaster {

	public final static int TYPE_BYTE		= DataBuffer.TYPE_BYTE;

	public final static int TYPE_DOUBLE		= DataBuffer.TYPE_DOUBLE;

	public final static int TYPE_FLOAT		= DataBuffer.TYPE_FLOAT;

	public final static int TYPE_INT		= DataBuffer.TYPE_INT;

	public final static int TYPE_SHORT		= DataBuffer.TYPE_SHORT;

	public final static int TYPE_UNDEFINED	= DataBuffer.TYPE_UNDEFINED;

	public final static int TYPE_USHORT		= DataBuffer.TYPE_USHORT;

	public int getWidth();

	public int getHeight();

	/**

	 * Devuelve el tipo de datos del raster

	 * @see java.awt.image.DataBuffer

	 */

	public int getDataType();

	/**

	 * Asigna tipo de datos del raster

	 * @see java.awt.image.DataBuffer

	 */

	public void setDataType(int dt);

	/**

	 * Inicializa el raster para querys (acceso a las celdas de manera

	 * individual.

	 * Esto permite a la implementaci?n del raster utilizar la resoluci?n

	 * base (en caso de raster multiresoluci?n, como los ECW o SID), y

	 * la gesti?n autom?tica de las vistas (en caso de rasters que no

	 * permiten acceso individual a los pixels, como los ECW).

	 */

	public void startQuerying();

	/**

	 * Obtiene el valor NODATA del raster.

	 * Este valor indica que el elemento no existe como dato en el

	 * raster.

	 * @return

	 */

	public Object getNoData();

	/**

	 * Obtiene el valor del raster en la coordenada que se le pasa.

	 * El valor ser? Double o Integer. dependiendo del tipo de

	 * raster.

	 * @param x	coordenada X

	 * @param y coordenada Y

	 * @return

	 */

	public Object getElem(int x, int y, int band);

	/**

	 * Obtiene el valor del raster en la coordenada que se le pasa.

	 * Devuelve un int.

	 * @param x	coordenada X

	 * @param y coordenada Y

	 * @return

	 */

	public int getElemInt(int x, int y, int band);

	public float getElemFloat(int x, int y, int band);

	/**

	 * Obtiene el valor del raster en la coordenada que se le pasa.

	 * Devuelve un double.

	 * @param x	coordenada X

	 * @param y coordenada Y

	 * @return

	 */

	public double getElemDouble(int x, int y, int band);

}