[R-sig-Geo] output gstat and maps

[Edzer J. Pebesma]

Mmm, yeah, the trouble is that R is not a GIS, and
GMT and ArcView are no statistics packages.

I think you can import shapefiles using packages
maptools, and perhaps Rmap (not on CRAN?).
I imaging you can then plot the values with

image(values.krige)

and add polygons to the plot with the plot
commands in maptools, using the add=TRUE
argument.

levelplot, which I like much more because it
draws a legend, can currently not do this
easily (i.e., without customizing panel functions);
has anyone looked into this?
--
Edzer

Marcelo Bruno wrote:

> Hi, members
>
> I'm Marcelo, oceanographer from Brazil
> and actually work with fishery acoustical data
> of brazilian shelf. I'm using R, GMT, ArcView (windows)
> in Linux redhat 9.0 (Athlon Xp 1.6Ghz, 20Gb 540mb mem).
> I run kriging analysis using gstat (R1.8.1) in
> fishery acoustic data:
>
> x     y     sa (backscatering acoustic parameter)
>
> the output is ordinary kriging surface.
> I need overlap these output to a map of
> Brazil coastal shelf and slope but i don't
> make this. I have a file (ascii) of Brazil coastal,
> and shapefile. Someone can helps me with this
> problem? I try in gimp, kpaint and others but i
> need a better solution to learn more in R.
>
> Thanks in advance!
> Marcelo
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at stat.math.ethz.ch
> https://www.stat.math.ethz.ch/mailman/listinfo/r-sig-geo