![]() ![]() |
Although this Federal Geographic Data Committee-compliant metadatafile is intended to document the data set in nonproprietary form,as well as in ArcInfo format, this metadata file may include someArcInfo-specific terminology.
Golden, Nadine E. , and Cochrane, Guy R. , 2007, chanwest_back.txt - Sidescan sonar backscatter image for Channel West in the Nearshore Benthic Habitat Mapping Project S. California map series.: Open-File Report 2007-1271, U.S. Geological Survey, Santa Cruz, CA.This is part of the following larger work.Online Links:
Cochrane, Guy R. , Golden, Nadine, Dartnell, Pete, Schroeder, Donna, Conrad, Jamie, and Stevenson, Andy, 2007, Seafloor Mapping and Benthic Habitat GIS for Southern California, Volume III.
This is a Raster data set.It contains the following raster data types:
(831) 427-4754 (voice)
(831) 427-4748 (FAX)
gcochrane@usgs.gov
These data are intended for science researchers, students, policymakers, and the general public. The data can beused with geographic information systems (GIS) softwareto display geologic and oceanographic information.
**** Ascii to ArcRaster & Project Pyramids (Python Script)**** # --------------------------------------------------------------------------- # ascii_to_arcraster_project_pyramids_stand_alone.py # Created on: Tuesday Oct 04 2006 # Nadine Golden # --------------------------------------------------------------------------- # 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.AddToolbox("E:/ArcGIS/ArcToolbox/Toolboxes/Data Management Tools.tbx") gp.AddToolbox("E:/ArcGIS/ArcToolbox/Toolboxes/Conversion Tools.tbx") files = glob.glob("F:/test/grids/*asc") for file in files: (basename, ext) = os.path.splitext(file) ingrid = file outgrid = basename + "tmp" os.system("scii_to_arcraster_project_pyramids_stand_alone.py %s %s" % (ingrid, outgrid)) # Process: ASCII to Raster... gp.ASCIIToRaster_conversion(ingrid, outgrid, "INTEGER") ## # Process: Define Projection... ## gp.DefineProjection_management(outgrid, "PROJCS['NAD_1983_UTM_Zone_11N',GEOGCS['GCS_North_American_1983',DATUM['D_North_American_1983', SPHEROID['GRS_1980',6378137.0,298.257222101]], PRIMEM['Greenwich',0.0], UNIT['Degree',0.0174532925199433]], PROJECTION['Transverse_Mercator'], PARAMETER['False_Easting',500000.0], PARAMETER['False_Northing',0.0], PARAMETER['Central_Meridian',-117.0], PARAMETER['Scale_Factor',0.9996], PARAMETER['Latitude_Of_Origin',0.0], UNIT['Meter',1.0]]") # Process: Project Raster... finalgrid = basename gp.ProjectRaster_management(outgrid, finalgrid, "PROJCS['WGS_1984_UTM_Zone_11N',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',-117.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(finalgrid) # Clean up the mess gp.Delete(outgrid) **** Rasters Statistics for Normalizing Data (Python Script)**** # --------------------------------------------------------------------------- # RasterStats.py # To Calculate the mean value of a raster # Created on: Thursday Oct 20, 2006 # Nadine Golden # --------------------------------------------------------------------------- # Import system modules import sys, string, os, win32com.client # Create the Geoprocessor object gp = win32com.client.Dispatch("esriGeoprocessing.GpDispatch.1") # Check out any necessary licenses gp.CheckOutExtension("spatial") # Load required toolboxes... gp.AddToolbox("E:/ArcGIS/ArcToolbox/Toolboxes/Spatial Analyst Tools.tbx") gp.AddToolbox("E:/ArcGIS/RasterStatsTool/RasterStatistics.tbx") #Directory of grids gp.Workspace = "F:\\test" rasterlist = gp.ListRasters() name = rasterlist.Next() while name: print name # Process: RasterBand Statistics... mean = gp.RasterBandStatistics(name, "MEAN") print mean # Process: Minus... (indir,basename) = os.path.split(name) outgrid = os.path.join("F:\\test\\Normalized", basename) gp.Minus_sa(name, mean, outgrid) print outgrid name = rasterlist.Next()On the order of 10 meters.
Please see the Methods section of <http://pubs.usgs.gov/of/2007/1271/> forinformation about omissions, selection criteria, generalization,definitions used, and other rules used to derive the data set.
Logical Consistency untested.
Are there legal restrictions on access or use of the data?
- Access_Constraints: None
- Use_Constraints:
- Public domain data from the U.S. government is freely redistributablewith proper metadata and source attribution.Please recognize the U.S. Geological Survey (USGS)as the source of this information.
(831) 427-4754 (voice)
(831) 427-4748 (FAX)
gcochrane@usgs.gov
Although this data set has been used by the USGS,no warranty, expressed or implied, is made by the USGSas 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 USGSin the use of these data or related materials.Any use of trade, product, or firm names isfor descriptive purposes only and does not implyendorsement by the U.S. Government.
(831) 427-4730 (voice)
(831) 427-4748 (FAX)
ngolden@usgs.gov