Application: gvSIG desktop

<project name="jlidarlib" default="clean-orig" basedir=".">

    <description>

        Genera jDielmoOpenLiDARlib.jar en /lib y una distribución

    </description>

	<!-- set global properties for this build -->

	<property name="src" location="src"/>

	<property name="build" location="bin"/>

	<property name="lib"  location="lib"/>

	<property name="doc"  location="doc"/>

	<property name="jarName" value="jDielmoOpenLiDARlib"/>

	<property name="sig_lidar_driver" location="../extDielmoOpenLidar/lib"/>

	<property name="dist"  location="dist"/>

	<property name="libDir" location="../_fwAndami/gvSIG/extensiones/com.iver.cit.gvsig/lib" />

	<property name="andamiLibs" location="../_fwAndami/lib" />

	<property name="compile-classpath" value="${libDir}/fmap.jar:${andamiLibs}/log4j-1.2.8.jar"/>

	<property name="JavaSourceVersion" value="1.4"/>

	<property name="JavaTargetVersion" value="1.4"/>

	<property name="encoding" value="ISO_8859_1"/>

  <target name="init">

    <!-- Create the time stamp -->

    <tstamp/>

	<echo>

		Compiling ${ant.project.name}...</echo>

  </target>

  <target name="make1" depends="init"

        description="generate jDielmoOpenLiDARlib.jar in /lib and in the gvSIG Driver" >

    <jar jarfile="${sig_lidar_driver}/${jarName}.jar">

      <fileset dir="${build}" />

    </jar>

  	<jar jarfile="${lib}/${jarName}.jar" basedir="${build}" includes="**/**"/>

  </target>

	<target name="clean-orig" depends="dist"

			description="clean dist directory" >

		<!-- Clean the distribution directory -->

		<delete dir="${dist}/src"/>

		<delete dir="${dist}/doc"/>

		<delete dir="${dist}/lib"/>

	</target>

  <target name="dist" depends="make1"

	description="generate the distribution" >

    <!-- Remove previous distribution directory -->

    <delete dir="${dist}"/>

    <!-- Create the distribution directory structure -->

    <mkdir dir="${dist}"/>

    <mkdir dir="${dist}/src"/>

    <mkdir dir="${dist}/doc"/>

    <mkdir dir="${dist}/lib"/>

    <!-- Copy necessary distribution files to dist dir -->

    <copy todir="${dist}/src">

      <fileset dir="${src}"/>

    </copy>

    <copy todir="${dist}/doc">

      <fileset dir="${doc}">

      </fileset>

    </copy>

  	<copy todir="${dist}/lib">

      <fileset dir="${lib}"/>

    </copy>

    <!-- Zip distribution -->

    <zip destfile="${dist}/jDielmoOpenLiDARlib.zip"

         basedir="${dist}"

         update="true"

    />

  </target>

	<target name="batch-build"

			description="compile the sources, create the jar file"

			depends="init,compile,create-jar,move-to-extDielmoOpenLidar">

	</target>

	<target name="compile" description="compile the source" >

		<!-- Compile the Java code from ${src} to ${build} -->

		<mkdir dir="${build}" />

		<javac  srcdir="${src}"

			destdir="${build}"

			classpath="${compile-classpath}"

			source="${JavaSourceVersion}"

			target="${JavaTargetVersion}"			

			debug="${debug}"

			debuglevel="${debuglevel}"

			excludes="com/iver/cit/gvsig/fmap/**"

			encoding="${encoding}"/>

	</target>

	<target name="create-jar"

			description="Creates the plugin jar">

		<mkdir dir="${dist}"/>

	    <jar 	jarfile="${dist}/${jarName}.jar"

	    		basedir="${build}"/>

	</target>

	<target name="move-to-extDielmoOpenLidar">

	  	<!-- jar file must be in gvSIG/lib directory -->

	  	<!--<copy file="${without_src}/${plugin}.jar" todir="${lib-dir}"/>-->

	    <move todir="${sig_lidar_driver}" file="${dist}/${jarName}.jar"/>

	</target>

	<target name="clean"

			description="clean up">

		<!-- Delete the ${build} and ${dist} directory trees -->

		<delete dir="${build}" />

		<delete dir="${dist}" />

	</target>

</project>

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

<classpath>

	<classpathentry kind="src" path="src"/>

	<classpathentry kind="con" path="org.eclipse.jdt.launching.JRE_CONTAINER"/>

	<classpathentry combineaccessrules="false" kind="src" path="/libIverUtiles"/>

	<classpathentry combineaccessrules="false" kind="src" path="/libGDBMS"/>

	<classpathentry combineaccessrules="false" kind="src" path="/libFMap"/>

	<classpathentry kind="output" path="bin"/>

