private void button2_Click(object sender, EventArgs e)
{
foreach (var i in listBox1.Items)
{
String f = (string)i;
string f0 = System.IO.Path.GetDirectoryName(f) + System.IO.Path.DirectorySeparatorChar + System.IO.Path.GetFileNameWithoutExtension(f) + "_corr.sgy";
SEGYFile s0 = new SEGYFile();
if (s0.Open(f) == 0)
{
break;
}
SEGYFile s1 = new SEGYFile();
if (s1.Open(f0) == 0)
{
s0.Close();
break;
}
s0.ReadFileHeader();
s1.Write(s0.FileHeader.Copy());
while (s0.ReadNextTrace())
{
SEGYTrace t = s0.currentTrace.Copy();
short waterDepthMsec = t.TraceHeader.delayRecordingTimeMsec;
t.TraceHeader.delayRecordingTimeMsec = 0;
//t.TraceHeader.waterDepthAtSource = (int)waterDepthMsec;
s1.Write(t);
}
s0.Close();
s1.Close();
}
}
private void button2_Click(object sender, EventArgs e)
{
foreach (var i in listBox1.Items)
{
String f = (string)i;
string f0 = System.IO.Path.GetDirectoryName(f) + System.IO.Path.DirectorySeparatorChar + System.IO.Path.GetFileNameWithoutExtension(f) + "_corr.sgy";
SEGYFile s0 = new SEGYFile();
if (s0.Open(f) == 0)
{
break;
}
SEGYFile s1 = new SEGYFile();
if (s1.Open(f0) == 0)
{
s0.Close();
break;
}
s0.ReadFileHeader();
var fh = s0.FileHeader.Copy();
fh.SetFileHeader(0, 39, "C39 File processed by FixOBSSEGY Sept/2019");
s1.Write(fh);
while (s0.ReadNextTrace())
{
if (s0.currentTrace.TraceHeader.numberOfSamplesInTrace < 1000 || s0.currentTrace.TraceHeader.numberOfSamplesInTrace > 60000)
{
break;
}
if (s0.currentTrace.timeTracedRecorded < new DateTime(1950, 1, 1))
{
break;
}
SEGYTrace t = s0.currentTrace.Copy();
if (checkBoxSwitch.Checked)
{
var s = t.TraceHeader.sourceCoordinateX;
t.TraceHeader.sourceCoordinateX = t.TraceHeader.sourceCoordinateY;
t.TraceHeader.sourceCoordinateY = s;
s = t.TraceHeader.groupCoordinateX;
t.TraceHeader.groupCoordinateX = t.TraceHeader.groupCoordinateY;
t.TraceHeader.groupCoordinateY = s;
}
if (checkBoxlldec.Checked)
{
t.TraceHeader.coordinateUnits = 3;
}
if (checkBoxllsa.Checked)
{
t.TraceHeader.coordinateUnits = 2;
}
s1.Write(t);
}
s0.Close();
s1.Close();
}
}
private void runToolStripMenuItem_Click(object sender, EventArgs e)
{
foreach (string f in listBox1.Items)
{
fseg = new SEGYFile();
fseg.Open(f);
fseg.ReadFileHeader();
SEGYFile fn = new SEGYFile();
string fnn = f.Replace(".sgy", "_clipped.sgy");
//create new filename
if (System.IO.File.Exists(fnn))
{
System.IO.File.Delete(fnn);
}
fn.Open(fnn);
fn.FileHeader = fseg.FileHeader.Copy();
fn.Write(fn.FileHeader);
while (fseg.ReadNextTrace())
{
if (fseg.currentTrace.TraceHeader.yearDataRecorded != 2001 || fseg.currentTrace.TraceHeader.dayOfYear > 366)
{
int i = 0;
break;
}
// set 1140 thru 1222 to zero in trace block ( las 82 entries )
byte[] buff = fseg.currentTrace.TraceData.TraceDataBuffer;
// since 4 byte words
for (int i = buff.Length - 82 * 4; i < buff.Length; i++)
{
buff[i] = 0;
}
fn.Write(fseg.currentTrace);
}
fn.Close();
fseg.Close();
}
}
private void button2_Click(object sender, EventArgs e)
{
foreach (var i in listBox1.Items)
{
String f = (string)i;
string f0 = System.IO.Path.GetDirectoryName(f) + System.IO.Path.DirectorySeparatorChar + System.IO.Path.GetFileNameWithoutExtension(f) + "_corr.sgy";
if (File.Exists(f0))
{
File.Delete(f0);
}
SEGYFile s0 = new SEGYFile();
if (s0.Open(f) == 0)
{
break;
}
SEGYFile s1 = new SEGYFile();
if (s1.Open(f0) == 0)
{
s0.Close();
break;
}
s0.ReadFileHeader();
s1.Write(s0.FileHeader.Copy());
int count = 0;
toolStripStatusLabel1.Text = "Writing " + f0 + " : ";
while (s0.ReadNextTrace())
{
SEGYTrace t = s0.currentTrace.Copy(); // seismic data
if (t.TraceHeader.traceIdentificationCode != 1)
{
continue;
}
// find water depth of source in meters
var srcDepth = Convert.ToDouble(t.TraceHeader.sourceDepthBelowSurface);
var fac = t.TraceHeader.scalarForAllElevationsAndDepths;
if (fac < 0)
{
srcDepth /= -fac;
}
else if (fac > 0)
{
srcDepth *= fac;
}
// calculate the corresponding time delay in seconds
var twt = 2 * srcDepth / 1500; // assume vw = 1500 m/s
var twtmsec = twt * 1e3;
t.TraceHeader.delayRecordingTimeMsec = Convert.ToInt16(Math.Floor(twtmsec));
var remainder = twtmsec - t.TraceHeader.delayRecordingTimeMsec;
// corresponding number of samples to pad trace for fractional delay times
var nsamp = Convert.ToInt32(remainder * 1e3 / t.TraceHeader.sampleIntervalUsec);
var data = t.TraceData.Data;
// mute start of trace
if (this.checkBox1.Checked)
{
int mutedelay = Convert.ToInt32(Convert.ToDouble(this.numericUpDown1.Value) * 1e3 / Convert.ToDouble(t.TraceHeader.sampleIntervalUsec));
for (int ii = 0; ii < mutedelay; ii++)
{
data[ii] = 0;
}
}
var datanew = new double[data.Length]; // should be initializes to zero
Array.Copy(data, 0, datanew, nsamp, data.Length - nsamp);
t.TraceData.Data = datanew;
s1.Write(t);
count++;
if (count % 100 == 0)
{
toolStripStatusLabel2.Text = count.ToString();
Application.DoEvents();
}
}
s0.Close();
s1.Close();
}
toolStripStatusLabel1.Text = "."; toolStripStatusLabel2.Text = ".";
}
