geosolutions-it / digital-twin-toolbox

Open source solution for inspecting and generating 3D Tiles for urban environments
GNU General Public License v3.0
46 stars 7 forks source link

Point cloud tutorial: Couldn't reproject point #24

Open bertt opened 4 months ago

bertt commented 4 months ago

I'm trying to follow tutorial 'Processing of point cloud data, vertical datum shift and colorization' https://github.com/geosolutions-it/digital-twin-toolbox/wiki/%5BTutorial%5D-Processing-of-point-cloud-data,-vertical-datum-shift-and-colorization, but get an error about projections after press 'Start processing':

digital_twin_toolbox          | micromamba run -n pdal_env pdal info static/data/1916.las --metadata
digital_twin_toolbox          | Horizontal WKT:<br><br>"PROJCS[\"NAD83 / Maryland\",GEOGCS[\"NAD83\",DATUM[\"North_American_Datum_1983\",SPHEROID[\"GRS 1980\",6378137,298.257222101,AUTHORITY[\"EPSG\",\"7019\"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY[\"EPSG\",\"6269\"]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],UNIT[\"degree\",0.0174532925199433,AUTHORITY[\"EPSG\",\"9122\"]],AUTHORITY[\"EPSG\",\"4269\"]],PROJECTION[\"Lambert_Conformal_Conic_2SP\"],PARAMETER[\"standard_parallel_1\",39.45],PARAMETER[\"standard_parallel_2\",38.3],PARAMETER[\"latitude_of_origin\",37.6666666666667],PARAMETER[\"central_meridian\",-77],PARAMETER[\"false_easting\",400000],PARAMETER[\"false_northing\",0],UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]],AXIS[\"X\",EAST],AXIS[\"Y\",NORTH],AUTHORITY[\"EPSG\",\"26985\"]]"
digital_twin_toolbox          | Vertical WKT:<br><br>"VERT_CS[\"NAVD88 height\",VERT_DATUM[\"North American Vertical Datum 1988\",2005,AUTHORITY[\"EPSG\",\"5103\"]],UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]],AXIS[\"Gravity-related height\",UP],AUTHORITY[\"EPSG\",\"5703\"]]"
digital_twin_toolbox          | [
digital_twin_toolbox          |   { filename: 'static/data/1916.las', type: 'readers.las' },
digital_twin_toolbox          |   { type: 'filters.colorization', raster: 'static/data/1916.tif' },
digital_twin_toolbox          |   {
digital_twin_toolbox          |     type: 'filters.reprojection',
digital_twin_toolbox          |     in_srs: '+init=EPSG:26985 +geoidgrids=./vdatum/egm08_25.gtx',
digital_twin_toolbox          |     out_srs: 'EPSG:26985+4269',
digital_twin_toolbox          |     error_on_failure: true
digital_twin_toolbox          |   },
digital_twin_toolbox          |   {
digital_twin_toolbox          |     type: 'writers.las',
digital_twin_toolbox          |     a_srs: 'EPSG:26985+4269',
digital_twin_toolbox          |     filename: './static/data/1916_processed.las',
digital_twin_toolbox          |     compression: true
digital_twin_toolbox          |   }
digital_twin_toolbox          | ]
digital_twin_toolbox          | Init point cloud processing...
digital_twin_toolbox          | micromamba run -n pdal_env pdal pipeline ./tmp/pipeline-80dd4010-32ce-11ef-b451-2b84e210b9b6.json
digital_twin_toolbox          | (pdal pipeline filters.reprojection Error) GDAL failure (1) No operation matching criteria found for coordinate
digital_twin_toolbox          |
digital_twin_toolbox          | PDAL: filters.reprojection: Couldn't reproject point with X/Y/Z coordinates of (395915.850000, 136883.980000, 41.470000).
digital_twin_toolbox          |
digital_twin_toolbox          |
digital_twin_toolbox          | Process completed!

was able to get it working using EPSG:26985+5703 as 'From CRS':

digital_twin_toolbox          | [
digital_twin_toolbox          |   { filename: 'static/data/1916.las', type: 'readers.las' },
digital_twin_toolbox          |   { type: 'filters.colorization', raster: 'static/data/1916.tif' },
digital_twin_toolbox          |   {
digital_twin_toolbox          |     type: 'filters.reprojection',
digital_twin_toolbox          |     in_srs: 'EPSG:26985+5703',
digital_twin_toolbox          |     out_srs: 'EPSG:26985+4269',
digital_twin_toolbox          |     error_on_failure: true
digital_twin_toolbox          |   },
digital_twin_toolbox          |   {
digital_twin_toolbox          |     type: 'writers.las',
digital_twin_toolbox          |     a_srs: 'EPSG:26985+4269',
digital_twin_toolbox          |     filename: './static/data/1916_processed.las',
digital_twin_toolbox          |     compression: true
digital_twin_toolbox          |   }
digital_twin_toolbox          | ]
digital_twin_toolbox          | Init point cloud processing...
digital_twin_toolbox          | micromamba run -n pdal_env pdal pipeline ./tmp/pipeline-83d86230-32cf-11ef-b451-2b84e210b9b6.json
digital_twin_toolbox          | Process completed!
digital_twin_toolbox          | Read point cloud statistic...
digital_twin_toolbox          | micromamba run -n pdal_env pdal info static/data/1916_processed.las --stats
digital_twin_toolbox          | Read point cloud metadata...
digital_twin_toolbox          | micromamba run -n pdal_env pdal info static/data/1916_processed.las --metadata
digital_twin_toolbox          | micromamba run -n pdal_env pdal pipeline ./tmp/pipeline-ce249de0-32cf-11ef-b451-2b84e210b9b6.json
tdipisa commented 2 months ago

Dear @bertt sorry for the late feedback. We have recently made significant changes and updates to the master branch of the project for which you may want to check again if the problem you previously identified still persists. We are looking forward to hearing from you.