I made two concentric buffers with 90-100km radius, from point shapefiles (hydroeletric dams here: http://maps.csr.ufmg.br/) cropped from Legal Amazon shape (http://www.dpi.inpe.br/amb_data/Shapefiles/UF_AmLeg_LLwgs84.zip), and use Difference (geoprocessing tool) with QGIS, resulting in:
After that, I used the rasterize of the raster package in the R, and resulted in these two buffers without the intersect clipping:
shp_100=readOGR(dsn="D:/Alunos/Abilio/buffer_2005",layer="90_100") rst_template <- raster(res= 0.01, crs = projection(des15AmL_1km), ext = extent(c(-74, -40, -19, 6))) buffer_100=rasterize(shp_100, rst_template, field="Nome")
Has anyone ever experienced this? Anyone know which shape feature might be resulting in this?