gshhsmapdata-polstere: GSHHS Low Resolution Coastlines in Polar Stereographic map...

In steingod/R-mipolsat: Useful functions for handling locally processed satellite data

Description

GSHHS - A Global Self-consistent, Hierarchical, High-resolution Shoreline Database. GSHHS is a high-resolution shoreline data set amalgamated from two data bases in the public domain. The data have undergone extensive processing and are free of internal inconsistencies such as erratic points and crossing segments. The shorelines are constructed entirely from hierarchically arranged closed polygons. The data can be used to simplify data searches and data selections, or to study the statistical characteristics of shorelines and land-masses. It comes with access software and routines to facilitate decimation based on a standard line-reduction algorithm.

The dataset provided herein have been extracted using the software provided by the National Geophysical Data Center of NOAA Satellite and Information Service. The original dataset contains lakes, islands within lakes etc. This dataset contain only the coastlines.

The dataset has been implemented in R using the getRgshhsMap function of maptools.

This dataset is a Proj.4 modified version of gshhsmapdata to a Polar Stereographic map projection correct at 60 degrees North.

Usage

    data(gshhsmapdata-polstere)

Format

A data.frame containing two elements, latitude (lat) and longitude (lon) in decimal degrees.

Source

http://www.ngdc.noaa.gov/mgg/shorelines/gshhs.html

GSHHS is developed and maintained by:

Dr. Paul Wessel, SOEST, University of Hawai'i, Honlulu, HI. [e-mail], and

Dr. Walter H. F. Smith, NOAA Laboratory for Satellite Altimetry, National Oceanographic Data Center, Silver Spring, MD. [e-mail]

References

Wessel, P., and W. H. F. Smith, A Global Self-consistent, Hierarchical, High-resolution Shoreline Database, J. Geophys. Res., 101, #B4, pp. 8741-8743, 1996.

Examples

    # Plot the ordinary map directly.
    plot(gshhsmapdata$lon,gshhsmapdata$lat,type="l")

    # To create a map projection use package proj4
    library(proj4)
    tmp <- project(list(gshhsmapdata$lon-180,gshhsmapdata$lat),"+proj=merc")
    plot(tmp$x,tmp$y,type="l")

steingod/R-mipolsat documentation built on May 30, 2019, 2:31 p.m.