<FONT face="Default Sans Serif,Verdana,Arial,Helvetica,sans-serif" size=2><P style="MARGIN: 0px">Hi,</P><P style="MARGIN: 0px">&nbsp;</P><P style="MARGIN: 0px">First, thanks for all help so far...but I still have some problems.<BR>Here is the code:<BR>cumark &lt;- read.table("k:/test/helamarkcu.txt", header=T, sep = ";")</P><P style="MARGIN: 0px">memory.limit(size = 1897)&nbsp;&nbsp;&nbsp; <BR>attach(cumark)<BR>vew &lt;- as.vector(ew)<BR>vns &lt;- as.vector(ns)<BR>vcu &lt;- as.vector(cu)</P><P style="MARGIN: 0px">fcumark &lt;- data.frame(vew, vns, vcu)</P><P style="MARGIN: 0px">coordinates(fcumark) &lt;- ~vew+vns&nbsp; </P><P style="MARGIN: 0px">L1 &lt;-round((max(vew)- min(vew))/100,0)<BR>L2 &lt;-round((max(vns)- min(vns))/100,0)</P><P style="MARGIN: 0px"><BR>fcumark.interp.old &lt;- interp.old(vew, vns, vcu,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; xo= seq(min(vew), max(vew), length = L1),<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; yo= seq(min(vns), max(vns), length = L2),<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ncp = 0,&nbsp; <BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; extrap= FALSE,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; duplicate = "mean",<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; dupfun = NULL)<BR>Then I get this;<BR>Error: cannot allocate vector of size 351.8 Mb</P><P style="MARGIN: 0px">When I tested with a smaller file (~2000 rows) I come to the end without problems.</P><P style="MARGIN: 0px">Rest of the code..</P><P style="MARGIN: 0px">fcumark.interp.old.grid &lt;- as.SpatialGridDataFrame.im(as.im(fcumark.interp.old))</P><P style="MARGIN: 0px">##&nbsp; spplot(fcumark.interp.old.grid, col.regions=bpy.colors(), sp.layout=list("sp.points", pch="+", fcumark))</P><P style="MARGIN: 0px">slot(slot(fcumark.interp.old.grid, "grid"), "cellsize") &lt;- rep(mean(slot(slot(fcumark.interp.old.grid,<BR>&nbsp;&nbsp; "grid"), "cellsize")), 2)</P><P style="MARGIN: 0px">write.asciigrid(fcumark.interp.old.grid[1], "k:/test/helamarkcu.asc", na.value = -9999)<BR>## End run</P><P style="MARGIN: 0px">helamarkcu.txt is a plain ascii file, looks like...<BR>ew;ns;cu<BR>1344628;6138662;24.0<BR>1340104;6138880;17.0<BR>1342246;6138946;21.0<BR>1347259;6140911;19.0<BR>1339968;6140932;15.0<BR>1349906;6140987;21.0<BR>...<BR>..<BR>.<BR>(38155 rows)<BR>Have anybody a solution? Writing result from interp.old direct to disk/file on the "fly"?</P><P style="MARGIN: 0px">Running R 2.7.2&nbsp;and WinXP</P><P style="MARGIN: 0px">Regards,<BR>Mikael</P></FONT>