Java读取Segy文件并利用vtk显示的例子

package org.sinoprobe.segy; /** * Segy道 * 包括道头和道数据 * 参照Segy 2002年的标准编写 * 下面是2002年标准中道头数据的说明 * * * 表3 道头 240字节道头 字节 描述 1-4 测线中道顺序号――若一条测线有若干SEG Y文件号数递增。强烈推荐所有类型的数据使用 5-8 SEG Y文件中道顺序号――每个文件以道顺序1开始。 9-12 野外原始记录号。强烈推荐所有类型的数据使用 13-16 野外原始记录的道号。强烈推荐所有类型的数据使用 17-20 震源点号――当在相同有效地表位置多于一个记录时使用。建议在道头197-202字节定义新的条目用于炮点号。 21-24 道集号（即CDP，CMP，CRP等） 25-28 道集的道数――每个道集从道号1开始 29-30 道识别码： -1=其他 0＝未知 1＝地震数据 2＝死道 3＝哑道 4＝时断 5＝井口 6＝扫描 7＝定时 8＝水断 9＝近场枪信号 10＝远场枪信号 11＝地震压力传感器 12＝多分量地震传感器――垂直分量 13＝多分量地震传感器――横向分量 14＝多分量地震传感器――纵向分量 15＝旋转多分量地震传感器――垂直分量 16＝旋转多分量地震传感器――切向分量 17＝旋转多分量地震传感器――径向分量 18＝可控源反应质量 19＝可控源底盘 20＝可控源估计地面力 21＝可控源参考 22＝时间速度对 23…N＝选用，（最大N＝32767） 强烈推荐所有类型的数据使用 31-32 产生该道的垂直叠加道数。（1是一道，2是两道求和，…） 33-34 产生该道的水平叠加道数。（1是一道，2是两道求和，…） 35-36 数据用途： 1＝生产 2＝试验 37-40 从震源中心点到检波器组中心的距离（若与炮激发线方向相反取负） 41-44 检波器组高程（所有基准以上高程为正，以下为负） 道头中69-70字节的因子应用于这些数值。在二进制文件头3255-3256字节指定单位为英尺或米。垂直基准应通过位置数据文本段定义（见D-1节）。 45-48 震源地表高程 49-52 震源距地表深度（正数） 53-56 检波器组基准高程 57-60 震源基准高程 61-64 震源水深 65-68 检波器组水深 69-70 应用于所有在道头41-68字节给定的真实高程和深度的因子。因子＝1，＋10，＋100，或＋1000。若为正，因子为乘数；若为负，因子为除数。 71-72 应用于所有在道头73-88字节和181-188字节给定的真实坐标值的因子。因子＝1，＋10，＋100，或＋1000。若为正，因子为乘数；若为负，因子为除数。 73-76 震源坐标――X 参考坐标系应通过扩展头位置数据文本段识别（见D-1节）。 若坐标单位为弧度秒、小数度或度/分/秒(DMS)，X值代表经度，Y值代表纬度。正值代表格林威治子午线以东或赤道以北，负值代表南或西。 77-80 震源坐标――Y 81-84 检波器组坐标――X 85-88 检波器组坐标――Y 89-90 坐标单位： 1＝长度（米或英尺） 2＝弧度秒 3＝小数度 4＝度，分，秒（DMS） 注意：为编码±DDDMMSS 89-90字节等于±DDD*104+MM*102+SS，71-72字节设置为1；为编码±DDDMMSS 89-90字节等于±DDD*106+MM*104+SS*102，71-72字节设置为-100。 91-92 风化层速度（如二进制文件头3255-3256字节指明的ft/s或m/s） 93-94 风化层下速度（如二进制文件头3255-3256字节指明的ft/s或m/s） 接表3 道头 95-96 震源处井口时间（毫秒） 毫秒表示的时间需应用道头215-216字节指定的因子。 97-98 检波器组处井口时间（毫秒） 99-100 震源的静校正量（毫秒） 101-102 检波器组的校正量（毫秒） 103-104 应用的总静校正量（毫秒）（如没有应用静校正量为零） 105-106 延迟时间A――以毫秒表示的240字节道识别头的结束和时间断点之间的时间。当时间断点出现在头之后，该值为正；当时间断点出现在头之前，该值为负。时间断点是最初脉冲，它由辅助道记录或由其他记录系统指定。 107-108 延迟时间B――以毫秒表示的时间断点到能量源起爆时间之间的时间。可正可负。 109-110 记录延迟时间――以毫秒表示的能量源起爆时间到数据采样开始记录之间的时间。在SEG Y修订版0中本条用来表示深水作业，如果数据记录不从0时间开始。该条可为负值以适应负的起始时间（即数据记录在零时间之前，假设静校正量应用于数据道的结果）。若某非零值（正或负）记录在该条，它造成的影响的注释应出现在原文文件头。 111-112 起始切除时间（毫秒） 113-114 终止切除时间（毫秒） 115-1167 该道采样点数。强烈推荐所有类型的数据使用 117-1187 该道采样间隔（微秒）。 一道记录的字节数必须和写在道头中的采样点数一致。对所有记录介质都重要，但对正确处理磁盘文件中的SEG Y数据尤为关键（见附录C）。 若二进制文件头中的3503-3504字节设置了固定长度道标志，SEG Y文件每道的采样间隔和采样点数必须与二进制文件头所记录的值一致。若没有设置固定长道标志，采样间隔和采样点数可能每道变化。 强烈推荐所有类型的数据使用 119-120 野外仪器增益类型： 1＝固定 2＝二进制 3＝浮点 4…N＝选用 121-122 仪器增益常数（分贝） 123-124 仪器初始增益（分贝） 125-126 相关： 1＝无 2＝有 127-128 起始扫描频率（赫兹） 129-130 终止扫描频率（赫兹） 131-132 扫描长度（毫秒） 133-134 扫描类型： 1＝线性 2＝抛物线 3＝指数 4＝其他 135-136 扫描道斜坡起始长度（毫秒） 137-138 扫描道斜坡终止长度（毫秒） 139-140 斜坡类型： 1＝线性 2＝cos2 3＝其他 141-142 假频滤波频率（赫兹），若使用 143-144 假频滤波坡度（分贝/倍频程） 145-146 陷波频率（赫兹），若使用 147-148 陷波坡度（分贝/倍频程） 149-150 低截频率（赫兹），若使用 151-152 高截频率（赫兹），若使用 153-154 低截坡度（分贝/倍频程） 155-156 高截坡度（分贝/倍频程） 157-158 数据记录的年――1975年标准没有说清应该记成2位还是4位还是两者都用。除SGE Y修订版0外，年份需记录成完整的4位罗马历法年（即2001年应记录为200110（7D116））。 159-160 日（以格林尼治标准时间和通用协调时间为基准的公元日） 161-162 时（24小时制） 163-164 分 165-166 秒 167-168 时间基准码： 1＝当地 2＝格林尼治标准时间 3＝其他，应在扩展原文文件头的用户定义文本段解释 4＝通用协调时间 169-170 道加权因子――最小有效位数定义为2-N伏。（N＝0，1，…，32767） 171-172 滚动开关位置1的检波器组号 173-174 野外原始记录中道号1的检波器组号 175-176 野外原始记录中最后一道的检波器组号 177-178 间隔大小（滚动时甩掉的总检波器组数） 179-180 相对测线斜坡起始或终止点的移动 1＝下（或后） 2＝上（或前） 181-184 该道的道集（CDP）位置X坐标（应用道头71-72字节的因子）。参考坐标系应通过扩展头位置数据文本段识别（见D-1节） 185-188 该道的道集（CDP）位置Y坐标（应用道头71-72字节的因子）。参考坐标系应通过扩展头位置数据文本段识别（见D-1节） 189-192 对于3-D叠后数据，本字段用来填纵向线号（In-line）。若每个SEG Y文件记录一条纵向线，文件中所有道的该值应相同，并且同样的值将记录在二进制文件头的3205-3206字节中。 193-196 对于3-D叠后数据，本字段用来填横向线号（Cross-line）。它应与道头21-24字节