summary.bathy = function(object, ...){
round(min(as.numeric(colnames(object))),2) -> lat.min
round(max(as.numeric(colnames(object))),2) -> lat.max
round(min(as.numeric(rownames(object))),2) -> lon.min
round(max(as.numeric(rownames(object))),2) -> lon.max
lon.max2 <- ifelse(lon.max > 180, lon.max-360, lon.max)
flag.l1 <- ifelse(lon.min < 0, "W", "E")
flag.l2 <- ifelse(lon.max2 < 0, "W", "E")
flag.l3 <- ifelse(lat.min < 0, "S", "N")
flag.l4 <- ifelse(lat.max < 0, "S", "N")
one.minute = 0.016667
as.numeric(rownames(object))[2] - as.numeric(rownames(object))[1] -> cell.size.centroid
round(cell.size.centroid / one.minute, 2) -> cell.size.minute
message(paste("Bathymetric data of class 'bathy', with",dim(object)[1],"rows and",dim(object)[2],"columns"))
message(paste("Latitudinal range: ", lat.min," to ", lat.max, " (",abs(lat.min)," ",flag.l3," to ",abs(lat.max)," ",flag.l4,")",sep=""))
message(paste("Longitudinal range: ", lon.min," to ", lon.max, " (",abs(lon.min)," ",flag.l1," to ",abs(lon.max2)," ",flag.l2,")",sep=""))
message(paste("Cell size:",cell.size.minute,"minute(s)"))
message("")
message("Depth statistics:")
print(summary(as.vector(object), ...))
message("")
message("First 5 columns and rows of the bathymetric matrix:")
object[1:5, 1:5]
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.