aadm / geophysical_notes Goto Github PK

Dear Great man.

Indeed great effort.... !!
Have you tried this wedge modeling. ? for me, some issues with agilegeo interlinked packages. even when they are available . like wavelet .. attributes. If you can comment.

greetings

Khawar

Dear Respected Fellows.

Good Morning. I hope you are fine. First of all, very nice work you have done with python. I am impressed.
I am facing an issue and want to request for help. I am trying to read the segy dataset. No matter it is your data of mine. it gives the error as below

\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\

ValueError Traceback (most recent call last)
~/.local/lib/python3.5/site-packages/segpy/header.py in set(self, instance, value)
299 try:
--> 300 self._instance_data[instance] = self._named_field._value_type(value)
301

~/.local/lib/python3.5/site-packages/segpy/field_types.py in class_new(cls, *args, **kwargs)
13 cls.name, instance,
---> 14 cls.MINIMUM, cls.MAXIMUM))
15 return instance

ValueError: NNInt16 value 40975 outside range 0 to 32767

The above exception was the direct cause of the following exception:

ValueError Traceback (most recent call last)
in ()
1 filename='reproduced2.sgy'
2 f=open(filename, 'rb')
----> 3 segy0 = create_reader(f)
4
5 ntraces = segy0.num_traces()

~/.local/lib/python3.5/site-packages/segpy/reader.py in create_reader(fh, encoding, binary_reel_header_format, trace_header_format, endian, progress, cache_directory, dimensionality)
151 if reader is None:
152 reader = _make_reader(fh, encoding, binary_reel_header_format, trace_header_format,
--> 153 endian, progress_callback, dimensionality)
154 if cache_file_path is not None:
155 _save_reader_to_cache(reader, cache_file_path)

~/.local/lib/python3.5/site-packages/segpy/reader.py in _make_reader(fh, encoding, binary_reel_header_format, trace_header_format, endian, progress, dimensionality)
256 encoding = ASCII
257 textual_reel_header = read_textual_reel_header(fh, encoding)
--> 258 binary_reel_header = read_binary_reel_header(fh, binary_reel_header_format, endian=endian)
259 extended_textual_header = read_extended_textual_headers(fh, binary_reel_header, encoding)
260 bps = bytes_per_sample(binary_reel_header)

~/.local/lib/python3.5/site-packages/segpy/toolkit.py in read_binary_reel_header(fh, binary_reel_header_format, endian)
180 header_packer = make_header_packer(binary_reel_header_format, endian)
181 buffer = fh.read(binary_reel_header_format.LENGTH_IN_BYTES)
--> 182 reel_header = header_packer.unpack(buffer)
183 return reel_header
184

~/.local/lib/python3.5/site-packages/segpy/packer.py in unpack(self, buffer)
198 str(e).capitalize())) from e
199 else:
--> 200 return self._unpack(values)
201
202 @AbstractMethod

~/.local/lib/python3.5/site-packages/segpy/packer.py in _unpack(self, values)
222
223 def _unpack(self, values):
--> 224 return self._header_format_class(*values)
225
226

~/.local/lib/python3.5/site-packages/segpy/header.py in init(self, *args, **kwargs)
29 """
30 for keyword, arg in zip(self.ordered_field_names(), args):
---> 31 setattr(self, keyword, arg)
32
33 for keyword, arg in kwargs.items():

~/.local/lib/python3.5/site-packages/segpy/header.py in set(self, instance, value)
303 except ValueError as e:
304 raise ValueError("Assigned value {!r} for {} attribute must be convertible to {}: {}"
--> 305 .format(value, self._name, self._named_field._value_type.name, e)) from e
306
307 def delete(self, instance):

ValueError: Assigned value 40975 for sample_interval attribute must be convertible to NNInt16: NNInt16 value 40975 outside range 0 to 32767
\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\
my file is little endian based. Please at least help me to be able to run your script to reproduce the results.
Please if you can help me resolving this issue would be great help for me.

Looking forward

Regards

Khawar
Khawar