
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]]

