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" ?
[,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"
[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.