Application: gvSIG desktop

	/**

	 * Obtiene una ventana de datos de la imagen a partir de coordenadas reales. 

	 * No aplica supersampleo ni subsampleo sino que devuelve una matriz de igual tama?o a los

	 * pixeles de disco. 

	 * @param x Posici?n X superior izquierda

	 * @param y Posici?n Y superior izquierda

	 * @param w Ancho en coordenadas reales

	 * @param h Alto en coordenadas reales

	 * @param adjustToExtent Flag que dice si el extent solicitado debe ajustarse al extent del raster o no.

	 * @param bandList

	 * @return Buffer de datos

	 */

	public IBuffer getWindowRaster(double ulx, double uly, double lrx, double lry) throws InvalidSetViewException {		

		Extent selectedExtent = new Extent(ulx, uly, lrx, lry);

		//Leemos pixels completos aunque el valor obtenido sea decimal. Esto se consigue redondeando

		//por arriba el m?s alto y por abajo el menor y luego restandolos

		Point2D p1 = ((RasterDataset)files.get(0)).worldToRaster(new Point2D.Double(ulx, uly));

		Point2D p2 = ((RasterDataset)files.get(0)).worldToRaster(new Point2D.Double(lrx, lry));

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

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

		//Ajustamos por si nos hemos salido del raster

		while(((int)(p1.getX() + width)) > getWidth()[0])

			width --;

		while(((int)(p1.getY() + height)) > getHeight()[0])

			height --;

		if(p1.getX() < 0)

			p1.setLocation(0, p1.getY());

		if(p1.getY() < 0)

			p1.setLocation(p1.getX(), 0);

		if(p2.getX() > getWidth()[0])

			p2.setLocation(getWidth()[0], p2.getY());

		if(p1.getY() > getHeight()[0])

			p1.setLocation(p2.getX(), getHeight()[0]);

		int mallocNBands = 0;

		if(bandList.getDrawableBands() != null)

			mallocNBands = bandList.getDrawableBands().length;

		else

			mallocNBands = bandList.getDrawableBandsCount();

		IBuffer raster = RasterBuffer.getBuffer(getDataType()[0], width, height, mallocNBands, false);

		//Caso en el que un buffer no se puedan conmutar las bandas (t?picamente RasterReadOnlyHugeBuffer)

		if(!raster.isBandSwitchable()) {

			RasterBuffer rb = RasterBuffer.getBuffer(getDataType()[0], width, height, getDataset(0).getBandCount(), false);

			if(rb instanceof RasterReadOnlyHugeBuffer) {

				try {

					((RasterReadOnlyHugeBuffer)rb).setBufferParams(getDataset(0).getFName(), selectedExtent);

				} catch (FileNotExistsException e) {

					//Esto no debe darse ya que se comprueba al hacer el open

				} catch (NotSupportedExtensionException e) {

					//Esto no debe darse ya que se comprueba al hacer el open

				}

				return rb;

			}

		}

		for(int iBand = 0; iBand < raster.getBandCount(); iBand ++)

			raster.assignBandToNotValid(iBand);

		//Reservamos memoria para los buffers por dataset

		IBuffer[] bufferDatasets = mallocBuffersDatasets(bandList, width, height);

		//Si hemos redondeado los pixeles de la petici?n (p1 y p2) por arriba y por abajo deberemos calcular un extent mayor 

		//equivalente a los pixeles redondeados.

		Point2D wc1 = ((RasterDataset)files.get(0)).rasterToWorld(new Point2D.Double(Math.floor(p1.getX()), Math.floor(p1.getY())));

		Point2D wc2 = ((RasterDataset)files.get(0)).rasterToWorld(new Point2D.Double(Math.ceil(p2.getX()), Math.ceil(p2.getY())));

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

			bufferDatasets[i] = ((RasterDataset)files.get(i)).getWindowRaster(wc1.getX(), wc1.getY(), wc2.getX(), wc2.getY(), bandList, bufferDatasets[i]);

		//Mezclamos los buffers de cada dataset en un solo buffer

		mergeBuffers(raster, bufferDatasets, bandList);

		return raster;

	}

		double lrx = (((int)(ext.getULX() - wRaster)) == ((int)ext.getLRX())) ? (ulx - w) : (ulx + w);

		double lry = (((int)(ext.getULY() - hRaster)) == ((int)ext.getLRY())) ? (uly - h) : (uly + h);