[其他]地震处理软件开发 GeoSong [复制链接]

建议你写一个长期规划。
基础工作一定要做好，尤其是数据管理这块。版本号从0.2开始算了，现在就是1.2版本，就有点太像某些国产软件了。

qt前景堪忧，gtk在开源社区更受支持

GeoSong系统里面中间数据都是用SEG-Y格式，因此需要把这个格式讲解一下先，懂的人就可以不用看。

一个SEG-Y格式的数据整体可以分为两个部分：文件卷头和数据。
头分为EBCDIC header和二进制数区。EBCDIC header在物探局处理中心常称作C卡，内容是3200个文本字符，用来对数据做一些说明。二进制数区则是400个字节的一些控制信息，其中比较有用的信息有
3213~3214 字节?每个记录的数据道数(每炮道数或总道数)。
3217~3218 字节?采样间隔(μs)。
3221~3222 字节?样点数/每道(道长)。
3225~3226 字节?数据样值格式码1-浮点；
3255~3256 字节?计量系统：1-米， 2-英尺。
3261~3262* 字节?文件中的道数(总道数)。
3269~3270* 字节?数据域(性质)：0-时域，1-振幅，2-相位谱
所以整个卷头部分共3600字节。


数据部分以道（trace）为单位，每道又由两个部分组成既道头和数据。道头包含240个字节信息，数据一般为n个样点×4字节浮点。

这么算来一个数据的内容就是  
3200字节EBCDIC + 400字节文件头 + {trace1, trace2, … {240字节道头+{sample1,sample2,…{4字节浮点数}…,samplem}}…,tracen }。

如果你要写一个读取SEG-Y数据的程序，你需要注意的东西有：
1，    在流文件和某些硬件条件下，一个4字节数的低位存储在前高位在后，然而其它大多数情况下（例如IBM PC兼容环境）都是刚好相反的。所以如果我们要读一个4字节浮点数表示的样点振幅，我们需要把这四个字节倒个顺序先。如果我们需要读2字节道头的第9-10字节，那么我们实际要把第11-12字节取出来再倒个个儿读；
2，    SEG-Y格式的浮点数支持IBM浮点数和IEEE浮点数两周标准。虽然IBM浮点是默认标准，但是IEEE在PC中应用的更加广泛，我是喜欢把文件保存成IEEE的，这样每次计算就不用多一次转格式的过程了。具体数据是哪个格式的，文件头400字节的信息就有用了，读取指定字节的内容，它标注了这个信息。

了解了SEG-Y格式并且注意以上两点，读取一个SEG-Y数据就不是问题了。

好了，程序员的问题解决了，下面处理员最关心的就是道头了。SEG-Y对于道头没有一致的标准，并且约定俗成的用法过时了，有些要表示的道头没有涉及到。
1-4* 一条测线中的道顺序号。如果一条测线有若干卷带，顺序号连续递增。
  5-8 在本卷磁带中的道顺序号。每卷带的道顺序号从1开始。
  9-12* 原始的野外记录号。
  13-16* 在原始野外记录中的道号。
  17-20 震源点号(在同一个地面点有多于一个记录时使用)。
  21-24 CMP号。
  25-28 在CMP道集中的道号(在每个CMP道集中道号从1开始)。
  29-30* 道识别码：
1=地震数据；4=时断；7=记时；
2=死道；5=井口时间；8=水断；
3=DUMMY；6=扫描道；9…N=选择使用(N=32767)
  31-32 产生这一道的垂直叠加道数(1是一道； 2是两道相加；…)。
  33-34 产生这一道的水平叠加道数(1是一道； 2是两道叠加；…)。
  35-36 数据类型：1=生产；2=试验。