</classpath>

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

<projectDescription>

	<name>libDielmoOpenLiDAR</name>

	<comment></comment>

	<projects>

	</projects>

	<buildSpec>

		<buildCommand>

			<name>org.eclipse.jdt.core.javabuilder</name>

			<arguments>

			</arguments>

		</buildCommand>

	</buildSpec>

	<natures>

		<nature>org.eclipse.jdt.core.javanature</nature>

	</natures>

</projectDescription>

outOFIndex=Punto inexperado: fuera del indice

outOFIndex=unexpected point: out of bounds

/* DielmoOpenLiDAR

 *

 * Copyright (C) 2008 DIELMO 3D S.L. (DIELMO) and 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.

 *

 * For more information, contact:

 *

 * DIELMO 3D S.L.

 * Plaza Vicente Andr?s Estell?s 1 Bajo E

 * 46950 Xirivella, Valencia

 * SPAIN

 *   

 * +34 963137212

 * dielmo@dielmo.com

 * www.dielmo.com

 * 

 * or

 * 

 * Generalitat Valenciana

 * Conselleria d'Infraestructures i Transport

 * Av. Blasco Ib??ez, 50

 * 46010 VALENCIA

 * SPAIN

 *

 * +34 963862235

 * gvsig@gva.es

 * www.gvsig.gva.es

 */

/*

 * AUTHORS (In addition to DIELMO and CIT):

 *  

 */

package com.dielmo.lidar;

import java.io.File;

import java.io.FileInputStream;

import java.io.FileNotFoundException;

import java.io.IOException;

import java.io.InputStream;

import javax.swing.JOptionPane;

/**

 * Inicizliza Lidar points and header

 *

 * @author Oscar Garcia

 */

public class InicializeLidar {

