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: mdsum...@gmail.com
_______________________________________________
gdal-dev mailing list
gdal-dev@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.
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev
