Application: gvSIG desktop

/*

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

 * 

 * Copyright (C) 2004-5. 

 *

 * 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

 * of the License, or (at your option) any later version.

 *

 * This program is distributed in the hope that it will be useful,

 * but WITHOUT ANY WARRANTY; without even the implied warranty of

 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

 * GNU General Public License for more details.

 *

 * You should have received a copy of the GNU General Public License

 * along with this program; if not, write to the Free Software

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

 *

 * For more information, contact:

 * 

 * cresques@gmail.com

 */

package org.cresques.io;

import java.awt.Image;

import java.awt.Point;

import java.awt.geom.Point2D;

import java.awt.image.BufferedImage;

import java.awt.image.DataBuffer;

import java.io.IOException;

import java.util.Vector;

import org.cresques.cts.ICoordTrans;

import org.cresques.cts.IProjection;

import org.cresques.io.raster.RasterBuf;

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;

        private String ext = "";

        private String shortName = "";

        // Polilinea con extent

        public class Contour extends Vector {

                final private static long serialVersionUID = -3370601314380922368L;

                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, aBandNr = 4;

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

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

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

                        WITH_OVERVIEWS = false;

                width = getRasterXSize();

                height = getRasterYSize();

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

                shortName = getDriverShortName();

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

                        //String [] metadata = getMetadata();

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

                        try{

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

                         }catch(GdalException exc){

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

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

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

                                 esq.add(new Point2D.Double(width, 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));

                  }

        }

        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, bandA = 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(this.getRasterCount() == 4 && shortName.equals("PNG"))

                                        bandA = getRasterBand(aBandNr); 

                        }

                        if (currentOverview > 0) {

                                bandR = bandR.getOverview(currentOverview);

                            if (this.getRasterCount() > 1) {

                                        bandG = bandG.getOverview(currentOverview);

                                        bandB = bandB.getOverview(currentOverview);

                                        if(this.getRasterCount() == 4 && shortName.equals("PNG"))

                                                bandA = bandA.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 w = (int)(((double)currentViewWidth)*step);

        int x = (int) currentViewX;

        int y = (int) lastReadLine;

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

        GdalBuffer a = new GdalBuffer();

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

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

                w = bandR.getRasterBandXSize()-x;

        if(bandR.getRasterColorTable() != null){

                p = bandR.readRasterWithPalette(x, y, w, 1, w, 1, dataType);

                a.buffByte = p.buffAPalette;

                r = new GdalBuffer();

                r.buffByte = p.buffRPalette;

                g = new GdalBuffer();

                g.buffByte = p.buffGPalette;

                b = new GdalBuffer();

                b.buffByte = p.buffBPalette;

        }else{

                a.buffByte = new byte[w];

                        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 i=0;

                  float j =0F;

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

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

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

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

                                                  if(k<3)

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

                                                  else

                                                          line[i][3] = 0xff;

                                          }

                              }

                          }else { // Multiband

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

                                  GdalBuffer [] bands = {r,g,b};

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

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

                                              if(k<3)

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

                                              else

                                                      line[i][3] = 0xff;

                                      }

                                  }

                          }

                  }else if(dataType == GDT_Float32){

                          GdalBuffer [] bands = {r,g,b};

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

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

                                      //Truncamos el float perdiendo centimetros. Es necesario para visualizar

                                      if(k < 3)

                                              line[i][k] = (int)bands[0].buffFloat[(int) j];

                                      else

                                              line[i][3] = 0xff;

                              }

                        }

                  }

                return;

        }

        int readLineRGBA(int [] line) throws GdalException {

                int err = 0;

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

        int x = (int) currentViewX;

        int y = (int) lastReadLine;

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

        GdalBuffer a = new GdalBuffer();

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

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

                w = bandR.getRasterBandXSize()-x;

        if(bandR.getRasterColorTable() != null){

                p = bandR.readRasterWithPalette(x, y, w, 1, w, 1, dataType);

                a.buffByte = p.buffAPalette;

                r = new GdalBuffer();

                r.buffByte = p.buffRPalette;

                g = new GdalBuffer();

                g.buffByte = p.buffGPalette;

                b = new GdalBuffer();

                b.buffByte = p.buffBPalette;

        }else{

                if(getRasterCount() == 4 && shortName.equals("PNG")){

                        a = bandA.readRaster(x, y, w, 1, w, 1, GDT_Byte);

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

                }else{

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

                            a.buffByte = new byte[w];

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

                                    a.buffByte[i] = (byte)255;

                }

        }

                  lastReadLine += step;

                  int 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 & ((a.buffByte[(int)j])& 0xff) << 24) + ((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 & ((a.buffByte[(int)j])& 0xff) << 24) + ((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 & ((a.buffByte[(int)j])& 0xff) << 24) + 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 & ((a.buffByte[(int)j])& 0xff) << 24) | (((r.buffShort[(int) j] & 0xfff0) << 12) & 0xff0000 ) | 

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

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

                              }

                          }

                  }

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

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

        GdalBuffer a = new GdalBuffer();

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

                w = bandR.getRasterBandXSize()-x;

        if(bandR.getRasterColorTable() != null){

                p = bandR.readRasterWithPalette(x, y, w, h, w, h, GDT_Byte);

                a.buffByte = p.buffAPalette;

                r = new GdalBuffer();

                r.buffByte = p.buffRPalette;

                g = new GdalBuffer();

                g.buffByte = p.buffGPalette;

                b = new GdalBuffer();

                b.buffByte = p.buffBPalette;

        }else{

                if(getRasterCount() == 4 && shortName.equals("PNG")){

                        a = bandA.readRaster(x, y, w, h, w, h, GDT_Byte);

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

                            if (bandG != null)

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

                            if (bandB != null)

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

                }else{

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

                            if (bandG != null)

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

                            if (bandB != null)

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

                            a.buffByte = new byte[w];

                            for (int i = 0;i < w*h;i++)

                                    a.buffByte[i] = (byte)255;

                }

        }

                  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 & ((a.buffByte[(int)j])& 0xff) << 24) + ((r.buffByte[(int) k]) << 16) + ((g.buffByte[(int) k]) << 8) + b.buffByte[(int) k];

                          }

                  }

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

                }

        }

        void pintaPaleta() {

        }

        public int getBlockSize(){

                return this.getBlockSize();

        }

}

