search for: showwkt

Hi, ? I am trying to write a shapefile with projection. I have my data in a data.frame called try and consists in xy coordinates and a numerical attribute value z1. ? Libraries loaded are: sp, rgdal, raster, maptools ? head(try) ???????? x?????? ??? y??? ? ? ? ? ? ? ??? z1 1 610237.1 ???????? 3375751 ???????????? 8.221 2 610236.1 ???????? 3375750 ???????????? 8.153 3 610236.1 ????????

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.