source("~/lib/R/filledcontour.R")
library(akima)
library(maps)
library(maptools)

map <- read.shape("~/lib/R/SHP_MAP/japan_jdg.shp")

x <- seq(139,140,0.1)
y <- 35.2 + runif(11)
z <- 1:11

filled.contour(interp(x,y,z))

plot(map, fg="transparent", add=T)