gimli-org / gimli

Geophysical Inversion and Modeling Library :earth_africa:
https://www.pygimli.org
Other
344 stars 131 forks source link

Removing columns for ERT data #686

Closed IsabellaLorenzen closed 1 month ago

IsabellaLorenzen commented 1 month ago

Problem description

Hallo!

I am unsure if you can call them columns but when i import my ERT data i get the following output:

Data: Sensors: 181 data: 4768, nonzero entries: ['a', 'b', 'ip', 'ip1', 'ip2', 'ip3', 'ip4', 'ip5', 'ip6', 'ip7', 'ip8', 'ip9', 'm', 'n', 'r', 'valid']

My IP decay curves shows very erratic behavior in the first IP-window, so i wish to remove the "columns" 'ip' and 'ip1' (they are the same)

Your environment

            OS : Windows
        CPU(s) : 8
       Machine : AMD64
  Architecture : 64bit
           RAM : 7.8 GiB
   Environment : IPython

Steps to reproduce

I am trying to remove the "columns" by using this code:

import pygimli as pg

dataOriginal_M.remove(dataOriginal_M['ip'])

Expected behavior

Removing of the "Column" 'ip'

Actual behavior

ArgumentError: Python argument types in
    DataContainerERT.remove(DataContainerERT, RVector)
did not match C++ signature:
    remove(GIMLI::DataContainerERT {lvalue}, GIMLI::Vector<bool> bvec)
    remove(DataContainerERT_wrapper {lvalue}, GIMLI::Vector<unsigned long long> idx)
    remove(GIMLI::DataContainerERT {lvalue}, GIMLI::Vector<unsigned long long> idx)

If possible, please add one or more labels to your issue, e.g. if you expect that your issue is rather a question than a problem with the code, please add the label "question".

halbmy commented 1 month ago

Actually, looking at dir(data) there seems to be no way of removing a data field (implemented in the C++ core).

What you can do is to prevent it from saving, e.g. data.save("myfile.dat", "a b m n u i ip4 ip5")

IsabellaLorenzen commented 1 month ago

Alright, thank you :)