	public static int testLASFormat(File name) {

		File f = name;

		InputStream input;

		int offset = 0,numRead = 0;

		char[] FileSignatureLASF = new char[4];

		char  versionMajor, versionMinor, format;

		byte[] cabecera = new byte[105]; // lleemos hasta el byte que nos indica la version de LAS en la que estamos

		try {

	    	input = new FileInputStream(f);

			while (offset < cabecera.length && (numRead = input.read(cabecera, offset, cabecera.length-offset) ) >= 0) {

			     offset += numRead;

			}	

		    FileSignatureLASF[0] = (char)(cabecera[0] & 0xFF);

		    FileSignatureLASF[1] = (char)(cabecera[1] & 0xFF);

		    FileSignatureLASF[2] = (char)(cabecera[2] & 0xFF);

		    FileSignatureLASF[3] = (char)(cabecera[3] & 0xFF);

		    if(FileSignatureLASF[0] == 'L'&& FileSignatureLASF[1] == 'A'&& FileSignatureLASF[2] == 'S' && FileSignatureLASF[3] == 'F') {

			    if (offset < cabecera.length) {

			    	JOptionPane.showMessageDialog(null, "Bad input format");

				}

			    versionMajor = (char)(cabecera[24] & 0xFF);

			    versionMinor = (char)(cabecera[25] & 0xFF);

			    format = (char)(cabecera[104] & 0xFF);

		    	if(versionMajor == 1 && versionMinor == 0) {

		    		if(format == 0)

		    			return LidarHeader.LAS10F0;

		    		else if(format == 1)

		    			return LidarHeader.LAS10F1;

		    	} else if(versionMajor == 1 && versionMinor == 1) {

		    		if(format == 0)

		    			return LidarHeader.LAS11F0;

		    		else if(format == 1)

		    			return LidarHeader.LAS11F1;

		    	} else if (versionMajor == 1 && versionMinor == 2){

		    		if(format == 0)

		    			return LidarHeader.LAS12F0;

		    		else if(format == 1)

		    			return LidarHeader.LAS12F1;

		    		else if(format == 2)

		    			return LidarHeader.LAS12F2;

		    		else if(format == 3)

		    			return LidarHeader.LAS12F3;

		    	} else {

		    		throw new UnexpectedLidarHeaderException("Header not expected");

		    	}

		    }

		} catch(UnexpectedLidarHeaderException e) {

			e.printStackTrace();			

		} catch (FileNotFoundException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		} catch (IOException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

		return LidarHeader.UNEXPECTED;

	}

	public static int testBINFormat(File name) {

		File f = name;

		InputStream input;

		int offset = 0,numRead = 0;

		byte[] cabecera = new byte[8];

		try {

	    	input = new FileInputStream(f);

			while (offset < cabecera.length && (numRead = input.read(cabecera, offset, cabecera.length-offset) ) >= 0) {

			     offset += numRead;

			}

		    int HdrVersionBin; 

	    	HdrVersionBin = ByteUtilities.arr2Int(cabecera, 4);

	    	if(HdrVersionBin == 20020715) {

	    		return LidarHeader.BIN20020715;

	    	}

	    	else if( HdrVersionBin == 20010712) {

	    		return LidarHeader.BIN20010712;

	    	} //  ( HdrVersionBin == 20010129 || HdrVersionBin == 970404)	

		} catch (FileNotFoundException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		} catch (IOException e) {

			// TODO Auto-generated catch block

			e.printStackTrace();

		}

		return LidarHeader.UNEXPECTED;

	}

	public static int testFormat(File name) {

		int f = testLASFormat(name);

		if(f != LidarHeader.UNEXPECTED) {

			return f;

		}

		f = testBINFormat(name);

		if(f != LidarHeader.UNEXPECTED) {

			return f;

		}

		if (f == LidarHeader.UNEXPECTED) {

	    	JOptionPane.showMessageDialog(null, "Unsupported format");

		}

		return LidarHeader.UNEXPECTED;

	}

	public static LidarPoint InizializeLidarPoint(File file) {

		boolean binColor, binTime;

		int type = testFormat(file);

		switch(type) {

			case LidarHeader.LAS10F0:

				return new LASPoint10F0();

			case LidarHeader.LAS10F1:

				return new LASPoint10F1();

			case LidarHeader.LAS11F0:

				return new LASPoint11F0();

			case LidarHeader.LAS11F1:

				return new LASPoint11F1();

			case LidarHeader.LAS12F0:

				return new LASPoint12F0();

			case LidarHeader.LAS12F1:

				return new LASPoint12F1();

			case LidarHeader.LAS12F2:

				return new LASPoint12F2();

			case LidarHeader.LAS12F3:

				return new LASPoint12F3();

			case LidarHeader.BIN20010712:

				BINHeader headerBin2001 = new BINHeader(file);

				headerBin2001.readLidarHeader();

				binColor = false;

				binTime = false;

				if(headerBin2001.getColor()>0)

					binColor = true;

				if(headerBin2001.getTime()>0)

					binTime = true;

				return new BINPoint2001(binColor, binTime);

			case LidarHeader.BIN20020715:

				BINHeader headerBin2002 = new BINHeader(file);

				headerBin2002.readLidarHeader();

				binColor = false;

				binTime = false;

				if(headerBin2002.getColor()>0)

					binColor = true;

				if(headerBin2002.getTime()>0)

					binTime = true;

				return new BINPoint2002(binColor, binTime);

		}

		return null;

	}

	public static LidarHeader InizializeLidarHeader(File file) {

		int type = testFormat(file);

		switch(type) {

			case LidarHeader.LAS10F0:

			case LidarHeader.LAS10F1:

				LASHeader_V10 header10 = new LASHeader_V10(file);

				header10.readLidarHeader();

				return header10;

			case LidarHeader.LAS11F0:

			case LidarHeader.LAS11F1:

				LASHeader_V11 header11 = new LASHeader_V11(file);

				header11.readLidarHeader();

				return header11;

			case LidarHeader.LAS12F0:

			case LidarHeader.LAS12F1:

			case LidarHeader.LAS12F2:

			case LidarHeader.LAS12F3:

				LASHeader_V12 header12 = new LASHeader_V12(file);

				header12.readLidarHeader();

				return header12;

			case LidarHeader.BIN20010712:

				BINHeader headerBin2001 = new BINHeader(file);

				headerBin2001.readLidarHeader();

				return headerBin2001;

			case LidarHeader.BIN20020715:

				BINHeader headerBin2002 = new BINHeader(file);

				headerBin2002.readLidarHeader();	

				return headerBin2002;

		}

		return null;

	}

}

/* DielmoOpenLiDAR

 *

 * Copyright (C) 2008 DIELMO 3D S.L. (DIELMO) and 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.

 *

 * For more information, contact:

 *

 * DIELMO 3D S.L.

 * Plaza Vicente Andr?s Estell?s 1 Bajo E

 * 46950 Xirivella, Valencia

 * SPAIN

 *   

 * +34 963137212

 * dielmo@dielmo.com

 * www.dielmo.com

 * 

 * or

 * 

 * Generalitat Valenciana

 * Conselleria d'Infraestructures i Transport

 * Av. Blasco Ib??ez, 50

 * 46010 VALENCIA

 * SPAIN

 *

 * +34 963862235

 * gvsig@gva.es

 * www.gvsig.gva.es

 */

/*

 * AUTHORS (In addition to DIELMO and CIT):

 *  

 */

package com.dielmo.lidar;

public class UnexpectedPointException extends Exception {

