library(RSEIS)
data(GH, package='RSEIS')
sta = GH$stafile
g1 = GH$pickfile
proj = GEOmap::setPROJ(type=2, LAT0 =median(sta$lat) , LON0 = median(sta$lon))
grcol = grey(seq(from=0.3, to=0.95, length=50))
contPFarrivals(g1, sta, proj=proj,cont=TRUE, POINTS=TRUE,
image=TRUE , col=grcol, phase="P",
add=FALSE )
Run the code above in your browser using DataLab