sl.plot.points <-
function (plot.init.res,lon,lat,col="black",pch=1,cex=1,ignore.visibility=FALSE,return.xy=FALSE) {
projection = plot.init.res$projection
if (projection == "platon" || projection == "3D") {
if (return.xy) {warning("'return.xy=TRUE' not used for multi-part projections")}
if (!is.null(plot.init.res$transform.function)) {
lonlat.trans = plot.init.res$transform.function(lon,lat)
lon = lonlat.trans$lon
lat = lonlat.trans$lat
}
npir = 1
repeat {
pir = plot.init.res[[npir]]
if (!is.list(pir)) {return()}
if (is.null(pir$projection)) {return()}
sl.plot.points(pir,lon,lat,col,pch,cex,ignore.visibility)
npir = npir + 1
}
}
if (!is.null(plot.init.res$transform.function)) {
lonlat.trans = plot.init.res$transform.function(lon,lat)
lon = lonlat.trans$lon
lat = lonlat.trans$lat
}
vsr.res = sl.vis.shift.rot(plot.init.res,lon,lat)
x = vsr.res$x
y = vsr.res$y
rot.lon = vsr.res$rot.lon
rot.lat = vsr.res$rot.lat
visible = vsr.res$visible
if (ignore.visibility) {
visible[] = TRUE
}
xshift = plot.init.res$xshift
yshift = plot.init.res$yshift
np = length(lon)
if (sum(visible) < np) {
if (length(col) > 1) {
if (length(col) < np) {
col = rep(col,np%/%length(col)+1)
}
col = col[(1:np)[visible]]
}
if (length(pch) > 1) {
if (length(pch) < np) {
pch = rep(pch,np%/%length(pch)+1)
}
pch = pch[(1:np)[visible]]
}
if (length(cex) > 1) {
if (length(cex) < np) {
cex = rep(cex,np%/%length(cex)+1)
}
cex = cex[(1:np)[visible]]
}
x = x[visible]
y = y[visible]
}
if (projection == "mollweide") {
xy = sl.proj.mollweide(lon = x, lat = y)
x = xy$x
y = xy$y
}
points(x=x+xshift,y=y+yshift,col=col,pch=pch,cex=cex)
if (return.xy) {return(list(x=x+xshift,y=y+yshift))}
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.