R包RSEIS-请举个例子吗?

时间:2019-04-29 05:26:28

标签: r

我有一个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软件包的人-一个如何显示数据部分或计划的示例将很有帮助。

0 个答案:

没有答案