/**

 * @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("bmp", GdalFile.class);

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

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

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

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

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

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

        }

        public GdalFile(IProjection proj, String fName){

                super(proj, fName);

                extent = new Extent();

                try {

                        file = new GdalNative(fName);

                        load();

                        showOnOpen();

                        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;

                }

                switch(file.getDataType()){

                case 1:setDataType(DataBuffer.TYPE_BYTE);break;//GDT_BYTE 

                case 2://GDT_UInt16 

                case 3:setDataType(DataBuffer.TYPE_SHORT);break;//GDT_Int16        

                case 4://GDT_UInt32

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

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

                case 7:setDataType(DataBuffer.TYPE_DOUBLE);break;//GDT_Float64

                case 8:setDataType(DataBuffer.TYPE_UNDEFINED);break;//GDT_CInt16

                case 9:setDataType(DataBuffer.TYPE_UNDEFINED);break;//GDT_CInt32

                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.

         */

        public GeoFile load() {

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

                return this;

        }

        /**

         * Cierra el fichero de imagen

         */

        public void close() {

                try {

                        if(file != null){

                                file.close();

                                file = null;

                        }

                } catch (GdalException e) {

                        // TODO Auto-generated catch block

                        e.printStackTrace();

                }

        }

        /**

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

         */

        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;

        }

        /**

         * Asigna el extent de la vista actual

         */

        public void setView(Extent e) { 

                v = new Extent(e); 

        }

        /**

         * Obtiene extent de la vista actual

         */

        public Extent getView() { 

                return v; 

        }

        /**

         * Obtiene la anchura del fichero

         */

        public int getWidth() {        

                return file.width; 

        }

        /**

         * Obtiene la altura del fichero

         */

        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#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);

                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*/];

                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 RasterBuf getRaster(int width, int height, ICoordTrans rp) {

                double dFileAspect, dWindowAspect;

                int line;

                RasterBuf 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 = new RasterBuf(DataBuffer.TYPE_INT, width, height, 4, new Point(0,0));

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

                                file.readLine(raster.getLineInt(line));

                        }

                } catch (Exception e) {

                        // TODO Auto-generated catch block

                        e.printStackTrace();

                }

                return raster;

        }

        /**

         * Asigna al objeto Image los valores con los dato de la imagen contenidos en el 

         * vector de enteros.

         * @param image        imagen con los datos actuales

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

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

         * @param w        Ancho de la imagen

         * @param h        Alto de la imagen

         * @param rgbArray        vector que contiene la banda que se va a sustituir

         * @param offset        desplazamiento

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

         */

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

        }

        /**

         * Asigna al objeto Image la mezcla entre los valores que ya tiene y los valores 

         * 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 image        imagen con los datos actuales

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

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

         * @param w        Ancho de la imagen

         * @param h        Alto de la imagen

         * @param rgbArray        vector que contiene la banda que se va a sustituir

         * @param offset        desplazamiento

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

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

         */

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

        }

        /**

         * Asigna al objeto Image la mezcla entre los valores que ya tiene y los valores 

         * 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 image        imagen con los datos actuales

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

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

         * @param w        Ancho de la imagen

         * @param h        Alto de la imagen

         * @param rgbArray        vector que contiene la banda que se va a sustituir

         * @param offset        desplazamiento

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

         * @param origBand        Banda origen del GeoRasterFile

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

         */

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

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

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

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

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

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

                                line[i] = (line[i] & 0x0000ffff) | (rgbArray[i] & 0xffff0000);

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

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

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

                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 == 0 && destBandFlag == GeoRasterFile.BLUE_BAND)

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

                                line[i] = (line[i] & 0xffffff00) | ((rgbArray[i] & 0x00ff0000) >> 16);

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

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

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

                                line[i] = (line[i] & 0xff00ffff) | ((rgbArray[i] & 0x000000ff) << 16);

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

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

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

                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, int)

         */

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

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

                if(width<=0)width=1;

                if(height<=0)height=1;

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

                try {

                        setBand(RED_BAND,   rBandNr);

                        setBand(GREEN_BAND, gBandNr);

                        setBand(BLUE_BAND,  bBandNr);

                        file.setAlpha(getAlpha());

                        if(img!=null){

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

                                        file.readLineRGBA(pRGBArray);

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

                                }

                                return img;

                        }else{

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

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

                                        file.readLineRGBA(pRGBArray);

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

                                }

                                return image;

                        }

                } 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;

        }

        /**

         * Devuelve los datos de una ventana solicitada

         * @param ulX        coordenada X superior izda.

         * @param ulY        coordenada Y superior derecha.

         * @param sizeX        tama?o en X de la ventana.

         * @param sizeY tama?o en Y de la ventana.

         * @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

         *@return N?mero de zona 

         */

        public int getUTM(){

                return 0;        

        }

        /**

         * Obtiene el sistema de coordenadas geograficas

         * @return Sistema de coordenadas geogr?ficas

         */

        public String getGeogCS(){

                return new String("");        

        }

        /**

         * Devuelve el tama?o de bloque

         * @return Tama?o de bloque

         */

        public int getBlockSize(){

     //TODO Nacho: Implementar getBlockSize de EcwFile        

          return file.getBlockSize();

        }

}