我有一个3D地震数据立方体,可以使用R进行实验。我发现R包RSEIS似乎可以加载SEG-Y文件,可以使用...
require(RSEIS)
fn <-'My.sgy'
dt <- read1segy(fn)
大文件(3G文件)被加载到具有两个字段“ HEAD”和“ Amp”的列表中,我将其解释为标题和振幅。 HEAD包含以下90个字段,我希望其中包含有关3D网格的信息-但我不确定如何解释这些字段-每个字段都有一个整数值
> names(dt$HEAD)
[1] "lineSeq" "reelSeq" "event_number"
[4] "channel_number" "energySourcePt" "cdpEns"
[7] "traceInEnsemble" "traceID" "vertSum"
[10] "horSum" "dataUse" "sourceToRecDist"
[13] "recElevation" "sourceSurfaceElevation" "sourceDepth"
[16] "datumElevRec" "datumElevSource" "sourceWaterDepth"
[19] "recWaterDepth" "elevationScale" "coordScale"
[22] "sourceLongOrX" "sourceLatOrY" "recLongOrX"
[25] "recLatOrY" "coordUnits" "weatheringVelocity"
[28] "subWeatheringVelocity" "sourceUpholeTime" "recUpholeTime"
[31] "sourceStaticCor" "recStaticCor" "totalStatic"
[34] "lagTimeA" "lagTimeB" "delay"
[37] "muteStart" "muteEnd" "sampleLength"
[40] "deltaSample" "gainType" "gainConst"
[43] "initialGain" "correlated" "sweepStart"
[46] "sweepEnd" "sweepLength" "sweepType"
[49] "sweepTaperAtStart" "sweepTaperAtEnd" "taperType"
[52] "aliasFreq" "aliasSlope" "notchFreq"
[55] "notchSlope" "lowCutFreq" "hiCutFreq"
[58] "lowCutSlope" "hiCutSlope" "year"
[61] "day" "hour" "minute"
[64] "second" "timeBasisCode" "traceWeightingFactor"
[67] "phoneRollPos1" "phoneFirstTrace" "phoneLastTrace"
[70] "gapSize" "taperOvertravel" "station_name"
[73] "sensor_serial" "channel_name" "totalStaticHi"
[76] "samp_rate" "data_form" "m_secs"
[79] "trigyear" "trigday" "trighour"
[82] "trigminute" "trigsecond" "trigmills"
[85] "scale_fac" "inst_no" "not_to_be_used"
[88] "num_samps" "max" "min"
“ AMP”字段包含7.81亿个大整数值,我希望需要使用标头值之一将其转换回十进制精度
$ amp : int [1:781504584] 1089794243 -1042038332 -672603328 -700417339 -19876061
对于那些可能熟悉SEG_Y格式的人来说,有两个问题
对于那些可能熟悉RSEIS软件包的人-一个如何显示数据部分或计划的示例将很有帮助。