<div dir="ltr"><div>I am producing DEMs from LiDAR data using gdal_grid. All processes compete fine, but the resulting DEM does not load into QGIS 2.18 on my Ubuntu 16.04 machine but they load fine on my Windows machines (same version of QGIS). When I look at layer properties on Ubuntu machine, the overviews that are produced are not recognized and when I try to rebuild within QGIS, I get an error that I cannot produce overviews for this file type -- it's a Geotif. <br></div><div><br></div><div>DEM is produced using the following:</div><div><br></div><div>gdal_grid -zfield field_3 -a nearest:radius1=20:radius2=20 -a_srs epsg:26916 -txe 607900.7    648882.51 -tye 5322415.07 5349297.05 -outsize 40981.81 26881.98 -co "BIGTIFF=YES" -co "COMPRESS=LZW" --config -a_nodata 0 -l ground2 ground2.vrt ground2_grid_1m_nearest_radius20.tif --config GDAL_NUM_THREADS ALL_CPUS --config GDAL_CACHEMAX 2000</div><div><br></div><div>I then turn 0s into NoData with Gdal_calc</div><div><br></div><div>gdal_calc.py --co="COMPRESS=LZW" -A ground2_grid_1m_nearest_radius20.tif --outfile=ground2_grid_1m_nearest_radius20_nodata0.tif --calc="A*(A>0)" --NoDataValue=0</div><div><br></div><div>Then I build overviews with gdaladdo</div><div><br></div><div>gdaladdo -clean --config COMPRESS_OVERVIEW LZW --config INTERLEAVE_OVERVIEW PIXEL -r average ground2_grid_1m_nearest_radius20_nodata0.tif 2 4 8 16 32 64 128 256 512 1024</div><div><br></div><div>Gdalinfo says....</div><div><br></div><div>Driver: GTiff/GeoTIFF<br>Files: ground2_grid_1m_nearest_radius20_nodata0.tif<br>Size is 40981, 26881<br>Coordinate System is:<br>PROJCS["NAD83 / UTM zone 16N",<br>    GEOGCS["NAD83",<br>        DATUM["North_American_Datum_1983",<br>            SPHEROID["GRS 1980",6378137,298.257222101,<br>                AUTHORITY["EPSG","7019"]],<br>            TOWGS84[0,0,0,0,0,0,0],<br>            AUTHORITY["EPSG","6269"]],<br>        PRIMEM["Greenwich",0,<br>            AUTHORITY["EPSG","8901"]],<br>        UNIT["degree",0.0174532925199433,<br>            AUTHORITY["EPSG","9122"]],<br>        AUTHORITY["EPSG","4269"]],<br>    PROJECTION["Transverse_Mercator"],<br>    PARAMETER["latitude_of_origin",0],<br>    PARAMETER["central_meridian",-87],<br>    PARAMETER["scale_factor",0.9996],<br>    PARAMETER["false_easting",500000],<br>    PARAMETER["false_northing",0],<br>    UNIT["metre",1,<br>        AUTHORITY["EPSG","9001"]],<br>    AXIS["Easting",EAST],<br>    AXIS["Northing",NORTH],<br>    AUTHORITY["EPSG","26916"]]<br>Origin = (607900.699999999953434,5322415.070000000298023)<br>Pixel Size = (1.000019765257072,1.000036456977029)<br>Metadata:<br>  AREA_OR_POINT=Area<br>Image Structure Metadata:<br>  COMPRESSION=LZW<br>  INTERLEAVE=BAND<br>Corner Coordinates:<br>Upper Left  (  607900.700, 5322415.070) ( 85d33' 8.00"W, 48d 2'45.22"N)<br>Lower Left  (  607900.700, 5349297.050) ( 85d32'43.44"W, 48d17'15.62"N)<br>Upper Right (  648882.510, 5322415.070) ( 85d 0' 9.50"W, 48d 2'15.55"N)<br>Lower Right (  648882.510, 5349297.050) ( 84d59'35.63"W, 48d16'45.69"N)<br>Center      (  628391.605, 5335856.060) ( 85d16'24.08"W, 48d 9'46.71"N)<br>Band 1 Block=40981x1 Type=Float64, ColorInterp=Gray<br>  NoData Value=0<br>  Overviews: 20491x13441, 10246x6721, 5123x3361, 2562x1681, 1281x841, 641x421, 321x211, 161x106, 81x53, 41x27<br></div><div><br></div><div>So as you can see, overviews are produced and the raster is comrpessed with LZW. But... it will not load in QGIS and causes it to crash (on Ubuntu). Loads great on Widows... (boo!)</div><div><br></div><div>Thoughts? <br></div><div><br></div><div>Cliff<br></div></div>