# [R-sig-Geo] Polygon dissolution or union issue

Sébastien Durand v8extra at gmail.com
Sun Dec 5 07:00:50 CET 2010

```Hello,

Could somebody find me the solution to this issue!

I cannot manage to merge or dissolve all the polygon  generated in the object "pol"

Here is my example:

library(akima)
library(maptools)
library(raster)

data(akima)
# define de dimension of grid
dimgrid=50

# Interpolate to regular grid
ak.li=interp(akima\$x, akima\$y, akima\$z, xo=seq(min(akima\$x), max(akima\$x), length = dimgrid),yo=seq(min(akima\$y), max(akima\$y), length = dimgrid),linear = TRUE, extrap=FALSE, duplicate = "mean")

# Show interpolation
image(ak.li)
points(akima)
with(akima, text(x, y, formatC(z,dig=2), adj = -0.1))

# Create a raster with the data
r <- raster(ak.li)
plot(r)
pol <- rasterToPolygons(r, fun=function(x){x>20 & x<25})

r <- raster(nrow=length(ak.li\$y), ncol=length(ak.li\$x), crs=NA)
tmp=ak.li\$z
tmp=apply(tmp, 1, rev)
r[] <- tmp

# Show the raster
plot(r)

# Convert the data found between ]20 and 25[ into polygons
pol <- rasterToPolygons(r, fun=function(x){x>20 & x<25})

# Show the polygons

# Attempt to merge similar adjacent into one polygon
union = unionSpatialPolygons(pol, ID=rep(1,
times=length(pol at polygons)))

# Show