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:

enter image description here

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")

enter image description here

Has anyone ever experienced this? Anyone know which shape feature might be resulting in this?

  • 1
    You are also missing one buffer top right. Add your R code to the question – BERA Mar 6 at 14:11
  • BERA, It's just the color that's white ... I've already looked at this buffer and it's present in the raster with all its values. Anyway, thank you! – Abilio André Mar 6 at 16:19

Your Answer

By clicking "Post Your Answer", you acknowledge that you have read our updated terms of service, privacy policy and cookie policy, and that your continued use of the website is subject to these policies.

Browse other questions tagged or ask your own question.