diff --git a/ogstools/meshlib/region/region.py b/ogstools/meshlib/region/region.py
index 1ae8660dbcee93abcacf7c40d0b200b0c7a157a5..994daff515d921baf756ea0fc7cf23efc37d7511 100644
--- a/ogstools/meshlib/region/region.py
+++ b/ogstools/meshlib/region/region.py
@@ -88,12 +88,10 @@ def to_boundary(
         filter_condition(surface_mesh["Normals"])
     ]
 
-    surface_mesh.rename_array("vtkOriginalPointIds", "BULK_NODE_ID")
-    surface_mesh.rename_array("vtkOriginalCellIds", "BULK_ELEMENT_ID")
-    surface_mesh.cell_data.remove("Normals")
     specific_cells = surface_mesh.extract_cells(ids)
-    specific_cells.cell_data.remove("vtkOriginalCellIds")
-    specific_cells.point_data.remove("vtkOriginalPointIds")
+    specific_cells.rename_array("vtkOriginalPointIds", "BULK_NODE_ID")
+    specific_cells.rename_array("vtkOriginalCellIds", "BULK_ELEMENT_ID")
+    specific_cells.cell_data.remove("Normals")
     return specific_cells
 
 
diff --git a/pyproject.toml b/pyproject.toml
index fb4f613297d8fffb97ed82a256642604a3265896..e06dee527a93c7221fae07f14b8d8363e33b65ee 100644
--- a/pyproject.toml
+++ b/pyproject.toml
@@ -16,7 +16,7 @@ dependencies = [
   "meshio==5.3.4",
   "matplotlib>=3.7.1",
   "nbconvert>=7.9.2",
-  "pyvista[all,trame]>=0.39.1,<0.43.0",
+  "pyvista[all,trame]>=0.39.1",
   "rich>=13.4.2",
   "scipy>=1.10.1",
   "Pint>=0.22",