R пакет RSEIS - пример, пожалуйста?

Question

У меня есть трехмерный сейсмический куб данных, который я бы экспериментировал с использованием R. Я нашел R-пакет RSEIS, который, кажется, загружает файл SEG-Y, хорошо, используя ...

require(RSEIS)
fn <-'My.sgy'
dt <- read1segy(fn)

Большой(Файл 3G) загружается в список с двумя полями «HEAD» и «Amp», которые я интерпретирую как заголовок и амплитуда.HEAD имеет 90 полей, как показано ниже, и я предполагаю, что где-то будет содержаться информация о трехмерной сетке - но я не уверен, как их интерпретировать - у каждого есть одно целое значение

> 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»"содержит 781 миллион больших целочисленных значений, которые, как я ожидаю, необходимо преобразовать обратно в десятичную точность, используя одно из значений заголовка

$ amp : int [1:781504584] 1089794243 -1042038332 -672603328 -700417339 -19876061

Так что есть два вопроса для тех, кто может быть знаком с форматом SEG_Y

1. Какие значения заголовка мне следует использовать для интерпретации источника, расстояния и количества узлов для этих трехмерных сейсмических данных?
2. Какое значение заголовка контролирует считывание точности "amp"поля амплитуд?

Для тех, кто может быть знаком с пакетом RSEIS - пример того, как отобразить сечение или план данных, будет полезен.