Application: gvSIG desktop » gvSIG raster

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

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

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

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

import org.gvsig.fmap.dal.coverage.store.props.HistogramComputer;

import org.gvsig.fmap.dal.exception.DataException;

import org.gvsig.fmap.dal.exception.ValidateDataParametersException;

			ArrayList comboSource = getHistogramPanel().getComboSource();

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

				String name = (String) ((ArrayList) comboSource.get(i)).get(1);

					((ArrayList) comboSource.get(i)).remove(0);

				}

					((ArrayList) comboSource.get(i)).remove(0);

					((ArrayList) comboSource.get(i)).add(0, ((FLyrRaster) lyr).getDataStore().getHistogramComputer());

				}

 * <code>HistogramDialog</code>. Creaci?n de la ventana de histograma para gvSIG.

 *

 * @version 20/03/2007

 */

public class HistogramDialog extends JPanel implements IWindow, ButtonsPanelListener {

	private static final long serialVersionUID = 7362459094802955247L;

 * 

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

 */

public interface FLyrRaster extends FLayer {

	/**

	 * Gets the uniform resource identifier

	 * @param t

	 */

	public void setNoDataTransparent(boolean t);

}

	static private IConfiguration   configuration                  = new DefaultLayerConfiguration();

	protected int                   zoomLevel                      =  1;

	public boolean                  recalcLevel                    = true;

	private static GeometryManager  geomManager   	               = GeometryLocator.getGeometryManager();

	private static final Logger     logger                         = LoggerFactory.getLogger(DefaultFLyrRaster.class);

	protected FileUtils             fileUtil                       = RasterLocator.getManager().getFileUtils();

	protected RasterUtils           rasterUtil                     = RasterLocator.getManager().getRasterUtils();

	protected CRSUtils              crsUtil                        = RasterLocator.getManager().getCRSUtils();

	protected MathUtils             mathUtil                       = RasterLocator.getManager().getMathUtils();

	public class RasterTaskStatus extends BaseTaskStatus {

		Cancellable c = null;

	 * 		com.iver.utiles.swing.threads.Cancellable)

	 */

	public void draw(BufferedImage image, Graphics2D g, ViewPort vp, Cancellable cancel, double scale) throws ReadException {

		TaskEventManager task = rManager.getRasterTask();

		task.setEvent(null);

	 * @param y Posici?n y

	 * @return valor de pixel

	 */

	/*

	 * (non-Javadoc)

		if (/*!isOpen() ||*/ !isVisible() || !isWithinScale(scale))

			return;

		if (!mustTilePrint)

			draw(null, g, viewPort, cancel,scale);

		else {

			Tiling tiles = new Tiling(maxTilePrintWidth, maxTilePrintHeight, g.getClipBounds());

			tiles.setAffineTransform((AffineTransform) viewPort.getAffineTransform().clone());

				// Parte que dibuja

				try {

					ViewPort vp = tiles.getTileViewPort(viewPort, tileNr);

					draw(null, g, vp, cancel, scale);

				} catch (NoninvertibleTransformException e) {

					throw new ReadException("Error en la transformaci?n.", e);

				}

		}

	}

	public void _print(Graphics2D g, ViewPort viewPort, Cancellable cancel,

			}

			if (this.isRemoveRasterFlag()) {

				if (getDataStore() != null) {

					String[] files = getFileName().clone();

			pReal.getY() < this.getMaxY()) {

			px = transformPoint(pReal);

		}

		info.addField("File", normalizeAsXMLTag(getName()), 0);

		info.addField("View Point", "[" + p.getX() + " , " + p.getY() + "]", 1);

		info.addField("World Point", "[" + mathUtil.format(pReal.getX(), 3) + " , " + mathUtil.format(pReal.getY(), 3) + "]", 2);

			info.addField("Pixel Point", "Out", 3);

		else

			info.addField("Pixel Point", "[" + (int) px.getX() + ",  " + (int) px.getY() + "]", 3);

		double[] cmyk = conv.RGBtoCMYK(rgb[1] & 0xff, rgb[2] & 0xff, rgb[3] & 0xff, 1D);

		info.addField("CMYK", "[" + mathUtil.format(cmyk[0], 4) + ",  " + mathUtil.format(cmyk[1], 4) + ",  " + mathUtil.format(cmyk[2], 4) + ",  " + mathUtil.format(cmyk[3], 4) + "]", 5);

		double[] hsl = conv.RGBtoHSL(rgb[1] & 0xff, rgb[2] & 0xff, rgb[3] & 0xff);

		hsl[0] = (int)(255.0 * hsl[0] / 360.0 + 0.5);

		hsl[2] = (int) (hsl[2] * 255. + 0.5);

		hsl[1] = (int) (hsl[1] * 255. + 0.5);

		String data = "[";

		try {

			if (px != null) {

	}

	/**

	 *

	 * @return ROIs asociadas a la capa raster.

	 */

		} 

	}

	/*

	 * (non-Javadoc)

	 * @see java.lang.Object#finalize()

		}

		strategy                       = null;

		configuration                  = null;

		fileUtil                       = null;

		rasterUtil                     = null;

		crsUtil                        = null;

		render                         = null;

	}

}

	public int getLastAlphaBandNumber();

	/**

	 * Informa de si el raster tiene tabla de color asociada o no.

	 * @return true si tiene tabla de color y false si no la tiene.

	 */

	 * @throws InvalidSetViewException

	 * @throws ArrayIndexOutOfBoundsException

	 */

		throws RasterDriverException, InvalidSetViewException, ProcessInterruptedException;

	/**

	public void drawThread(Graphics2D g, ViewPortData vp);

	/**

	 * Sets the view structures to draw

	 * @param g

	 * @param vp

	 */

	public void refreshHistogram();

}

