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

Moritz Lennert mlennert at club.worldonline.be
Sun Nov 18 18:08:39 EST 2007


Tomas Lanczos wrote:
> Hello,
> 
> thank You Boris, it is something new for me, how can I integrate octave 
> to grass? is there a tutorial for that?

No, I think you have to export your map with r.out.mat and then use the 
result in matlab.

However, maybe r3.cross.rast can do the trick for you directly in GRASS ?

Moritz

> 
> thank You
> 
> Tomas
> 
> 
> Boris Avdeev wrote:
>> 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
>>>
>>
>>
> 
> _______________________________________________
> grass-user mailing list
> grass-user at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-user



More information about the grass-user mailing list