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 = "."; }