简单的segy地震数据读取


swpu 物探系 杨连刚
#include
               #i nclude
               using namespace std;

               int main()
               {
                     int i,j,ns,ntr; 
                     float data;
                     ifstream infile;
                     infile.open("E:\Flei\test.sgy",ios::binary);
                     ntr = 201;\\Segy文件总道数
                     ns = 501;\\Segy文件的采样点
                     infile.seekg(3600,ios::beg); \\跳过3600字节的卷头
                     for(i=0;i
                           infile.seekg(2 40,ios::cur); \\跳过240字节的道头
                           for(j=0;j
                                 infile.read((char*)&data,sizeof(float));
                                 \\读取每个采样点的数据,格式为float类型
                                 cout<
                     }
                 }
                 infile.close();
                 return 0;
                 }

你可能感兴趣的:(地震编程)