Application: gvSIG desktop » gvSIG raster

<?xml version="1.0" encoding="UTF-8"?>

<project xmlns="http://maven.apache.org/POM/4.0.0" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://maven.apache.org/POM/4.0.0 http://maven.apache.org/xsd/maven-4.0.0.xsd">

	<modelVersion>4.0.0</modelVersion>

	<artifactId>org.gvsig.raster.gdal.app</artifactId>

	<packaging>pom</packaging>

	<name>${project.artifactId}</name>

	<parent>

		<groupId>org.gvsig</groupId>

		<artifactId>org.gvsig.raster.gdal</artifactId>

		<version>2.2.16</version>

	</parent>

	<modules>

		<module>org.gvsig.raster.gdal.app.mainplugin</module>

	</modules>

</project>

<?xml version="1.0" encoding="ISO-8859-1"?>

<project xmlns="http://maven.apache.org/POM/4.0.0" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://maven.apache.org/POM/4.0.0 http://maven.apache.org/maven-v4_0_0.xsd">

    <modelVersion>4.0.0</modelVersion>

    <artifactId>org.gvsig.raster.gdal.app.mainplugin</artifactId>

    <packaging>jar</packaging>

    <name>${project.artifactId}</name>

    <description>GDAL raster data provider for gvSIG: jpeg, tiff, png, gif.</description>

    <parent>

        <groupId>org.gvsig</groupId>

        <artifactId>org.gvsig.raster.gdal.app</artifactId>

        <version>2.2.16</version>

    </parent>

    <dependencies>

        <dependency>

            <groupId>org.gvsig</groupId>

            <artifactId>org.gvsig.andami</artifactId>

            <scope>compile</scope>

        </dependency>

        <dependency>

            <groupId>org.gvsig</groupId>

            <artifactId>org.gvsig.raster.gdal.io</artifactId>

            <scope>runtime</scope>

        </dependency>

    </dependencies>

    <properties>

        <svn.tags.folder>https://devel.gvsig.org/svn/gvsig-raster/org.gvsig.raster.gdal/tags/</svn.tags.folder>

        <gvsig.package.info.name>Formats: GDAL raster file formats support</gvsig.package.info.name>     

        <gvsig.package.info.state>testing</gvsig.package.info.state>

        <gvsig.package.info.official>true</gvsig.package.info.official>

        <gvsig.package.info.operatingSystem>all</gvsig.package.info.operatingSystem>

        <gvsig.package.info.architecture>all</gvsig.package.info.architecture>

        <gvsig.package.info.categories>Raster, Formats</gvsig.package.info.categories>

        <gvsig.package.info.poolURL>https://devel.gvsig.org/download/projects/Raster/pool/</gvsig.package.info.poolURL>

        <gvsig.package.info.codealias>org.gvsig.raster.gdal.app</gvsig.package.info.codealias>

        <gvsig.package.info.description>${project.parent.description}</gvsig.package.info.description>

        <gvsig.package.info.dependencies>required: org.gvsig.gdal.app.mainplugin -ge 1, required: org.gvsig.raster.tilecache.app -ge 2</gvsig.package.info.dependencies>

        <gvsig.install.plugin>true</gvsig.install.plugin>

        <gvsig.install.plugin.package>true</gvsig.install.plugin.package>

    </properties>

    <build>

        <plugins>

            <plugin>

                <!-- Skip compilation tests -->

                <groupId>org.apache.maven.plugins</groupId>

                <artifactId>maven-compiler-plugin</artifactId>

                <executions>

                    <execution>

                        <id>default-testCompile</id>

                        <phase>process-test-sources</phase>

                        <goals>

                            <goal>testCompile</goal>

                        </goals>

                        <configuration>

                            <skip>true</skip>

                        </configuration>

                    </execution>

                </executions>

            </plugin>

            <plugin>

                <!-- Skip test execution -->

                <groupId>org.apache.maven.plugins</groupId>

                <artifactId>maven-surefire-plugin</artifactId>

                <configuration>

                    <skipTests>true</skipTests>

                </configuration>

            </plugin>

        </plugins>

    </build>

</project>

#Fri Oct 16 19:09:02 CEST 2015

buildNumber=68

<!--

    gvSIG. Desktop Geographic Information System.

    Copyright (C) 2007-2013 gvSIG Association.

    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 3

    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.

    For any additional information, do not hesitate to contact us

    at info AT gvsig.com, or visit our website www.gvsig.com.

-->

<assembly>

    <id>gvsig-plugin-package</id>

    <formats>

        <format>zip</format>

    </formats>

    <baseDirectory>${project.artifactId}</baseDirectory>

    <includeBaseDirectory>true</includeBaseDirectory>

    <files>

        <file>

            <source>target/${project.artifactId}-${project.version}.jar</source>

            <outputDirectory>lib</outputDirectory>

        </file>

        <file>

            <source>target/package.info</source>

        </file>

    </files>

    <fileSets>

        <fileSet>

            <directory>src/main/resources-plugin</directory>

            <outputDirectory>.</outputDirectory>

        </fileSet>

    </fileSets>

    <dependencySets>

        <dependencySet>

            <useProjectArtifact>false</useProjectArtifact>

            <useTransitiveDependencies>false</useTransitiveDependencies>

            <outputDirectory>lib</outputDirectory>

            <includes>

                <include>org.gvsig:org.gvsig.raster.gdal.io:jar</include>

            </includes>

        </dependencySet>

    </dependencySets>