	/**

	 * serial ID

	 */

	private static final long serialVersionUID = 207114730916728752L;

	/**

	 * Crea un nuevo unExpected Point.

	 *

	 * @param s string that identify the cause

	 */

	public UnexpectedPointException(String string){

		super(string);

	}

}

/* DielmoOpenLiDAR

 *

 * Copyright (C) 2008 DIELMO 3D S.L. (DIELMO) and 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.

 *

 * For more information, contact:

 *

 * DIELMO 3D S.L.

 * Plaza Vicente Andr?s Estell?s 1 Bajo E

 * 46950 Xirivella, Valencia

 * SPAIN

 *   

 * +34 963137212

 * dielmo@dielmo.com

 * www.dielmo.com

 * 

 * or

 * 

 * Generalitat Valenciana

 * Conselleria d'Infraestructures i Transport

 * Av. Blasco Ib??ez, 50

 * 46010 VALENCIA

 * SPAIN

 *

 * +34 963862235

 * gvsig@gva.es

 * www.gvsig.gva.es

 */

/*

 * AUTHORS (In addition to DIELMO and CIT):

 *  

 */

package com.dielmo.lidar;

import java.io.InputStream;

import java.nio.ByteBuffer;

public abstract class LASVariableLengthRecord_1X {	

	/**

	 * data after the variable length record.

	 * 

	 * Make sure that maximum size is 65536 bytes

	 */

	protected byte[] data;

	/**

	 * The userID field is ASCII character data that identifies the

	 * user which created the variable length record.

	 */

	protected char[] userID = new char[16];

	/**

	 * Record ID is dependent upon the User ID. There can be 0 to 65535

	 * record Ids for every UserID. The ASPRS standard will manage its own record Ids (User Ids owned by the

	 * specification), otherwise record Ids will be managed by the owner of the given User ID. So each	 

	 * User ID is allowed to assign 0 to 65535 record Ids as they wish. Publicizing the meaning of a

	 * given record ID will be left to the owner of the given User ID. Unknown User ID/Record ID

	 * combinations should be ignored.

	 */

	protected int recordID; 

	/**

	 * The record length is the number of bytes for the record after the

	 * end of the standard part of the header.

	 */

	protected int RecordLengthAfterHeader;

	/**

	 * Optional null terminated text description of the data.

	 */

	protected char[] description = new char[32]; 

	/**

	 * Default constructor, without arguments.

	 * Initializes all components to zero.

	 */ 

	public LASVariableLengthRecord_1X() {

		int i;

		for(i=0;i<16;i++)

			userID[i] = 0;

		recordID = 0;

		RecordLengthAfterHeader=0;

		for(i=0;i<32;i++)

			description[i] = 0;

	}

	/**

	 * Return the userID field as string that identifies the

	 * user which created the variable length record.

	 * 

	 * @return user ID

	 */

	public String getUserID() {	

		String s;

		char[] G = new char[16];

		int i;

		for(i=0;i<16;i++)

		{

			if(userID[i]!=0)

				G[i] = userID[i];

			else

				G[i] = ' ';

		}

		s = new String(G);

		return s;

	}

	/**

	 * set the userID field as string that identifies the

	 * user which created the variable length record.

	 * 

	 * Only accept 16 characters

	 * 

	 * @param u new userID

	 */

	public void setUserID(String u) {	

		int i;

		for(i=0;i<16 && i<u.length() ;i++){

			userID[i] = u.charAt(i);

		}

		// put to zero.

		while(i<16){

			userID[i] = 0;

			i++;

		}

	}

	/**

	 * Return record ID.

	 * 

	 * Record ID is dependent upon the User ID. There can be 0 to 65535 

	 * record Ids for every UserID. The ASPRS standard will manage its own

	 * record Ids (User Ids owned by the specification), otherwise record

	 * Ids will be managed by the owner of the given User ID. So each User

	 * ID is allowed to assign 0 to 65535 record Ids as they wish. 

	 * Publicizing the meaning of a given record ID will be left to the

	 * owner of the given User ID. Unknown User ID/Record ID combinations

	 * should be ignored. 

	 *  

	 * @return record ID

	 */

