[R-sig-Geo] Error In proj4string

[Alexandre Villers]

Good evening !

You should either read a bit more about S4 sp object (in particular,
if you want to combine object that have exactly the same structure) or
use the stack function of the raster package
(st.bio<-stack("bio1.asc", bio2.asc", ..., "bio11.asc")).
You can't just cbind() Spatial Objects (someone might correct this
statment if not true). SpatialGridDataFrames as the ones produced with
readGDAL() have a data slot where you can store different vectors of
values.
With some googling, you should find what you are looking for.

HTH

Alex


2014-10-20 20:21 GMT+02:00 alannie <agg18 at pitt.edu>:
> Hi list!
>
> I have been trying to run some code that I have run tons of times before,
> only to get a warning message:
>
> ##add in ASCII raster maps bio1-11
>> bio1 <- readGDAL("bio1.asc")
> bio1.asc has GDAL driver AAIGrid
> and has 1486 rows and 1347 columns
>> bio2 <- readGDAL("bio2.asc")
> bio2.asc has GDAL driver AAIGrid
> and has 1486 rows and 1347 columns
>> bio3 <- readGDAL("bio3.asc")
> bio3.asc has GDAL driver AAIGrid
> and has 1486 rows and 1347 columns
>> bio4 <- readGDAL("bio4.asc")
> bio4.asc has GDAL driver AAIGrid
> and has 1486 rows and 1347 columns
>> bio5 <- readGDAL("bio5.asc")
> bio5.asc has GDAL driver AAIGrid
> and has 1486 rows and 1347 columns
>> bio6 <- readGDAL("bio6.asc")
> bio6.asc has GDAL driver AAIGrid
> and has 1486 rows and 1347 columns
>> bio7 <- readGDAL("bio7.asc")
> bio7.asc has GDAL driver AAIGrid
> and has 1486 rows and 1347 columns
>> bio8 <- readGDAL("bio8.asc")
> bio8.asc has GDAL driver AAIGrid
> and has 1486 rows and 1347 columns
>> bio9 <- readGDAL("bio9.asc")
> bio9.asc has GDAL driver AAIGrid
> and has 1486 rows and 1347 columns
>> bio10 <- readGDAL("bio10.asc")
> bio10.asc has GDAL driver AAIGrid
> and has 1486 rows and 1347 columns
>> bio11 <- readGDAL("bio11.asc")
> bio11.asc has GDAL driver AAIGrid
> and has 1486 rows and 1347 columns
> ###combine raster maps
>> temp <- cbind("bio1", "bio2", "bio3",
>> "bio4","bio5","bio6","bio7","bio8","bio9","bio10","bio11")
>> names(temp) <- c("bio1", "bio2", "bio3",
>> "bio4","bio5","bio6","bio7","bio8","bio9","bio10","bio11")
> ###read in point file lin.shp
>> lin <- readShapePoints("lin.shp")
> ###set proj4strings to the same proj4string: *this is where the issue is*
>> proj4string(lin)= proj4string(temp)
> *Error in (function (classes, fdef, mtable)  :
>   unable to find an inherited method for function ‘proj4string’ for
> signature ‘"matrix"’*
> ###what is in "temp" ?
>> temp
>      [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8]   [,9]   [,10]
> [,11]
> [1,] "bio1" "bio2" "bio3" "bio4" "bio5" "bio6" "bio7" "bio8" "bio9" "bio10"
> "bio11"
> attr(,"names")
>  [1] "bio1"  "bio2"  "bio3"  "bio4"  "bio5"  "bio6"  "bio7"  "bio8"  "bio9"
> "bio10" "bio11"
>> class(temp)
> [1] "character"
> ###end
>
> I am not sure why this issue is occurring.
> I have checked all the projections - all files are in the same projection.
> I have restarted R - nothing gets better.
> I have tried data for a different organism (it is the same data, but
> belonging to a different species) - this does not work.
> I have the following packages loaded:
> library(maptools)
> library(adehabitatHS)
> library(SDMTools)
> library(spatial)
> library(sp)
> library(rgdal)
> library(dismo)
> library(shapefiles)
>
>
> Thanks, List!
>
>
>
>
>
>
>
>
>
> --
> View this message in context: http://r-sig-geo.2731867.n2.nabble.com/Error-In-proj4string-tp7587305.html
> Sent from the R-sig-geo mailing list archive at Nabble.com.
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo



-- 
Alexandre Villers, PhD.
Postdoctoral researcher
Team Agripop
CEBC-CNRS, UMR 7372
79360 Beauvoir sur Niort, FR