</assembly>

/* 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.gdal.app;

import org.gvsig.andami.plugins.Extension;

/**

 * Extension to add Gdal library support

 *

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

 */

public class GdalExtension extends Extension {

	public void execute(String actionCommand) {

	}

	public void initialize() {

	}

	public boolean isEnabled() {

		return false;

	}

	public boolean isVisible() {

		return false;

	}

}

<?xml version="1.0" encoding="ISO-8859-1"?>

<plugin-config>

    <depends plugin-name="org.gvsig.gdal.app.mainplugin" />

    <depends plugin-name="org.gvsig.raster.mainplugin" />

    <depends plugin-name="org.gvsig.raster.tilecache.app" optional="true" />

    <libraries library-dir="lib"/>

    <resourceBundle name="text"/>

    <extensions>

        <extension class-name="org.gvsig.raster.gdal.app.GdalExtension"

                           description=""

                           active="true"

                           priority="1">

        </extension>

    </extensions>

</plugin-config>

<?xml version="1.0" encoding="UTF-8"?>

<project xmlns="http://maven.apache.org/POM/4.0.0" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://maven.apache.org/POM/4.0.0 http://maven.apache.org/xsd/maven-4.0.0.xsd">

	<modelVersion>4.0.0</modelVersion>

	<artifactId>org.gvsig.raster.gdal</artifactId>

	<packaging>pom</packaging>

	<version>2.2.16</version>

	<name>${project.artifactId}</name>

	<description>GDAL raster data provider for gvSIG</description>

	<inceptionYear>2011</inceptionYear>

	<parent>

		<groupId>org.gvsig</groupId>

		<artifactId>org.gvsig.desktop</artifactId>

		<version>2.0.105</version>

	</parent>

	<properties>

		<!-- El plugin versions:use-latest-versions falla con scope import -->

		<!-- asi que toca usar el versions:update-properties que si que funciona -->

		<org.gvsig.gdal.version>1.0.4</org.gvsig.gdal.version>

		<org.gvsig.raster.version>2.2.18</org.gvsig.raster.version>

	</properties>

	<repositories>

		<repository>

			<id>gvsig-public-http-repository</id>

			<name>gvSIG maven public HTTP repository</name>

			<url>http://devel.gvsig.org/m2repo/j2se</url>

			<releases>

				<enabled>true</enabled>

				<updatePolicy>daily</updatePolicy>

				<checksumPolicy>warn</checksumPolicy>

			</releases>

			<snapshots>

				<enabled>true</enabled>

				<updatePolicy>daily</updatePolicy>

				<checksumPolicy>warn</checksumPolicy>

			</snapshots>

		</repository>

	</repositories>

	<scm>

		<connection>scm:svn:https://devel.gvsig.org/svn/gvsig-raster/org.gvsig.raster.gdal/tags/org.gvsig.raster.gdal-2.2.16</connection>

		<developerConnection>scm:svn:https://devel.gvsig.org/svn/gvsig-raster/org.gvsig.raster.gdal/tags/org.gvsig.raster.gdal-2.2.16</developerConnection>

		<url>https://devel.gvsig.org/redmine/projects/gvsig-raster/repository/show/org.gvsig.raster.gdal/tags/org.gvsig.raster.gdal-2.2.16</url>

	</scm>

	<build>

		<plugins>

			<plugin>

				<groupId>org.apache.maven.plugins</groupId>

				<artifactId>maven-release-plugin</artifactId>

				<configuration>

					<tagBase>https://devel.gvsig.org/svn/gvsig-raster/org.gvsig.raster.gdal/tags/</tagBase>

					<goals>deploy</goals>

				</configuration>

			</plugin>

		</plugins>

	</build>

	<developers>

		<developer>

			<id>nbrodin</id>

			<name>Nacho Brodin</name>

			<email>nachobrodin@gmail.com</email>

			<roles>

				<role>Architect</role>

				<role>Developer</role>

			</roles>

		</developer>

	</developers>

	<dependencyManagement>

		<dependencies>

			<dependency>

				<groupId>org.gvsig</groupId>

				<artifactId>org.gvsig.raster</artifactId>

				<version>${org.gvsig.raster.version}</version>

				<type>pom</type>

				<scope>import</scope>

			</dependency>

			<dependency>

				<groupId>org.gvsig</groupId>

				<artifactId>org.gvsig.gdal</artifactId>

				<version>${org.gvsig.gdal.version}</version>

				<type>pom</type>

				<scope>import</scope>

			</dependency>

                        <dependency>

				<groupId>org.gvsig</groupId>

				<artifactId>org.gvsig.raster.gdal.app.common</artifactId>

				<version>2.2.14-SNAPSHOT</version>

				<type>zip</type>

			</dependency>

			<dependency>

				<groupId>org.gvsig</groupId>

				<artifactId>org.gvsig.raster.gdal.io</artifactId>

				<version>2.2.16</version>

			</dependency>

		</dependencies>

	</dependencyManagement>

	<modules>

		<module>org.gvsig.raster.gdal.io</module>

		<module>org.gvsig.raster.gdal.app</module>

	</modules>

