Using Matlab (or octave): n=1;index=zeros(prod(dim(rast2d)),3) for [i j]=1:dim(rast2d,[1 2]) k=round(rast2d(i,j)) index(n,:)=[i j k] end answer=rast3d(index)
Boris On Nov 17, 2007 6:14 PM, Tomas Lanczos <[EMAIL PROTECTED]> wrote: > Hello everybody, > > I have a 3D raster and several 2D elevation rasters. Is it possible to > perform somehow (within GRASS, or outside) intersection betwwen the 2D > and 3D rasters and get 2D raster with values of the 3D rasters where > they are intersected? > > Tomas > GRASS 6.3 CVS, 2.2, Ubuntu 7.10 Gutsy > > > _______________________________________________ > grass-user mailing list > [email protected] > http://lists.osgeo.org/mailman/listinfo/grass-user > _______________________________________________ grass-user mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-user
