Hi Hamish, thanks for the hints. I've already tried using r.region, but it seemed to me it just managed boundaries. How could it be used to shift xy?<br><br><br><div><span class="gmail_quote">2008/1/25, Hamish <<a href="mailto:hamish_b@yahoo.com">
hamish_b@yahoo.com</a>>:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">G. Allegri wrote:<br>> I don't have a raster that fits my region.
<br><br>this will do it:<br> r.mapcalc 'regionrast = 1'<br><br><br>> The unsolved problem is rotation (and shift): making affine<br>> tranformations using just a x-shift and y-shift parameter, and a<br>> rotation angle having set a point as origin (or just around the
<br>> "center" of the raster).<br><br>some ideas.<br><br>simple translation without rotation: r.region<br><br>90 deg rotation / flipping:<br> r.out.mat -> Matlab or GNU Octave.<br> ML>> array_rot = orig_array' ;
<br> % flip array left-right or up-down<br> ML>> array_flip = fliplr( flipud( orig_array ) ) ;<br> r.in.mat<br><br>for rotation by angle you could do more complicated stuff in Matlab.<br><br><br>In GRASS i.rectify may be an option.
<br><br><br>Hamish<br><br><br><br> ____________________________________________________________________________________<br>Never miss a thing. Make Yahoo your home page.<br><a href="http://www.yahoo.com/r/hs">http://www.yahoo.com/r/hs
</a><br><br></blockquote></div><br>