Gustaf, Thanks, I think I found the bug. It happens when there is no
overlap at all between polygons and raster. I think it is fixed in raster
2.0-39 (you can install from here
https://r-forge.r-project.org/R/?group_id=294). Robert
On Tue, Dec 18, 2012 at 10:35 AM, Gustaf Granath
<gustaf.gran...@ebc.uu.se>wrote:
> Robert,
>
> Thanks for looking at this. Unfortunately, updating R and Raster did not
> help. Here is info you requested. Hope that helps!
>
> Gustaf
>
> > tt<-rasterize(SpPDF,tes,field="idd")
>
> Error in writeValues(out, fun(n), tr$row[i]) :
> error in evaluating the argument 'v' in selecting a method for function
> 'writeValues': Error in fun(n) : unused argument(s) (n)
> > traceback()
> 6: writeValues(out, fun(n), tr$row[i])
> 5: init(rstr, function() NA)
> 4: .polygonsToRaster(x, y, field = field, fun = fun, background =
> background,
> mask = mask, update = update, updateValue = updateValue,
> filename = filename, getCover = getCover, silent = silent,
> ...)
> 3: .local(x, y, ...)
> 2: rasterize(SpPDF, tes, field = "idd")
> 1: rasterize(SpPDF, tes, field = "idd")
> >
>
> > show(SpPDF)
> class : SpatialPolygonsDataFrame
> nfeatures : 58824
> extent : -11.08333, 31.91753, 33.91667, 71.91743 (xmin, xmax, ymin,
> ymax)
> coord. ref. : +proj=longlat +datum=WGS84
> nvariables : 92
> names : X.Lon, Lat, JulT00, P1900, NH1900, N1900, JulT10,
> P1910, NH1910, N1910, JulT20, P1920, NH1920, N1920, JulT30, ...
> min values : -10.5, 34.8333, 0.3, 260, 0.103, 0.0163, 0,
> 261, 0.103, 0.0184, 0.1, 245, 0.118, 0.0215, 0.3, ...
> max values : 31.8333, 71, 27.2, 3267, 14.837, 1.5572, 27.2,
> 3226, 11.927, 1.2981, 27.8, 3303, 12.945, 1.4166, 27.5, ...
> Warning messages:
> 1: closing unused connection 4
> (C:/Users/veGustaf/AppData/Local/Temp/R_raster_tmp/veGustaf/raster_tmp_2012-12-18_192108_37929.gri)
>
> 2: closing unused connection 3
> (C:/Users/veGustaf/AppData/Local/Temp/R_raster_tmp/veGustaf/raster_tmp_2012-12-18_191936_11753.gri)
>
> >
>
> > show(tes)
>
> class : RasterLayer
> dimensions : 12768, 11050, 141086400 (nrow, ncol, ncell)
> resolution : 250, 250 (x, y)
> extent : 2782250, 5544750, 2220250, 5412250 (xmin, xmax, ymin, ymax)
> coord. ref. : +proj=longlat +datum=WGS84
>
>
> > sessionInfo()
> R version 2.15.2 (2012-10-26)
> Platform: i386-w64-mingw32/i386 (32-bit)
>
> locale:
> [1] LC_COLLATE=Swedish_Sweden.1252 LC_CTYPE=Swedish_Sweden.1252
> [3] LC_MONETARY=Swedish_Sweden.1252 LC_NUMERIC=C
> [5] LC_TIME=Swedish_Sweden.1252
>
> attached base packages:
> [1] stats graphics grDevices utils datasets methods base
>
> other attached packages:
> [1] raster_2.0-31 sp_1.0-2
>
> loaded via a namespace (and not attached):
> [1] grid_2.15.2 lattice_0.20-10 tools_2.15.2
> >
>
>
> On 2012-12-17 07:32, Robert J. Hijmans wrote:
>
> Gustaf,
>
> When reporting that something does not work, please include the result
> of your
>
> sessionInfo()
>
> Just after the error occurs.
>
> traceback()
>
> can also be helpful. In this case it would appear that that you are not
> using the current version of raster. If that is the case, please
> update.packages() first (perhaps after getting the most recent version of
> R). Some more info on the SpatialPolygonsDataFrame that causes the error
> would also be helpful (use the show function after raster is loaded).
>
>
> Robert
>
>
> On Wed, Dec 12, 2012 at 11:55 PM, Gustaf Granath <gustaf.gran...@ebc.uu.se
> > wrote:
>
>> Hi,
>> I have used the rasterize function to get values from a
>> spatialPolyDataFrame variable (@data$variable) to raster:
>> rasterize(SPDFobj, rasterObj, field="variable")
>>
>> For a specific case a get a, to me, strange error message .
>> "Error in writeValues(out, fun(n), tr$row[i]) :
>> error in evaluating the argument 'v' in selecting a method for function
>> 'writeValues': Error in fun(n) : unused argument(s) (n)"
>>
>> It works with one SpatialPolyDataFrame that Im working with, but not with
>> another one and I cannot figure out why. Any ideas of what is going on?
>>
>> Unfortunately, I have not been able to reproduce the error with simulated
>> data (and my data files are huge) but this is essentially what I have done
>> (but the code below works - actually the raster layer in this code works
>> with my SPDFobj):
>> p1 <- rbind(c(-180,-20), c(-140,55), c(10, 0), c(-140,-60), c(-180,-20))
>> p2 <- rbind(c(-10,0), c(140,60), c(160,0), c(140,-55), c(-10,0))
>> p3 <- rbind(c(-125,0), c(0,60), c(40,5), c(15,-45), c(-125,0))
>> pols <- SpatialPolygons( list( Polygons(list(Polygon(p1),
>> Polygon(hole)), 1),
>> Polygons(list(Polygon(p2)), 2),
>> Polygons(list(Polygon(p3)), 3)))
>> pols <-SpatialPolygonsDataFrame(pols, data.frame("W"=1:3))
>> r <- raster(ncol=80, nrow=90)
>> test<-rasterize(pols, r, field="W")
>>
>> Here is an overview of the objects I am using.
>>
>> > summary(SPDFobj)
>> Object of class SpatialPolygonsDataFrame
>> Coordinates:
>> min max
>> x -11.08333 31.91753
>> y 33.91667 71.91743
>> Is projected: FALSE
>> proj4string : [+proj=longlat +datum=WGS84]
>> Data attributes:
>> variable
>> ...
>> ...
>>
>> > rasterObj_orig
>> class : RasterLayer
>> dimensions : 12768, 11050, 141086400 (nrow, ncol, ncell)
>> resolution : 250, 250 (x, y)
>> extent : 2782250, 5544750, 2220250, 5412250 (xmin, xmax, ymin, ymax)
>> coord. ref. : +proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000
>> +ellps=GRS80 +units=m +no_defs
>> data source : C:\Users\veGustaf\Desktop\phd\Where is Ndep neg\bogs.grd
>> names : layer
>> values : 36, 36 (min, max)
>>
>> > rasterObj #projection fixed etc through an extent object, raster -->
>> extent --> raster
>> class : RasterLayer
>> dimensions : 12768, 11050, 141086400 (nrow, ncol, ncell)
>> resolution : 250, 250 (x, y)
>> extent : 2782250, 5544750, 2220250, 5412250 (xmin, xmax, ymin, ymax)
>> coord. ref. : +proj=longlat +datum=WGS84
>>
>>
>> Cheers,
>>
>> Gustaf
>>
>> --
>> Gustaf Granath (PhD)
>>
>> Dept of Plant Ecology and Evolution
>> Uppsala University
>>
>> _______________________________________________
>> R-sig-Geo mailing list
>> R-sig-Geo@r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>>
>
>
> --
> Gustaf Granath (PhD)
> Dept of Plant Ecology and Evolution
> Uppsala University
>
>
>
[[alternative HTML version deleted]]
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
