MG4 Lidar Driver contribution

[mrosen]

Per earlier discussion, attached find an archive that integrates the ability to read MG4 Lidar point cloud data and expose it as a Raster. It depends on the current, freely-available-though-not-open-source MG4 Lidar SDK v1.1 (​http://www.lizardtech.com/developer/members/

I have tested the build procedure against gdal 7.1.

[warmerdam]

Source code, and sample data available at:

​http://home.gdal.org/tmp/GDAL_MG4Lidar_Src.zip

Thanks Mike, I'll have someone start integration in the coming weeks.

[warmerdam]

Ivan,

Could you work on integrating this please?

[ilucena]

OK Frank.

[ilucena]

Driver added on revision r19355

http://trac.osgeo.org/gdal/changeset/19355

Examples of use:

% gdalinfo --format MG4Lidar
Format Details:
  Short Name: MG4Lidar
  Long Name: MrSID Generation 4 / Lidar (.sid)
  Extension: view
  Help Topic: frmt_mrsid_lidar.html
  Creation Datatypes: Float64

% gdalinfo Tetons_200k.view
Driver: MG4Lidar/MrSID Generation 4 / Lidar (.sid)
Files: Tetons_200k.view
Size is 457, 431
Coordinate System is:
PROJCS["NAD83 / UTM zone 12N",
    GEOGCS["NAD83",
        DATUM["North_American_Datum_1983",
            SPHEROID["GRS 1980",6378137,298.2572221010002,
                AUTHORITY["EPSG","7019"]],
            AUTHORITY["EPSG","6269"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4269"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",-111],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","26912"]]
Origin = (504489.919999999980000,4795848.389999999700000)
Pixel Size = (3.078227571115974,-3.078259860787739)
Corner Coordinates:
Upper Left  (  504489.920, 4795848.390)
Lower Left  (  504489.920, 4794521.660)
Upper Right (  505896.670, 4795848.390)
Lower Right (  505896.670, 4794521.660)
Center      (  505193.295, 4795185.025)
Band 1 Block=457x431 Type=Float32, ColorInterp=Undefined
  Description = Mean of Z (filtered by Classification)
  Minimum=2402.990, Maximum=2935.860, Mean=2681.271, StdDev=125.189
  NoData Value=3.4028234663852886e+038
  Overviews: 229x216

% gdal_translate Tetons_200k.view Tetons_200k.tif
Input file size is 457, 431
0...10...20...30...40...50...60...70...80...90...100 - done.

% gdalinfo Tetons_200k.tif
Driver: GTiff/GeoTIFF
Files: \data\Tetons_200k.tif
Size is 457, 431
Coordinate System is:
PROJCS["NAD83 / UTM zone 12N",
    GEOGCS["NAD83",
        DATUM["North_American_Datum_1983",
            SPHEROID["GRS 1980",6378137,298.2572221010002,
                AUTHORITY["EPSG","7019"]],
            AUTHORITY["EPSG","6269"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4269"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",-111],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","26912"]]
Origin = (504489.919999999980000,4795848.389999999700000)
Pixel Size = (3.078227571115974,-3.078259860787739)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  504489.920, 4795848.390)
Lower Left  (  504489.920, 4794521.660)
Upper Right (  505896.670, 4795848.390)
Lower Right (  505896.670, 4794521.660)
Center      (  505193.295, 4795185.025)
Band 1 Block=457x4 Type=Float32, ColorInterp=Gray
  NoData Value=3.4028234663852886e+038
  Metadata:
    LAYER_TYPE=athematic

Everything looks fine from what I can see. I will try in Linux now.