</project>

org.gvsig.raster.gdal.io.DefaultGdalIOLibrary

/* 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.addo;

import java.io.File;

import org.gdal.gdal.Dataset;

import org.gdal.gdal.gdal;

import org.gvsig.jgdal.JNIBase;

/**

 * Clase para la construccion de overviews de un raster.

 *

 * 15-nov-2007

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

 */

public class Jaddo extends JNIBase implements IOverviewIncrement{

	public final static int NEAREST = 0;

	public final static int AVERAGE = 1;

	public final static int AVERAGE_MP = 2;

	public final static int AVERAGE_MAGPHASE = 3;

	public final static int MODE = 4;

	private IOverviewIncrement  incrementListener = null;

	private int value = 0;

//	private native int buildOverviewsNative(int resamplingAlg, String file, int[] overviews);

	/**

	 * Construccion de overviews

	 */

	public void buildOverviews(int resamplingAlg, String file, int[] overviews)

		throws BuildingOverviewsException, WritingException {

		File f = new File(file);

		if(!f.exists() || !f.isFile())

			throw new BuildingOverviewsException("File does not exist.");

		if(!f.canWrite())

			throw new WritingException("File is not writeable");

		Dataset data = gdal.Open(file);

		if(data.BuildOverviews(overviews) >= 0)

			throw new BuildingOverviewsException("Problems building overviews");

	}

	/**

	 * Devuelve el porcentaje del incremento.

	 * @return int

	 */

	public int getPercent() {

		return value;

	}

	/**

	 * Asigna el porcentaje de incremento de la construccion de overview.

	 * Esto se hace automaticamente desde el callback que asigna el porcentaje.

	 */

	public void setPercent(int value) {

		this.value = value;

		if(incrementListener != null)

			incrementListener.setPercent(value);

	}

	/**

	 * Asigna el listener para la asignacion del incremento

	 * @param incr IOverviewIncrement

	 */

	public void setIncrementListener(IOverviewIncrement incr) {

		this.incrementListener = incr;

	}

}

/* 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.addo;

/**

 * Es generada no se puede escribir las overviews

 * 

 * @author Nacho Brodin (nachobrodin@gmail.com).<BR> Equipo de desarrollo gvSIG.<BR> http://www.gvsig.gva.es

 * @version 0.0

 * @link http://www.gvsig.org

 */

public class WritingException extends Exception {

	public WritingException(String msg){

		super(msg);

	}

}

/* 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.addo;

/**

 * Interfaz para la asignacion y recuperacion del incremento

 * de la tarea.

 *

 * 17-nov-2007

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

 */

public interface IOverviewIncrement {

	/**

	 * Devuelve el porcentaje del incremento

	 * @return int

	 */

	public int getPercent();

	/**

	 * Asigna el porcentaje de incremento de la construccion de overview.

	 * Esto se hace automaticamente desde el callback que asigna el porcentaje. 

	 */

	public void setPercent(int value);

}

/* 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.addo;

/**

 * Es generada cuando hay algun problema en la construccion de 

 * las overviews y la funcion nativa devuelve 0

 * 

 * @author Nacho Brodin (nachobrodin@gmail.com).<BR> Equipo de desarrollo gvSIG.<BR> http://www.gvsig.gva.es

 * @version 0.0

 * @link http://www.gvsig.org

 */

public class BuildingOverviewsException extends Exception {

	public BuildingOverviewsException(String msg) {

		super(msg);

	}

}

/* 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.gdal.overview;

import org.gvsig.addo.BuildingOverviewsException;

import org.gvsig.addo.IOverviewIncrement;

import org.gvsig.addo.Jaddo;

import org.gvsig.addo.WritingException;

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

import org.gvsig.fmap.dal.coverage.process.BaseIncrementableTask;

import org.gvsig.fmap.dal.coverage.process.overview.OverviewBuilder;

/**

 * Servicios ofrecidos por el proceso de construcci?n de overviews. Esta

 * implementaci?n usa gdal para construir las overviews.

 * 

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

 */

public class GdalOverviewBuilder implements OverviewBuilder {

	private Jaddo        build         = new Jaddo();

	private int          value         = 0;

	/**

	 * Bindings between an IOverviewIncrement and a BaseIncrementableTask

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

	 */

	class ProcessBinding implements IOverviewIncrement {

		private BaseIncrementableTask  task = null;

		public ProcessBinding(BaseIncrementableTask task) {

			this.task = task;	

		}

		public int getPercent() {

			return task.getPercent();

		}

		public void setPercent(int value) {

			task.setPercent(value);

		}

	}

	/*

	 * (non-Javadoc)

	 * @see org.gvsig.fmap.dal.coverage.process.overview.OverviewBuilder#buildOverviews(int, java.lang.String, int[])

	 */

	public void buildOverviews(int resamplingAlg, String file, int[] overviews)

			throws OverviewException {

		try {

			build.buildOverviews(resamplingAlg, file, overviews);

		} catch (BuildingOverviewsException e) {

			throw new OverviewException("Error building overviews", e);

		} catch (WritingException e) {

			throw new OverviewException("Error writing overviews", e);

		}

	}

