Currently, coalescence data are saved as obsy Stream objects output to miniSEED. Value

I have added the following snippet to the default scan class: <div class="snippet-

np.int32 overflow about quakemigrate HOT 2 CLOSED

quakemigrate commented on May 29, 2024

np.int32 overflow

from quakemigrate.

Comments (2)

hemmelig commented on May 29, 2024

I have added the following snippet to the default scan class:

c_precision = 1e6
precision = 1e6

and

        data["COA"] = coa.select(station="COA")[0].data / self.c_precision
        data["COA_N"] = coa.select(station="COA_N")[0].data / self.c_precision
        data["X"] = coa.select(station="X")[0].data / self.precision
        data["Y"] = coa.select(station="Y")[0].data / self.precision
        data["Z"] = coa.select(station="Z")[0].data / self.precision

        st = Stream(Trace(data=(np.array(data["COA"]) * self.c_precision).astype(np.int32),
                          header=stats_COA))
        st += Stream(Trace(data=(np.array(data["COA_N"]) * self.c_precision).astype(np.int32),
                           header=stats_COA_N))
        st += Stream(Trace(data=(np.array(data["X"]) * self.precision).astype(np.int32),
                           header=stats_X))
        st += Stream(Trace(data=(np.array(data["Y"]) * self.precision).astype(np.int32),
                           header=stats_Y))
        st += Stream(Trace(data=(np.array(data["Z"]) * self.precision).astype(np.int32),
                           header=stats_Z))

This is similar to the way in which obspy handles UTCDateTime precision (class variable that can be overwritten). Need to think of a way to double check the input from the user and tell them if this may cause them to exceed the capacity of np.int32 (as well as potentially giving them the option to change it to np.int64). Currently it's looking like it'll have to be by extending the signature of the SeisPlot class to take the precision and save dtype.

from quakemigrate.

hemmelig commented on May 29, 2024

Final solution described in commit that closes this issue.

from quakemigrate.

Recommend Projects