[GRASS-user] 2D - 3D raster intersection

[Boris Avdeev]

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 <lanczos at t-zones.sk> 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
> grass-user at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-user
>