	/*

	 * (non-Javadoc)

	 * @see org.gvsig.fmap.dal.coverage.process.overview.OverviewBuilder#getPercent()

	 */

	public int getPercent() {

		return value;

	}

	/*

	 * (non-Javadoc)

	 * @see org.gvsig.fmap.dal.coverage.process.overview.OverviewBuilder#setIncrementListener(org.gvsig.addo.IOverviewIncrement)

	 */

	public void setIncrementListener(BaseIncrementableTask incr) {

		build.setIncrementListener(new ProcessBinding(incr));

	}

	/*

	 * (non-Javadoc)

	 * @see org.gvsig.fmap.dal.coverage.process.overview.OverviewBuilder#setPercent(int)

	 */

	public void setPercent(int value) {

		this.value = value;

	}

}

/* 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.gdal.io;

import org.gvsig.fmap.dal.DALFileLibrary;

import org.gvsig.fmap.dal.DALLibrary;

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

import org.gvsig.raster.gdal.overview.GdalOverviewBuilder;

import org.gvsig.raster.gdal.util.DefaultCRSUtils;

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

import org.gvsig.tools.ToolsLibrary;

import org.gvsig.tools.library.AbstractLibrary;

import org.gvsig.tools.library.Library;

import org.gvsig.tools.library.LibraryException;

/**

 *

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

 */

public class DefaultGdalIOLibrary extends AbstractLibrary {	

	public DefaultGdalIOLibrary() {

		super(DefaultGdalIOLibrary.class,Library.TYPE.IMPL);

		require(ToolsLibrary.class);

		require(DALLibrary.class);

		require(DALFileLibrary.class);

	}

	@Override

	protected void doInitialize() throws LibraryException {

		//RasterLibrary.wakeUp();

	}

	@Override

	protected void doPostInitialize() throws LibraryException {

		//AbstractRasterFileDataParameters.registerDynClass();

		AbstractNewRasterStoreParameters.registerDynClass();

		// Registro de los drivers de lectura

		GdalProvider.register();

		// Registro de los drivers de escritura

		GdalWriter.register();

		JpegWriter.register();

		PngWriter.register();

		GdalDataParameters.registerDynClass();

		//Register GdalOverviewBuilder

		RasterLocator.getManager().getProviderServices().registerOverviewBuilderImplementation(

				"GDAL",

				new GdalOverviewBuilder());

		//Register CRSUtils

		RasterLocator.getManager().getProviderServices().registerCRSUtilImplementation(

						new DefaultCRSUtils());

	}

}

	/* 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.gdal.io;

import java.awt.Color;

import java.awt.Rectangle;

import java.awt.geom.AffineTransform;

import java.awt.geom.NoninvertibleTransformException;

import java.awt.geom.Point2D;

import java.io.IOException;

import java.util.ArrayList;

import java.util.List;

import org.gdal.gdal.Dataset;

import org.gdal.gdal.gdal;

import org.gdal.gdal.gdalJNI;

import org.gdal.ogr.ogrJNI;

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

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.ColorItem;

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

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

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

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

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

import org.gvsig.fmap.dal.coverage.util.FileUtils;

import org.gvsig.jgdal.Gdal;

import org.gvsig.jgdal.GdalBuffer;

import org.gvsig.jgdal.GdalColorEntry;

import org.gvsig.jgdal.GdalColorTable;

import org.gvsig.jgdal.GdalDataset;

import org.gvsig.jgdal.GdalException;

import org.gvsig.jgdal.GdalRasterBand;

import org.gvsig.jgdal.GeoTransform;

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

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

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

import org.gvsig.raster.impl.process.RasterTask;

import org.gvsig.raster.impl.process.RasterTaskQueue;

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

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

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

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

import org.gvsig.tools.dispose.Disposable;

import org.gvsig.tools.task.TaskStatus;

/**

 * Soporte 'nativo' para ficheros desde GDAL.

 * 

 * @author Luis W. Sevilla (sevilla_lui@gva.es)

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

 */

public class GdalNative extends Gdal implements Disposable {

//	private boolean isInitialized = false;

	private String                       fileName                = null;

	private String                       shortName               = "";

	public 	GeoTransform                 trans                   = null;

	public int                           width                   = 0, height = 0;

	public double                        originX                 = 0D, originY = 0D;

	public String                        version                 = "";

	protected int                        rBandNr                 = 1, gBandNr = 2, bBandNr = 3, aBandNr = 4;

	private int[]                        dataType                = null;

	DataStoreMetadata                    metadata                = null;

	protected boolean                    georeferenced           = true;

	/**

	 * Vectores que contiene los desplazamientos de un pixel cuando hay supersampling.

	 * , es decir el n?mero de pixels de pantalla que tiene un pixel de imagen. Como todos

	 * los pixeles no tienen el mismo ancho y alto ha de meterse en un array y no puede ser

	 * una variable. Adem?s hay que tener en cuenta que el primer y ?ltimo pixel son de 

	 * distinto tama?o que el resto.   

	 */

	public int[]                              stepArrayX             = null;

	public int[]                              stepArrayY             = null;

	protected GdalRasterBand[]                gdalBands              = null;