	public int getRecordID() {	

		return recordID;

	}

	/**

	 * set record ID.

	 * 

	 * Record ID is dependent upon the User ID. There can be 0 to 65535 

	 * record Ids for every UserID. The ASPRS standard will manage its own

	 * record Ids (User Ids owned by the specification), otherwise record

	 * Ids will be managed by the owner of the given User ID. So each User

	 * ID is allowed to assign 0 to 65535 record Ids as they wish. 

	 * Publicizing the meaning of a given record ID will be left to the

	 * owner of the given User ID. Unknown User ID/Record ID combinations

	 * should be ignored. 

	 *  

	 * @param r new record ID

	 */

	public void setRecordID( int r) {	

		try{

			if(r>=0 && r<=LidarHeader.UNSIGNED_SHORT_MAX)

				recordID=r;

			else

				throw new OutOfRangeLidarException("Out of range of record ID (Variable Length Record)");

		} catch(OutOfRangeLidarException e) {

			e.printStackTrace();

		}

	}

	/**

	 * Return record length after header. This is the number of bytes

	 * for the record after the end of the standard part of the header.

	 *  

	 * @return record length after header

	 */

	public int getRecordLengthAfterHeader() {	

		return RecordLengthAfterHeader;

	}

	/**

	 * set record length after header. This is the number of bytes

	 * for the record after the end of the standard part of the header.

	 *  

	 * @param r new record length after header

	 */

	public void setRecordLengthAfterHeader(int r) {	

		try{

			if(r>=0 && r<=LidarHeader.UNSIGNED_SHORT_MAX)

				RecordLengthAfterHeader = r;

			else

				throw new OutOfRangeLidarException("Out of range of record length after header (Variable Length Record)");

		} catch(OutOfRangeLidarException e) {

			e.printStackTrace();

		}

	}

	/**

	 * Return the text description of the data.

	 * 

	 * @return description

	 */

	public String getDescription() {	

		String s;

		char[] D = new char[32];

		int i;

		for(i=0;i<32;i++) {

			if(description[i]!=0)

				D[i] = description[i];

			else

				D[i] = ' ';

		}

		s = new String(D);

		return s;

	}

	/**

	 * set the text description of the data.

	 * Only accept 32 characters

	 * 

	 * @param d new description

	 */

	public void setDescription(String d) {	

		int i;

		for(i=0;i<32 && i<d.length();i++) {

			description[i] = d.charAt(i);

		}

		// put to zero.

		while(i<32){

			description[i] = 0;

			i++;

		}

	}

	/**

	 * Read the collection of bytes of variable length record 

	 * 

	 * @param input input file to read

	 * @param VarLengthRecord buffer into which the data is read 

	 * @return true if success else return false 

	 */

	protected abstract boolean readVarLegthRecordData(InputStream input);

	/**

	 * Read the Variable length record header of LAS file

	 * 

	 * @param input input file to read

	 * @return true if success else return false 

	 */

	public abstract boolean readVarLegthRecord(InputStream input);

	/**

	 * Write the Variable length record header of LAS file

	 * 

	 * @param bb byte buffer to write

	 * @return true if success else return false 

	 */

	public abstract boolean writeVarLegthRecord(ByteBuffer bb);

	/**

	 * return data after the variable length record.

	 * 

	 * maximum size is 65536 bytes

	 * 

	 * @return data after variable length record.

	 */

	public byte[] getData() {

		return data;

	}

	/**

	 * set data after the variable length record.

	 * 

	 * Make sure that maximum size is 65536 bytes

	 * 

	 * @return data after variable length record.

	 */

	public void setData(byte[] data) {

		this.data = data;

	}

}

/* DielmoOpenLiDAR

 *

 * Copyright (C) 2008 DIELMO 3D S.L. (DIELMO) and 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.

 *

 * For more information, contact:

 *

 * DIELMO 3D S.L.

 * Plaza Vicente Andr?s Estell?s 1 Bajo E

 * 46950 Xirivella, Valencia

 * SPAIN

 *   

 * +34 963137212

 * dielmo@dielmo.com

 * www.dielmo.com

 * 

 * or

 * 

 * Generalitat Valenciana

 * Conselleria d'Infraestructures i Transport

 * Av. Blasco Ib??ez, 50

 * 46010 VALENCIA

 * SPAIN

 *

 * +34 963862235

 * gvsig@gva.es

 * www.gvsig.gva.es

 */

/*

 * AUTHORS (In addition to DIELMO and CIT):

 *  

 */

package com.dielmo.lidar;

/**

 * Static class to work with a given vector of bytes.

 *

 * @author Oscar Garcia

 */

public class ByteUtilities {

