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!
------------------------------------------------------------------------
If you reply to this email, your message will be added to the
discussion below:
http://r-sig-geo.2731867.n2.nabble.com/Error-In-proj4string-tp7587305.html
To unsubscribe from R-sig-geo, click here
<http://r-sig-geo.2731867.n2.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=2731867&code=bC5oaWxnZXJAa3UuZGV8MjczMTg2N3wtODM0Mzg5NjIw>.
NAML
<http://r-sig-geo.2731867.n2.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml>