R: Read one SEGY/SAC file

read1segy {RSEIS} | R Documentation

## Read one SEGY/SAC file

### Description

Read one SEGY/SAC file

### Usage

    read1segy(fname, Iendian = 1, HEADONLY = FALSE, BIGLONG = FALSE)
    read1sac(fname, Iendian = 1, HEADONLY = FALSE, BIGLONG = FALSE )

### Arguments

`fname` | character, file name
`Iendian` | Endian of the input file name
`HEADONLY` | logical, TRUE=return only header (default=FALSE)
`BIGLONG` | logical, indicating whether long is 8 or 4 bytes.

### Details

Segy format files are in integer format. The time series ususally represents counts recorded in a data acquisition system. The header includes meta-data and other identifying information.

SAC data is stored as floats, typically volts.

### Value

list of header and times series

### Note

The Endian-ness of the input files is set by the system that created them. If the read1segy or read1sac does not make sense, try a different endian or BIGLONG setting.

### Author(s)

Jonathan M. Lees<jonathan.lees@unc.edu>

### See Also

write1sac, write1segy, sac2rseis, segy2rseis, prepSEIS

### Examples

    ## Not run:

    fnames = list.files(path="/Users/lees/Site/Buf2014/ALL_SEGY",
                                         pattern=NULL, full.names=TRUE )

    for(i in 1:length(fnames))
    {
    fn1  = fnames[i]

    bed =  read1segy(fn1   , Iendian = 1 , HEADONLY=TRUE, BIGLONG=FALSE)
    scaly[i] =  bed$HEAD$scale_fac
    print(paste(i, fn1, bed$HEAD$station_name , bed$HEAD$channel_name,
           bed$HEAD$scale_fac, bed$HEAD$gainConst)  )

    }

    ## End(Not run)

[Package RSEIS version 3.7-4 Index]