NOAA Electronic Navigational Chart [ENC] -- Bathymetry Grid

Metadata also available as

Metadata:


Identification_Information:
Citation:
Citation_Information:
Originator: U.S. EPA - Mid-Continent Ecology Division-Duluth
Publication_Date: Unpublished Material
Title: NOAA Electronic Navigational Chart [ENC] -- Bathymetry Grid
Geospatial_Data_Presentation_Form: raster digital data
Online_Linkage:
M:\base_data\water\bathymetry\noaa_enc\lake\superior\data\bathy_grids\
Description:
Abstract:
Bathymetric integer grid generated from interpolating original point sounding and depth line features acquired from NOAA Electronic Navigational Charts [ENCs]. ENC data is available for various geographic regions from the Office of Coast Survey [<http://chartmaker.ncd.noaa.gov/>]>. Lake Superior data was identified, downloaded, and the ENC data imported into geodatabase format.
The narrow coastal zone corresponding to the depth area outlined by depth contours was interpolated from DepthLineFeature contours. The rest of the lake bathymetry data came from SoundingFeature point interpolations.
RELATIONSHIP BETWEEN SCALE AND SPATIAL RESOLUTION
Vector to raster interpolation requires that certain assumptions are made regarding an appropriate raster cell size corresponding to a vector map scale. A literature review was performed in order to try and identify some sort of industry standard, logic, or methodology, which could be used to relate vector scale to raster spatial resolution. The literature review revealed a few common relationships between cell sizes and original map scales, but the conclusion to this literature review is there is no definite industry standard with defined rules. Some of the sources found along the way are as follows:
 <http://srmwww.gov.bc.ca/gis/gisscale.html>
 <http://cbc.rs-gis.amnh.org/guides/basic_concepts/scale_res.html>
 <http://www.profc.udec.cl/~gabriel/tutoriales/giswb/vol2/cp1/t1-2.gif>
 <http://listes.ulaval.ca/cgi-bin/wa?A2=ind0503&L=applied-gis-rs&F=&S=&P=3331>
The goal is to get the most information from the sounding features as possible, but also not to use an output cell size which assumes more accuracy than is actually present in the data, as the cell size should not be smaller than the uncertainty in the vector linework. In order to determine what output cell size was most appropriate for the production of bathymetry grids, the vector map scale [i.e., CSCL or compilation scale given as a subfield of the Data Set Parameter [DSPM] field] of each NOAA ENC would be identified. And because with a vector data model resolution is tied to the scale of the data, the cell size of the vector data set can be calculated using this scale. Sources say that the smallest line that can be drawn on a map which is still recognizable is about .5 millimeters. On a 1:24,000 scale map, 1 millimeter equals 24,000 millimeters on the ground. Since you can only draw something as small as .5 mm, the resolution of the data set is 12,000 mm or 12 meters. Using this logic, the output cell size of the bathymetry grids was selected.
This methodology was proposed via email to Adeline Wong [Adeline.Wong@noaa.gov] with the NOAA National Ocean Service [NOS] IN Silver Spring, MD. Ms. Wong provided the following statement regarding my proposed methodology, "Your approach is sensible with regard to ENC soundings, which are often downsampled and scaled by chart." Ms. Wong was also asked about how to determine the accuracy of the ENC data. During the same literature review, it was discovered in section 7.3.4, Data Set Accuracy Record Structure, of "IHO Transfer Standard for Digital Hydrographic Data, Edition 3.1 - November 2000" [i.e., 31Main.pdf], there is mention of a field table named SDAC or Data Set Accuracy. However, such a table has never been found with any of the ENCs downloaded at MED. It appears as though the information in this table would answer a lot of questions regarding ENC accuracy. Ms. Wong responded by saying, "Your question on the accuracy of charted features is on the cutting edge of research. We're currently developing methods to attribute accuracy (SDAC field) for ENCs!"
The compilation scale of all 19 NOAA ENC charts available for Lake Superior ranged from 1:15,000 to 1:600,000, or 8 meters to 300 meters.
Lake NOAA Chart NOAA Cell CSCL Scale Resolution [m] Superior 14965 US4MI70M 1:120,000 60 Superior 14963 US4MI76M 1:120,000 60 Superior 14962 US4MI77M 1:120,000 60 Superior 14966 US4MN11M 1:120,000 60 Superior 14967 US4MN21M 1:120,000 60 Superior 14968 US4MN22M 1:120,000 60 Superior 14970 US5MI75M 1:15,000 8 Superior 14975 US5WI11M 1:15,000 8 Superior 14974 US5WI22M 1:15,000 8 Superior 14964 US4MI71M 1:180,000 90 Superior 14971 US5MI72M 1:30,000 15 Superior 14972 US5MI73M 1:30,000 15 Superior 14969 US5MI74M 1:30,000 15 Superior 14882 US5MI61M 1:40,000 20 Superior 14883 US5MI62M 1:40,000 20 Superior 14884 US5MI63M 1:40,000 20 Superior 14976 US5MI78M 1:40,000 20 Superior 14973 US4WI21M 1:60,000 30 Superior 14961 US3MI79M 1:600,000 300
Purpose:
Use bathymetry grids to calculate bathymetric zonal statistics for GLEI embayment sites.
Supplemental_Information:
Metadata produced by CSC under the Fair II contract 68-W-02-032 Task Order 2024.
Time_Period_of_Content:
Time_Period_Information:
Range_of_Dates/Times:
Beginning_Date: 20010711
Ending_Date: Present
Currentness_Reference: publication date
Status:
Progress: Complete
Maintenance_and_Update_Frequency: As needed
Spatial_Domain:
Bounding_Coordinates:
West_Bounding_Coordinate: -92.295368
East_Bounding_Coordinate: -91.978650
North_Bounding_Coordinate: 46.799942
South_Bounding_Coordinate: 46.637911
Keywords:
Theme:
Theme_Keyword_Thesaurus: MED Theme Thesaurus
Theme_Keyword: Depth
Theme_Keyword: Bathymetry
Theme_Keyword: Sounding Features
Theme_Keyword: NOAA Electronic Navigational Chart
Place:
Place_Keyword_Thesaurus: MED Place Thesaurus
Place_Keyword: Lake Superior
Access_Constraints: None.
Use_Constraints:
None. Acknowledgement of the originating agencies would be appreciated in products derived from these data.
Point_of_Contact:
Contact_Information:
Contact_Organization_Primary:
Contact_Organization: U.S. EPA - Mid-Continent Ecology Division-Duluth
Contact_Person: Jack Kelly
Contact_Address:
Address_Type: mailing and physical address
Address: 6201 Congdon Blvd
City: Duluth
State_or_Province: MN
Postal_Code: 55804
Country: USA
Contact_Voice_Telephone: 218-529-5000
Contact_Facsimile_Telephone: 218-529-5003
Security_Information:
Security_Classification_System: None.
Security_Classification: Unclassified
Native_Data_Set_Environment:
Microsoft Windows XP Version 5.1 (Build 2600) Service Pack 1; ESRI ArcCatalog 9.0.0.535

Data_Quality_Information:
Attribute_Accuracy:
Attribute_Accuracy_Report:
The NOAA ENC® Direct to GIS data is derived from NOAA ENC®'s which conforms to the International Hydrographic Organization Transfer Standard for Digital Hydrographic Data, Special Publication S-57, ENC Product Specification Edition 2.0.
Logical_Consistency_Report:
The NOAA ENC® Direct data is derived from NOAA ENC®'s which conforms to the International Hydrographic Organization Transfer Standard for Digital Hydrographic Data, Special Publication S-57, ENC Product Specification Edition 2.0. NOAA ENC® Direct to GIS data has added generalized nautical chart symbology to aid in the visualazition of this data.
Completeness_Report:
This dataset is complete.
NOAA ENC® Direct to GIS data represents current coverage of all large scale Harbor and Approach (1:150,000 to 1:50,000 or larger) ENC Charts. The NOAA ENC® Chart Suite is not scheduled to be completed until 2007. As new NOAA ENC®s and NOAA ENC editions are created they will be added to the NOAA ENC Direct to GIS data on the first of every month.
Positional_Accuracy:
Horizontal_Positional_Accuracy:
Horizontal_Positional_Accuracy_Report:
The NOAA ENC® Direct to GIS data is derived from data that conforms to the International Hydrographic Organization Transfer Standard for Digital Hydrographic Data, Special Publication S-57, ENC Product Specification Edition 2.0.
Vertical_Positional_Accuracy:
Vertical_Positional_Accuracy_Report:
The NOAA ENC® Direct to GIS data is dereived from data that conforms to the International Hydrographic Organization Transfer Standard for Digital Hydrographic Data, Special Publication S-57, ENC Product Specification Edition 2.0.
Lineage:
Process_Step:
Process_Description:
METHODOLOGY
Technical Procedures
This section describes the pre-processing, processing, and post-processing steps used to process NOAA ENC bathymetry grids and to calculate embayment bathymetric statistics. Steps 1-9 in the technical process description below have to be implemented separately for each individual ENC database. Steps 10-13 generate the final bathymetry statistics for all embayments of interest.
1. In ArcMap, set the Data Frame coordinate system to the USA Contiguous Albers Equal Area Conic USGS projection. The Data Frame Properties can be accessed from the View menu in ArcMap and the coordinate system specified from the Coordinate System tab.
2. Add the DepthAreaFeature, DepthLineFeature, and SoundingFeature data sets from the identified NOAA ENC Geodatabases [GDB] into ArcMap.
3. Once added into ArcMap, export each of the three GDB features classes into shapefiles by right-clicking on the individual feature, selecting Data, and then selecting Export Data. Make sure to check the radio button next to "Use the same Coordinate System as the data frame" so that all shapefiles will be exported in the Albers projection defined above.
4. In ArcMap, activate the Spatial Analyst Extension from the Tools menu. Then, select Options from under the Spatial Analyst toolbar dropdown list. From the General tab, check the radio button next to "Analysis output will be saved in the same coordinate system as the active data frame," from the Extent tab, select the shapefile representing the DepthAreaFeature as the Analysis Extent, and from the Cell Size tab, specify your desired output cell size.
5. Next, create two mask grids from the DepthAreaFeature polygon shapefile by using:
a]. all of depth area polygons contained within the shapefile, and
b]. by selecting the largest inside polygon[s] and then switching the selection to the remaining coastal polygons only.
To create the two mask grids, make your a]. and b]. selection of features, select Convert from the Spatial Analyst toolbar dropdown list, and finally Features to Raster. Remember to select the DepthAreaFeature polygon shapefile as the input features, DEPTH_TYPE as the field, your chosen output cell size, and a desired location and name for your output raster. Name the two masks appropriately so the a]. ENTIRE and the b]. COASTAL depths area masks are easy to distinguish from one another.
6. Once you have created your masks, access and use the Topo to Raster interpolation tool from the Spatial Analysis Toolbox from ArcToolbox and interpolate the DepthLineFeature, and SoundingFeature shapefiles [created in step 3] into two different floating-point bathymetry grids:
a]. Bathymetry Grid 1 - For the sounding features shapefile, specify DEPTH for the Field value and PointElevation for the Type value. Specify an output surface raster name [i.e., SOUNDING90] and cell size, select NO_ENFORCE for drainage enforcement and SPOT for the primary type of input data.
b]. Bathymetry Grid 2 - For the depth line features shapefile, specify VALDCO [value of depth contour] for the Field value and Contour for the Type value. Specify an output surface raster name [i.e., CONTOUR90] and cell size, select NO_ENFORCE for drainage enforcement and CONTOUR for the primary type of input data.
7. After the soundings and contour interpolations are complete, select Options from the Spatial Analyst toolbar dropdown list. From the General tab, specify the ENTIRE mask created in step 5a as the Analysis Mask. Next, select Raster Calculator from the Spatial Analyst toolbar dropdown list and convert the floating-point bathymetry grid of sounding features created in step 6a to an integer grid with vertical units in centimeters using the following expression: BATHY_SOUND = INT [100 * SOUNDING90]. Repeat the same steps to create an integer grid of the coastal contours using the COASTAL mask created in step 5b as the Analysis Mask and the following expression: BATHY_CONTOUR = INT [100 * CONTOUR90].
8. Now, merge the BATHY_ SOUND and the BATHY_CONTOUR grids together using the ENTIRE mask created in step 5a as the Analysis Mask, and the following expression in the Raster Calculator: ENC_BATHY = MERGE [BATHY_CONTOUR, BATHY_ SOUND]. Be sure to put the BATHY_CONTOUR grid first in the expression so that its values will take precedence over the values of BATHY_ SOUND. Note: "ENC" in each bathymetry grid name is replaced by the appropriate ENC dataset cell name [i.e. MN21, MI78, etc.].
The vertical units of each final output ENC_BATHY grid are in centimeters and represent bathymetry for the spatial extent of the NOAA ENC Geodatabase DepthLineFeature [contours] and SoundingFeature [soundings] data sets. Note: Using the raster command MOSAIC instead of MERGE will produce different results along the shorelines. Depth values in coastal areas become altered as the mosaic algorithm attempts to reconcile values within overlapping portions of coastal and deep water grids.
9. If desired, you can now delete all intermediate vector and raster datasets. If you would like to merge two or more ENC bathymetry grids together, please following the remaining steps below. If you do not wish to merge grids together, please skip step 10.
10. Select Options from the Spatial Analyst toolbar dropdown list. From the General tab, specify the Analysis Mask as <None> and from the Extent tab specify "Maximum of Inputs" for the Analysis Extent. Finally, merge all desired ENC bathymetry grids together using the following expression: MERGED_BATHY = MERGE [ENC_BATHY1, ENC_BATHY2, . ENC_BATHYn].
11. A limited number of cells with negative values could be found in bathymetry grids because of an interpolation result of parallel contours labeled by zero depth values. Therefore, it is necessary to convert all negative bathymetry grid values to a value of '0' so that future bathymetric calculations will not be affected. The following expression can be used in the Spatial Analyst Raster Calculator to accomplish this task: FINAL_BATHY = CON [MERGED_BATHY < 0, 0, MERGED_BATHY].
12. Once the FINAL_BATHY grid has been created, access and use the Zonal Statistics as Tables tool from the Spatial Analysis Toolbox from ArcToolbox to calculate zonal statistics for the embayments of interest. If there is more than one embayment coverage for which zonal statistics are needed, this step will need to be repeated for each additional embayment coverage that exits. Scripts calc_embay_grid.py and clac_embay_shape in M:\med_projects\glgr\wo500227\output\scripts\ directory were used to automate this step. The Input feature zone data will be the polygon coverage or shapefile containing embayment geometry, the Zone Field should be a numeric field that has unique values for each embayment and can later be used to join the zonal statistics results back to the polygon features, the Input Value Raster should be the FINAL_BATHY grid, and lastly just specify and output location and DBF table name.
13. Open the output DBF table and add a new field which will hold the water volume calculation values. Give the field a descriptive name and specify the Type as DOUBLE accepting the default values for precision and scale. Provided that depth values/vertical units are in centimeters and each raster cell area is 90 x 90 = 8,100 square meters, calculate the water volume in cubic meters for each embayment using the mean depth and cell count as: [[MEAN] / 100] * [[COUNT] * 8100]. An alternative expression for calculating water volume is: [[SUM] / 100] * 8100. These expressions should be modified by using a multiplier of 900 instead of 8,100 if 30 meter grids are used for water volume calculations.
Bathymetry grids have centimeters as vertical units and represents bathymetry for an area covered by ENC sounding and contour data.
End of technical procedures.
Process_Date: May 2005
Process_Time: unknown
Process_Contact:
Contact_Information:
Contact_Organization_Primary:
Contact_Organization: EPA - Mid-Continent Ecology Division-Duluth
Contact_Person: Jack Kelly
Contact_Address:
Address_Type: mailing and physical address
Address: 6201 Congdon Blvd.
City: Duluth
State_or_Province: MN
Postal_Code: 55804
Country: USA
Contact_Voice_Telephone: 218-529-5000
Contact_Facsimile_Telephone: 218-529-5003
Process_Step:
Process_Description:
The NOAA ENC® database has been built from a combination of charted information as well as original "source" information. NOAA has compiled critical features such as channel limits, aids to navigation, and obstructions from the original documents that were used to put the feature on the paper chart. The objective is to use the most accurate information for features that are critical to the safety of navigation. NOAA uses a number of sources in compiling NOAA ENC®S including U.S. Army Corps of Engineers surveys, drawings, and permits, U.S. Coast Guard Local Notices to Mariner, National Imagery and Mapping Agency Notices to Mariners, NOAA hydrographic surveys, and the largest scale paper chart of an area. ENC® Direct to GIS data was created by transforming Approach and Harbor NOAA ENC® version 2 cells to ESRI's shapefile format using Safe Software's Feature Manipulation Engine.
Process_Date: 200305 - present
Process_Contact:
Contact_Information:
Contact_Organization_Primary:
Contact_Organization:
National Oceanic and Atmospheric Administration (NOAA), National Ocean Service (NOS), Office of Coast Survey (OCS)
Contact_Position: NOAA ENC® Direct to GIS Technical Contact
Contact_Address:
Address_Type: mailing and physical address
Address: 1315 East West Highway N/CS2
City: Silver Spring
State_or_Province: Maryland
Postal_Code: 20910
Country: USA
Contact_Voice_Telephone: 301-713-2645
Contact_Facsimile_Telephone: 301-713-4501
Contact_Electronic_Mail_Address: enc.gis@noaa.gov
Hours_of_Service: 0700-1500 EST
Contact_Instructions: E-mail is preferred
Process_Step:
Process_Description: Dataset copied.
Source_Used_Citation_Abbreviation:
E:\s_drive\GLEI-GIS\Base\Bathymetry\Lake Superior - EPA\bathymetry\noaa\enc\lake\superior\data\bathy_grids\us5wi11_bathy
Process_Date: 20090803
Process_Time: 15444800

Spatial_Data_Organization_Information:
Direct_Spatial_Reference_Method: Raster
Raster_Object_Information:
Raster_Object_Type: Grid Cell
Row_Count: 2125
Column_Count: 2954
Vertical_Count: 1

Spatial_Reference_Information:
Horizontal_Coordinate_System_Definition:
Planar:
Map_Projection:
Map_Projection_Name: Albers Conical Equal Area
Albers_Conical_Equal_Area:
Standard_Parallel: 29.500000
Standard_Parallel: 45.500000
Longitude_of_Central_Meridian: -96.000000
Latitude_of_Projection_Origin: 23.000000
False_Easting: 0.000000
False_Northing: 0.000000
Planar_Coordinate_Information:
Planar_Coordinate_Encoding_Method: coordinate pair
Coordinate_Representation:
Abscissa_Resolution: 0.000005
Ordinate_Resolution: 0.000005
Planar_Distance_Units: meters
Geodetic_Model:
Horizontal_Datum_Name: North American Datum of 1983
Ellipsoid_Name: Geodetic Reference System 80
Semi-major_Axis: 6378137.000000
Denominator_of_Flattening_Ratio: 298.257222

Entity_and_Attribute_Information:
Detailed_Description:
Entity_Type:
Entity_Type_Label: NOAA ENC Bathymetry Grid
Attribute:
Attribute_Label: Rowid
Attribute_Definition: Internal feature number.
Attribute_Definition_Source: ESRI
Attribute_Domain_Values:
Unrepresentable_Domain:
Sequential unique whole numbers that are automatically generated.
Attribute:
Attribute_Label: Value
Attribute_Definition: Bathymetry depth values. Vertical units are in centimeters.
Attribute_Definition_Source: USEPA-MED
Attribute:
Attribute_Label: Count