37-40 炮检距(如果是相反向激发为负值)。
41-44 接收点高程。高于海平而的高程为正，低于海平面为负。
45-48 炮点的地面高程。
49-52 炮点低于地面的深度(正数)(井深)。
53-56 接收点的基准面高程。
57-60 炮点的基准面高程。
61-64 炮点的水深。
65-68 接收点的水深。
69-70 对41-68字节中的所有高程和深度应用了此因子给出真值。
比例因子=1，±10，±100，±1000或者±10000。如果为正，乘以因子； 如果为负，则除以因子。
18-2 71-72 对73-88字节中的所有坐标应用了此因子给出真值。比例因子=1，±10，±100，±1000或者±10000。如果为正，乘以因子；如果为负，则除以因子(在GRISYS中为10)。
73-76 炮点坐标--X |- 如果坐标单位是弧度的秒，X值代表径度，
77-80 炮点坐标--Y | Y值代表纬度。正值代表格林
81-84 检波点坐标--X | 威治子午线东或者赤道北的秒数。负
85-88 检波点坐标--Y |- 值则为西或者南的秒数
89-90 坐标单位； 1=长度(米或者英尺)； 2=弧度的秒。
91-92 风化层速度。
93-94 降速层速度。
95-96 震源处的井口时间。
97-98 接收点处的井口时间。
99-100 炮点的静校正。
101-102 接收点的静校正。
103-104 应用的总静校正量(如果没有应用静校正为零)。
105-106 延迟时间-A，以 ms 表示。240字节的道标识的结束和时间信号之间的时间。如果时间信号出现在道头结束之前为正。如果时间信号出现在道头结束之后为负。时间信号就是起始脉冲，它记录在辅助道上或者由记录系统指定。
107-108 时间延迟-B，以 ms 表示。为时间信号和能量起爆之间的时间。可正可负。
109-110 时间延迟时间，以 ms 表示。能量源的起爆时间和开始记录数据样点之间的时间(深水时，数据记录不从时间零开始。)
111-112 起始切除时间。
113-114 结束切除时间。
115-116* 本道的采样点数。
117-118* 本道的采样间隔，以 ms 表示。
119-120 野外仪器的增益类型：
1=固定增益；2=二进制增益；3=浮点增益；4…N=选择使用。

引用第11楼dedream于2011-03-06 17:50发表的  :
qt前景堪忧，gtk在开源社区更受支持

nokia收购了QT，这会儿又跟MS合作下一代手机OS了，估计是不乐观。不过我先用着吧，挺好的

引用第12楼silverwing于2011-03-06 18:28发表的  :
GeoSong系统里面中间数据都是用SEG-Y格式，因此需要把这个格式讲解一下先，懂的人就可以不用看。
一个SEG-Y格式的数据整体可以分为两个部分：文件卷头和数据。
头分为EBCDIC header和二进制数区。EBCDIC header在物探局处理中心常称作C卡，内容是3200个文本字符，用来对数据做一些说明。二进制数区则是400个字节的一些控制信息，其中比较有用的信息有
3213~3214 字节?每个记录的数据道数(每炮道数或总道数)。
.......

SEGY格式与IEEE格式的转换还是需要一定时间的，处理任务数据量巨大啊。原始数据可以直接从SEGY中读取，但中间成果写成一个内部格式或者IEEE格式的数据应该能提高效率。

只有默认的SEGY标准格式，不合适。即使是Y格式，实际的道头和文件头也是各种各样
道头和文件头，可以给一个标准的作为默认的，还应该可以让用户指定这些头的格式和含义，程序按用户定义去读数据。
以前写过一个文件转换程序，就可读多种格式的文件（有头无头，非标准头等）、IBM－IEEE互转。
MPICH可以在单台机器上运行多个结点，实现相互间的通信。
建议LZ有时间看看软件工程学之类的书，做好系统的规划设计，不然以后模块、数据多了以后，再回头来改，工作量很大
祝LZ成功

支持开发啊 不过相对而言 数据处理软件的开发实在是太难了 一个人很难走到最后的 开发的东西必须要有特色才行啊 否则的话根本支持不到最后 而这特色看你的描述如果是能做一个很易用的借口那也行 呵呵 支持

支持原创！！！！！！

很好，学习下~欢迎探讨~