cmutel / pandarus

GIS toolkit for regionalized life cycle assessment

Pandarus

Pandarus is a GIS software toolkit for regionalized life cycle assessment. It is designed to work with brightway LCA framework, brightway2-regional, and Constructive Geometries.

For more information, including installation instructions, see the online documentation, including a usage example. Pandarus is only compatible with >= Python 3.8.

Pandarus provides the following functionality:

• Overlay two vector datasets, calculating the areas of each combination of features using the Mollweide projection.
• Calculate the area of the geometric difference (the areas present in one input file but not in the other) of one vector dataset with another vector dataset.
• Calculate statistics such as min, mean, and max when overlaying a raster dataset with a vector dataset.
• Normalize raster datasets, including use of compatible nodata values
• Vectorization of raster datasets

Requirements

• appdirs
• fiona
• numpy
• pyproj
• rasterio
• rasterstats
• Rtree
• shapely

Development

Pandarus is developed by Chris Mutel as part of his work at the Technology Assessment Group at the Paul Scherrer Institut and previously at the Ecological Systems Design group at ETH Zürich.

Source code is available on GitHub.

Contributing

Contributions are very welcome. To learn more, see the Contributor Guide.

License

Distributed under the terms of the BSD-2-Clause, pyilcd is free and open source software.

Issues

If you encounter any problems, please file an issue along with a detailed description.

pandarus's Issues

RemoteError: Can't reach https://pandarus.brightway.dev/

Hi,
Seems that pandarus link is not working.

RemoteError Traceback (most recent call last)
Cell In[16], line 1
----> 1 bwr.remote.intersection("world", "watersheds-hh")
2 bwr.remote.intersection("world", "watersheds-eq-sw-certain")
3 bwr.remote.intersection("world", "watersheds-eq-sw-all")

File ~/.conda/envs/bwcourse23/lib/python3.11/site-packages/bw2regional/pandarus_remote.py:73, in check_alive(wrapped, instance, args, kwargs)
70 @wrapt.decorator
71 def check_alive(wrapped, instance, args, kwargs):
72 if not instance.alive:
---> 73 raise RemoteError("Can't reach {}".format(instance.url))
74 return wrapped(*args, **kwargs)

RemoteError: Can't reach https://pandarus.brightway.dev/

Please, I would appreciate any help you can give me with this issue :)

Intersecting point and vector raises error

File ~/miniconda3/envs/pandarus_remote/lib/python3.9/site-packages/pandarus/calculate.py:258, in intersect(first_fp, first_field, second_fp, second_field, first_kwargs, second_kwargs, dirpath, cpus, driver, compress, log_dir)
    251     with fiona.open(
    252             fiona_fp, 'w',
    253             crs=WGS84,
    254             driver=driver,
    255             schema=schema,
    256         ) as sink:
    257         for f in as_features(data):
--> 258             sink.write(f)
    260 data_fp = json_exporter(
    261     {
    262         'data': [(k[0], k[1], v['measure']) for k, v in data.items()],
   (...)
    271     compress=compress
    272 )
    274 return fiona_fp, data_fp

File ~/miniconda3/envs/pandarus_remote/lib/python3.9/site-packages/fiona/collection.py:367, in Collection.write(self, record)
    365 def write(self, record):
    366     """Stages a record for writing to disk."""
--> 367     self.writerecords([record])

File ~/miniconda3/envs/pandarus_remote/lib/python3.9/site-packages/fiona/collection.py:361, in Collection.writerecords(self, records)
    359 if self.mode not in ('a', 'w'):
    360     raise IOError("collection not open for writing")
--> 361 self.session.writerecs(records, self)
    362 self._len = self.session.get_length()
    363 self._bounds = None

File fiona/ogrext.pyx:1285, in fiona.ogrext.WritingSession.writerecs()

GeometryTypeValidationError: Record's geometry type does not match collection schema's geometry type: 'MultiPoint' != 'MultiPolygon'

Proj.is_latlong is deprecated

The fonction Proj.is_latlong that is used in pandarus.projection.project line 47 is deprecated in the package Pyproj >= 2.2.0. It has been replaced by a Proj.crs.is_geographic property.