	private double                            lastReadLine           = -1;

	private int                               overviewWidth          = -1;

	private int                               overviewHeight         = -1;

	private int                               currentViewWidth       = -1;

	private int                               currentViewHeight      = -1;

	private double                            currentViewX           = 0D;

	private double                            viewportScaleX         = 0D;

	private double                            viewportScaleY         = 0D;

	private double                            stepX                  = 0D;

	private double                            stepY                  = 0D;

	public boolean                            isSupersampling        = false;

	private boolean                           open                   = false;

	/**

	 * Estado de transparencia del raster.

	 */

	protected DataStoreTransparency           fileTransparency       = null;

	protected DataStoreColorTable             palette                = null;

	protected DataStoreColorInterpretation    colorInterpr           = null;

	protected AffineTransform                 ownTransformation      = null;

	protected AffineTransform                 externalTransformation = new AffineTransform();

	public static int getGdalTypeFromRasterBufType(int rasterBufType) {

		switch (rasterBufType) {

			case Buffer.TYPE_BYTE: return Gdal.GDT_Byte;

			case Buffer.TYPE_USHORT: return Gdal.GDT_UInt16;

			case Buffer.TYPE_SHORT: return Gdal.GDT_Int16;

			case Buffer.TYPE_INT: return Gdal.GDT_Int32;

			case Buffer.TYPE_FLOAT: return Gdal.GDT_Float32;

			case Buffer.TYPE_DOUBLE: return Gdal.GDT_Float64;

			case Buffer.TYPE_UNDEFINED: return Gdal.GDT_Unknown;

			case Buffer.TYPE_IMAGE: return Gdal.GDT_Byte;

		}

		return Gdal.GDT_Unknown;

	}

	/**

	 * Conversi?n de los tipos de datos de gdal a los tipos de datos de RasterBuf

	 * @param gdalType Tipo de dato de gdal

	 * @return Tipo de dato de RasterBuf

	 */

	public static int getRasterBufTypeFromGdalType(int gdalType) {

		switch (gdalType) {

			case 1:// Eight bit unsigned integer GDT_Byte = 1

				return Buffer.TYPE_BYTE;

			case 3:// Sixteen bit signed integer GDT_Int16 = 3,

				return Buffer.TYPE_SHORT;

			case 2:// Sixteen bit unsigned integer GDT_UInt16 = 2

				//return RasterBuffer.TYPE_USHORT;

				return Buffer.TYPE_SHORT; //Apa?o para usar los tipos de datos que soportamos

			case 5:// Thirty two bit signed integer GDT_Int32 = 5

				return Buffer.TYPE_INT;

			case 6:// Thirty two bit floating point GDT_Float32 = 6

				return Buffer.TYPE_FLOAT;

			case 7:// Sixty four bit floating point GDT_Float64 = 7

				return Buffer.TYPE_DOUBLE;

				// TODO:Estos tipos de datos no podemos gestionarlos. Habria que definir

				// el tipo complejo y usar el tipo long que de momento no se gasta.

			case 4:// Thirty two bit unsigned integer GDT_UInt32 = 4,

				return Buffer.TYPE_INT;

				//return RasterBuffer.TYPE_UNDEFINED; // Deberia devolver un Long

			case 8:// Complex Int16 GDT_CInt16 = 8

			case 9:// Complex Int32 GDT_CInt32 = 9

			case 10:// Complex Float32 GDT_CFloat32 = 10

			case 11:// Complex Float64 GDT_CFloat64 = 11

				return Buffer.TYPE_UNDEFINED;

		}

		return Buffer.TYPE_UNDEFINED;

	}

	/**

	 * Overview usada en el ?ltimo setView

	 */

	int currentOverview = -1;

	public GdalNative(String fName) throws GdalException, IOException {

		super();

		init(fName);

	}

	private static long initializeGdal(String fName) throws GdalException {

		gdalJNI.AllRegister();

		ogrJNI.RegisterAll();

		Dataset data = gdal.Open(fName, 1);

		if (data == null)

			throw new GdalException("Error en la apertura del fichero. El fichero no tiene un formato v?lido.");

		return GdalDataset.getCPtr(data);

	}

//	public GdalNative(long cPtr, String fName) throws GdalException, IOException {

//		super(cPtr);

//		isInitialized = true;

//		init(fName);

//	}

