Nothing
targetLL <-
function(sta, rdist=100)
{
### get center for a projection
mlat = median(sta$lat)
mlon = median(sta$lon)
proj = setPROJ(2, LAT0 =mlat, LON0 =mlon )
print(c(mlat, mlon))
Jlat = range(sta$lat)
print(Jlat)
Jlon = range(sta$lon)
print(Jlon)
labs = c("LL", "LR", "UL", "UR" )
TARG = XY.GLOB(c(-rdist, rdist, -rdist, rdist ) , c(-rdist, -rdist, rdist, rdist ), proj)
A = cbind( TARG$lon, TARG$lat )
blon = TARG$lon
blon[TARG$lon>180] = blon-360
B = cbind( blon , TARG$lat )
rownames(A)<-labs
rownames(B)<-labs
write.table(A, quote = FALSE,col.names =FALSE )
cat("\n")
write.table(B, quote = FALSE, col.names =FALSE)
invisible(list(A=A, B=B, mlat=mlat, mlon=mlon, Jlat =Jlat, Jlon =Jlon, proj=proj))
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.