[BUG] Notebook examples - AttributeError: 'StratigraphicColumn' object has no attribute 'keys'

[gpirot]

Describe your issue

Working with LoopStructural v1.6.16, I get the AttributeError: 'StratigraphicColumn' object has no attribute 'keys' when running the provided example notebooks (1b. Implicit surface modelling, 1c. Visualising models, 1d. Using Stratigraphic Columns, 1j. Exporting models)

Minimal reproducing code example

https://loop3d.org/LoopStructural/_auto_examples/1_basic/plot_2_surface_modelling.html 

https://loop3d.org/LoopStructural/_auto_examples/1_basic/plot_3_model_visualisation.html

https://loop3d.org/LoopStructural/_auto_examples/1_basic/plot_4_using_stratigraphic_column.html

https://loop3d.org/LoopStructural/_auto_examples/1_basic/plot_7_exporting.html

Error message

ERROR FOR 1b.
---------------------------------------------------------------------------
AttributeError Traceback (most recent call last)
Cell In[8], line 2
 1 viewer = Loop3DView(model)
----> 2 viewer.plot_block_model(cmap="tab20")
 3 viewer.display()

File ~\Anaconda3\envs\loop2025\lib\site-packages\loopstructuralvisualisation\_3d_viewer.py:293, in Loop3DView.plot_block_model(self, cmap, model, pyvista_kwargs, scalar_bar, slicer, threshold, name)
 291 name = 'block_model'
 292 name = self.increment_name(name) # , 'block_model')
--> 293 block, codes = model.get_block_model()
 294 block = block.vtk()
 295 block.set_active_scalars('stratigraphy')

File ~\Anaconda3\envs\loop2025\lib\site-packages\LoopStructural\modelling\core\geological_model.py:1755, in GeologicalModel.get_block_model(self, name)
 1752 def get_block_model(self, name='block model'):
 1753 grid = self.bounding_box.structured_grid(name=name)
-> 1755 grid.cell_properties['stratigraphy'] = self.evaluate_model(
 1756 self.rescale(self.bounding_box.cell_centres())
 1757 )
 1758 return grid, self.stratigraphic_ids()

File ~\Anaconda3\envs\loop2025\lib\site-packages\LoopStructural\modelling\core\geological_model.py:1497, in GeologicalModel.evaluate_model(self, xyz, scale)
 1495 logger.warning("No stratigraphic column defined")
 1496 return strat_id
-> 1497 for group in reversed(self.stratigraphic_column.keys()):
 1498 if group == "faults":
 1499 continue

AttributeError: 'StratigraphicColumn' object has no attribute 'keys'



ERROR FOR 1c.
---------------------------------------------------------------------------
AttributeError Traceback (most recent call last)
Cell In[3], line 10
 6 viewer.plot_surface(strati, value=vals, cmap="prism", paint_with=strati)
 9 viewer.plot_scalar_field(strati, cmap="prism")
---> 10 print(viewer._build_stratigraphic_cmap(model))
 11 viewer.plot_block_model()
 12 # Add the data addgrad/addvalue arguments are optional

File ~\Anaconda3\envs\loop2025\lib\site-packages\loopstructuralvisualisation\_3d_viewer.py:371, in Loop3DView._build_stratigraphic_cmap(self, model)
 369 boundaries = []
 370 data = []
--> 371 for g in model.stratigraphic_column.keys():
 372 if g == "faults":
 373 continue

AttributeError: 'StratigraphicColumn' object has no attribute 'keys'

ERROR FOR 1d.
---------------------------------------------------------------------------
AttributeError Traceback (most recent call last)
Cell In[3], line 4
 1 model.set_stratigraphic_column(stratigraphic_column)
 3 viewer = Loop3DView(model)
----> 4 viewer.plot_block_model()
 5 viewer.display()

File ~\Anaconda3\envs\loop2025\lib\site-packages\loopstructuralvisualisation\_3d_viewer.py:293, in Loop3DView.plot_block_model(self, cmap, model, pyvista_kwargs, scalar_bar, slicer, threshold, name)
 291 name = 'block_model'
 292 name = self.increment_name(name) # , 'block_model')
--> 293 block, codes = model.get_block_model()
 294 block = block.vtk()
 295 block.set_active_scalars('stratigraphy')

File ~\Anaconda3\envs\loop2025\lib\site-packages\LoopStructural\modelling\core\geological_model.py:1755, in GeologicalModel.get_block_model(self, name)
 1752 def get_block_model(self, name='block model'):
 1753 grid = self.bounding_box.structured_grid(name=name)
-> 1755 grid.cell_properties['stratigraphy'] = self.evaluate_model(
 1756 self.rescale(self.bounding_box.cell_centres())
 1757 )
 1758 return grid, self.stratigraphic_ids()

File ~\Anaconda3\envs\loop2025\lib\site-packages\LoopStructural\modelling\core\geological_model.py:1497, in GeologicalModel.evaluate_model(self, xyz, scale)
 1495 logger.warning("No stratigraphic column defined")
 1496 return strat_id
-> 1497 for group in reversed(self.stratigraphic_column.keys()):
 1498 if group == "faults":
 1499 continue

AttributeError: 'StratigraphicColumn' object has no attribute 'keys'


ERROR FOR 1j.
---------------------------------------------------------------------------
AttributeError Traceback (most recent call last)
Cell In[7], line 1
----> 1 model.save('model.geoh5')

File ~\Anaconda3\envs\loop2025\lib\site-packages\LoopStructural\modelling\core\geological_model.py:1788, in GeologicalModel.save(self, filename, block_model, stratigraphic_surfaces, fault_surfaces, stratigraphic_data, fault_data)
 1786 s.save(f'{parent}/{name}_{s.name}{extension}')
 1787 if block_model:
-> 1788 grid, ids = self.get_block_model()
 1789 if extension == ".geoh5" or extension == '.omf':
 1790 grid.save(filename)

File ~\Anaconda3\envs\loop2025\lib\site-packages\LoopStructural\modelling\core\geological_model.py:1755, in GeologicalModel.get_block_model(self, name)
 1752 def get_block_model(self, name='block model'):
 1753 grid = self.bounding_box.structured_grid(name=name)
-> 1755 grid.cell_properties['stratigraphy'] = self.evaluate_model(
 1756 self.rescale(self.bounding_box.cell_centres())
 1757 )
 1758 return grid, self.stratigraphic_ids()

File ~\Anaconda3\envs\loop2025\lib\site-packages\LoopStructural\modelling\core\geological_model.py:1497, in GeologicalModel.evaluate_model(self, xyz, scale)
 1495 logger.warning("No stratigraphic column defined")
 1496 return strat_id
-> 1497 for group in reversed(self.stratigraphic_column.keys()):
 1498 if group == "faults":
 1499 continue

AttributeError: 'StratigraphicColumn' object has no attribute 'keys'

[gpirot]

After upgrading to 1.6.18, I still have the same AttributeError: 'StratigraphicColumn' object has no attribute 'keys'
in the example provided: 1_basic/plot_7_exporting (but not for the other examples)