[pdal] las to geotif

[Javier Jimenez Shaw]

I have seen point clouds in WGS84, and GeoTIFFs also in WGS84. Does it mean
that a point cloud in geographic coordinates is conceptually wrong? I find
it strange, but is it wrong? Why does PDAL care?

On Sat, 8 Jun 2024 at 03:17, Sergey Yefimov via pdal <pdal at lists.osgeo.org>
wrote:

> It helps. Thank you very much.
> ------------------------------
> *From:* Adam Steer <adam.d.steer at gmail.com>
> *Sent:* Friday, June 7, 2024 7:22 PM
> *To:* Sergey Yefimov <syefimov at ptfs.com>
> *Cc:* pdal at lists.osgeo.org <pdal at lists.osgeo.org>
> *Subject:* Re: [pdal] las to geotif
>
> THIS EMAIL IS FROM AN * UNTRUSTED* EMAIL PROVIDER.
>
> Hi Sergey - try reprojecting your source las file into a projected (in
> metres) system and running again.
>
> PDAL/GDAL pass projections through. Your input data are in degrees, your
> output resolution is interpreted as the same unit, so resolution 1.0 is
> interpreted as '1 degree'.
>
> Hope that helps,
>
> Adam
>
> --
> Dr Adam Steer
> https://spatialised.net
> Available for consulting/open for headhunting
>
> On Sat, Jun 8, 2024, 07:33 Sergey Yefimov via pdal <pdal at lists.osgeo.org>
> wrote:
>
> Hi,
> I try to convert las to geotiff using this pipeline
>
> "pipeline":[
> {"type" : "readers.las",
> "filename":"C:\\KVData\\source\\las\\test.las"
> },
> {
> "type":"writers.gdal",
> "filename":"C:/KVData/source/tiff/ttest.tif",
> "resolution":1.0,
> "output_type":"all",
> "gdaldriver":"GTiff"
> }
>     ]
>
> After conversion invalid 2 KB tiff file created. las info:
>
> {
>   "file_size": 137422487,
>   "filename": "C:\\KVData\\source\\las\\test.las",
>   "now": "2024-06-07T17:09:19-0400",
>   "pdal_version": "2.7.1 (git-version: Release)",
>   "reader": "readers.las",
>   "stats":
>   {
>     "bbox":
>     {
>       "EPSG:4326":
>       {
>         "bbox":
>         {
>           "maxx": -75.606374,
>           "maxy": 35.9592838,
>           "maxz": 29.71,
>           "minx": -75.6350255,
>           "miny": 35.942234,
>           "minz": -11.88
>         },
>         "boundary": { "type": "Polygon", "coordinates": [ [ [ -75.6350255,
> 35.942234, -11.88 ], [ -75.6350255, 35.9592838, -11.88 ], [ -75.606374,
> 35.9592838, 29.71 ], [ -75.606374, 35.942234, 29.71 ], [ -75.6350255,
> 35.942234, -11.88 ] ] ] }
>       }, ..
> "statistic":
>     [
>       {
>         "average": -75.62226764,
>         "count": 4907929,
>         "maximum": -75.606374,
>         "minimum": -75.6350255,
>         "name": "X",
>         "position": 0,
>         "stddev": 0.005315505055,
>         "variance": 2.825459399e-05
>       },
>       {
>         "average": 35.95136764,
>         "count": 4907929,
>         "maximum": 35.9592838,
>         "minimum": 35.942234,
>         "name": "Y",
>         "position": 1,
>         "stddev": 0.002706953159,
>         "variance": 7.327595404e-06
>       },..
> }
>
>
>
>
>
> *|** P **|** T **|** F **|** S **|*
>
> Empowering Content
>
> 1801 Research Blvd
>
> Suite 310
>
> Rockville, MD  20850
>
> 301 654 8088 Ext 131
>
> www.ptfs.com
>
> syefimov at ptfs.com <udoh at ptfs.com>
> _______________________________________________
> pdal mailing list
> pdal at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/pdal
>
> _______________________________________________
> pdal mailing list
> pdal at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/pdal
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/pdal/attachments/20240608/f726709f/attachment.htm>