[R-sig-Geo] Error In proj4string

[Tom Philippi]

Without more information, we can't give you a definitive answer.

Try removing the quotes from your cbind.
temp <- cbind(bio1, bio2, bio3,
           bio4, bio5, bio6, bio7, bio8, bio9, bio10, bio11)

Another possibility is that one or more of your 11 grids has a different
topology, although that might throw an error at cbind rather than produce a
different class of object.

You might consider the raster package, as I think you are attempting to
stack the grid values into a raster brick or raster stack.

I hope that this helps...

Tom 2


On Mon, Oct 20, 2014 at 11:21 AM, alannie <agg18 at pitt.edu> 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")
> > 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
>

	[[alternative HTML version deleted]]