	private void init(String fName) throws GdalException, IOException {

		gdalJNI.AllRegister();

		ogrJNI.RegisterAll();

		fileName = fName;

		open(fName, GA_ReadOnly);

		open = true;

		if (getDataset() == null)

			throw new GdalException("Error en la apertura del fichero. El fichero no tiene un formato v?lido.");

//		if(isInitialized){

//		ext = RasterUtilities.getExtensionFromFileName(fName);

		width = getRasterXSize();

		height = getRasterYSize();

		int[] dt = new int[getRasterCount()];

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

			dt[i] = this.getRasterBand(i + 1).getRasterDataType();

		setDataType(dt);

		shortName = getDriverShortName();

		colorInterpr = new DataStoreColorInterpretation(getRasterCount());

		fileTransparency = new DataStoreTransparency(colorInterpr);

		metadata = new DataStoreMetadata(getMetadata(), colorInterpr);

		// Asignamos la interpretaci?n de color leida por gdal a cada banda. Esto

		// nos sirve para saber que banda de la imagen va asignada a cada banda de

		// visualizaci?n (ARGB)

		metadata.initNoDataByBand(getRasterCount());

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

			GdalRasterBand rb = getRasterBand(i + 1);

			String colorInt = getColorInterpretationName(rb.getRasterColorInterpretation());

			metadata.setNoDataEnabled(rb.existsNoDataValue());

			if(rb.existsNoDataValue()) {

				metadata.setNoDataValue(i, rb.getRasterNoDataValue());

				metadata.setNoDataEnabled(rb.existsNoDataValue());

			}

			colorInterpr.setColorInterpValue(i, colorInt);

			if (colorInt.equals("Alpha"))

				fileTransparency.setTransparencyBand(i);

			if (rb.getRasterColorTable() != null && palette == null) {

				palette = new DataStoreColorTable(gdalColorTable2ColorItems(rb.getRasterColorTable()), false);

//				fileTransparency.setTransparencyRangeList(palette.getTransparencyRange());

			}

		}

		fileTransparency.setTransparencyByPixelFromMetadata(metadata);

