Application: gvSIG desktop » gvSIG raster

/* gvSIG. Geographic Information System of the Valencian Government

 *

 * Copyright (C) 2007-2008 Infrastructures and Transports Department

 * of the Valencian Government (CIT)

 *

 * 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., 51 Franklin Street, Fifth Floor, Boston,

 * MA  02110-1301, USA.

 *

 */

package org.gvsig.raster.netcdf.io;

import java.awt.geom.AffineTransform;

import java.awt.geom.Point2D;

import java.awt.geom.Rectangle2D;

import java.io.IOException;

import java.util.List;

import org.gvsig.fmap.dal.DALFileLocator;

import org.gvsig.fmap.dal.DALLocator;

import org.gvsig.fmap.dal.DataStore;

import org.gvsig.fmap.dal.coverage.RasterLocator;

import org.gvsig.fmap.dal.coverage.dataset.Buffer;

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

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

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

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

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

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

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

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

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

import org.gvsig.fmap.dal.coverage.store.parameter.RasterFileStoreParameters;

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

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

import org.gvsig.fmap.dal.spi.DataManagerProviderServices;

import org.gvsig.fmap.dal.spi.DataStoreProviderServices;

import org.gvsig.metadata.MetadataLocator;

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

import org.gvsig.raster.impl.DefaultRasterManager;

import org.gvsig.raster.impl.buffer.SpiRasterQuery;

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

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

import org.gvsig.raster.impl.provider.AbstractRasterProvider;

import org.gvsig.raster.impl.provider.RasterProvider;

import org.gvsig.raster.impl.provider.tile.FileTileServer;

import org.gvsig.raster.impl.store.DefaultRasterStore;

import org.gvsig.raster.impl.store.DefaultStoreFactory;

import org.gvsig.raster.impl.store.properties.DataStoreColorInterpretation;

import org.gvsig.raster.impl.store.properties.DataStoreTransparency;

import org.gvsig.tools.ToolsLocator;

import org.slf4j.Logger;

import org.slf4j.LoggerFactory;

import ucar.ma2.Array;

import ucar.ma2.DataType;

import ucar.ma2.InvalidRangeException;

import ucar.ma2.Range;

import ucar.nc2.dt.GridDatatype;

import ucar.nc2.dt.grid.GridDataset;

import ucar.unidata.geoloc.ProjectionRect;

/**

 * Data provider for NetCDF files with raster data

 *

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

 */

public class NetCDFProvider extends AbstractRasterProvider {

        public static String                     NAME                     = "NetCDF Raster";

        public static String                     DESCRIPTION              = "NetCDF Raster file";

        public final String                      METADATA_DEFINITION_NAME = NAME;

        private Extent                           viewRequest              = null;

        private TileServer                       tileServer               = null;

        private boolean                          open                     = false;

    private GridDataset                      gridNetCDF               = null;

    private List<GridDatatype>               gridList                 = null;

    private GridDatatype                     selectedGridDataType     = null;

    private DataStoreTransparency            fileTransparency         = null;

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

    protected static String[]                formatList               = null;

        public static void register() {

                DataManagerProviderServices dataman = (DataManagerProviderServices) DALLocator.getDataManager();

                RasterLocator.getManager().getProviderServices().registerFileProvidersTiled(NetCDFProvider.class);

                registerFormats();

                if (dataman != null && !dataman.getStoreProviders().contains(NAME)) {

                        dataman.registerStoreProvider(NAME,

                                        NetCDFProvider.class, NetCDFDataParametersImpl.class);

                }

                /*if (!dataman.getExplorerProviders().contains(NetCDFFilesystemServerExplorer.NAME)) {

                        dataman.registerExplorerProvider(NetCDFFilesystemServerExplorer.NAME, NetCDFFilesystemServerExplorer.class, NetCDFServerExplorerParameters.class);

                }*/

                if(DALFileLocator.getFilesystemServerExplorerManager() != null)

                        DALFileLocator.getFilesystemServerExplorerManager().registerProvider(

                                        NAME, DESCRIPTION,

                                        NetCDFFilesystemServerExplorer.class);

                dataman.registerStoreFactory(NAME, DefaultStoreFactory.class);

        }

        private static void registerFormats() {

                formatList      = new String[] {"nc", "nc4"};

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

                        RasterLocator.getManager().getProviderServices().addFormat(formatList[i], NetCDFProvider.class);

        }

        public String[] getFormatList() {

                return formatList;

        }

        /**

         * Returns true if the extension is supported and false if doesn't

         * @param ext

         * @return

         */

        public boolean isExtensionSupported(String ext) {

                if(ext.indexOf(".") != -1)

                        ext = ext.substring(ext.lastIndexOf(".") + 1, ext.length());

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

                        if(formatList[i].compareTo(ext) == 0)

                                return true;

                }

                return false;

        }

        public NetCDFProvider() {

        }

        /**

         * Opens the dataset.

         * @param proj Projection

         * @param fName File name

         * @throws NotSupportedExtensionException

         */

        public NetCDFProvider(String params) throws NotSupportedExtensionException, OpenException {

                super(params);

                if(params instanceof String) {

                        NetCDFDataParameters p = new NetCDFDataParametersImpl();

                        p.setURI((String)params);

                        super.init(p, null, ToolsLocator.getDynObjectManager()

                                        .createDynObject(

                                                        MetadataLocator.getMetadataManager().getDefinition(

                                                                        DataStore.METADATA_DEFINITION_NAME)));

                        init(p, null);

                }

        }

        public NetCDFProvider (NetCDFDataParameters params,

                        DataStoreProviderServices storeServices) throws NotSupportedExtensionException, OpenException {

                super(params, storeServices, ToolsLocator.getDynObjectManager()

                                .createDynObject(

                                                MetadataLocator.getMetadataManager().getDefinition(

                                                                DataStore.METADATA_DEFINITION_NAME)));

                init(params, storeServices);

        }

        /**

         * Build file references

         * @param proj Projection

         * @param param Load parameters

         * @throws NotSupportedExtensionException

         */

        public void init (NetCDFDataParameters params,

                        DataStoreProviderServices storeServices) throws NotSupportedExtensionException, OpenException {

                if(((RasterFileStoreParameters)params).getFile().exists()) {

                        String fileName = ((RasterFileStoreParameters)params).getFile().getAbsolutePath();

            try { 

                    gridNetCDF = GridDataset.open(fileName);

                    gridList = gridNetCDF.getGrids();

                    if(gridList.size() == 0)

                            throw new OpenException("There is not a grid variable", null);

                    selectedGridDataType = gridList.get(0);

                //netCDFFile = NetcdfFile.open(((RasterFileStoreParameters)params).getFile().getAbsolutePath());

            } catch (IOException e) {

                throw new OpenException("Imposible to read the file", e);

            }  

            /*List<Variable> variableList = netCDFFile.getVariables();

            Iterator<Variable> it = variableList.iterator();

            while(it.hasNext()) {

                    Variable var = it.next();

                    System.out.println("===>>" + var.getName());

            }*/

                        setParam(storeServices, params);

                        reloadMetadataFromGrid();

                        noData = new DefaultNoData(Double.NaN, Double.NaN, fileName);

                        load();

                } else

                        setParam(storeServices, params);

                super.init();

                selectSubdataset(getId(0, 0, 0));

                try {

                        loadFromRmf(getRmfBlocksManager());

                } catch (ParsingException e) {

                        //No lee desde rmf

                }

                open = true;

        }

        /**

         * Reloads metadata using the selected grid

         */

        private void reloadMetadataFromGrid() {

                //wktProjection = null;

                //CrsWkt crs = new CrsWkt(wktProjection);

                //IProjection proj = CRSFactory.getCRS("EPSG:23030");

                /*LatLonRect gcs = selectedGridDataType.getCoordinateSystem().getLatLonBoundingBox();

                getColorInterpretation();

                double scaleX = gcs.getWidth() / selectedGridDataType.getXDimension().getLength();

                double scaleY = gcs.getHeight() / selectedGridDataType.getYDimension().getLength();

                ownTransformation = new AffineTransform(

                                scaleX, 0, 

                                0, -scaleY, 

                                gcs.getLonMin(), 

                                gcs.getLatMax());*/

                ProjectionRect pRect = selectedGridDataType.getCoordinateSystem().getBoundingBox();

                double scaleX = pRect.getWidth() / selectedGridDataType.getXDimension().getLength();

                double scaleY = pRect.getHeight() / selectedGridDataType.getYDimension().getLength();

                ownTransformation = new AffineTransform(

                                scaleX, 0, 

                                0, -scaleY, 

                                pRect.getMinX(), 

                                pRect.getMaxY());

                externalTransformation = (AffineTransform)ownTransformation.clone();

                bandCount = 1; //One variable is always shown

                setDataType();

        }

        /**

         * @param dataType The dataType to set.

         */

        private void setDataType() {

                DataType dt = selectedGridDataType.getDataType();

                int type = Buffer.TYPE_UNDEFINED;

                if(dt.name() == DataType.BYTE.name()) {

                        type = Buffer.TYPE_BYTE;

                }

                if(dt.name() == DataType.SHORT.name()) {

                        type = Buffer.TYPE_SHORT;

                }

                if(dt.name() == DataType.INT.name()) {

                        type = Buffer.TYPE_INT;

                }

                if(dt.name() == DataType.DOUBLE.name()) {

                        type = Buffer.TYPE_DOUBLE;

                }

                if(dt.name() == DataType.FLOAT.name()) {

                        type = Buffer.TYPE_FLOAT;

                }

                if(dt.name() == DataType.LONG.name()) {

                        type = Buffer.TYPE_DOUBLE;

                }

                int[] dtype = new int[getBandCount()];

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

                        dtype[i] = type;

                setDataType(dtype);

        }

        public RasterProvider load() {

                return this;

        }

        public boolean isOpen() {

                return open;

        }

        public void close() {

                try {

                        gridNetCDF.close();

                } catch (IOException e) {

                        logger.error("Error closing file", e);

                }

        }

        public String translateFileName(String fileName) {

                return fileName;

        }

        /**

         * Asigna el extent de la vista actual. existe un fichero .rmf debemos hacer una transformaci�n

         * de la vista asignada ya que la petici�n viene en coordenadas del fichero .rmf y la vista (v)

         * ha de estar en coordenadas del fichero.

         */

        public void setView(Extent e) {

                viewRequest = new ExtentImpl(e);

        }

        public Extent getView() {

                return viewRequest;

        }

        public double getWidth() {

                return selectedGridDataType.getXDimension().getLength();

        }

        public double getHeight() {

                return selectedGridDataType.getYDimension().getLength();

        }

        public Object readBlock(int pos, int blockHeight, double scale)

                throws InvalidSetViewException, FileNotOpenException, RasterDriverException, ProcessInterruptedException {

                NetCDFDataParameters p = (NetCDFDataParameters)param;

                if(pos < 0)

                        throw new InvalidSetViewException("Request out of grid");

                if((pos + blockHeight) > getHeight())

                        blockHeight = Math.abs(((int)getHeight()) - pos);

                Buffer buf = DefaultRasterManager.getInstance().createBuffer(getDataType()[0], (int)getWidth(), blockHeight, 1, true);

                try {

                        int time = p.getFieldTime();

                        int level = getLevelValue();

                        Range rangeY = new Range(pos, pos + blockHeight - 1, 1);

                        Range rangeX = new Range(0, (int)(getWidth() - 1), 1);

                        GridDatatype dt = selectedGridDataType.makeSubset(null, null, null, null /*getTime(1), getLevel(1)*/, rangeY, rangeX);

                        Array values = dt.readDataSlice(time, level, -1, -1);

                        return arrayValuesToBuffer(values, buf, rangeX.length(), rangeY.length(), null);

                } catch (IOException e) {

                        throw new RasterDriverException("Error reading a slice", e);

                } catch (InvalidRangeException e) {

                        throw new RasterDriverException("Error reading a slice", e);

                }

        }

        public Object getData(int x, int y, int band)throws InvalidSetViewException, FileNotOpenException, RasterDriverException {

                NetCDFDataParameters p = (NetCDFDataParameters)param;

                if(x < 0 || y < 0 || x >= getWidth() || y >= getHeight())

                        throw new InvalidSetViewException("Request out of grid");

                try {

                        int strideX = 1;

                        int strideY = 1;

                        Range rangeY = new Range((int)(getHeight() - y), (int)(getHeight() - y), strideY);

                        Range rangeX = new Range(x, x, strideX);

                        selectSubdataset();

                        int time = p.getFieldTime();

                        int level = getLevelValue();

                        GridDatatype dt = selectedGridDataType.makeSubset(null, null, null, null, rangeY, rangeX);

                        Array values = dt.readDataSlice(time, level, -1, -1);

                        Object data = null;

                        if(getDataType()[0] == Buffer.TYPE_BYTE) {

                                data = new java.lang.Integer(values.getByte(0));

                        }

                        if(getDataType()[0] == Buffer.TYPE_SHORT) {

                                data = new java.lang.Integer(values.getShort(0));

                        }

                        if(getDataType()[0] == Buffer.TYPE_INT) {

                                data = new java.lang.Integer(values.getInt(0));

                        }

                        if(getDataType()[0] == Buffer.TYPE_FLOAT) {

                                data = new java.lang.Float(values.getFloat(0));

                        }

                        if(getDataType()[0] == Buffer.TYPE_DOUBLE) {

                                data = new java.lang.Double(values.getDouble(0));

                        }

                        return data;

                } catch (IOException e) {

                        throw new RasterDriverException("Error reading a slice", e);

                } catch (InvalidRangeException e) {

                        throw new RasterDriverException("Error reading a slice", e);

                }

        }

        @Override

        public void loadBuffer(SpiRasterQuery query)

                        throws ProcessInterruptedException, RasterDriverException {

                NetCDFDataParameters p = (NetCDFDataParameters)param;

                setView(query.getAdjustedRequestBoundingBox());

                Point2D ul = new Point2D.Double(viewRequest.getULX(), viewRequest.getULY());

                Point2D lr = new Point2D.Double(viewRequest.getLRX(), viewRequest.getLRY());

                ul = worldToRaster(ul);

                lr = worldToRaster(lr);

                ul.setLocation(ul.getX() < 0 ? 0 : ul.getX(), ul.getY() < 0 ? 0 : ul.getY());

                lr.setLocation(lr.getX() < 0 ? 0 : lr.getX(), lr.getY() < 0 ? 0 : lr.getY());

                ul.setLocation(ul.getX() >= getWidth() ? getWidth() - 1 : ul.getX(), ul.getY() >= getHeight() ? getHeight() - 1 : ul.getY());

                lr.setLocation(lr.getX() >= getWidth() ? getWidth() - 1 : lr.getX(), lr.getY() >= getHeight() ? getHeight() - 1 : lr.getY());

                adjustPoints(ul, lr);

                //int width = Math.abs(((int)lr.getX()) - ((int)ul.getX()));

                //int height = Math.abs(((int)lr.getY()) - ((int)ul.getY()));

                try {

                        int strideX = 1;

                        int strideY = 1;

                        Range rangeY = new Range((int)(getHeight() - lr.getY()), (int)(getHeight() - ul.getY() - 1), strideY <= 0 ? 1 : strideY);

                        Range rangeX = new Range((int)ul.getX(), (int)(lr.getX() - 1), strideX <= 0 ? 1 : strideX);

                        selectSubdataset();

                        int time = p.getFieldTime();

                        int level = getLevelValue();

                        GridDatatype dt = selectedGridDataType.makeSubset(null, null, null, null /*getTime(strideX), getLevel(strideX)*/, rangeY, rangeX);

                        Array values = dt.readDataSlice(time, level, -1, -1);

                        Rectangle2D r = query.getBufferForProviders().getDataExtent();

                        Buffer rasterBuf = arrayValuesToBuffer(

                                        values, 

                                        query.getBufferForProviders(), 

                                        rangeX.length(), 

                                        rangeY.length(), 

                                        query.getBandList());

                        rasterBuf.setDataExtent(r);

                        query.setBufferResult(rasterBuf);

                } catch (IOException e) {

                        throw new RasterDriverException("Error reading a slice", e);

                } catch (InvalidRangeException e) {

                        throw new RasterDriverException("Error reading a slice", e);

                }

        }

        /*public Buffer getWindow(Extent extent, BandList bandList, Buffer rasterBuf, TaskStatus status) 

                throws ProcessInterruptedException, RasterDriverException {

                NetCDFDataParameters p = (NetCDFDataParameters)param;

                setView(extent);

                Point2D ul = new Point2D.Double(viewRequest.getULX(), viewRequest.getULY());

                Point2D lr = new Point2D.Double(viewRequest.getLRX(), viewRequest.getLRY());

                ul = worldToRaster(ul);

                lr = worldToRaster(lr);

                ul.setLocation(ul.getX() < 0 ? 0 : ul.getX(), ul.getY() < 0 ? 0 : ul.getY());

                lr.setLocation(lr.getX() < 0 ? 0 : lr.getX(), lr.getY() < 0 ? 0 : lr.getY());

                ul.setLocation(ul.getX() > getWidth() ? getWidth() : ul.getX(), ul.getY() > getHeight() ? getHeight() : ul.getY());

                lr.setLocation(lr.getX() > getWidth() ? getWidth() : lr.getX(), lr.getY() > getHeight() ? getHeight() : lr.getY());

                adjustPoints(ul, lr);

                try {

                        int strideX = 1;

                        int strideY = 1;

                        Range rangeY = new Range((int)(getHeight() - lr.getY()), (int)(getHeight() - ul.getY() - 1), strideY);

                        Range rangeX = new Range((int)ul.getX(), (int)(lr.getX() - 1), strideX);

                        selectSubdataset();

                        int time = p.getFieldTime();

                        int level = getLevelValue();

                        GridDatatype dt = selectedGridDataType.makeSubset(null, null, null, null , rangeY, rangeX);

                        Array values = dt.readDataSlice(time, level, -1, -1);

                        rasterBuf = arrayValuesToBuffer(values, rasterBuf, rangeX.length(), rangeY.length(), bandList);

                } catch (IOException e) {

                        throw new RasterDriverException("Error reading a slice", e);

                } catch (InvalidRangeException e) {

                        throw new RasterDriverException("Error reading a slice", e);

                }

                return rasterBuf;

        }*/

        /*public Buffer getWindow(Extent extent, int bufWidth, int bufHeight, 

                        BandList bandList, Buffer rasterBuf, boolean adjustToExtent, TaskStatus status) throws ProcessInterruptedException, RasterDriverException {

                NetCDFDataParameters p = (NetCDFDataParameters)param;

                setView(extent);

                Point2D ul = new Point2D.Double(viewRequest.getULX(), viewRequest.getULY());

                Point2D lr = new Point2D.Double(viewRequest.getLRX(), viewRequest.getLRY());

                ul = worldToRaster(ul);

                lr = worldToRaster(lr);

                ul.setLocation(ul.getX() < 0 ? 0 : ul.getX(), ul.getY() < 0 ? 0 : ul.getY());

                lr.setLocation(lr.getX() < 0 ? 0 : lr.getX(), lr.getY() < 0 ? 0 : lr.getY());

                ul.setLocation(ul.getX() >= getWidth() ? getWidth() - 1 : ul.getX(), ul.getY() >= getHeight() ? getHeight() - 1 : ul.getY());

                lr.setLocation(lr.getX() >= getWidth() ? getWidth() - 1 : lr.getX(), lr.getY() >= getHeight() ? getHeight() - 1 : lr.getY());

                adjustPoints(ul, lr);

                int width = Math.abs(((int)lr.getX()) - ((int)ul.getX()));

                int height = Math.abs(((int)lr.getY()) - ((int)ul.getY()));

                try {

                        int strideX = width / rasterBuf.getWidth();

                        int strideY = height / rasterBuf.getHeight();

                        Range rangeY = new Range((int)(getHeight() - lr.getY()), (int)(getHeight() - ul.getY() - 1), strideY <= 0 ? 1 : strideY);

                        Range rangeX = new Range((int)ul.getX(), (int)(lr.getX() - 1), strideX <= 0 ? 1 : strideX);

                        selectSubdataset();

                        int time = p.getFieldTime();

                        int level = getLevelValue();

                        GridDatatype dt = selectedGridDataType.makeSubset(null, null, null, null , rangeY, rangeX);

                        Array values = dt.readDataSlice(time, level, -1, -1);

                        rasterBuf = arrayValuesToBuffer(values, rasterBuf, rangeX.length(), rangeY.length(), bandList);

                } catch (IOException e) {

                        throw new RasterDriverException("Error reading a slice", e);

                } catch (InvalidRangeException e) {

                        throw new RasterDriverException("Error reading a slice", e);

                }

                return rasterBuf;

        }*/

        /*public Buffer getWindow(int x, int y, int w, int h, 

                        BandList bandList, Buffer rasterBuf, TaskStatus status) throws ProcessInterruptedException, RasterDriverException {

                try {

                        NetCDFDataParameters p = (NetCDFDataParameters)param;

                        int strideX = 1;

                        int strideY = 1;

                        Range rangeY = new Range((int)(h - y), (int)(h - y - 1), strideY);

                        Range rangeX = new Range((int)x, (int)(x + w - 1), strideX);

                        selectSubdataset();

                        int time = p.getFieldTime();

                        int level = getLevelValue();

                        GridDatatype dt = selectedGridDataType.makeSubset(null, null, null, null, rangeY, rangeX);

                        Array values = dt.readDataSlice(time, level, -1, -1);

                        rasterBuf = arrayValuesToBuffer(values, rasterBuf, rangeX.length(), rangeY.length(), bandList);

                } catch (IOException e) {

                        throw new RasterDriverException("Error reading a slice", e);

                } catch (InvalidRangeException e) {

                        throw new RasterDriverException("Error reading a slice", e);

                }

                return rasterBuf;

        }*/

        /**

         * Gets the Range of the selected time or null if there is not a selected time

         * @param strideX

         * @return

         * @throws InvalidRangeException

         */

        @SuppressWarnings("unused")

        private Range getTime(int strideX) throws InvalidRangeException {

                NetCDFDataParameters p = (NetCDFDataParameters)param;

                int time = p.getFieldTime();

                return new Range(time, time, strideX);

        }

        /*private int getTimeValue() {

                int time = 0;

                if(param.hasDynValue(NetCDFDataParameters.FIELD_TIME))

                        time = ((Integer)param.getDynValue(NetCDFDataParameters.FIELD_TIME)).intValue();

                return time;

        }*/

        /**

         * Gets the Range of the selected level or null if there is not a selected level

         * @param strideX

         * @return

         * @throws InvalidRangeException

         */

        @SuppressWarnings("unused")

        private Range getLevel(int strideX) throws InvalidRangeException {

                int level = getLevelValue();

                return new Range(level, level, strideX);

        }

        private int getLevelValue() {

                int level = 0;

                if(param.hasDynValue(NetCDFDataParameters.FIELD_LEVEL) && param.getDynValue(NetCDFDataParameters.FIELD_LEVEL) != null)

                        level = ((Integer)param.getDynValue(NetCDFDataParameters.FIELD_LEVEL)).intValue();

                return level;

        }

        /**

         * Selects the GridDataType using the selected variable in the parameters

         */

        public void selectSubdataset() {

                if(param.hasDynValue(NetCDFDataParameters.FIELD_VARIABLE)) {

                        NetCDFDataParameters p = (NetCDFDataParameters)param;

                        String variable = (String)param.getDynValue(NetCDFDataParameters.FIELD_VARIABLE);

                        if(variable != null) {

                                for (int j = 0; j < gridList.size(); j++) {

                                        if(gridList.get(j).getName().compareTo(variable) == 0) {

                                                selectedGridDataType = gridList.get(j);

                                                reloadMetadataFromGrid();

                                                super.selectSubdataset(getId(j, getLevelValue(), p.getFieldTime()));

                                        }

                                }

                        }

                }

        }

        /**

         * Gets the identifier of a subdataset

         * @param grid

         * @param level

         * @param time

         * @return

         */

        private String getId(int grid, int level, int time) {

                return grid + "-" + level + "-" +  time;

        }

        /**

         * Loads a Buffer object from an netCDF Array. 

         * @param values

         * @param rasterBuf

         * @param w

         * @param h

         * @return

         */

        private Buffer arrayValuesToBuffer(Array values, Buffer rasterBuf, int w, int h, BandList bandList) {

                Buffer buf = null;

                boolean resampling = false;

                if((rasterBuf.getWidth() * rasterBuf.getHeight()) != values.getSize()) {

                        buf = DefaultRasterManager.getInstance().createBuffer(getDataType()[0], w, h, rasterBuf.getBandCount(), true);

                        resampling = true;

                } else

                        buf = rasterBuf;

                int[] drawableBands = bandList != null ? bandList.getBufferBandToDraw(getURIOfFirstProvider(), 0) : new int[1];

                if(drawableBands == null || (drawableBands.length == 1 && drawableBands[0] == -1))

                        return rasterBuf;

                if(getDataType()[0] == Buffer.TYPE_BYTE) {

                        for (int i = 0; i < values.getSize(); i++) {

                                int[] rc = getColAndRow(i, buf);

                                for (int iBands = 0; iBands < drawableBands.length; iBands++) 

                                        buf.setElem(rc[0], rc[1], drawableBands[iBands], values.getByte(i));

                        }

                }

                if(getDataType()[0] == Buffer.TYPE_SHORT) {

                        for (int i = 0; i < values.getSize(); i++) {

                                int[] rc = getColAndRow(i, buf);

                                for (int iBands = 0; iBands < drawableBands.length; iBands++) 

                                        buf.setElem(rc[0], rc[1], drawableBands[iBands], values.getShort(i));

                        }

                }

                if(getDataType()[0] == Buffer.TYPE_INT) {

                        for (int i = 0; i < values.getSize(); i++) {

                                int[] rc = getColAndRow(i, buf);

                                for (int iBands = 0; iBands < drawableBands.length; iBands++) 

                                        buf.setElem(rc[0], rc[1], drawableBands[iBands], values.getInt(i));

                        }

                }

                if(getDataType()[0] == Buffer.TYPE_FLOAT) {

                        for (int i = 0; i < values.getSize(); i++) {

                                int[] rc = getColAndRow(i, buf);

                                for (int iBands = 0; iBands < drawableBands.length; iBands++) 

                                        buf.setElem(rc[0], rc[1], drawableBands[iBands], values.getFloat(i));

                        }

                }

                if(getDataType()[0] == Buffer.TYPE_DOUBLE) {

                        for (int i = 0; i < values.getSize(); i++) {

                                int[] rc = getColAndRow(i, buf);

                                for (int iBands = 0; iBands < drawableBands.length; iBands++) 

                                        buf.setElem(rc[0], rc[1], drawableBands[iBands], values.getDouble(i));

                        }

                }

                if(resampling) {

                        try {

                                Buffer result = null;

                                result = buf.getAdjustedWindow(rasterBuf.getWidth(), rasterBuf.getHeight(), Buffer.INTERPOLATION_NearestNeighbour);

                                if(result != buf)

                                        buf.dispose();

                                return result;

                        } catch (ProcessInterruptedException e) {

                                return buf;

                        }

                }

                return buf;

        }

        /**

         * Calculates the row and column number for the position i in the

         * array of data

         * @param i

         * @param buf

         * @return

         */

        private int[] getColAndRow(int i, Buffer buf) {

                int auxRow = (int)(i / buf.getWidth());

                return new int[]{ 

                                (int)buf.getHeight() - auxRow - 1,

                                i - (auxRow * buf.getWidth())

                                };

        }

        private void adjustPoints(Point2D ul, Point2D lr) {

                double a = (ul.getX() - (int)ul.getX());

                double b = (ul.getY() - (int)ul.getY());

                ul.setLocation(        (a > 0.95 || a < 0.05) ? Math.round(ul.getX()) : ul.getX(), 

                                                (b > 0.95 || b < 0.05) ? Math.round(ul.getY()) : ul.getY());

                lr.setLocation(        (a > 0.95 || a < 0.05) ? Math.round(lr.getX()) : lr.getX(), 

                                                (b > 0.95 || b < 0.05) ? Math.round(lr.getY()) : lr.getY());

        }

        public int getBlockSize(){

                return 0;

        }

        public ColorInterpretation getColorInterpretation() {

                if(super.getColorInterpretation() == null) {

                        setColorInterpretation(new DataStoreColorInterpretation(new String[]{ColorInterpretation.GRAY_BAND}));

                }

                return super.getColorInterpretation();

        }

        public DataStoreTransparency getTransparency() {

                if(fileTransparency == null)

                        fileTransparency = new DataStoreTransparency(getColorInterpretation());

                return fileTransparency;

        }

        /**

         * Informa de si el driver ha supersampleado en el �ltimo dibujado. Es el driver el que colocar�

         * el valor de esta variable cada vez que dibuja.

         * @return true si se ha supersampleado y false si no se ha hecho.

         */

        public boolean isSupersampling() {

                return false;

        }

        public void setAffineTransform(AffineTransform t){

                super.setAffineTransform(t);

        }

        public int getOverviewCount(int band) throws BandAccessException, RasterDriverException {

                if(band >= getBandCount())

                        throw new BandAccessException("Wrong band");

                return 0;

        }

        public int getOverviewWidth(int band, int overview) throws BandAccessException, RasterDriverException {

                if (band >= getBandCount())

                        throw new BandAccessException("Wrong band");

                return 0;

        }

        public int getOverviewHeight(int band, int overview) throws BandAccessException, RasterDriverException {

                if (band >= getBandCount())

                        throw new BandAccessException("Wrong band");

                return 0;

        }

        public boolean isOverviewsSupported() {

                return true;

        }

        public boolean isReproyectable() {

                return true;

        }

        public String getProviderName() {

                return NAME;

        }

        public void setStatus(RasterProvider provider) {

                if(provider instanceof NetCDFProvider) {

                        //Not implemented yet

                }

        }

        public boolean isTimeSupported() {

                return true;

        }

        public TileServer getTileServer() {

                if(tileServer == null) {

                        DefaultRasterStore store = new DefaultRasterStore();

                        store.setProvider(this);

                        tileServer = new FileTileServer(store);

                }

                return tileServer;

        }

}