las2su.sh

su 的源代码库

#!/bin/sh
# Convert a LAS file to SU format
# Author: Werner M. Heigl
# Usage: ./las2su.sh lasfile ndc nsamp start dz units
#	 All parameters are required until a program is written
#	 that extracts them from the LAS header.

set -x

### Parameter declaration
lasfile=$1	# name of LAS file
ndc=$2		# no. of data channels in LAS file
nsamp=$3	# no. of samples, usually (START-STOP)/DZ +1
start=$4	# start depth
dz=$5		# sampling interval
unit=$6		# depth units: m or ft


### separate LAS file into header and data channels
lasstrip <$lasfile

### convert from ascii to binary floats
a2b <las_data.asc n1=$ndc outpar=/dev/tty >las_data.bin

### transpose from columns to rows
transp <las_data.bin n1=$ndc n2=$nsamp verbose=1 |

### add trace headers
suaddhead ns=$nsamp | segyclean |

### set trace header fields
# need to define trid for logs
# to be added in future: get this information from LAS header
sushw key=trid,d1,f1,ntr a=300,$dz,$start,$ndc | #>las.su

### remove depth channel
# this requires a modified suwind
suwind key=tracl reject=1 >las.su

### clean up
rm ./las_data.bin ./las_data.asc

surange <las.su

exit 0