[R-sig-Geo] asciigrid R

Michael Sumner mdsumner at utas.edu.au
Mon Dec 3 21:54:47 CET 2007


Try this, you were nulling each column as you read it, and you didn't 
let us know how "expanatory" was created in the first place:

explanatory = read.asciigrid("silt.asc")
names(explanatory) <- "silt"

explanatory$sand = read.asciigrid("sandasc.asc")$band1

explanatory$nd = read.asciigrid("nd.asc")$band1

explanatory$bnd3 = read.asciigrid("bnd3.asc")$band1


Jose Funes wrote:
> ASCIIGRID
>
> I am using read.asciigrid function from R to import 4 grid maps from
> ArcGIS to a single grid in R. However, I am having troubles to make it
> work, when I finish to import the 4 grid to R the final grid appear as
> NULL. Below I did include the code, I will appreciate your help to
> solve this problem.I tried to register to GTAT group but my email is
> returned.
>
> Thanks,
>
> R code:
>
> Four layers were stacked:
> explanatory$silt = read.asciigrid("silt.asc")$band1
> explanatory$silt = explanatory$band1
>
> explanatory$sand = read.asciigrid("sandasc.asc")$band1
> explanatory$sand = explanatory$band1
>
> explanatory$nd = read.asciigrid("nd.asc")$band1
> explanatory$nd = explanatory$band1
>
> explanatory$bnd3 = read.asciigrid("bnd3.asc")$band1
> explanatory$bnd3 = explanatory$band1
>
>
> Deleting double band:
> explanatory$band1 = NULL
>
>  str(explanatory)
>  NULL
>
> Jose Funes
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
>
>




More information about the R-sig-Geo mailing list