fixCoastwrap: Correct the Wrapping problem

In GEOmap: Topographic and Geologic Mapping

Correct the Wrapping problem

Description

Correct wrapping for GEOmaps

Usage

fixCoastwrap(Z, maxdis = 100)

Arguments

Z	list of x, y
maxdis	maximum distance for differences

Details

Based on mapswrap program

Value

List:

x	x-coordinates (longitudes)
y	y-coordinates (latitudes)

Author(s)

Jonathan M. Lees<[e-mail]>

Examples

data(coastmap)
SEL = which(coastmap$STROKES$nam=="AFROASIA")

A = ExcludeGEOmap(coastmap, SEL, INOUT="in")

plot(A$POINTS$lon, A$POINTS$lat, type='n')

points(A$POINTS$lon, A$POINTS$lat, pch='.')

######  note that the map wraps around.

B = fixCoastwrap(list(x=A$POINTS$lon, y=A$POINTS$lat), 100)
 which(is.na(B$x))



lines(B)


polygon(B, col=rgb(.8,1, .8))

GEOmap documentation built on Sept. 1, 2023, 5:09 p.m.