Hi,<br><br>Howard asked me to write up use case for using liblas in gdal.&nbsp; Here&#39;s some background - we currently use gdal to read coverages and write them to different formats, read metadata if its in the header, and to subset and resample images that we use in our WMS and WCS.<br>
<br>The use case for liblas and gdal is that gdal would be able to read, write, and subset las files in its native format.&nbsp; Gdal should be able to also read and write metadata from a las file.&nbsp; Gdal should be be able to convert a las file to gridded coverage such as geotiff to support a preview of the data.<br>
<br>Here some specifics on how gdal should behave for our use case:<br>
<br>
The program should do the following:<br>
<br>- provides metadata about the file, eg:<br>
<br>
Source=E:/data/Coverages/02pkbm2b.i42 <br>
&nbsp;&nbsp; Driver=NITF <br>
&nbsp;&nbsp; PixelSize=(1536,1536) <br>
&nbsp;&nbsp; Srs=EPSG:4326 <br>
&nbsp;&nbsp; Box_UpperLeft=(68.948223,34.526854) <br>
&nbsp;&nbsp; Box_LowerLeft=(68.948223,34.513043) <br>
&nbsp;&nbsp; Box_UpperRight=(68.966497,34.526854) <br>
&nbsp;&nbsp; Box_LowerRight=(68.966497,34.513043) <br>
&nbsp;&nbsp; NBands=1 <br>
&nbsp;&nbsp; DataType=Byte <br>
&nbsp;&nbsp; HasMetadata=false 129 <br>
<br>
or<br>
<br>
Source=E:/Data/Coverages/Imagery.tif <br>
&nbsp;&nbsp; Driver=GeoTIFF <br>
&nbsp;&nbsp; PixelSize=(6000,6000) <br>
&nbsp;&nbsp; Srs=EPSG:4326 <br>
&nbsp;&nbsp; GCP_UpperLeft=(1,1);(94.905,6.03351,0) <br>
&nbsp;&nbsp; GCP_UpperRight=(6000,1);(95.5752,5.88181,0) <br>
&nbsp;&nbsp; GCP_LowerLeft=(6000,6000);(95.4525,5.35223,0) <br>
&nbsp;&nbsp; GCP_LowerRight=(1,6000);(94.7828,5.5035,0) <br>
&nbsp;&nbsp; NBands=4 <br>
&nbsp;&nbsp; DataType=Byte,Byte,Byte,Byte <br>
&nbsp;&nbsp; HasMetadata=true <br>
&nbsp;&nbsp; TIFFTAG_IMAGEDESCRIPTION=XS3 XS2 XS1 SWIR <br>
&nbsp;&nbsp; TIFFTAG_DATETIME=2005:01:06 11:03:20<br>
<br>
note that SRS can be in either EPSG or WKT<br>
<br>
- able to write all or part of the file and write in the native file format<br>
<br>
The software is able to subset a file using a bounding box. When writing in the native format, interpolation is not allowed.<br>
<br>
- able to write all or part of the file to a DEM<br>
<br>
The software is able to subset a file using a bounding box. When writing to a DEM, interpolation is not allowed.<br>
<br>
able to write all or part of the file to a geotiff<br>
<br>
The software is able to subset a file using a bounding box. When writing to a geotiff, interpolation is allowed.<br>
<br><br>sophia<br>