joferkington / scipy2014-3d_seismic Goto Github PK

Dear Joe Kington,

Thank you very much for sharing your code, experience and knowledge.
Your code and presentation are really helpful,useful and interesting for me.
I have successfully loaded and viewed the 3D volume, however, when I try to load horizon, I get the the following error:
I will be very grateful to you for any assistance or hint.

Mohammed

EOFError Traceback (most recent call last)
in ()
4 import scipy.ndimage
5
----> 6 hor = geoprobe.horizon('C:\Users\Desktop\Python\scipy2014-3d_seismic-master\scipy2014-3d_seismic-master\data\Horizons\channels_smoothed.hzn')
7 vol = geoprobe.volume('C:\Users\Desktop\Python\scipy2014-3d_seismic-master\scipy2014-3d_seismic-master\data\Volumes\example.vol')
8 data = vol.load()

C:\Users\AppData\Local\Enthought\Canopy32\User\lib\site-packages\geoprobe-0.3-py2.7.egg\geoprobe\horizon.pyc in init(self, _args, *_kwargs)
77 # and make a horizon object by reading from disk
78 if (len(args) == 1) and isinstance(args[0], basestring):
---> 79 self._readHorizon(args[0])
80
81 # Otherwise, pass the args on to _make_horizon_from_data for

C:\Users\AppData\Local\Enthought\Canopy32\User\lib\site-packages\geoprobe-0.3-py2.7.egg\geoprobe\horizon.pyc in _readHorizon(self, filename)
106 ' horizon')
107
--> 108 self.data = self._file.readAll()
109
110 # Surface and line attributes

C:\Users\AppData\Local\Enthought\Canopy32\User\lib\site-packages\geoprobe-0.3-py2.7.egg\geoprobe\horizon.pyc in readAll(self)
451 # Read lines section
452 line_info = [None]
--> 453 self.numlines = self.readSectionHeader()
454 for i in xrange(self.numlines):
455 line_info.append(self.readLineHeader())

C:\Users\AppData\Local\Enthought\Canopy32\User\lib\site-packages\geoprobe-0.3-py2.7.egg\geoprobe\horizon.pyc in readSectionHeader(self)
423
424 def readSectionHeader(self):
--> 425 return self.readBinary(self._sectionHdrFmt)
426
427 def readLineHeader(self):

C:\Users\AppData\Local\Enthought\Canopy32\User\lib\site-packages\geoprobe-0.3-py2.7.egg\geoprobe\common.pyc in readBinary(self, fmt)
65 # Reading beyond the end of the file just returns ''
66 if len(data) != size:
---> 67 raise EOFError('End of file reached')
68 data = struct.unpack(fmt, data)
69

EOFError: End of file reached

Dear Mr. Kington,

The code you have provided seems to be very useful for visualization of seismic data.

I have a SEGY file. Can you please tell how to convert segy file to vol file to be used with geoprobe library?

Regards,
Mohammad Farzanullah