mercredi 28 janvier 2015

Georeferencing using Liblas error


I am trying to write georeferencing information to a '.las' file using liblas and c++. Here is my code:



ofstream ofs;
ofs.open("test.las", ios::out | ios::binary);
liblas::Header header;
liblas::Writer writer(ofs, header);

liblas::Point point(&header);
point.SetCoordinates(1.0, 2.0, 3.0);

writer.WritePoint(point);
header.SetPointRecordsCount(1);

liblas::SpatialReference srs;
srs.SetFromUserInput("EPSG:4326");
header.SetSRS(srs);
writer.SetHeader(header);
writer.WriteHeader();


But, when I check the file using lasinfo tool I get:



Spatial Reference: None


If I check proj4, the information is stored, so I don't know why lasinfo doesn't print it:



Code : std::cout << srs.GetProj4() << std::endl;
Output: +proj=longlat +datum=WGS84 +no_defs


However, when I use lasinfo tool with others las files not written by me I can see the WKT and geotiff output.


What am I doing wrong?





Aucun commentaire:

Enregistrer un commentaire