Terrain Tiler fails to product tileset

I’m trying to create a terrain tileset using v4.5.2 of the terrain-tiler.exe. Messages emitted by the tiler are not helpful in understanding the problem. The data is 1 meter lidar derived GeoTIFFs for Pennsylvania (275 GB zip files).

Command line:

C:\...\Software\Cesium\Cesium-ion-3D-Tiling-Pipeline-Windows-4.5.2\bin\terrain-tiler.exe --verbose --input-list H:\Tilesets\FileList.txt --height-reference MeanSeaLevel --output H:\Tilesets\PAMAP_LiDAR.terraindb > H:\Tilesets\PAMAP_verbose.txt 2>&1

Sampling of messages (duplicates removed for brevity):

e[0mIdentifying contiguous north-up raster sourcese[0m

e[0mStitching together contiguous north-up raster sourcese[0m

e[35mWarning: Warning 1: Too large for int: nane[0m
...
e[35mWarning: More than 1000 errors or warnings have been reported. No more will be reported from now.e[0m

e[0m13390 total raster sources.e[0m

e[33m0% complete.e[0m
e[0m(1 of 81) Adding 1.vrt to tileset. (Using 6352 sources: "38001640PAS_dem.tif", 

e[35mWarning: Raster contains invalid elevations that exceed max/min height defined on planet.e[0m

Sample dataset gdalinfo:

e[me[32me]9;8;"USERNAME"e\@e]9;8;"COMPUTERNAME"e\ e[92mJ:\PAMAP_LiDAR\NORTH\17002190PAN_deme[90m
e[90m>e[m e]9;12e\gdalinfo.exe 17002190PAN_dem.tif 
Driver: GTiff/GeoTIFF
Files: 17002190PAN_dem.tif
       17002190PAN_dem.aux
       17002190PAN_dem.rrd
       17002190PAN_dem.tif.aux.xml
Size is 3125, 3125
Coordinate System is:
PROJCRS["NAD83 / Pennsylvania North",
    BASEGEOGCRS["NAD83",
        DATUM["North American Datum 1983",
            ELLIPSOID["GRS 1980",6378137,298.257222101,
                LENGTHUNIT["metre",1]]],
        PRIMEM["Greenwich",0,
            ANGLEUNIT["degree",0.0174532925199433]],
        ID["EPSG",4269]],
    CONVERSION["unknown",
        METHOD["Lambert Conic Conformal (2SP)",
            ID["EPSG",9802]],
        PARAMETER["Latitude of false origin",40.1666666666667,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8821]],
        PARAMETER["Longitude of false origin",-77.75,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8822]],
        PARAMETER["Latitude of 1st standard parallel",41.95,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8823]],
        PARAMETER["Latitude of 2nd standard parallel",40.8833333333333,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8824]],
        PARAMETER["Easting at false origin",1968500,
            LENGTHUNIT["us_survey_feet",0.304800609601219],
            ID["EPSG",8826]],
        PARAMETER["Northing at false origin",0,
            LENGTHUNIT["us_survey_feet",0.304800609601219],
            ID["EPSG",8827]]],
    CS[Cartesian,2],
        AXIS["easting (X)",east,
            ORDER[1],
            LENGTHUNIT["us_survey_feet",0.304800609601219]],
        AXIS["northing (Y)",north,
            ORDER[2],
            LENGTHUNIT["us_survey_feet",0.304800609601219]]]
Data axis to CRS axis mapping: 1,2
Origin = (2190000.000000000000000,170000.000000000000000)
Pixel Size = (3.200000000000000,-3.200000000000000)
Metadata:
  AREA_OR_POINT=Area
  TIFFTAG_RESOLUTIONUNIT=1 (unitless)
  TIFFTAG_SOFTWARE=IMAGINE TIFF Support
Copyright 1991 - 1999 by ERDAS, Inc. All Rights Reserved
@(#)$RCSfile: etif.c $ $Revision: 1.10.1.9.1.9.2.11 $ $Date: 2004/09/15 18:42:01EDT $
  TIFFTAG_XRESOLUTION=1
  TIFFTAG_YRESOLUTION=1
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  ( 2190000.000,  170000.000) ( 76d57' 7.31"W, 40d37'49.60"N)
Lower Left  ( 2190000.000,  160000.000) ( 76d57' 8.50"W, 40d36'10.80"N)
Upper Right ( 2200000.000,  170000.000) ( 76d54'57.63"W, 40d37'48.67"N)
Lower Right ( 2200000.000,  160000.000) ( 76d54'58.87"W, 40d36' 9.87"N)
Center      ( 2195000.000,  165000.000) ( 76d56' 3.08"W, 40d36'59.74"N)
Band 1 Block=64x64 Type=Float32, ColorInterp=Gray
  Min=382.690 Max=812.590 
  Minimum=382.690, Maximum=812.590, Mean=456.418, StdDev=89.927
  Overviews: 1563x1563, 782x782, 391x391, 196x196, 98x98, 49x49, 25x25
  Metadata:
    LAYER_TYPE=athematic
    STATISTICS_COVARIANCES=8086.79774158226
    STATISTICS_MAXIMUM=812.59002685547
    STATISTICS_MEAN=456.41761162347
    STATISTICS_MEDIAN=402.83999633789
    STATISTICS_MINIMUM=382.69000244141
    STATISTICS_MODE=384.79998779297
    STATISTICS_SKIPFACTORX=1
    STATISTICS_SKIPFACTORY=1
    STATISTICS_STDDEV=89.926624208753

Hello,

Thank you for providing those details. We have a theory as to what’s going wrong, but in order to confirm it we would need access to the source data to test internally. Would you be able to share it with us, through Google Drive or some other means?

Also, did the tiling job fail, or did it print those warnings but carry on?

Thanks,
Matt

It carried on and started to create a terraindb but those are all of the messages.

Sharing is difficult because the data is on an air-gapped network. Let me see if I can copy it to another drive so that it can be shared.

BTW, what’s your theory? Share privately if you like.

The data is available via FTP from ftp.pasda.psu.edu using anonymous login. You’ll need to copy the files from the directories in /pub/pasda/pamap/pamap_lidar/cycle1/DEM/{North,South}. They also provide mosaics for the counties which I’m preparing to tile today or tomorrow depending on how long it takes to unzip all of the files. I’m concerned with how the tiler will deal with overlapping datasets and possible NODATA values though.