Dear list, I´m getting the following error when trying to convert a spatial lines object into a raster using the rasterize function (package raster):
> my.line An object of class "SpatialLines" Slot "lines": [[1]] An object of class "Lines" Slot "Lines": [[1]] An object of class "Line" Slot "coords": lon lat [1,] -22.02833 63.40222 [2,] -22.00111 63.40500 [3,] -21.98444 63.40722 [4,] -21.95722 63.40861 [5,] -21.93944 63.40861 [6,] -21.91306 63.40583 [7,] -21.88556 63.40194 [8,] -21.86833 63.39750 [9,] -21.84139 63.39222 [10,] -21.82556 63.38722 [11,] -21.79972 63.39306 [12,] -21.80583 63.39083 [13,] -21.84028 63.39556 [14,] -21.86222 63.40194 [15,] -21.88694 63.40528 [16,] -21.92389 63.40500 [17,] -21.94972 63.40417 [18,] -21.97500 63.40389 [19,] -21.99444 63.40194 [20,] -22.01750 63.40444 Slot "ID": [1] "a" Slot "bbox": min max x -22.02833 -21.79972 y 63.38722 63.40861 Slot "proj4string": CRS arguments: +proj=longlat +ellps=WGS84 > my.raster class : RasterLayer filename : nrow : 10 ncol : 52 ncell : 520 min value : NA max value : NA projection : +proj=longlat +ellps=WGS84 extent : -22.04, -21.785, 63.375, 63.42 (xmin, xmax, ymin, ymax) resolution : 0.004903846, 0.0045 (x, y) > rasterize(my.line,my.raster) Error in if (outcol < 1) next : missing value where TRUE/FALSE needed Is this a bug in the raster function? Or I am missing something? Many thanks, Julián Julián Mariano Burgos Hafrannsóknastofnunin/Marine Research Institute Skúlagata 4, 121 Reykjavík, Iceland Sími/Telephone : +354-5752037 Bréfsími/Telefax: +354-5752001 Netfang/Email: jul...@hafro.is, jmbur...@uw.edu [[alternative HTML version deleted]]
_______________________________________________ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo