Jose,
This looks like a bug. Enclosed polygons can get overwritten by the hole in
the enclosing polygons. More later. Thanks for the useful example.
Robert

On Fri, Sep 7, 2012 at 6:55 AM, Jose Tanago <josetanago@gmail.com> wrote:

> Thanks roman
>
> I am newbie in this kind of posting, so i dont know exactly what you
> mean with self-confined reproductible example.
>
> I have created a super small subset of the data as example (few
> polygons and a super small raster with same extent) as well as the
> rasterized output. all in a dropbox folder of 40kb in total.
>
> the lines of my code I am running are:
>
> >segmshp<-readShapePoly("Segm2rasterize.shp")
> >ima<-raster("imagecroped.tiff")
> >rastsegm<-rasterize(segmshp, ima, field='ID', fun=max,
> filename='segmRast.tif')
>
> After running the rasterize() I got this message in the console:
>
> > Found 11 region(s) and 14 polygon(s)
>
> A priori I identified two polygons that are not rasterized: (ID 636)
> and (ID 163)
>
> I am sending to you the dropbox link here:
> https://www.dropbox.com/sh/snrfz1p5tlaok5i/4z_OGOJ4Y-
>
> Please let me know if that works fine for you.
>
> Thank you vey much
> Jose
>
>
> 2012/9/7 Roman Luštrik <roman.lustrik@gmail.com>:
> > Can you provide a small, self-confined reproducible example?
> >
> > Cheers,
> > Roman
> >
> >
> > On Fri, Sep 7, 2012 at 1:57 PM, Jose Tanago <josetanago@gmail.com>
> wrote:
> >>
> >> Dear R-geo experts
> >> I have experienced a weird problem rasterizing polygons with the
> >> function "rasterize()" of Raster package.
> >> Polygons that falls within the boundaries of another polygon are not
> >> rasterized for an unknown reason. The output leave the inner pixels
> >> (the ones overlaying the ignored polygons) without data
> >>
> >> The ignored polygons are completely surrounded by another polygon but
> >> the outer polygon do not spatially overlay (just enclosed). All the
> >> polygons were created by segmenting a raster with the same pixel size
> >> and exactly equal spatial structure. So, the concerned empty pixels
> >> centers are perfectly covered by the inner polygon (the one ignored),
> >> while the outer polygon do not cover at all the concerned pixels.
> >>
> >> Do you have any idea what is wrong and how to overcome it?
> >> (I have rasterized the same vector with ARGIS and those polygons
> >> ignored in rasterize() are well rasterized by the ARCGIS tool)
> >> thanks in advance for any hint.
> >> Best
> >> Jose
> >>
> >> _______________________________________________
> >> R-sig-Geo mailing list
> >> R-sig-Geo@r-project.org
> >> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
> >
> >
> >
> >
> > --
> > In God we trust, all others bring data.
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>

	[[alternative HTML version deleted]]

