gimli-org / example-data

Collection of example data for the use in tutorials & examples
3 stars 5 forks source link

how to create .sgt data fomat for single shot and several geophones? #2

Closed ad1995-geek closed 1 year ago

halbmy commented 3 years ago

The format is easy, it first describes the sensor (both shot and geophone for traveltime) positions and then the data. Assume a single shot with 4 receivers, each spaced 5m with an offset of 10m:

5 # 1 shot and 4 receivers
# x z
-10 0
0 0
5 0
10 0
15 0
# s g t
1 2 0.2
1 3 0.3
1 4 0.5
1 5 0.6

You can create this in Python

import numpy as np
import pygimli as pg

data = pg.DataContainer()
data.registerSensorIndex("s")
data.registerSensorIndex("g")
data.createSensor([-10, 0]) # shot
for i in range(4):
    data.createSensor([i*5, 0]) # receiver

data.resize()
data["s"] = np.zeros(4)
data["g"] = np.arange(1, 5)
data.save("out.sgt")
ad1995-geek commented 3 years ago

Thank you for the reply and explanation.