	// operaciones de vectores de bytes a tipos de java

	/**

	 * From an array of bytes (byte[]), takes over from the position given as start

	 * and gives value to int variable using the following 4 bytes.

	 * 

	 * @param arr array of bytes

	 * @param start initial position

	 */

	public static int arr2Int (byte[] arr, int start) {

		int i = 0;

		int length = 4;

		int cnt = 0;

		byte[] tmp = new byte[length];

		for (i = start; i < (start + length); i++) {

			tmp[cnt] = arr[i];

			//System.out.println(java.lang.Byte.toString(arr[i]) + " " + i);

			cnt++;

		}

		int accum = 0;

		i = 0;

		for ( int shiftBy = 0; shiftBy < 32; shiftBy += 8 ) {

			accum |= ( (int)( tmp[i] & 0xff ) ) << shiftBy;

			i++;

		}

		return accum;

	}

	/**

	 * From an array of bytes (byte[]), takes over from the position given as start

	 * and gives value to double variable using the following 8 bytes.

	 * 

	 * @param arr array of bytes

	 * @param start initial position

	 */

	public static double arr2Double (byte[] arr, int start) {

		int i = 0;

		int length = 8;

		int cnt = 0;

		byte[] tmp = new byte[length];

		for (i = start; i < (start + length); i++) {

			tmp[cnt] = arr[i];

			//System.out.println(java.lang.Byte.toString(arr[i]) + " " + i);

			cnt++;

		}

		long accum = 0;

		i = 0;

		for ( int shiftBy = 0; shiftBy < 64; shiftBy += 8 ) {

			accum |= ( (long)( tmp[i] & 0xff ) ) << shiftBy;

			i++;

		}

		return Double.longBitsToDouble(accum);

	}

	/**

	 * From an array of bytes (byte[]), takes over from the position given as start

	 * and gives value to long variable using the following 4 bytes.

	 * 

	 * @param arr array of bytes

	 * @param start initial position

	 */

	public static long arr2UnsignedInt (byte[] arr, int start) {

		int i = 0;

		int length = 4;

		int cnt = 0;

		byte[] tmp = new byte[length];

		for (i = start; i < (start + length); i++) {

			tmp[cnt] = arr[i];

			//System.out.println(java.lang.Byte.toString(arr[i]) + " " + i);

			cnt++;

		}

		long accum = 0;

		i = 0;

		for ( int shiftBy = 0; shiftBy < 32; shiftBy += 8 ) {

			accum |= ( (long)( tmp[i] & 0xff ) ) << shiftBy;

			i++;

		}

		return accum;

	}

	/**

	 * From an array of bytes (byte[]), takes over from the position given as start

	 * and gives value to long variable using the following 8 bytes.

	 * 

	 * @param arr array of bytes

	 * @param start initial position

	 */

	public static long arr2Long(byte[] arr, int start) {

		int i = 0;

		int length = 8;

		int cnt = 0;

		byte[] tmp = new byte[length];

		for (i = start; i < (start + length); i++) {

			tmp[cnt] = arr[i];

			//System.out.println(java.lang.Byte.toString(arr[i]) + " " + i);

			cnt++;

		}

		long accum = 0;

		i = 0;

		for ( int shiftBy = 0; shiftBy < 64; shiftBy += 8 ) {

			accum |= ( (long)( tmp[i] & 0xff ) ) << shiftBy;

			i++;

		}

		return accum;

	}

	/**

	 * From an array of bytes (byte[]), takes over from the position given as start

	 * and gives value to int variable using the following 2 bytes.

	 * 

	 * @param arr array of bytes

	 * @param start initial position

	 */

	public static int arr2Unsignedshort (byte[] arr, int start) {

		int i = 0;

		int length = 2;

		int cnt = 0;

		byte[] tmp = new byte[length];

		for (i = start; i < (start + length); i++) {

			tmp[cnt] = arr[i];

			//System.out.println(java.lang.Byte.toString(arr[i]) + " " + i);

			cnt++;

		}

		int accum = 0;

		i = 0;

		for ( int shiftBy = 0; shiftBy < 16; shiftBy += 8 ) {

			accum |= ( (int)( tmp[i] & 0xff ) ) << shiftBy;

			i++;

		}

		return accum;

	}	

	/**

	 * Write the bytes of "var" into new byte array.

	 *

	 * @param var the int to encode

	 * @param arrayBytes The byte array to store into.

	 * @param startIndex index data to begin write.

	 */

