[gdal-dev] convert raster to raw table ?

[Even Rouault]

Michael,

doesn't do exactly what you want, but should be easy to modify to 
satisfy your requirements: 
https://github.com/OSGeo/gdal/blob/master/apps/gdal2ogr.c

Even

Le 12/12/2024 à 06:40, Michael Sumner via gdal-dev a écrit :
> Hi, I'd like to be able to create a tabular version efficiently 
> directly from a raster. I can do this in scripting but is there a 
> utility- or api- way to do this?
>
> I literally just want a Parquet file with a column for each band from 
> a raster, with raw pixel values, no index, no coordinates, no rat or 
> anything like that. I have very common workflows for this where the 
> georeferencing (bbox, shape, crs) is just used manually.
>
> If this doesn't exist a barebones outline would be amazing (!), and 
> (I) could add features like column names, cell index or row/col index, 
> coordinates, sparse representation for missing values etc etc.
>
> An imagined -info output would look like this for a 2x3 raster with 2 
> bands and values 1:6 in band 1, and values 7:12 in band 2 described  
> below.
>
> thanks so much, Mike
>
>
>
>
> gdalinfo file.tif
> Driver: GTiff/GeoTIFF
> Files: file.tif
> Size is 3, 2
> Origin = (0.000000000000000,2.000000000000000)
> Pixel Size = (1.000000000000000,-1.000000000000000)
> Image Structure Metadata:
>   COMPRESSION=LZW
>   INTERLEAVE=PIXEL
> Corner Coordinates:
> Upper Left  (   0.0000000,   2.0000000)
> Lower Left  (   0.0000000,   0.0000000)
> Upper Right (   3.0000000,   2.0000000)
> Lower Right (   3.0000000,   0.0000000)
> Center      (   1.5000000,   1.0000000)
> Band 1 Block=3x2 Type=UInt32, ColorInterp=Gray
>   Description = a
>   Min=1.000 Max=6.000
>   Minimum=1.000, Maximum=6.000, Mean=-9999.000, StdDev=-9999.000
>   NoData Value=4294967295
>   Metadata:
>     STATISTICS_MINIMUM=1
>     STATISTICS_MAXIMUM=6
>     STATISTICS_MEAN=-9999
>     STATISTICS_STDDEV=-9999
> Band 2 Block=3x2 Type=UInt32, ColorInterp=Undefined
>   Description = b
>   Min=7.000 Max=12.000
>   Minimum=7.000, Maximum=12.000, Mean=-9999.000, StdDev=-9999.000
>   NoData Value=4294967295
>   Metadata:
>     STATISTICS_MINIMUM=7
>     STATISTICS_MAXIMUM=12
>     STATISTICS_MEAN=-9999
>     STATISTICS_STDDEV=-9999
>
>
>
> INFO: Open of `file.parquet'
>       using driver `Parquet' successful.
>
> Layer name: file
> Geometry: None
> Feature Count: 6
> Layer SRS WKT:
> (unknown)
> a: Integer (0.0)
> b: Integer (0.0)
> OGRFeature(file):0
>   a (Integer) = 1
>   b (Integer) = 7
>
> OGRFeature(file):1
>   a (Integer) = 3
>   b (Integer) = 9
>
> OGRFeature(file):2
>   a (Integer) = 5
>   b (Integer) = 11
>
> OGRFeature(file):3
>   a (Integer) = 2
>   b (Integer) = 8
>
> OGRFeature(file):4
>   a (Integer) = 4
>   b (Integer) = 10
>
> OGRFeature(file):5
>   a (Integer) = 6
>   b (Integer) = 12
>
>
> -- 
> Michael Sumner
> Research Software Engineer
> Australian Antarctic Division
> Hobart, Australia
> e-mail: mdsumner at gmail.com
>
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/gdal-dev

-- 
http://www.spatialys.com
My software is free, but my time generally not.
Butcher of all kinds of standards, open or closed formats. At the end, this is just about bytes.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20241212/8543367a/attachment.htm>