chancen_bath.html

Metadata also available as

Metadata:


Identification_Information:
Citation:
Citation_Information:
Originator: Nadine E. Golden
Originator: Guy R. Cochrane
Publication_Date: 2007
Title:
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

Geospatial_Data_Presentation_Form:
Larger_Work_Citation Citation_Information Series_Information Series_Name Publication_Information Publication_Place Online_Linkage <https://pubs.usgs.gov/of/2007/1271/>

Description:
Abstract:
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.


Purpose:
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.

Supplemental_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.


Time_Period_of_Content:
Time_Period_Information:
Range_of_Dates/Times:
Beginning_Date: 20050808
Ending_Date: 20050827
Currentness_Reference: ground condition


Status:
Progress: Complete
Maintenance_and_Update_Frequency: As needed



Spatial_Domain:
Bounding_Coordinates:
West_Bounding_Coordinate: -119.643628
East_Bounding_Coordinate: -119.471725
North_Bounding_Coordinate: 34.390598
South_Bounding_Coordinate: 34.320467


Keywords:
Theme:
Theme_Keyword_Thesaurus: None
Theme_Keyword: bathymetry
Theme_Keyword: sidescan sonar
Theme_Keyword: remote-sensing
Theme_Keyword: MIPS
Theme_Keyword: U.S. Geological Survey
Theme_Keyword: USGS
Theme_Keyword: Coastal and Marine Geology Program
Theme_Keyword: CMGP
Theme:
Theme_Keyword_Thesaurus: ISO 19115 Topic Category
Theme_Keyword: geoscientificinformation
Theme_Keyword: imagereyBaseMapsEarthCover
Theme_Keyword: oceans
Theme_Keyword: elevation
Place:
Place_Keyword_Thesaurus: None
Place_Keyword: Channel Islands National Marine Sanctuary
Place:
Place_Keyword_Thesaurus: Geographic Names Information System (GNIS)
Place_Keyword: Santa Barbara
Place_Keyword: Ventura
Place_Keyword: Santa Barbara Channel
Place_Keyword: Channel Islands
Place_Keyword: California



Access_Constraints: None

Use_Constraints:
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.

Point_of_Contact:
Contact_Information:
Contact_Person_Primary:
Contact_Person: Guy R. Cochrane
Contact_Organization: U.S. Geological Survey, Coastal and Marine Geology
Contact_Position: Geophysicist
Contact_Address:
Address_Type: mailing and physical address
Address: USGS Pacific Science Center, 400 Natural Bridges Drive
City: Santa Cruz
State_or_Province: CA
Postal_Code: 95060-5792
Country: USA
Contact_Voice_Telephone: (831) 427-4754
Contact_Facsimile_Telephone: (831) 427-4748
Contact_Electronic_Mail_Address: gcochrane@usgs.gov


Browse_Graphic:
Browse_Graphic_File_Name: images/chancen_bath.jpg
Browse_Graphic_File_Description:
Image showing coverage of Channel Central bathymetry data ArcInfo format.
Browse_Graphic_File_Type: JPEG



Native_Data_Set_Environment:
SunOS, 5.7, sun4u UNIX Windows XP Professional, v. 2002, service pack 2 SEA Grid Processor v. 2.051 ArcGIS version 9.1

Cross_Reference:
Citation_Information:
Originator: U.S. Geological Survey, Coastal and Marine Geology Program
Publication_Date: 2007
Title: InfoBank
Online_Linkage:
<http://walrus.wr.usgs.gov/infobank/s/s105sc/html/s-1-05-sc.meta.html>

Cross_Reference:
Citation_Information:
Originator:
Guy Cochrane, U.S. Geological Survey, Coastal and Marine Geology Program
Publication_Date: 200510
Title:
Sonar Survey of Sea-Floor Habitats Southeast of Santa Barbara, California in Sound Waves Monthly Newsletter, Coastal Science & Research News from Across the USGS
Online_Linkage:
<http://soundwaves.usgs.gov/2005/10/fieldwork2.html>


Data_Quality_Information:

Logical_Consistency_Report: Logical Consistency untested.

Completeness_Report:
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.
Positional_Accuracy:
Horizontal_Positional_Accuracy:
Horizontal_Positional_Accuracy_Report: On the order of 2 meters.

Lineage:
Process_Step:
Process_Description:
Started with text files that were the output from A 234 kHz Interferometric Submetix Swath Bathy Sonar System.
Process_Date: 20060702
Process_Step:
Process_Description:
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

Process_Date: 20060722
Process_Step:
Process_Description:
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)

Process_Date: 20060722
Process_Step:
Process_Description:
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.
Process_Date: 20060810
Process_Step:
Process_Description:
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.
Process_Date: 20060812

Process_Step:
Process_Description:
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.
Process_Date: 2010



Spatial_Data_Organization_Information:
Direct_Spatial_Reference_Method: Raster
Raster_Object_Information:
Raster_Object_Type: Pixel
Row_Count: 15638
Column_Count: 7390


Spatial_Reference_Information:
Horizontal_Coordinate_System_Definition:
Planar:
Map_Projection:
Map_Projection_Name: WGS_1984_UTM_Zone_11N
Polar_Stereographic:
Straight-Vertical_Longitude_from_Pole: -117
Standard_Parallel: 0.999600
False_Easting: 500000.000
False_Northing: 0.000
Planar_Coordinate_Information:
Planar_Coordinate_Encoding_Method: coordinate pair



Distribution_Information:
Distributor:
Contact_Information:
Contact_Organization_Primary:
Contact_Organization: U.S. Geological Survey, Coastal and Marine Geology
Contact_Person: Guy R. Cochrane
Contact_Position: Geophysicist
Contact_Address:
Address_Type: mailing and physical address
Address: USGS Pacific Science Center, 400 Natural Bridges Drive
City: Santa Cruz
State_or_Province: CA
Postal_Code: 95060-5792
Country: USA
Contact_Voice_Telephone: (831) 427-4754
Contact_Facsimile_Telephone: (831) 427-4748
Contact_Electronic_Mail_Address: gcochrane@usgs.gov

Distribution_Liability:
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.

Metadata_Reference_Information:
Metadata_Date: 2010
Metadata_Contact:
Contact_Information:
Contact_Organization_Primary:
Contact_Organization: U.S. Geological Survey, Coastal and Marine Geology
Contact_Person: Peter Dartnell
Contact_Position: Physical Scientist
Contact_Address:
Address_Type: mailing and physical address
Address: USGS 345 Middlefield Rd
City: Menlo Park
State_or_Province: CA
Postal_Code: 94025
Country: USA
Contact_Voice_Telephone: (650)329-5460
Contact_Electronic_Mail_Address: pdartnell@usgs.gov
Metadata_Standard_Name: FGDC Content Standards for Digital Geospatial Metadata
Metadata_Standard_Version: FGDC-STD-001-1998

Generated by mp version 2.9.8 on Mon Jun 28 14:05:26 2010