NAME
hlsplot - map one, two or three input files to a single
output file using an HLS color scheme.
SYNOPSIS
hlsplot [ -Hfile_h ] [ -Lfile_l ] [ -Sfile_s ] [ -Ofile_out
] [ -hlsfile_hls ] [ -uhminuhmin ] [ -uhmaxuhmax ] [
-ulminulmin ] [ -ulmaxulmax ] [ -usminusmin ] [ -usmaxusmax
] [ -chminchmin ] [ -chmaxchmax ] [ -clminclmin ] [ -clmax-
clmax ] [ -csmincsmin ] [ -csmaxcsmax ] [ -tstarttstart ] [
-tendtend ] [ -starttracestarttrace ] [ -endtraceendtrace ]
[ -startrecstartrec ] [ -endrecendrec ] [ -V ] [ -? ]
DESCRIPTION
hlsplot reads in one, two or three USP format files and maps
them to HLS color space.
hlsplot gets all its parameters from command line arguments.
These arguments specify one, two or three input files, the
output files and data limits used for scaling.
-H file_h
Enter the input data set name or file to be mapped to
hue immediately after typing -H.
-L file_l
Enter the input data set name or file to be mapped to
lightness immediately after typing -L.
-S file_s
Enter the input data set name or file to be mapped to
saturation immediately after typing -S.
-O file_out
Enter -O followed by the output file name to produce a
USP format multicomponent or multiatribute display. The
values in file_out have been chosen such that (when
converted to an 8 bit integer) the left most 6 bits
correspond to a valid seisworks, iesx, promax or xsd
color table. This color table corresponds to the HLS
color model described below and is generated using
using program hls2rgb on file_hls.
-hls file_hls
Enter -hls followed by the hls table file name to out-
put an ascii flat file containing the hue, lightness
and saturation of each sample contained in the output
file_out file. This file will be input to program
hls2rgb to generate the RGB (red, green, blue) color
lookup table needed to load into programs xsd, iesx,
promax traceDisplay or seisworks in order to properly
display file_out. See NOTES 1 and 2 below. (Default
file_hls = hls.table)
-tstart tstart
Enter the beginning of the analysis window in ms. The
output record will be tend-tstart ms long. (Default =
first sample, or 0 ms).
-tend tend
Enter the end of the analysis window in ms. The output
record will be tend-tstart ms long. (Default = last
sample of the trace)
-starttrace starttrace
Enter the first output trace to be processed (Default:
starttrace=1, the first trace found in each seismic
record).
-endtrace endtrace
Enter the last output trace to be processed (Default:
endtrace= the last trace to be found in each seismic
record).
-startrec startrec
Enter the first output record to be processed (Default:
startrec=1, the first record of the data set).
-endrec endrec
Enter the last output line to be processed (Default:
endrec= the last lin e of the data set).
-chmin chmin
Enter the minimum hue in degrees to be used in mapping
file_h. See Note 1 below. (Default: chmin==-180.)
-chmax chmax
Enter the maximum hue in degrees to be used in mapping
file_h. See Note 1 below. (Default: chmax==+180.)
-clmin clmin
Enter the minimum percent lightness to be used in map-
ping file_l. (Default: clmin==10.)
-clmax clmax
Enter the maximum percent lightness to be used in map-
ping file_l. (Default: clmax==100.)
-csmin csmin
Enter the minimum percent saturation to be used in map-
ping file_s. (Default: csmin==0.)
-csmax csmax
Enter the maximum percent saturation to be used in map-
ping file_s. (Default: csmax==100.)
-uhmin uhmin
Enter the value of file_h that corresponds to the
minimum hue chmin . (Default: minimum value present in
the data).
-uhmax uhmax
Enter the value of file_h that corresponds to the max-
imum hue chmax . (Default: maximum value present in the
data).
-ulmin ulmin
Enter the value of file_l that corresponds to the
minimum lightness clmin . (Default: minimum value
present in the data).
-ulmax ulmax
Enter the value of file_l that corresponds to the max-
imum lightness clmax . (Default: maximum value present
in the data).
-usmin usmin
Enter the value of file_s that corresponds to the
minimum saturation csmin . (Default: minimum value
present in the data).
-usmax usmax
Enter the value of file_s that corresponds to the max-
imum saturation csmax . (Default: maximum value present
in the data).
-V Enter the command line argument '-V' to get additional
printout.
-? Enter the command line argument '-?' to get online
help. The program terminates
after the help screen is printed.
Multiattribute Display Implementation:
The values of file_h, file_l and file_s are mapped to
the HLS (hue, lightness, saturation) color model as follows:
The hues are pure, or 100% saturated colors, and correspond to the following 1994 96 crayon Crayola standard:
hue Crayola Color
0 blue
30 plum
60 magenta
90 salmon
120 red
150 orange-red
180 yellow
210 lime-green
240 green
270 forest-green
300 cyan
330 cerulean
360 blue
Partial 50% saturation corresponds to 'dirtier' or 'muddier' colors:
@ hue @ Crayola Color (hue)
0 cadet blue
45 fuscia
90 maroon
135 sepia
180 gold
225 olive
270 sea green
315 steel blue
360 cadet blue
0% saturation corresponds to no color pigments:
hue Crayola Color (hue)
0 grey
90 grey
180 grey
270 grey
360 grey
Low values of lightness correspond to 'dark' colors.
Intermediate values of lightness correspond to 'deep' colors.
High values of lightness correspond to 'pastel' colors.
The user is referred to Marfurt et al. (1995) for a more detailed description.
NOTE 1: Counterclockwise mapping of hue.
If we wish to map low values of file_h to magenta (h=60) and
high values of hue to red (h=120) with intermediate values
being mapped to blue, cyan, green and yellow, we need to
unwrap the color wheel and define clmin=420 and clmax=120.
Setting clmin=60 will result in a clockwise interpolation
and intermediate values being mapped to salmon (eg. h=90).
NOTE 2: Loading multicomponent color displays
To generate the RGB color table corresponding to the
file_hls file for use
in program xsd the user should type in :
hls2rgb -Nhls.table -X -V -Ocolor_bar.xsd
Once in xsd, the user loads in file_out (or
file_color_legend) with the Fixed Scaling option chosen, and
with Scalar=.0025, Offset= 0%. The user then loads in the
color table out_color_bar.xsd.
NOTE 3: Loading multicomponent color displays
To generate the RGB color table corresponding to the
file_hls file for use
in Landmark's program SeisWorks 3D the user should first
type in :
hls2rgb -Nhls.table -L -V -Oout_color_bar.clm
Within SeisWorks, the user
1) opens a 'Seismic View' window
2) selects the 'Seismic' option from the control bar, which results in a pull do
wn menu,
3) selects 'parameters' on the pull down menu, which opens the "Seismic Display
Parameters"
dialog box.
In this dialog box, the user
1) Under the 'Seismic Files' heading selects the 3dv format solid angle file lo
aded from fBstolm3dv.
(eg out01) as both the Vertical (Inline) and Vertical (Crossline) file to be
displayed.
2) Under the 'Seismic Display Scales' heading, sets the 'Variable Density Interp
olation Method' = 'constant',
3) Under the 'Dynamic Scale and Clip' heading, selects the 'No Scaling' option,
and
sets the "largest unclipped Amplitude"=127, and
4) Clicks the 'OK' button to apply and remove the dialog box.
In the 'Seismic View' window, the user selects the color control bar ikon, and o
pens the .clm color file 'out_color_bar'
created above.
NOTE 4: Loading multicomponent color displays
To generate the RGB color table corresponding to the
file_hls file for use
in Advance's PROMAX plotting software traceDisplay the user
should first type in under USP:
hls2rgb -Nhls.table -P -V -Oout_color_bar.rgb
Within promax's traceDisplay, under 'controls' the user clicks 'trace display',
then within the pop-up window:
1) clicks display type='variable density'
2) clicks colormap='color'
3) clicks interpolation='nearest sample'
4) clicks 'ok'
Next, under 'controls' the user clicks 'trace scaling', then within the pop-up window:
1) clicks polarity='normal'
2) clicks scaling='range limited'
3) clicks estimate from data='no'
4) types minimum='-128'
5) types maximum='+127'
6) clicks 'ok'.
Finally, under 'controls' the user clicks 'edit colormap', then within the pop-up window:
1) clicks 'file'
2) clicks 'open'
3) types in the path to the above generated color map,
4) clicks 'out_color_bar.rgb' or the appropriate name given to program hls2rgb.
5) clicks 'ok'.
EXAMPLE 1:
run a data cube through the D1 spectral decomposition algo-
rithm tune3d, extracting the frequency having the peak
amplitude (file peakf_file) and the peak amplitude itself
(file apeakf_file) using routine tune3d. Scale the amplitude
cube to a known range of values. Combine both these attri-
butes using hlsplot.
#
# define named pipes:
#
/etc/mknod peakf_file p
/etc/mknod apeakf_file p
/etc/mknod apeakf_file_scaled p
#launch program hlsplot in the background.
# Note that we will assign a maximum lightness to those values 150% of the rms value (1.50*205=300.).
# the limits on peak frequency are better known, lieing between 5 and 70 Hz.
#
hlsplot -L apeakf_file_scaled -H peakf_file -O hlsplot.freq.amp
-hls hls_table -nh 18 -nl 11
-clmin 0 -clmax 60. -ulmin 0. -ulmax 300.
-chmin 60 -chmax 360 -uhmin 5 -uhmax 70 &
#
#launch scaling program gasp in the background
#
# scale the amplitude of each line of the 3-D survey to fall between 0 and 2047.
# if the rms value is set to be 10% of 2047, then this corresponds to 205.
# use record constant scaling to allow us to use named pipes.
#
gasp -Napeakf_file -Oapeakf_file_scaled -rec -s 10. &
#
#
#launch program tune3d in the foreground.
#
tune3d -N seismic.cube -fmin 5 -fmax 70 -peakf peakf_file -apeakf apeakf_file
#
# remove named pipes
#
"rm" -f peakf_file apeakf_file apeakf_file_scaled
#
See Also:
asig3d, asig1d, correl3d, princ3d, semb3d, hls2rgb asig3d
REFERENCES:
Marfurt, K.J., Kirlin, R.L., Farmer, S. F. and Bahorich,
M.S, 1995, 3D seismic attributes using a running window sem-
blance technique. Amoco Geos. Tech. Bull. F94-G-55.
Bucher, R.H., Marfurt, K.J. and Stanley, T.D, 1988, U.S.
Patent No. 4,970,699: M ethod for color mapping geophysical
data.
Partyka, G.A., Gridley, J., Mims, C.V., and Lopez, J.A..,
1995, The tuning cube: Interpretational aspects of spectral
decomposition: Amoco Geos. Tech. Bull. F95-G-46.
CONTRACT AGREEMENT
This product is brought to you by Research Agreement D96-
2548 (The seismic coherency cube) as a 2nd quarter deliver-
able entitled "Development and geostatistical calibration of
new thin bed analysis algorithms". Thank you for your sup-
port.
AUTHOR
Kurt. J. Marfurt (E&PTG, Tulsa, OK, USA. 1st Quarter,1996).
COPYRIGHT
copyright 2001, Amoco Production Company
All Rights Reserved
an affiliate of BP America Inc.
Man(1) output converted with
man2html