Error In proj4string
Alannie, you get an error, not a warning, and I doubt this code has ever worked before. Read on:
On 10/20/2014 08:21 PM, alannie wrote:
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")
You are cbind-ing length-one character vectors, but you probably want to: temp <- cbind(bio1, bio2, bio3) # etc i.e. cbind-ing the SpatialGridDataFrame objects. Indeed,
library(sp) methods(cbind)
[1] cbind.data.frame cbind.SpatialGridDataFrame* [3] cbind.ts* sp provides a method for this. After this, the rest should work.
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
Edzer Pebesma, Co-Editor-in-Chief Computers & Geosciences Institute for Geoinformatics (ifgi), University of M?nster Heisenbergstra?e 2, 48149 M?nster, Germany. Phone: +49 251 83 33081 http://ifgi.uni-muenster.de GPG key ID 0xAC227795 -------------- next part -------------- A non-text attachment was scrubbed... Name: signature.asc Type: application/pgp-signature Size: 490 bytes Desc: OpenPGP digital signature URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20141020/ae324244/attachment.bin>