Is ground level varying, over your domain? If yes, the new set of voxels may no longer be a regular grid in the new CRS; you could represent them as a SpatialPointsDataFrame. If no, you could simply modify obj@grid@cellcentre.offset[3]

On 02/13/2013 10:05 PM, Nathalie Morin wrote:
Hello,



I would like to change or assign a new value of the z coordinates
intrinsically in a 3 dimensional SpatialGridDataFrame from absolute height
in UTM to relative height above the ground, in order to remove slope. So
basically, I would like to move the current position of the voxels according
to a new vector of coordinates, but by keeping the same data attribute
contained in each of the original voxels.



I tried creating a new SpatialGridDataFrame from (i) a new GridTopology with
z coordinates in relative height and from (ii) the same vector of attribute,
but of course this does not actually move the cubes accordingly in the 3
dimensional space. It just changes the z value of the voxel from UTM to
relative height above the ground. Same thing when I try doing this with a
regular 3D Array.



Does anybody have any hints in how to achieve this with (or without) sp ?

Cheers,

Nathalie



--

Nathalie MORIN

Doctorante en télédétection forestière (Ph.D. student in remote sensing
applied to forestry)

Centre d’application et de recherches en télédétection (CARTEL)

Département de géomatique appliquée

Université de Sherbrooke

Sherbrooke (QC) Canada J1K 2R1

Tel. : + 1 (819)  821-8000 #62962

e-mail :  <mailto:nathalie.m.mo...@usherbrooke.ca>
nathalie.m.mo...@usherbrooke.ca

Page Web :  <http://pages.usherbrooke.ca/cartel/morin_nathalie.php>
http://pages.usherbrooke.ca/cartel/morin_nathalie.php




        [[alternative HTML version deleted]]



_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo


--
Edzer Pebesma
Institute for Geoinformatics (ifgi), University of Münster
Weseler Straße 253, 48151 Münster, Germany. Phone: +49 251
8333081, Fax: +49 251 8339763  http://ifgi.uni-muenster.de
http://www.52north.org/geostatistics      e.pebe...@wwu.de

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo

Reply via email to