		try {

			trans = getGeoTransform();

			boolean isCorrect = false;

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

				if (trans.adfgeotransform[i] != 0)

					isCorrect = true;

			if (!isCorrect)

				throw new GdalException("");

			double psX = trans.adfgeotransform[1];

			double psY = trans.adfgeotransform[5];

			double rotX = trans.adfgeotransform[4];

			double rotY = trans.adfgeotransform[2];

			double offX = trans.adfgeotransform[0];

			double offY = trans.adfgeotransform[3];

			ownTransformation = new AffineTransform(psX, rotX, rotY, psY, offX, offY);

					//trans.adfgeotransform[1], trans.adfgeotransform[4], trans.adfgeotransform[2], trans.adfgeotransform[5], trans.adfgeotransform[0], trans.adfgeotransform[3]);

			externalTransformation = (AffineTransform) ownTransformation.clone();

			overviewWidth = width;

			overviewHeight = height;

			this.georeferenced = true;

		} catch (GdalException exc) {

			// Transformaci�n para ficheros sin georreferenciaci�n. Se invierte la Y

			// ya que las WC decrecen de

			// arriba a abajo y los pixeles crecen de arriba a abajo

			ownTransformation = new AffineTransform(1, 0, 0, -1, 0, height);

			externalTransformation = (AffineTransform) ownTransformation.clone();

			overviewWidth = width;

			overviewHeight = height;

			this.georeferenced = false;

		}

//		}

	}

	/**

	 * Returns true if this provider is open and false if don't

	 * @return

	 */

	public boolean isOpen() {

		return open;

	}

	/**

	 * Obtiene el flag que informa de si el raster tiene valor no data o no.

	 * Consultar� todas las bandas del mismo y si alguna tiene valor no data

	 * devuelve true sino devolver� false.

	 * @return true si tiene valor no data y false si no lo tiene

	 * @throws GdalException

	 */

	public boolean existsNoDataValue() throws GdalException {

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

			GdalRasterBand rb = getRasterBand(i + 1);

			if (rb.existsNoDataValue())

				return true;

		}

		return false;

	}

	/**

	 * Obtiene el flag que informa de si el raster tiene valor no data o no

	 * en una banda concreta.

	 * @return true si tiene valor no data en esa banda y false si no lo tiene

	 * @param band Posici�n de la banda a consultar (0..n)

	 * @throws GdalException

	 */

	public boolean existsNoDataValue(int band) throws GdalException {

		GdalRasterBand rb = getRasterBand(band + 1);

		return rb.existsNoDataValue();

	}

	/**

	 * Gets nodata value

	 * @return

	 */

	public NoData getNoDataValue() {

		Number value = null;

		int type = getRasterBufTypeFromGdalType(getDataType()[0]);

		if (metadata != null && metadata.isNoDataEnabled() && metadata.getNoDataValue().length > 0) {

			switch (type) {

			case Buffer.TYPE_BYTE:

				if (metadata == null || metadata.getNoDataValue().length == 0)

					value = new Byte(RasterLibrary.defaultByteNoDataValue);

				else

					value = new Byte((byte)metadata.getNoDataValue()[0]);

				break;

			case Buffer.TYPE_SHORT:

				if (metadata == null || metadata.getNoDataValue().length == 0)

					value = new Short(RasterLibrary.defaultShortNoDataValue);

				else

					value = new Short((short)metadata.getNoDataValue()[0]);

				break;

			case Buffer.TYPE_INT:

				if (metadata == null || metadata.getNoDataValue().length == 0)

					value = new Integer((int)RasterLibrary.defaultIntegerNoDataValue);

				else

					value = new Integer((int)metadata.getNoDataValue()[0]);

				break;

			case Buffer.TYPE_FLOAT:

				if (metadata == null || metadata.getNoDataValue().length == 0)

					value = new Float(RasterLibrary.defaultFloatNoDataValue);

				else

					value = new Float(metadata.getNoDataValue()[0]);

				break;

			case Buffer.TYPE_DOUBLE:

				if (metadata == null || metadata.getNoDataValue().length == 0)

					value = new Double(RasterLibrary.defaultFloatNoDataValue);

				else

					value = new Double(metadata.getNoDataValue()[0]);

				break;

			}

		}

		NoData nodata = new DefaultNoData(value, value, fileName);

		nodata.setNoDataTransparent(false);

		return nodata;

	}

	/**

	 * Asigna el tipo de dato

	 * @param dt entero que representa el tipo de dato

	 */

	public void setDataType(int[] dt) { 

		dataType = dt; 

	}

	/**

	 * Obtiene el tipo de dato

	 * @return entero que representa el tipo de dato

	 */

	public int[] getDataType() { 

		return dataType; 

	}

	/**

	 * Gets the color interpretation

	 * @return

	 */

	public ColorInterpretation getColorInterpretation() { 

		return colorInterpr; 

	}

	/**

	 * Gets the color table

	 * @return

	 */

	public ColorTable getColorTable() {

		return palette;

	}

	/**

	 * Obtiene un punto 2D con las coordenadas del raster a partir de uno en coordenadas

	 * del punto real.

	 * Supone rasters no girados

	 * @param pt	punto en coordenadas del punto real

	 * @return	punto en coordenadas del raster

	 */

	public Point2D worldToRasterWithoutRot(Point2D pt) {

		Point2D p = new Point2D.Double();

		AffineTransform at = new AffineTransform(	externalTransformation.getScaleX(), 0, 

													0, externalTransformation.getScaleY(), 

													externalTransformation.getTranslateX(), externalTransformation.getTranslateY());

		try {

			at.inverseTransform(pt, p);

		} catch (NoninvertibleTransformException e) {

			return pt;

		}

		return p;

	}

	/**

	 * Obtiene un punto 2D con las coordenadas del raster a partir de uno en coordenadas

	 * del punto real.

	 * Supone rasters no girados

	 * @param pt	punto en coordenadas del punto real

	 * @return	punto en coordenadas del raster

	 */

	public Point2D worldToRaster(Point2D pt) {

		Point2D p = new Point2D.Double();

		try {

			externalTransformation.inverseTransform(pt, p);

		} catch (NoninvertibleTransformException e) {

			return pt;

		}

		return p;

	}

	/**

	 * Obtiene un punto del raster en coordenadas pixel a partir de un punto en coordenadas

	 * reales. 

	 * @param pt Punto en coordenadas reales

	 * @return Punto en coordenadas pixel.

	 */

	public Point2D rasterToWorld(Point2D pt) {

		Point2D p = new Point2D.Double();

		externalTransformation.transform(pt, p);

		return p;

	}

	/**

	 * Calcula el overview a usar. Hay que tener en cuenta que tenemos que tener calculadas las variables

	 * viewPortScale, currentFullWidth y currentFulHeight

	 * @param coordenada pixel expresada en double que indica la posici�n superior izquierda

	 * @throws GdalException

	 */

	private void calcOverview(Point2D tl, Point2D br) throws GdalException {

		gdalBands[0] = getRasterBand(1);

		currentOverview = -1;

		if (gdalBands[0].getOverviewCount() > 0) {

			GdalRasterBand ovb = null;

			for (int i = gdalBands[0].getOverviewCount() - 1; i > 0; i--) {

				ovb = gdalBands[0].getOverview(i);

				if (ovb.getRasterBandXSize() > getRasterXSize() * viewportScaleX) {

					currentOverview = i;

					viewportScaleX *= ((double) width / (double) ovb.getRasterBandXSize());

					viewportScaleY *= ((double) height / (double) ovb.getRasterBandYSize());

					stepX = 1D / viewportScaleX;

					stepY = 1D / viewportScaleY;

					overviewWidth = ovb.getRasterBandXSize();

					overviewHeight = ovb.getRasterBandYSize();

					currentViewX = Math.min(tl.getX(), br.getX());

					lastReadLine = Math.min(tl.getY(), br.getY());

					break;

				}

			}

		}

	}

	public void setView(double dWorldTLX, double dWorldTLY,

						double dWorldBRX, double dWorldBRY,

						int nWidth, int nHeight) throws GdalException {

		overviewWidth = width;

		overviewHeight = 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;

//		wcWidth = Math.abs(br.getX() - tl.getX());

		currentViewX = Math.min(tl.getX(), br.getX());

		viewportScaleX = (double) currentViewWidth / (br.getX() - tl.getX());

		viewportScaleY = (double) currentViewHeight / (br.getY() - tl.getY());

		stepX = 1D / viewportScaleX;

		stepY = 1D / viewportScaleY;

		lastReadLine = Math.min(tl.getY(), br.getY());

		//Para lectura del renderizado (ARGB). readWindow selecciona las bandas que necesita.

		// calcula el overview a usar

		gdalBands = new GdalRasterBand[4];

		calcOverview(tl, br);

	}

	/**

	 * Selecciona bandas y overview en el objeto GdalRasterBand[] para el n�mero de bandas solicitado.

	 * @param nbands N�mero de bandas solicitado.

	 * @throws GdalException

	 */

	public void selectGdalBands(int nbands) throws GdalException {

		gdalBands = new GdalRasterBand[nbands];

		// Selecciona las bandas y los overviews necesarios

		gdalBands[0] = getRasterBand(1);

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

			gdalBands[i] = gdalBands[0];

		assignDataTypeFromGdalRasterBands(gdalBands);

//		setDataType(gdalBands[0].getRasterDataType());

		for (int i = 2; i <= nbands; i++) {

			if (getRasterCount() >= i) {

				gdalBands[i - 1] = getRasterBand(i);

				for (int j = i; j < nbands; j++)

					gdalBands[j] = gdalBands[i - 1];

			}

		}

		if (currentOverview > 0) {

			gdalBands[0] = gdalBands[0].getOverview(currentOverview);

			for (int i = 2; i <= nbands; i++) {

				if (getRasterCount() >= i)

					gdalBands[i - 1] = gdalBands[i - 1].getOverview(currentOverview);

			}

		}

	}

	int lastY = -1;

	/**

	 * Lee una l�nea de bytes

	 * @param line Buffer donde se cargan los datos

	 * @param initOffset Desplazamiento inicial desde el margen inzquierdo. Esto es necesario para cuando

	 * se supersamplea ya que cada pixel de imagen ocupa muchos pixeles de pantalla y puede empezar a dibujarse

	 * por la izquierda a mitad de pixel

	 * @param gdalBuffer Buffer con la l�nea de datos original

	 */

	private void readLine(byte[][] line, double initOffset, GdalBuffer[] gdalBuffer) {

		double j = 0D;

		int i = 0;

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

			for (i = 0, j = initOffset; i < currentViewWidth && j < gdalBuffer[0].getSize(); i++, j += stepX) {

				line[iBand][i] = gdalBuffer[iBand].buffByte[(int) j];

			}

		}

	}

	/**

	 * Lee una l�nea de shorts

	 * @param line Buffer donde se cargan los datos

	 * @param initOffset Desplazamiento inicial desde el margen inzquierdo. Esto es necesario para cuando

	 * se supersamplea ya que cada pixel de imagen ocupa muchos pixeles de pantalla y puede empezar a dibujarse

	 * por la izquierda a mitad de pixel

	 * @param gdalBuffer Buffer con la l�nea de datos original

	 */

	private void readLine(short[][] line, double initOffset, GdalBuffer[] gdalBuffer) {

		double j = 0D;

		int i = 0;

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

			for (i = 0, j = initOffset; i < currentViewWidth && j < gdalBuffer[0].getSize(); i++, j += stepX) {

				line[iBand][i] = (short) (gdalBuffer[iBand].buffShort[(int) j] & 0xffff);

			}

		}

	}

	/**

	 * Lee una l�nea de ints

	 * @param line Buffer donde se cargan los datos

	 * @param initOffset Desplazamiento inicial desde el margen inzquierdo. Esto es necesario para cuando

	 * se supersamplea ya que cada pixel de imagen ocupa muchos pixeles de pantalla y puede empezar a dibujarse

	 * por la izquierda a mitad de pixel

	 * @param gdalBuffer Buffer con la l�nea de datos original

	 */

	private void readLine(int[][] line, double initOffset, GdalBuffer[] gdalBuffer) {

		double j = 0D;

		int i = 0;

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

			for (i = 0, j = initOffset; i < currentViewWidth && j < gdalBuffer[0].getSize(); i++, j += stepX) {

				line[iBand][i] = (gdalBuffer[iBand].buffInt[(int) j] & 0xffffffff);

			}

		}

	}

	/**

	 * Lee una l�nea de float

	 * @param line Buffer donde se cargan los datos

	 * @param initOffset Desplazamiento inicial desde el margen izquierdo. Esto es necesario para cuando

	 * se supersamplea ya que cada pixel de imagen ocupa muchos pixeles de pantalla y puede empezar a dibujarse

	 * por la izquierda a mitad de pixel

	 * @param gdalBuffer Buffer con la l�nea de datos original

	 */

	private void readLine(float[][] line, double initOffset, GdalBuffer[] gdalBuffer) {

		double j = 0D;

		int i = 0;

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

			for (i = 0, j = initOffset; i < currentViewWidth && j < gdalBuffer[0].getSize(); i++, j += stepX) {

				line[iBand][i] = gdalBuffer[iBand].buffFloat[(int) j];

			}

		}

	}

	/**

	 * Lee una l�nea de doubles

	 * @param line Buffer donde se cargan los datos

	 * @param initOffset Desplazamiento inicial desde el margen inzquierdo. Esto es necesario para cuando

	 * se supersamplea ya que cada pixel de imagen ocupa muchos pixeles de pantalla y puede empezar a dibujarse

	 * por la izquierda a mitad de pixel

	 * @param gdalBuffer Buffer con la l�nea de datos original

	 */

	private void readLine(double[][] line, double initOffset, GdalBuffer[] gdalBuffer) {

		double j = 0D;

		int i = 0;

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

			for (i = 0, j = initOffset; i < currentViewWidth && j < gdalBuffer[0].getSize(); i++, j += stepX) {

				line[iBand][i] = gdalBuffer[iBand].buffDouble[(int) j];

			}

		}

	}

	/**

	 * Lee una l�nea completa del raster y devuelve un array del tipo correcto. Esta funci�n es util