Skip to content
Prev 156173 / 398506 Next

Blank certain areas of a contour plot

Hi,

I'm interpolating a list of syncronous accumulated precipitation 
observations collected over a number of raingauge stations sited over 
land, over a regular lat/lon grid using akima's interp().

Then, I plot and locate geographycally the resulting field with a 
filled.contour() and a call to map(). Everything is fine but I need to 
"blank" (in Golden Surfer nomenclature) the areas of my domain that lay 
over the sea, as I do not have observations there and structures are 
completely artificial.

Is there any function out there that can use the mapsdata database to 
"blank" the sea areas of my plot?
Otherwise, can anybody think of a shortcut to get a similar effect on 
the final map (not showing the interpolated precipitation areas over the 
sea?

Thanks a lot in advance,

V?ctor.