mdsumner
@mdsumner.bsky.social
290 followers 380 following 500 posts
maxing the pixels to glow meshes, R, Python, GDAL, virtual Zarr, coordinate systems, southern ocean research
Posts Media Videos Starter Packs
Reposted by mdsumner
mdsumner.bsky.social
It worked, replaced all saved RasterLayer objects with actively-bound exported names, first use after package load triggers load from tif to SpatRaster
Reposted by mdsumner
ctuise.bsky.social
"a lot like a protection racket for the gambling industry" is a very polite way of describing current state of gov't action.
#Auspol

Australia’s gambling harm is likely underreported – and authorities are still failing to act
theconversation.com/australias-g...
Australia’s gambling harm is likely underreported – and authorities are still failing to act
New research shows Australian gambling harms, including suicide, may have been systematically under reported.
theconversation.com
mdsumner.bsky.social
#rstats what do folks do when converting from raster to terra with saved data objects? active bindings is a possibility

github.com/ropensci/gra...

interested in any perspectives, I don't want to change the entire way the package works (I'm not the maintainer, just toying with ideas)
Reposted by mdsumner
kylebarron.dev
New blog post! kylebarron.dev/blog/fast-py...

This dives into performance pitfalls of Python bindings, why projects like Lonboard are so fast, and why I care about Apache Arrow so much. The post is high level and you don't need to know C/C++/Rust to understand it!
Reposted by mdsumner
eliocamp.mastodon.social.ap.brid.gy
The solution was to not use coord_sf(). I need to remember to not use sf in any project unless absolutely necessary.
Reposted by mdsumner
eriksdotter.bsky.social
in short:
it’s complicated and we’re fucked

#climatecrisis
bjerknes.uib.no
Rising sea levels are making extreme winter cold events MORE frequent and intense in East Asia. Yes, you read that right. New paper in Nature Communications reveals a fascinating connection. (1/4) 🧪🌊
Reposted by mdsumner
peterfox59.bsky.social
Great story on ABC #4Corners now on how online gambling companies operate, the damage they do to the community & how they keep the major political parties in the pockets.
Reposted by mdsumner
paleolimbot.bsky.social
Definitely depends what you're up to! People are definitely experimenting with serving features directly from Parquet files for visualization (e.g., medium.com/center-for-c... )
Optimal GeoParquet Partitioning Strategy
My quest for one-size-fits-all format
medium.com
mdsumner.bsky.social
is Scott Morrison in jail yet?

Impressive to outdo John Howard for evil
Reposted by mdsumner
govpritzker.illinois.gov
This evening, President Trump is ordering 400 members of the Texas National Guard for deployments to Illinois, Oregon, and other locations within the United States. No officials from the federal government called me directly to discuss or coordinate.
Reposted by mdsumner
aaliyahxx.bsky.social
Watch | Israeli occupation forces have released new drone footage from Gaza, revealing the mass destruction and the erasure of every means of life in the Strip.
Reposted by mdsumner
satojuniper.bsky.social
Amsterdam, we salute you!
#Auspol
pauleric70.bsky.social
They tried to erase Palestine, but the world has become Palestine.

The outrage is global. The world is rising against genocide.

Thank you Amsterdam.
Reposted by mdsumner
pp0196.bsky.social
The claim seems widely overstated because in the end it is about selecting the partitioning in order to push as much computation into storage layer as possible. Selecting the partitioning is mostly about access pattern just like db indexes
Reposted by mdsumner
funferall.bsky.social
Oh yes I'm using gdalwarp here.

User draws on wmts and it looks at a .vrt file to match which geotiffs it overlays. That way I don't have to display the tifs.

Thinking of adding more crs options for the output and choice between geotiff and jpg2000. Just prototyping this
Reposted by mdsumner
funferall.bsky.social
Got the gdal clip working. First time hearing about ALB, security groups, ECS, services, task definitions, cloudfront, a bunch of other bs. I think I can scale this (so says Claude 🤷‍♂️)
funferall.bsky.social
So deep in AWS land again. If this eventually works it will be a miracle
Reposted by mdsumner
paleolimbot.bsky.social
- FGB is great! Lookups with small ranges work well but full scans are much slower (particularly if zipped!)
- SedonaDB, DuckDB, and pyarrow with fsspec can all do partial Parquet reads from http (give geoarrow.org/data a try!)
Data – GeoArrow
geoarrow.org
mdsumner.bsky.social
I did a filter on the ibis example

I think you're talking about tooling problems
mdsumner.bsky.social
Nice! Could use gdalwarp (gdal raster reproject) to clip and write (and resize or change crs) in one step fwiw

(gdal raster reproject has a --bbox-crs arg, and this can be pipelined or written in json as a GDALG virtual dataset)
Reposted by mdsumner
funferall.bsky.social
Hacking together a Clip-Zip-Ship tool: draw custom polygon over wmts imagery service in web mercator, converts to UTM zone 19 coordinates and clips the corresponding geotiff (4band, 15cm, 16bit) and delivers to user.

Next: Select town as clipping boundary, give 100m buffer
draw custom polygon over wmts imagery service in web mercator, converts to UTM zone 19 and clips the corresponding geotiff (4band, 15cm, 16bit) and delivers to user
mdsumner.bsky.social
(open issue in nanoparquet was the unfinished bit there)