	public static void int2Arr(int var, byte[] arrayBytes, int startIndex) {

		int length = 4;

		if (arrayBytes != null && startIndex+length <= arrayBytes.length) {

			for (int j = startIndex; j < startIndex+length; j++) {

				arrayBytes[j] = (byte) var; // se copian los 8 primeros bits de la variable

				var >>= 8;

			}

		}

	}

	/**

	 * Write the bytes of "var" into new byte array.

	 *

	 * @param var the int that represents a unsigned short of 2 bytes to encode

	 * @param arrayBytes The byte array to store into.

	 * @param startIndex index data to begin write.

	 */

	public static void unsignedShort2Arr(int var, byte[] arrayBytes, int startIndex) {

		int length = 2;

		if (arrayBytes != null && startIndex+length <= arrayBytes.length) {

			for (int j = startIndex; j < startIndex+length; j++) {

				arrayBytes[j] = (byte) var; // se copian los 8 primeros bits de la variable

				var >>= 8;

			}

		}

	}

	/**

	 * Write the bytes of "var" into new byte array.

	 *

	 * @param var the long that represents a unsigned int of 4 bytes to encode

	 * @param arrayBytes The byte array to store into.

	 * @param startIndex index data to begin write.

	 */

	public static void unsignedInt2Arr(long var, byte[] arrayBytes, int startIndex) {

		int length = 4;

		if (arrayBytes != null && startIndex+length <= arrayBytes.length) {

			for (int j = startIndex; j < startIndex+length; j++) {

				arrayBytes[j] = (byte) var; // se copian los 8 primeros bits de la variable

				var >>= 8;

			}

		}

	}

	/**

	 * Write the bytes of "var" into new byte array.

	 *

	 * @param var the long of 8 bytes to encode

	 * return tmp The byte array to store into.

	 */

	public static void long2Arr(long var, byte[] arrayBytes, int startIndex) {

		int length = 8;

		if (arrayBytes != null && startIndex+length <= arrayBytes.length) {

			for (int j = startIndex; j < startIndex+length; j++) {

				arrayBytes[j] = (byte) var; // se copian los 8 primeros bits de la variable

				var >>= 8;

			}

		}

	}

	/**

	 * Write the bytes of "var" into new byte array.

	 *

	 * @param var the double of 8 bytes to encode

	 * return tmp The byte array.

	 */

	public static void double2Arr(double var, byte[] arrayBytes, int startIndex) {

		long bits = Double.doubleToLongBits(var);

		long2Arr(bits, arrayBytes, startIndex);

	}

}

/* DielmoOpenLiDAR

 *

 * Copyright (C) 2008 DIELMO 3D S.L. (DIELMO) and 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.

 *

 * For more information, contact:

 *

 * DIELMO 3D S.L.

 * Plaza Vicente Andr?s Estell?s 1 Bajo E

 * 46950 Xirivella, Valencia

 * SPAIN

 *   

 * +34 963137212

 * dielmo@dielmo.com

 * www.dielmo.com

 * 

 * or

 * 

 * Generalitat Valenciana

 * Conselleria d'Infraestructures i Transport

 * Av. Blasco Ib??ez, 50

 * 46010 VALENCIA

 * SPAIN

 *

 * +34 963862235

 * gvsig@gva.es

 * www.gvsig.gva.es

 */

/*

 * AUTHORS (In addition to DIELMO and CIT):

 *  

 */

package com.dielmo.lidar;

import java.io.File;

public abstract class LASHeader_1X implements LidarHeader{

	/**

	 * path of LAS file 

	 */

	protected File m_Fich;

	/**

	 * file Signature, must contain LASF

	 */

	protected char[] fileSignature = new char[4];

	/** 

	 * Global Unique IDentifier 1; unsigned int

	 */

	protected long GUID1 ; 

	/**

	 * Global Unique IDentifier 2; unsigned short

	 */

	protected int GUID2 ;

	/**

	 * Global Unique IDentifier 3; unsigned short

	 */

	protected int GUID3 ; 

	/**

	 * Global Unique IDentifier 4

	 */

	protected char[]GUID4 = new char[8];

	/**

	 * version format number mayor

	 */

	protected char versionMayor; 

	/**

	 * version format number minor

	 */

	protected char versionMinor;

	/**

	 * software generated by: Original Data, Merge, Modification, Extraction, Transformation or Other

	 */

	protected char[] systemIDentifier = new char[32]; 

	/**

	 * Who generate software

	 */

	protected char[] generatingSoftware = new char[32];

	/**

	 * LAS 1.0 Julian day

	 * LAS 1.1 day of year 1-365

	 */

