Example of masked areal interpolation

[martinfleis]

Do we have an example notebook for masked area interpolation? The one using raster. I thought I have seen it before and wanted to share it with a student but can't find anything now. @knaaptime do you have it in some research code available somewhere, at least?

[knaaptime]

dang. the logic is here, but that's all been wrapped into the function now and it doesnt look like i ever pushed up an updated example. Lame on my part. Let me find something for you

[knaaptime]

https://github.com/pysal/tobler/blob/main/notebooks/binary_dasymetric.ipynb (wong link earlier)

[martinfleis]

I've seen that one but was hoping to find a direct application of masked_area_interpolate

[knaaptime]

yeah i'll grab you one

[knaaptime]

not much explanation here (and it works a lot better with the interactive map), but here's a barebones example https://gist.github.com/knaaptime/7067008cac2837aa131793f0f3f1faca

[martinfleis]

You're the best! Thanks!

[knaaptime]

@martinfleis unrelated, but whats up with shapely sending those invalid value encountered in intersection warnings?

[martinfleis]

GEOS issue that in fact does not affect the outcome. Some of it was fixed in GEOS 3.12 but I'm not sure if all.