Support specifying coordinate reference system (CRS) of geometries

[Hofer-Julian]

Fixes #1254

Breaking changes:

• Adds new required TOML key "crs"
• Adds new required Model parameter "crs" in Ribasim-Python.

Note: we could avoid adding the TOML key by requiring Model.read to take crs as parameter. Not sure if that's an improvement though.

Upgrade instructions

Specify the CRS. For instance for WGS 84 latitude longitude, use crs = "EPSG:4326", for the Netherlands use crs = "EPSG:28992".

[visr]

Yes I think adding the TOML key is fine. I added some commits. Notably a test that "Basin / area" also gets the CRS, but that fails, so I think that still needs to be fixed.

If it is hard to fix it is fine as-is, but currently the Node tables only gets a CRS when it is assembled by node_table(). That means if you just want the GeoDataframe of the Basin nodes with model.basin.node you get a GeoDataframe with no CRS set.

[Hofer-Julian]

If it is hard to fix it is fine as-is, but currently the Node tables only gets a CRS when it is assembled by node_table(). That means if you just want the GeoDataframe of the Basin nodes with model.basin.node you get a GeoDataframe with no CRS set.

Yeah, I don't know how that would be done. When we run model.basin.add, we don't have access to the model state and therefore no access to the crs.

[visr]

Would it make sense to define model.set_crs() that sets the CRS of all spatial tables? Then we can call that on save, but also earlier if needed.

[Huite]

Indeed. As I mentioned in #1254, I would probably implement both set_crs (to just assign a crs) and to_crs (to reproject).

[Hofer-Julian]

Would it make sense to define model.set_crs() that sets the CRS of all spatial tables? Then we can call that on save, but also earlier if needed.

Sounds good to me.