Distribution_Information:
Distributor:
Contact_Information:
Contact_Organization_Primary:
Contact_Organization: U.S. EPA - Mid-Continent Ecology Division-Duluth
Contact_Person: Jack Kelly
Contact_Address:
Address_Type: mailing and physical address
Address: 6201 Congdon Blvd
City: Duluth
State_or_Province: MN
Postal_Code: 55804
Country: USA
Contact_Voice_Telephone: 218-529-5000
Contact_Facsimile_Telephone: 218-529-5003
Resource_Description: Downloadable Data
Distribution_Liability:
Although these data have been processed successfully on a computer system, no warranty expressed or implied is made by the EPA regarding the utility of the data on any other system, nor shall the act of distribution constitute any such warranty.
Standard_Order_Process:
Digital_Form:
Digital_Transfer_Information:
Transfer_Size: 0.031

Metadata_Reference_Information:
Metadata_Date: 20050513
Metadata_Contact:
Contact_Information:
Contact_Organization_Primary:
Contact_Organization: U.S. EPA - Mid-Continent Ecology Division-Duluth
Contact_Person: Jack Kelly
Contact_Address:
Address_Type: mailing and physical address
Address: 6201 Congdon Blvd
City: Duluth
State_or_Province: MN
Postal_Code: 55804
Country: USA
Contact_Voice_Telephone: 218-529-5000
Contact_Facsimile_Telephone: 218-529-5003
Metadata_Standard_Name: FGDC Content Standards for Digital Geospatial Metadata
Metadata_Standard_Version: FGDC-STD-001-1998
Metadata_Time_Convention: local time
Metadata_Access_Constraints: None.
Metadata_Use_Constraints:
None. Acknowledgement of the originating agencies would be appreciated in products derived from these data.
Metadata_Security_Information:
Metadata_Security_Classification_System: None.
Metadata_Security_Classification: Unclassified
Metadata_Extensions:
Online_Linkage: <http://www.esri.com/metadata/esriprof80.html>
Profile_Name: ESRI Metadata Profile

Generated by mp version 2.9.6 on Wed Aug 05 14:03:49 2009