| gvSIG desktop 1

/*

 * EcwFile.java

 */

package org.cresques.io;

import java.awt.Image;

import java.awt.image.BufferedImage;

import com.ermapper.ecw.JNCSFile;

import com.ermapper.util.JNCSDatasetPoint;

import org.cresques.geo.ReProjection;

import org.cresques.px.*;

/**

 * Filtro para raster. Permite eliminar el fondo;

 * 

 * @author "Luis W. Sevilla" <sevilla_lui@gva.es>

 */

class PixelFilter {

        boolean debug = false;

        int transparente = 0x10ffff00;

        int orColor = 0xff000000;

        int umbral = 0xf0f0f0;

        /**

         * Constructor para el Oleico

         */

        public PixelFilter() {

                transparente = 0x10ffff00;

                orColor = 0xff000000;

                umbral = 0xf0f0f0;

        }

        /**

         * Constructor.

         * 

         * @param trans

         * @param or

         * @param pxMax

         */

        public PixelFilter(int trans, int or, int pxMax) {

                transparente = trans;

                orColor = or;

                umbral = pxMax;

        }

        /**

         * Filtra una l?nea de pixeles.

         * 

         * @param pRGBArray

         */

        public void filterLine(int [] pRGBArray) {

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

                        if (debug)

                                System.out.print(""+i+":"+Integer.toHexString(pRGBArray[i])+",");

                        if (pRGBArray[i]  >= umbral)

                                pRGBArray[i] = transparente;

                        else

                                pRGBArray[i] |= orColor;

                }

                if (debug)

                        System.out.println("");

        }

}

/**

 * Soporte para los ficheros .ecw de ErMapper

 * 

 * @author "Luis W. Sevilla" <sevilla_lui@gva.es>

 */

public class EcwFile extends GeoRasterFile {

        JNCSFile file = null;

        boolean bErrorOnOpen = false;

        String errorMessage = null;

        Extent v = null;

        boolean doTransparency = false;

        PixelFilter tFilter = null;

        public EcwFile(String fname) {

                super(null, null); // TODO aqui van las proyecciones

                fname = DataSource.normalize(fname);

                super.setName(fname);

                extent = new Extent();

                try {

                        System.err.println("Abriendo "+fname);

                        file = new JNCSFile(fname, false);

                        load();

                } catch(Exception e) {

                  bErrorOnOpen = true;

                  errorMessage = e.getMessage();

//                  g.drawString(errorMessage, 0, 50);

                  System.err.println(errorMessage);

                  e.printStackTrace();

                }

        }

        public void setExtent(Extent e) { extent = new Extent(e); }

        public Extent getExtent() { return extent; }

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

        public Extent getView() { return v; }

        public void setTransparency(boolean t) {

                doTransparency = t;

                tFilter = new PixelFilter();

        }

        public boolean getTransparency() { return doTransparency; }

        /**

         * Carga un ECW.

         * 

         * @param fname

         */

        public GeoFile load() {

                double minX, minY, maxX, maxY;

                System.out.println("ECW size = ("+file.width+","+file.height+")\n"+

                        " inc=("+file.cellIncrementX+","+file.cellIncrementY+")\n"+

                        " datum='"+file.datum+"', proyeccion='"+file.projection+"'");

                minX = file.originX;

                maxY = file.originY;

                maxX = file.originX + (double)(file.width-1)*file.cellIncrementX;

                minY = file.originY + (double)(file.height-1)*file.cellIncrementY;

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

                        System.out.println(extent);

                return this;

        }

        public void close() {

                file.close(true);

                file = null;

        }

        /**

         * Obtiene un trozo de imagen (determinado por la vista y los par?metros.

         * 

         * @param width

         * @param height

         */

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

                // TODO reproyectar

                BufferedImage ecwImage = null;

                if (file != null) {

                  try {

                        double dFileAspect, dWindowAspect;

                        //double dWorldTLX, dWorldTLY, dWorldBRX, dWorldBRY;

                        int bandlist[];

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

                        }

                        // Create an image of the ecw file.

                        if (doTransparency)

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

                        else

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

                        pRGBArray = new int[width];

                        // Setup the view parameters for the ecw file.

                        bandlist = new int[file.numBands];

                        for (int i=0; i< file.numBands; i++) {

                          bandlist[i] = i;

                        }

                        // Set the view

                        try {

                                JNCSDatasetPoint ptMin = file.convertWorldToDataset(v.minX(), v.minY());

                                JNCSDatasetPoint ptMax = file.convertWorldToDataset(v.maxX(), v.maxY());

                                System.out.println("Dataset coords Width = "+(ptMax.x-ptMin.x)+", px width ="+width);

                                file.setView(file.numBands, bandlist, 

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

                                        width, height);

                        } catch (com.ermapper.ecw.JNCSInvalidSetViewException e) {

                                System.err.println(errorMessage);

                                e.printStackTrace();

                        }

                        // Read the scan lines

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

                          file.readLineRGBA(pRGBArray);

                          // Prueba de sustituci?n de color transparente

                          if (doTransparency) {

                                        if (line == 0) tFilter.debug = true;

                                          tFilter.filterLine(pRGBArray);

                                        tFilter.debug = false; 

                          }

                          ecwImage.setRGB(0, line, width, 1, pRGBArray, 0, width);

                        }

                  } catch(Exception e) {

                        bErrorOnOpen = true;

                        errorMessage = e.getMessage();

//                        g.drawString(errorMessage, 0, 50);

                        System.err.println(errorMessage);

                        e.printStackTrace();

                  }

                }

                return ecwImage;

        }

        /**

         * Cambia el raster de proyecci?n.

         */

        public void reProject(ReProjection rp) {

                // TODO metodo reProject pendiente de implementar

        }

}