aadm / geophysical_notes

Collection of geophysical notes in the form of IPython/Jupyter notebooks.
147 stars 75 forks source link

reading script #1

Closed Tortoise17 closed 6 years ago

Tortoise17 commented 6 years ago

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
Tortoise17 commented 6 years ago

Resolved. !!