r - R包RSEIS - 请举例？

Question

我有一个 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 格式的人来说有两个问题

1. 我应该使用哪些标题值来解释此 3D 地震数据的原点、间距和节点数？
2. 什么标题值控制“amp”幅度字段的精度读数？

对于那些可能熟悉 RSEIS 包的人 - 如何显示数据的部分或计划的示例会有所帮助。

标签： r