Nadine E. Golden
Guy R. Cochrane
2007
chancen_bath.txt - Bathymetric sidescan sonar data of Channel Central
in the Nearshore Benthic Habitat Mapping Project S. California map series.Edition: Version 2.0, 2010
Larger_Work_Citation Citation_Information
Series_Information Series_Name
Publication_Information Publication_Place
Online_Linkage https://pubs.usgs.gov/of/2007/1271/
Bathymetric sidescan sonar data of the nearshore seafloor (0 to 100 m water depths) of the Santa Barbara Channel area collected in 2005. A 234 kHz Interferometric Submetix Swath Bathy Sonar System; SEA Swath Processor, v. 2.05, SEA Grid Processor v. 2.05, was used for geophysical surveying. The 2005 survey was navigated with a CodaOctopus, Model F180, Differential Global Positioning System (DGPS). A KVH Industries Inc. azimuth digital gyro-compass provided ship headings with 0.5 degree accuracy. Navigation data were recorded using Yo-Nav version 1.19 (Gann, 1992). This is one of a collection of digital files of a geographic information system of spatially referenced data related to the USGS Coastal and Marine Geology Program Nearshore Benthic Habitat Mapping Project. See <http://walrus.wr.usgs.gov/nearshorehab> for more information.
These data are intended for science researchers, students, policy
makers, and the general public. The data can be
used with geographic information systems (GIS) software
to display geologic and oceanographic information.
Any use of trade, product, or firm names is for descriptive purposes
only and does not imply endorsement by the
U.S. Government.
Although this Federal Geographic Data Committee-compliant metadata
file is intended to document the data set in nonproprietary form,
as well as in ArcInfo format, this metadata file may include some
ArcInfo-specific terminology.
20050808
20050827
ground condition
As needed
-119.643628
-119.471725
34.390598
34.320467
None
bathymetry
sidescan sonar
remote-sensing
MIPS
U.S. Geological Survey
USGS
Coastal and Marine Geology Program
CMGP
ISO 19115 Topic Category
geoscientificinformation
imagereyBaseMapsEarthCover
oceans
elevation
None
Channel Islands National Marine Sanctuary
Geographic Names Information System (GNIS)
Santa Barbara
Ventura
Santa Barbara Channel
Channel Islands
California
None
Public domain data from the U.S. government is freely redistributable
with proper metadata and source attribution.
Please recognize the U.S. Geological Survey (USGS)
as the source of this information.
Guy R. Cochrane
U.S. Geological Survey, Coastal and Marine Geology
Geophysicist
mailing and physical address
USGS Pacific Science Center, 400 Natural Bridges Drive
Santa Cruz
CA
95060-5792
USA
(831) 427-4754
(831) 427-4748
gcochrane@usgs.gov
images/chancen_bath.jpg
Image showing coverage of
Channel Central bathymetry data ArcInfo format.
JPEG
SunOS, 5.7, sun4u UNIX
Windows XP Professional, v. 2002, service pack 2
SEA Grid Processor v. 2.051
ArcGIS version 9.1
U.S. Geological Survey, Coastal and Marine Geology Program
2007
InfoBank
http://walrus.wr.usgs.gov/infobank/s/s105sc/html/s-1-05-sc.meta.html
Guy Cochrane, U.S. Geological Survey, Coastal and Marine Geology Program
200510
Sonar Survey of Sea-Floor Habitats
Southeast of Santa Barbara, California
in Sound Waves Monthly Newsletter,
Coastal Science & Research News from Across the USGS
http://soundwaves.usgs.gov/2005/10/fieldwork2.html
Logical Consistency untested.
Please see the Methods section of https://pubs.usgs.gov/of/2007/1271/ for
information about omissions, selection criteria, generalization,
definitions used, and other rules used to derive the data set.
On the order of 2 meters.
Started with text files that were the output from A 234 kHz Interferometric Submetix
Swath Bathy Sonar System.
20060702
>Ran Python script, "Bathy Processor v1.5", designed to filter the bathymetry text file
>
> ****Bathy_Proc (Python Script)****
> #!/usr/bin/python
> # Process bathymetry for P1-06-MB
> # Once this stabalizes, I will code into either a batch file or something more robust
> import os, glob, sys, time
>
> # Configure Directories
> SCRIPT_DIR = "E:/subx_filters"
> TXT_DIR = "D:/S205SC/Bathy/TXT"
> GRID_DIR = "D:/S205SC/Bathy/asciigrids"
>
> # Trackline Filtering Parameters
> MAX_SWATH_MULTIPLE = 4.0 # Trim depth samples > MAX * Water Depth (trims far-field noise)
> MIN_SWATH_MULTIPLE = 0.0 # Trim depth samples < MIN * Water Depth (trims nadir noise)
>
> ROLL_BEAMS = 150 # Number of beams to average to eliminate roll artifacts
> NOISE_BEAMS = 300 # Number of beams to average to eliminate along-track noise
>
> # Output Grid Parameters
> GRID_RESOLUTION = 1 # Output grid resolution
> SEARCH_RADIUS = GRID_RESOLUTION + 5 # Search radius for samples while gridding
>
> # Get the files for processing
> files = glob.glob("%s/*txt" % TXT_DIR)
>
> # Get the total size of the input files to estimate processing progress
> bytes_total = sum([os.path.getsize(x) for x in files])
> bytes_sofar = 0
>
> for i, file in enumerate(files):
> time_start = time.time()
> (indir, infile) = os.path.split(file)
> (basename, ext) = os.path.splitext(infile)
>
> print
> print "-------------------------------------------------------------------------"
> print "Begin Processing line: %s (%d of %d)" % (basename, i+1, len(files))
> print
>
> if not os.path.exists(basename[0:9]):
> print "Step 1 of 8: converting file to binary format..."
> outfile = basename + ".bin"
> os.system('python %s/seaconvert.py -bnxzad -o %s %s' % (SCRIPT_DIR, outfile, file))
>
> print "Step 2 of 8: trimming beams..."
> infile = outfile
> outfile = basename + ".bin.flt"
> os.system('python -u %s/trim.py -d%f -n%f %s > %s' % (SCRIPT_DIR, MAX_SWATH_MULTIPLE, MIN_SWATH_MULTIPLE, infile, outfile))
> os.remove(infile)
>
> print "Step 3 of 8: removing roll artifacts..."
> infile = outfile
> outfile = basename + ".bin.flt.roll"
> os.system('python -u %s/deroll.py %s %d > %s' % (SCRIPT_DIR, infile, ROLL_BEAMS, outfile))
> os.remove(infile)
>
> print "Step 4 of 8: removing system noise..."
> infile = outfile
> outfile = basename + ".bin.flt.roll.noise"
> os.system('python -u %s/denoise.py %s %d > %s' % (SCRIPT_DIR, infile, NOISE_BEAMS, outfile))
> os.remove(infile)
>
> print "Step 5 of 8: calculate spatial extent of line..."
> infile = outfile
> outfile = basename + ".region"
> os.system('gmtset D_FORMAT %.2f')
> os.system('gmtconvert %s -bi6 -F1,2 -bo > temp1' % infile)
> os.system('minmax -bi2 -I%d temp1 > %s' % (GRID_RESOLUTION, outfile))
> os.remove('temp1')
>
> print "Step 6 of 8: gridding line..."
> Rfile = basename + ".region"
> infile = basename + ".bin.flt.roll.noise"
> outfile = basename + ".nc"
> region = open(Rfile).readline().strip()
> os.system('gmtconvert %s -bi6 -F1,2,3 -bo > temp1' % (infile))
> os.system('blockmean -I%d %s -C -bi3 -bo temp1 > temp2' % (GRID_RESOLUTION, region))
> os.remove('temp1')
> os.system('nearneighbor -G%s -I%d -bi3 -N4 -R -S%d temp2' % (outfile, GRID_RESOLUTION, SEARCH_RADIUS))
> os.remove('temp2')
> os.remove(infile)
> os.remove(Rfile)
>
> print "Step 7 of 8: convert line to ESRI ASCII Grid format..."
> infile = basename + ".nc"
> outfile = GRID_DIR + "/" + basename + ".asc"
> os.system('grdmath %s 100 MUL RINT = temp.nc' % (infile))
> os.system('grd2xyz temp.nc -E > %s' % (outfile))
> os.remove('temp.nc')
> os.remove(infile)
>
> print "Step 8 of 8: convert line to ESRI Raster format..."
> infile = GRID_DIR + "/" + basename + ".asc"
> outfile = GRID_DIR + "/" + basename
> #os.system('ascii2arc.py %s %s' % (infile, outfile))
> os.system("ascii_to_arcraster_project_pyramids.py %s %s" % (infile, outfile))
> os.remove(infile)
>
> # Calculate Estimated time remaining
> time_stop = time.time()
> if time_stop == time_start:
> time_stop = time_start + 1
> bytes_file = os.path.getsize(file)
> bytes_sofar = bytes_sofar + bytes_file
> bytes_permin = float(bytes_file) / (time_stop - time_start) * 60.0
>
> bytes_remaining = bytes_total - bytes_sofar
> time_remaining = bytes_remaining / bytes_permin
>
> percent_complete = float(bytes_sofar) / float(bytes_total) * 100
>
> # End this line
> print
> print "Line Completed. Completed: %.0f%%, Est. Time Remaining: %.0f minutes" % (percent_complete, time_remaining)
> print "-------------------------------------------------------------------------"
> print
>
20060722
>Ran Python script, "ascii2arc", designed to convert ASCII grids to ESRI ArcGrid fromat
>
> ****ascii2arc Script)****
> # ---------------------------------------------------------------------------
> # ascii_to_arcraster_project_pyramids.py
> # Created on: Mon Oct 02 2006 01:45:01 PM
> # Nadine Golden
> # Usage: ascii2arc.py <temp1_txt> <Output_raster>
> # ---------------------------------------------------------------------------
>
> # Import system modules
> import glob, sys, string, os, win32com.client
>
> # Create the Geoprocessor object
> gp = win32com.client.Dispatch("esriGeoprocessing.GpDispatch.1")
>
> # Load required toolboxes...
> gp.SetProduct("ArcInfo")
> gp.AddToolbox("E:/Mapping/ArcGIS/ArcToolbox/Toolboxes/Data Management Tools.tbx")
> gp.AddToolbox("E:/Mapping/ArcGIS/ArcToolbox/Toolboxes/Conversion Tools.tbx")
>
> # Get script arguments...
> if len(sys.argv) < 3:
> print >> sys.stderr, "ascii_to_arcraster_project_pyramids.py: error1: not enough arguments"
> sys.exit(1)
>
> # Convert ascii grid to arc grid
> ingrid = sys.argv[1]
> outgrid = sys.argv[2]
> (tail, head) = os.path.split(ingrid)
> (root, ext) = os.path.splitext(head)
>
> if len(root) > 13:
> print >> sys.stderr, "ascii2arc.py: error 1: Grid name exceeds 13 characters: %s" % root
> sys.exit(1)
>
> gp.Workspace = tail
> gp.ScratchWorkspace = tail
> #print "Workspace:", gp.Workspace
>
> tmpgrid = tail + "/temp"
>
> # Check to see if temp already exists and delete it if necessary
> if os.path.exists(tmpgrid):
> gp.Delete(tmpgrid)
> #print "ingrid1: %s outgrid1: %s" % (ingrid, tmpgrid)
> gp.ASCIIToRaster_conversion(ingrid, tmpgrid, "INTEGER")
>
> # Process: Project Raster...
> ingrid = tmpgrid
> outgrid = tail + "/" + root
> if os.path.exists(outgrid):
> gp.Delete(outgrid)
> gp.ProjectRaster_management(ingrid, outgrid,
> "PROJCS['WGS_1984_UTM_Zone_10N',
> GEOGCS['GCS_WGS_1984',
> DATUM['D_WGS_1984',
> SPHEROID['WGS_1984',6378137.0,298.257223563]],
> PRIMEM['Greenwich',0.0],
> UNIT['Degree',0.0174532925199433]],
> PROJECTION['Transverse_Mercator'],
> PARAMETER['False_Easting',500000.0],
> PARAMETER['False_Northing',0.0],
> PARAMETER['Central_Meridian',-123.0],
> PARAMETER['Scale_Factor',0.9996],
> PARAMETER['Latitude_Of_Origin',0.0],
> UNIT['Meter',1.0]];
> -10000 -10000 100000;0 100000;0 100000",
> "NEAREST", "1")
>
> # Process: Build Pyramids...
> gp.BuildPyramids_management(outgrid)
>
> # Clean up the mess
> if os.path.exists(tmpgrid):
> gp.Delete(tmpgrid)
>
20060722
In ArcGIS, imported lines output from "Bathy_Proc" according to track line number and mosaiced
adjacent bathymetry lines using the "blend" option on the ArcGIS Raster "Merge" Tool.
20060810
In ArcGIS, imported lines according to track line number and mosaiced adjacent bathymetry lines.
Mosaicing was done by manually drawing a mask around the best possible data for each line and
it's overlapping lines. Best possible data was determined subjectively and included the least
no data values as possible in the overlapping areas.
20060812
Projected chancen_bath file available at https://pubs.usgs.gov/of/2007/1271/ from UTM, zone 10N, WGS84 coordinates to UTM, zone 11N, WGS84 coordinates and added a value of 3.6 m to the bathymetry data to account for a vertical offset.
2010
Raster
Pixel
15638
7390
WGS_1984_UTM_Zone_11N
-117
0.999600
500000.000
0.000
coordinate pair
U.S. Geological Survey, Coastal and Marine Geology
Guy R. Cochrane
Geophysicist
mailing and physical address
USGS Pacific Science Center, 400 Natural Bridges Drive
Santa Cruz
CA
95060-5792
USA
(831) 427-4754
(831) 427-4748
gcochrane@usgs.gov
Although this data set has been used by the USGS,
no warranty, expressed or implied, is made by the USGS
as to the accuracy of the data and/or related materials.
The act of distribution shall not constitute any such warranty,
and no responsibility is assumed by the USGS
in the use of these data or related materials.
Any use of trade, product, or firm names is
for descriptive purposes only and does not imply
endorsement by the U.S. Government.
2010
U.S. Geological Survey, Coastal and Marine Geology
Peter Dartnell
Physical Scientist
mailing and physical address
USGS 345 Middlefield Rd
Menlo Park
CA
94025
USA
(650)329-5460
pdartnell@usgs.gov
FGDC Content Standards for Digital Geospatial Metadata
FGDC-STD-001-1998