import org.gvsig.raster.cache.tile.exception.TileGettingException;

import org.gvsig.raster.cache.tile.provider.TileListener;

import org.gvsig.raster.impl.DefaultRasterManager;

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

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

import org.gvsig.raster.impl.grid.GridImpl;

	 */

	private int[]            renderBands              = { 0, 1, 2 };

	/**

	 * Ultima transparencia aplicada en la visualizaci?n que es obtenida desde el

	 * grid

	 */

	private RasterFilterList filterList               = null;

	/**

	 * Ancho y alto del objeto Image en una petici?n de dibujado a un raster

	 * raster

		}

		if (dataStore != null) {

			// Asignamos la banda de transparencia si existe esta

			RasterQuery query = DefaultRasterManager.getInstance().createQuery();

			query.setSupersamplingLoadingBuffer(false); // Desactivamos el supersampleo en la carga del buffer.

			query.setDrawableBands(getRenderBands());

			query.setFrameWidth(0);

			query.setAreaOfInterest(vp.getExtent(), 

					(int)Math.round(vp.getWidth()), 

					(int)Math.round(vp.getHeight()), this);

			return;

	}

		throws RasterDriverException, InvalidSetViewException, ProcessInterruptedException {

		lastGraphics = g;

		lastViewPortData = vp;

		if(util.isOutside(vp.getExtent(), dataStore.getExtent())) {

			endReading();

		}

		Extent adjustedRotedRequest = request(vp, dataStore);

		if ((widthImage <= 0) || (heightImage <= 0)) {

			endReading();

		}

		if (dataStore == null) 

		Buffer buf = null; 

		double[] step = null;

					(int)Math.round(heightImage), 

					vp.getTime(), 

					renderBands, 

		} else {

			// Asignamos la banda de transparencia si existe esta

			RasterQuery query = DefaultRasterManager.getInstance().createQuery();

			query.setTaskStatus(taskStatus);

			query.setTime(vp.getTime());

			query.setSupersamplingLoadingBuffer(false); // Desactivamos el supersampleo en la carga del buffer.

			query.setAreaOfInterest(adjustedRotedRequest, (int)Math.round(widthImage), (int)Math.round(heightImage));

			}

			query.setAreaOfInterest(adjustedRotedRequest, (int)Math.round(widthImage), (int)Math.round(heightImage));

			query.setDrawableBands(getRenderBands());

			buf = dataStore.query(query);

			query.setSupersamplingLoadingBuffer(true);

			step = dataStore.getStep();

		}

			init();

		}

	}

	/*

	 * (non-Javadoc)

	 * @see org.gvsig.fmap.dal.coverage.grid.render.Render#draw(java.awt.Graphics2D, org.cresques.geo.ViewPortData)

	 */

		throws RasterDriverException, InvalidSetViewException, ProcessInterruptedException {

		grid = bufferPreprocessing(buf, lastTransparency);

		//Buffer filtrado para renderizar

		drawer.setSupersamplingOn(step); // Desplazamiento para supersampleo

		drawer.setOutputSize((int)Math.round(widthImage), (int)Math.round(heightImage)); // Ancho y alto del buffer

		drawer.setLastTransparency(getLastTransparency());

		Image geoImage = drawer.drawBufferOverImageObject(); // Acci?n de renderizado

		// Borramos el buffer de transparencia para que siempre se tenga que regenerar.

		//En el caso de no tenga rotaci?n y el tama?o de pixel sea positivo en X y negativo en Y no aplicamos ninguna

		//transformaci?n. Esto no es necesario hacerlo, sin ello se visualiza igual. Unicamente se hace porque de esta

			Point2D lastGraphicOffset = new Point2D.Double(adjustedRotedRequest.getULX(), adjustedRotedRequest.getULY());

			((DefaultViewPortData)vp).mat.transform(lastGraphicOffset, lastGraphicOffset);

			g.drawImage(geoImage, (int) Math.round(lastGraphicOffset.getX()), (int) Math.round(lastGraphicOffset.getY()), null);

		}

		/*

			LoggerFactory.getLogger(getClass()).debug("Transformation error", e);

		}

		// long t2 = new Date().getTime();

		// System.out.println("Renderizando Raster: " + ((t2 - t1) / 1000D) + ", secs.");

	}

			} catch (ProcessInterruptedException e2) {

				return;

			}

			lastTransparency.setAlphaBand(null);

			AffineTransform at = new AffineTransform();

		return out;

	}

	public Transparency getLastTransparency() {

		return lastTransparency;

	}

	/*

	 * (non-Javadoc)

	 * @see org.gvsig.fmap.dal.coverage.grid.render.Render#setFilterList(org.gvsig.fmap.dal.coverage.grid.RasterFilterList)

	 */

	public void setFilterList(RasterFilterList filterList) {

			grid.dispose();

		if (getFilterList() != null)

			getFilterList().dispose();

		try {

			finalize();

		} catch (Throwable e) {

        drawer                   = null;

		lastTransparency         = null;

		filterList               = null;

        ulPxRequest              = null;

        lrPxRequest              = null;

		lastGraphics             = null;

	public void dispose() {

		if(rasterBuf != null)

			rasterBuf.dispose();

		try {

			finalize();

		} catch (Throwable e) {

	protected void finalize() throws Throwable {

		step        = null;

		rasterBuf   = null;

		super.finalize();

	}

	protected int			            progressHistogram    = 0;

	private double[]                    limits               = null;

	private boolean                     refresh              = false;

	/**

	 * Constructor

		this.buffer = buffer;

	}

	/*

	 * (non-Javadoc)

	 * @see org.gvsig.fmap.dal.coverage.store.props.HistogramComputer#getBufferHistogram()

	private int            offset = 1;

	private int            bufWidth, bufHeight;

	/**

	 * Constructor. Asigna las variables de inicializaci?n y crea la estructura de 

	 */

	public RasterReadOnlyBuffer(int dataType, int width, int height, int nBand) {

		this.dataType = dataType;

	}

	/**

		this.maxX = maxX;

		this.maxY = maxY;

			bufWidth = width;

			bufHeight = height;

		init();

	}

	 */

	private void init() {

		//Calculo de los bits por p?gina

		while(h > 1) {

			h >>= 1;

			bitsPag ++;

		}

		//Calculo del desplazamiento

		//Creamos el bufferFactory asociado a la carga de p?ginas. Cada p?gina se fuerza

		//a que sea cargada en memoria.

		query.setMemoryBuffer(true);

		query.setDrawableBands(bandList.getDrawableBands());

		//double aux = ((double)height / (double)RasterLibrary.blockHeight);

		//heightLastPage = (int)((aux - (int)aux) * RasterLibrary.blockHeight);

		else {

		}

		drawableBands = new int[datasource.getBandCount()];

			drawableBands[i] = i;

	}

	 */

	public Extent[] calcExtentPages(DefaultRasterProvider dataset, int nPages) {

		Extent datasetExtent = dataset.getExtent();

		Extent[] ext = new Extent[nPages];

		double minX = datasetExtent.getMin().getX();

	private Buffer queryPage(int pag) throws RasterDriverException, ProcessInterruptedException, InvalidSetViewException {

		int w = maxX - minX + 1;

		int h = maxY - minY + 1;

		if(bufWidth == 0 && bufHeight == 0) {

			if(pag == nPages - 1)

			else 

		} else {

			if(pag == nPages - 1) {

			} else 

		}

		return datasource.query(query);

	}

		//TODO: FUNCIONALIDAD: Falta asignar lo necesario para la transparencia por selecci?n

		this.transparencyRanges = (ArrayList<TransparencyRange>) t.getTransparencyRange().clone();

		this.opacity = t.getOpacity();

		this.alphaBandNumber = t.alphaBandNumber;

		this.noData = (NoData)t.getNoData().clone();

		this.transparencyActive = t.transparencyActive;

	}

	 * @see org.gvsig.fmap.dal.coverage.store.props.Transparency#activeTransparency()

	 */

	public void activeTransparency() {

			(transparencyRanges.size() > 0) ||

			(opacity != 0xff))

			transparencyActive = true;

			mask.dispose();

		if (originalData != null)

			originalData.dispose();

	}

	/*

	private RasterProvider          provider  = null;

	private int                     percent   = 0;

	private boolean                 refresh   = false;

	/**

	 * Constructor

	public SimpleProviderHistogramComputer(RasterProvider provider) {

		this.provider = provider;

	}

	/**

	 * Obtiene el minimo valor de las estadisticas de un histograma.

		for (int block = 0; block < provider.getHeight(); block += h) {

			Object buf = null;

			try {

			} catch (InvalidSetViewException e) {

				// La vista se asigna autom?ticamente

			}

			if ((block + hB) > provider.getHeight())

				hB = Math.abs(((int)provider.getHeight()) - block);

			if(buf instanceof Buffer) {

				loadHistogramFromBuffer((Buffer)buf);

 */

package org.gvsig.raster.impl.store.properties;

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

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

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

	private boolean                     refresh     = false;

	private int                         percent     = 0;

	private int                         files       = 0;

	/**

	 * Constructor

		files = provider.getInternalProviderCount();

	}

	/**

	 * Obtiene el histograma. Pregunta a todos los datasets que forman el multidataset

	 * por su histograma y luego los fusiona.

					return null;

				try {

				} catch (FileNotOpenException e) {

					throw new HistogramException("");

				} catch (RasterDriverException e) {

		}

	}

}

		double[][] res = new double[5][];

		double[] average = new double[bandCount]; // Valor de pixel medio (Media)

		double[] middle = new double[bandCount]; // Mediana

		double[] nPixelsBand = new double[bandCount];// N?mero de pixels por banda

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

			if (bands[iBand]) {

				for (int i = beginPos; i <= endPos; i++) {

					// Calculo del n?mero de pixeles

					nPixelsBand[showBandCounter] += histogram[iBand][i].getValue();

					average[showBandCounter] += histogram[iBand][i].getValue() * ((histogram[iBand][i].getMax() + histogram[iBand][i].getMin())/2);

				}

				// Calculo de la media

				// Calculo de mediana

				double stopPos = nPixelsBand[showBandCounter] / 2;

				int i = beginPos;

				for (i = beginPos; aux <= stopPos; i++) {

					if (i >= histogram[iBand].length)

					}

				}

				middle[showBandCounter] = (histogram[iBand][i - 1].getMax() + histogram[iBand][i - 1].getMin()) / 2;

				showBandCounter++;

			}

		}

			if(lastblock > 0)

				nBlocks ++;

			int initYSrc = inputWindow[1];

			int initYBuffer = 0;

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

				if(lastblock > 0 && i == (nBlocks - 1)) {

					int[] newStepBuffer = new int[]{0, initYBuffer, stepBuffer[2], initYBuffer + lastblock};

					readData(buf, 

							bandList, 

							bufWidth, lastblock, 0, 0, newStepBuffer);

				} else {

					int[] newStepBuffer = new int[]{0, initYBuffer, stepBuffer[2], initYBuffer + buf.getBlockHeight()};

					readData(buf, 

							bandList, 

							bufWidth, buf.getBlockHeight(), 0, 0, newStepBuffer);

					initYBuffer += buf.getBlockHeight();

				}

			}