digitalcourage.social is one of the many independent Mastodon servers you can use to participate in the fediverse.
Diese Instanz wird betrieben von Digitalcourage e.V. für die Allgemeinheit. Damit wir das nachhaltig tun können, erheben wir einen jährlichen Vorausbeitrag von 1€/Monat per SEPA-Lastschrifteinzug.

Server stats:

828
active users

#gdal

6 posts6 participants0 posts today

okay #GDAL fans - how the hell am I supposed to work with jpeg2000 files and why on God's earth is the Copernicus Data Space Ecosystem using this format? I would much rather make requests to a european server but it looks like I'm expected to just download entire files and nothing is cloud optimised 😩 - really hoping I'm wrong...

I've been quiet about {vrtility} recently but I've made some nice steady progress with improved asynch read/write of raster files. A colleague wanted some help getting at the amazing Digital Earth Africa GeoMAD landsat catalogs so I put togther this gist - sharing incase it's useful to anyone else! with an already beautifully processed annual image composite - this is a classic and simple use case for VRT 🤓 . #rstats #rspatial #gdal gist.github.com/h-a-graham/461

Found a cool feature of ogr2ogr: you can use SQL like 'SELECT * EXCEPT (col1, col2) FROM table' with the -sql parameter to exclude a few fields from the conversion process. Handy when you've got lots of fields and don't want to list them all out in the SELECT.

Now I just wish Postgres had this functionality (DuckDB does!)

See gdal.org/en/stable/user/ogr_sq for details, and I found it from github.com/OSGeo/gdal/issues/8

gdal.orgOGR SQL dialect — GDAL documentation

#GDAL is there a way to specify an empty raster source entirely as a text string (compactly)?

I can use MEM and generate a datapointer description from an array in memory, but what about the equivalent of an empty cog? something like

"MEM:::EMPTY_COG,87,61,1,Float64,0/1/0/61/0/-1"

(I certainly don't recommend experimenting with arbitrary values for DATAPOINTER)

GDAL modernization webinar and FAQ: The #GDAL Project Steering Committee recently hosted a webinar highlighting the major #modernization of the GDAL command-line interface, featuring the unification of commands and improved parameter consistency. You can now catch up on all the...
spatialists.ch/posts/2025/06-0 #GIS #GISchat #geospatial #SwissGIS

Spatialists – geospatial newsGDAL modernization webinar and FAQ – Spatialists – geospatial news
More from Spatialists

Quite excited by the new artihmetic #GDAL band on-the-fly capabilities of github.com/OSGeo/gdal/pull/125

```python
with gdal.Open("rgb.tif") as ds:
R = ds.GetRasterBand(1)
G = ds.GetRasterBand(2)
B = ds.GetRasterBand(3)
gray = (0.299 * R + 0.587 * G + 0.114 * B).as_type(gdal.GDT_Byte)
gdal.GetDriverByName("GTiff").CreateCopy("gray.tif", gray)
```

GitHub[WIP] GDALRasterBand: define +, -, *, /, AsType(), gdal::min()/max() to perform band artihmetics from C++, C and SWIG (Python) by rouault · Pull Request #12507 · OSGeo/gdalBy rouault

Fun #GDAL contest. Craft the smallest file (in terms of byte size) that is recognized (opened by GDAL, using any official driver of your choice) with 2147483647 x 2147483647 raster size. My current record is 27 bytes. Any improvement over that?

Kicked off a HPC job, 64 CPUs, converting global 11 variable netcdf to single var cog, from and to the local object storage.

There's 120K files to be created, and 20 minutes or so in it's up to 37K done, each is 2880x1440, Int32 (zstd, predictor 2)

Blows my mind how fast stuff is when obstacles and complications are removed
...thanks to #GDAL
and #rstats {gdalraster}, and {mirai} via {purrr}