scallops <- read.table("scallops.txt", header=T)

hist(scallops$tcatch)

myscallops <- scallops
myscallops[,"lgcatch"] <- log(scallops$tcatch+1)
summary(myscallops$lgcatch)
hist(myscallops$lgcatch)

plot(myscallops$long, myscallops$lat)

library(akima)
int.scp <- interp(myscallops$long,myscallops$lat,myscallops$lgcatch)

lat.lim <- range(myscallops$lat)
lon.lim <- range(myscallops$lon)
image(int.scp, xlim=lon.lim, ylim=lat.lim, xlab="Longitude", ylab="Latitude")
contour(int.scp, add=T)

#persp(int.scp,xlim=range(myscallops$lon),ylim=range(myscallops$lat))