[Liblas-devel] Bug?: problem copying file with offset

Peder Axensten Peder.Axensten at slu.se
Sun Mar 20 11:57:31 EDT 2011


I have a las-file with non-zero offsets (attached below).

If I used the code below to copy this file to a las-file I "lose" the offset. 
However, if I copy to a compressed laz-file all is well.
If I add a line printing the coordinates, they look normal in both cases (but still lose their offset in the las-case).

The first plots in the las-file becomes:
-3366.35,749.70,-62.21
-3366.34,750.40,-60.36
-3366.46,751.90,-62.03

Whereas in the laz-file they (correctly) are: 
419059.89,6483734.30,106.79
419059.89,6483734.37,106.97
419059.88,6483734.52,106.81

Here is the code (compress_file(out) returns true if the file extension is 'laz'): 
	std::size_t copy (
		const std::vector< std::string > & in,
		const std::string & out
	) {
		std::size_t					count = 0;
		try {
			std::ofstream			ofs( out.c_str(), std::ios::out | std::ios::binary );
			liblas::Header			header;
			header.SetCompressed( compress_file( out ) );
			liblas::Writer			writer( ofs, header );
			
			for( std::size_t i=0; i<in.size(); i++ ) {
				std::ifstream			ifs( in[i].c_str(), std::ios::in | std::ios::binary );
				liblas::ReaderFactory	rf;
				liblas::Reader			reader = rf.CreateWithStream( ifs );
		
				while( reader.ReadNextPoint() ) {
					++count;
					writer.WritePoint( reader.GetPoint() );
				}
			}
		} catch (std::exception const& e) {
			std::cerr << "Error: " << e.what() << std::endl;
		}
		return( count );
	}


Here is the las-file:

-------------- next part --------------
A non-text attachment was scrubbed...
Name: offset_test_in.las
Type: application/octet-stream
Size: 5605 bytes
Desc: offset_test_in.las
Url : http://lists.osgeo.org/pipermail/liblas-devel/attachments/20110320/5068355a/offset_test_in.obj
-------------- next part --------------


Have I misunderstood something?
Any help appreciated!

Best regards,

Peder Axensten
Research engineer

Swedish University of Agricultural Sciences
Dept. of Forest Resource Management
Forest Remote Sensing
Skogmarksgränd
se-90183 Umeå
Sweden

Phone: +46-90-786.8500
http://www.srh.slu.se/rsl



More information about the Liblas-devel mailing list