	protected int day; 

	/**

	 * In LAS 1.0 year that the data was collected.

	 * In LAS 1.1 year on which this file was created.

	 */

	protected int year;

	/**

	 * Size, in bytes, of the header file 

	 */

	protected int hdrSize;

	/**

	 * Number of bytes from the beginning of the file

	 * to the first point record data (including the two 

	 * bytes for data start signature in case of LAS1.0).

	 */

	protected long offsetData;

	/**

	 * number of variable length records

	 */

	protected long numVarLengthRecord;

	/**

	 * point data record format type (0-1)

	 */

	protected char pointDataFormatID;

	/**

	 * size of the point data record

	 */

	protected int pointDataRecordLength;

	/**

	 * total number of points records within the file

	 */

	protected long numPointsRecord;

	/**

	 * array of the total point records per return

	 */

	protected long[] numPointsByReturn = new long[5];

	/**

	 * X Scale factor for scale corresponding X values within the point

	 * record. 

	 */

	protected double xScale;

	/**

	 * Y Scale factor for scale corresponding Y values within the point

	 * record. 

	 */

	protected double yScale;

	/**

	 * Z Scale factor for scale corresponding Z values within the point

	 * record. 

	 */

	protected double zScale;

	/**

	 * X Offset used to set the overall offset for the point records.

	 */

	protected double xOffset;

	/**

	 * Y Offset used to set the overall offset for the point records.

	 */

	protected double yOffset;

	/**

	 * Z Offset used to set the overall offset for the point records.

	 */

	protected double zOffset;

	/**

	 * The actual maximum X of coordinate extents of the file

	 */

	protected double maxX;

	/**

	 * The actual minimum X of coordinate extents of the file

	 */

	protected double minX;

	/**

	 * The actual maximum Y of coordinate extents of the file

	 */

	protected double maxY;

	/**

	 * The actual minimum Y of coordinate extents of the file

	 */

	protected double minY;

	/**

	 * The actual maximum Z of coordinate extents of the file

	 */

	protected double maxZ;

	/**

	 * The actual minimum Z of coordinate extents of the file

	 */

	protected double minZ;

	/**

	 * Default constructor, without arguments.

	 * Initializes all components to zero.

	 */ 

	public LASHeader_1X(File file) {

		m_Fich = file;

		int i;

		fileSignature[0] = 'L';

		fileSignature[1] = 'A';

		fileSignature[2] = 'S';

		fileSignature[3] = 'F';

		GUID1 = 0;

		GUID2 = 0;

		GUID3 = 0;

		for(i=0;i<8;i++)

			GUID4[i] = 0;

		versionMayor = 0;

		versionMinor = 0;

		for(i=0;i<32;i++)

			systemIDentifier[i] = 0;

		for(i=0;i<32;i++)

			generatingSoftware[i] = 0;

		generatingSoftware[0] = 'D';

		generatingSoftware[1] = 'I';

		generatingSoftware[2] = 'E';

		generatingSoftware[3] = 'L';

		generatingSoftware[4] = 'M';

		generatingSoftware[5] = '0';

		generatingSoftware[6] = ' ';

		generatingSoftware[7] = '3';

		generatingSoftware[8] = 'D';

		generatingSoftware[9] = ' ';

		generatingSoftware[10] = 'S';

		generatingSoftware[11] = '.';

		generatingSoftware[12] = 'L';

		generatingSoftware[13] = '.';

		day = 0;

		year = 0;

		hdrSize=227;

		offsetData=0;

		numVarLengthRecord=0;

		pointDataFormatID=0;

		pointDataRecordLength=0;

		numPointsRecord=0;

		for(i=0;i<5;i++)

			numPointsByReturn[i] = 0;

		xScale=0;

		yScale=0;

		zScale=0;

		xOffset=0;

		yOffset=0;

		zOffset=0;

		maxX=0;

		minX=0;

		maxY=0;

		minY=0;

		maxZ=0;

		minZ=0;

	}

	// GET MEHODS

	/**

	 * Return file signature. This must be always "LASF"

	 * 

	 * @return String with file signature

	 */

	public String getFileSignature() {

		String f = new String(fileSignature);

		// comprobamos los datos y si son nulos se ponen a cero 		

		return f;

	}

	/**

	 * Return Global unique identifier 1 (GUID1)

	 * 

	 * @return GUID1

	 */

	public long getGUID1() {

		return GUID1;

	}

	/**

	 * Return Global unique identifier 2 (GUID2)

	 * 

	 * @return GUID2

	 */