大型地震数据的读取和修改软件

对大型地震数据（200G左右的数据均可）读取以及修改道头和截取任意道，以及输出最大振幅等功能---已经测试过

#include
#include
#include
#include“SEGY_RW.h“
#include“write_txt.h“
#include“Max.h“

using namespace std;

void main（）
{
	float a[5]={1.0 2 3 4 -5};
	cout<
	string Input_Seismic_FileOutput_Seismic_File;
	int code;
	//functions
		cout<<“************Choose the function*************“<		cout<<“function codes: “<		cout<<“1. 重置道头“<			<<“2. 截取segy“<			<<“3. 输出单道地震数据ASCII“<		    <<“4. 输出各道最大振幅“<			
	    cout<<“Enter the code of function: “<		cin>>code;		
		cout<<“======Input Seismic Data Name======“<		cin>>Input_Seismic_File;
		switch（code）
		{
		case 1:  Output_Seismic_File=“hdr_“+Input_Seismic_File;
			break;
		case 2:  Output_Seismic_File=“trun_“+Input_Seismic_File;
			break;
     	case 3:  Output_Seismic_File=“trace_“+Input_Seismic_File;
			break;
		case 4:
			Output_Seismic_File=“Max_“+Input_Seismic_File;
			break;
	//	default:  cout<<“No function is choosed!“<		//	break;
		}
 
	SEGY_RW segy_rw（Input_Seismic_FileOutput_Seismic_File）;

	int TraceNumber=segy_rw.Return_TraceNumber（）;
	short SampleNumber=segy_rw.Return_SampleNumber（）;
	float SampleInterval=segy_rw.Return_SampleInterval（）;
	SampleInterval=SampleInterval;
	//Display the Seismic Data Property
	cout<<“>>========================<<“<	//cout<<“>>Seismic Data Name: “<	cout<<“>>Input Seismic Trace Number: “<	segy_rw.Read_TraceHeader（1）;
	int first=segy_rw.Return_Trace（）;
	segy_rw.Read_TraceHeader（TraceNumber）;
	int end=segy_rw.Return_Trace（）;
	cout<<“Trace Number:  from “<	cout<<“>>Input Seismic Sample Interval: “<	cout<<“>>Input Seismic Sample Number: “<	//cout<<“>>Output Q File: “<	cout<<“>>========================<<“<    float *one_trace=new float[SampleNumber];
		switch（code）
		{
		case 1:  
			
		case 2:  {	
		 
	       int s_trae_tra;
		   cout<<“输入起始道： “<		   cin>>s_tra;
		   s_tra=s_tra-（first-1）;
		   cout<<“输入终止道： “<		   cin>>e_tra;
		   e_tra=e_tra-（first-1）;
          for（int i=s_tra;i<=e_tra;i++）
	         {
		segy_rw.Read_TraceHeader（i）;
		segy_rw.Read_TraceData（one_tracei）;
		segy_rw.Write_TraceHeader（）;
		segy_rw.Write_TraceData（one_trace）;
		cout<<“trace: “<	         }	
				 }
           break;  

		case 3: {
			int trace;
			cout<<“输出地震道： “<			cin>>trace;
		    //segy_rw.Read_TraceHeader（trace）;
		    segy_rw.Read_TraceData（one_tracetrace）;
           
			ofstream Output_File;
	//Write the stratify data matrix
	        Output_File.open（Input_Seismic_File+“.txt“ios::out）;//不能有in
	         if （Output_File.bad（）） cout<<“bad!“;
	         else
	        {
	             for（int i=0;i	           {
	              Output_File.setf（ios::fi

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2012-12-29 20:43  wr\
     目录           0  2012-12-29 20:22  wr\Debug\
     文件      132608  2012-12-29 20:42  wr\Debug\wr.exe
     文件      980640  2012-12-29 20:42  wr\Debug\wr.ilk
     文件     1059840  2012-12-29 20:42  wr\Debug\wr.pdb
     目录           0  2012-12-29 16:33  wr\ipch\
     目录           0  2012-12-29 20:29  wr\ipch\wr-d38eb0c9\
     文件    17235968  2012-12-29 20:29  wr\ipch\wr-d38eb0c9\wr-5bb99bd5.ipch
     目录           0  2012-12-29 20:42  wr\wr\
     文件     5984256  2012-12-29 20:43  wr\wr.sdf
     文件         873  2012-12-21 11:25  wr\wr.sln
     文件       26112  2012-12-29 20:43  wr\wr.suo
     目录           0  2012-12-29 20:42  wr\wr\Debug\
     文件        2846  2012-12-29 20:42  wr\wr\Debug\cl.command.1.tlog
     文件       40186  2012-12-29 20:42  wr\wr\Debug\CL.read.1.tlog
     文件        3616  2012-12-29 20:42  wr\wr\Debug\CL.write.1.tlog
     文件           2  2012-12-29 20:42  wr\wr\Debug\link-cvtres.read.1.tlog
     文件           2  2012-12-29 20:42  wr\wr\Debug\link-cvtres.write.1.tlog
     文件           2  2012-12-29 20:42  wr\wr\Debug\link.11272-cvtres.read.1.tlog
     文件           2  2012-12-29 20:42  wr\wr\Debug\link.11272-cvtres.write.1.tlog
     文件           2  2012-12-29 20:42  wr\wr\Debug\link.11272.read.1.tlog
     文件           2  2012-12-29 20:42  wr\wr\Debug\link.11272.write.1.tlog
     文件           2  2012-12-29 20:42  wr\wr\Debug\link.12396-cvtres.read.1.tlog
     文件           2  2012-12-29 20:42  wr\wr\Debug\link.12396-cvtres.write.1.tlog
     文件           2  2012-12-29 20:42  wr\wr\Debug\link.12396.read.1.tlog
     文件           2  2012-12-29 20:42  wr\wr\Debug\link.12396.write.1.tlog
     文件           2  2012-12-29 20:42  wr\wr\Debug\link.7820-cvtres.read.1.tlog
     文件           2  2012-12-29 20:42  wr\wr\Debug\link.7820-cvtres.write.1.tlog
     文件           2  2012-12-29 20:42  wr\wr\Debug\link.7820.read.1.tlog
     文件           2  2012-12-29 20:42  wr\wr\Debug\link.7820.write.1.tlog
     文件        2128  2012-12-29 20:42  wr\wr\Debug\link.command.1.tlog
............此处省略32个文件信息