Hi,
Doing more search i've discovered package RGDAL that can write a geotiff file
with projection. I saved a geotiff file in UTM projection and if i read the
file back in R and check the projection seems that everything is OK. But if i
load the file in ArcGIS (ESRI product) i get the warning that the file is
missing spatial reference so it cannot be projected, but it is displayed
correctly. I guess somewhere i am doing a mistake when i define the projection
in R. My code follows:
data.grid <- read.csv(x, header=TRUE)
gridded(data.grid) = ~East.m.+North.m.
proj4string(data.grid) = CRS("+proj=tmerc +lat_0=0.00000000000
+lon_0=-81.00000000000 +k=0.99960000 +x_0=500000.000000 +y_0=0.0000000
+ellps=GRS80 +units=m")
tr <- "e:\\JELA_veg\\test_gtiff\\test.tif"
writeGDAL(data.grid["class.pca"], tr)
mg3 <- readGDAL(tr)
proj4string(mg3)
[1] " +proj=utm +zone=17 +ellps=GRS80 +units=m +no_defs"
I will really appreciate if anybody can point me in the right dirrection.
Thanks,
Monica
_________________________________________________________________
[[alternative HTML version deleted]]
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
