<div dir="ltr">And this is for a raster with around 5 milliion cells? Your results surprise me a bit, especially if the bottleneck is at the call to buffer()/distance().  Your best bet is to wait for a response from Robert Hijmans or someone else working on the package, they might have some other tricks up their sleeves. In any case, I agree, 5+ hours is extremely excessive.<div><br></div><div>Sincerely,</div><div>Forrest</div></div><br><div class="gmail_quote"><div dir="ltr">On Mon, Apr 11, 2016 at 6:57 PM Nelly Reduan <<a href="mailto:nell.redu@hotmail.fr">nell.redu@hotmail.fr</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">




<div dir="ltr">
<div style="font-size:12pt;color:#000000;background-color:#ffffff;font-family:Calibri,Arial,Helvetica,sans-serif">
<p>Thank you very much Forrest for your answer. I don't think I have memory problems (my computer has a 32 GB memory).
</p>
<p><br>
</p>
<p>Here is my code to draw the 100-m buffers with the function "buffer" (package raster):</p>
<p>r1 <- r ## r is the original raster and has a resolution of 10m<br>
r1[(r1[]!=2 & r1[]!=5)]=NA<br>
plot(r1)<br>
r2 <- buffer(r1, 10)</p>
<p><br>
</p>
<p>I also tested the function distance (package raster):</p>
<p>r1 <- r ## r is the original raster and has a resolution of 10m<br>
r1[(r1[]!=2 & r1[]!=5)]=NA<br>
plot(r1)<br>
r3 <- distance(r1)</p>
<p><br>
</p>
<p>In the two cases, I stopped the functions afer 5 hours as this was too long ! For the moment, I haven't find solutions.</p></div></div><div dir="ltr"><div style="font-size:12pt;color:#000000;background-color:#ffffff;font-family:Calibri,Arial,Helvetica,sans-serif">
<p><br>
</p>
<p></p>
<p>Thanks a lot for your time.</p>
<p>Have a nice day.</p>
<p>Nell</p>
<p></p>
<p><br>
</p>
<br>
<br>
</div></div><div dir="ltr"><div style="font-size:12pt;color:#000000;background-color:#ffffff;font-family:Calibri,Arial,Helvetica,sans-serif"><div style="color:rgb(0,0,0)">
<hr style="width:98%;display:inline-block">
<div dir="ltr"><font color="#000000" face="Calibri, sans-serif" style="font-size:11pt"><b>De :</b> Forrest Stevens <<a href="mailto:r-sig-geo@forreststevens.com" target="_blank">r-sig-geo@forreststevens.com</a>><br>
<b>Envoyé :</b> lundi 11 avril 2016 09:09<br>
<b>À :</b> Nelly Reduan; <a href="mailto:r-sig-geo@r-project.org" target="_blank">r-sig-geo@r-project.org</a><br>
<b>Objet :</b> Re: [R-sig-Geo] How to reduce the buffering time with the function "buffer" (package raster)</font>
<div> </div>
</div></div></div></div><div dir="ltr"><div style="font-size:12pt;color:#000000;background-color:#ffffff;font-family:Calibri,Arial,Helvetica,sans-serif"><div style="color:rgb(0,0,0)">
<div>
<div dir="ltr">Five million cells isn't all that many, how slow is too slow? Buffering a raster of about 10 million cells on my laptop takes on the order of 20 seconds or so for a binary raster. Is it possible that you're fighting memory problems?
<div><br>
</div>
<div>In the past when doing multiple ring buffers I've found it faster to calculate a distance-to raster first, and then apply multiple logical comparisons on the distance raster.  I don't know if this applies to your situation or not but it's one trick that
 might help.</div>
<div><br>
</div>
<div>Sincerely,</div>
<div>Forrest<br>
<br>
<div class="gmail_quote">
<div dir="ltr">On Mon, Apr 11, 2016 at 10:58 AM Nelly Reduan <<a title="mailto:nell.redu@hotmail.fr
Ctrl+Cliquer ou appuyer pour suivre le lien" href="mailto:nell.redu@hotmail.fr" target="_blank">nell.redu@hotmail.fr</a>> wrote:<br>
</div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;padding-left:1ex;border-left-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid">
<div dir="ltr">
<div style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;background-color:rgb(255,255,255)">
<p><span></span></p>
<p>Hello,</p>
<p><br>
</p>
<p>I would like to build 100-m buffers (representing ecological corridors) around land cover attributes in a raster layer and to assign a given value (i.e., value of 13) to buffer cells.</p>
<p><br>
</p>
<p>To do this, I'm using the function "buffer" (package raster) to draw the buffers around particular raster cells (i.e., cell values of 2 and 5). For example, in the image below, the buffers are represented in blue (cell values of 13) and the particular raster
 cells are represented in yellow (cell values of 2).</p>
<p><br>
</p>
<p>The problem is that the time to run the function "buffer" is very low because I have a raster of 5000000 cells. Does anyone know how I can reduce the buffering time ?</p>
<p><br>
</p>
<p>Thanks a lot for your time.</p>
<p>Have a nice day.</p>
<p>Nell</p>
<span></span><br>
<p></p>
<img style="max-width:100%" alt="Capture.PNG" src="cid:154060fc7216d0e7b411">
<p></p>
</div>
</div>
_______________________________________________<br>
R-sig-Geo mailing list<br>
<a href="mailto:R-sig-Geo@r-project.org" target="_blank">R-sig-Geo@r-project.org</a><br>
<a href="https://stat.ethz.ch/mailman/listinfo/r-sig-geo" rel="noreferrer" target="_blank">https://stat.ethz.ch/mailman/listinfo/r-sig-geo</a></blockquote>
</div>
</div>
</div>
</div>
</div></div></div></blockquote></div>