R/summary.bathy.R

In ericpante/marmap: Import, Plot and Analyze Bathymetric and Topographic Data

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]
}