ENH: enable reading topography from res2dinv data (#629)

gimli-org · Dec 13, 2023 · 48ce372 · 48ce372
1 parent 5f1033f
commit 48ce372
Showing 1 changed file with 10 additions and 0 deletions.
diff --git a/pygimli/physics/ert/importData.py b/pygimli/physics/ert/importData.py
@@ -276,6 +276,16 @@ def getNonEmptyRow(i, comment='#'):
         if hasIP:
             data.set('ip', dataBody[nn - 1])
 
+    istopo = int(getNonEmptyRow(it, comment=';'))
+    if istopo:
+        ntopo = int(getNonEmptyRow(it, comment=';'))
+        ap = data.additionalPoints()
+        for i in range(ntopo):
+            strs = getNonEmptyRow(it, comment=';').replace(',', ' ').split()
+            ap.push_back(pg.Pos([float(s) for s in strs]))
+
+        data.setAdditionalPoints(ap)
+
     data.sortSensorsX()
     data.sortSensorsIndex()
     if return_header: