<div dir="ltr"><div><div><div><div><br><br></div>Hi,<br><br></div>I tried the way you told. I created a new location in laea without setting the default boundary. Then from this mapset I used<br><br><span style="font-family:monospace,monospace">r.proj location<em>=input_location</em> mapset=<em>input_mapset</em> input=<em>input_raster name</em> output=</span><em><span style="font-family:monospace,monospace">output_raster<br><br></span></em><br>I get an error:<br><br><br><span style="color:rgb(0,0,255)">ERROR: Input raster map is outside current region</span><br><br></div>Any suggestions.<br><br></div>Uttam.<br><div><div><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Sep 22, 2015 at 11:59 PM, sajid pareeth <span dir="ltr"><<a href="mailto:spareeth@gmail.com" target="_blank">spareeth@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><div>Hi Uttam<br><br></div>If I understand you correctly you want to reproject a landsat image from sinusoidal to laea projection!! <br></div>You dont have to manually set the bounding box while reprojecting a raster, as GRASS will take care of it.<br><br></div><div>First you create a new location and mapset with desired output projection , which in this case is Lamberts Azimuthal Equal Area.<br></div><div>Then get into this mapset, use the following r.proj command.<br><span style="font-family:monospace,monospace">r.proj location<em>=input_location</em> mapset=<em>input_mapset</em> input=<em>input_raster name</em> output=</span><em><span style="font-family:monospace,monospace">output_raster</span><br></em></div><div><em><br></em></div><div>More details here: <a href="https://grass.osgeo.org/grass70/manuals/r.proj.html" target="_blank">https://grass.osgeo.org/grass70/manuals/r.proj.html</a><br><em><br></em></div><div>If you then need to set your working area to a particular bounding box, use g.region<br><br></div><div><span style="font-family:monospace,monospace">g.region raster=<em>output_raster</em></span></div><br></div>regards<br><br></div>Sajid <br></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5">On Tue, Sep 22, 2015 at 11:07 PM, Uttam Sinha <span dir="ltr"><<a href="mailto:uttamsinha09@gmail.com" target="_blank">uttamsinha09@gmail.com</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div dir="ltr"><div><div><div><div><div><div><br><br></div>Hi All,<br><br></div>I have a set of Landsat images in Sinusoidal Projection, WGS84 Datum. I know the north, south, east and west extent of the boundary of the image.<br><br>I want to reproject the image to Lamberts Azimuthal Equal Area, WGS84 Datum. I do not know what will be the coordinates of the north, south, east and west bounds in this projection.<br><br></div>1.) How do I convert upper left coordinates and bottom right coordinates of Sinusoidal Projection to Lamberts Azimuthal Equal Area to create a new destination Location?<br><br></div>2.) How do I reproject the images from Sinusoidal to Lamberts without georeferencing the image at source location (in sinusoidal projection)<br><br></div>Appreciate any reply.<span><font color="#888888"><br><br></font></span></div><span><font color="#888888">Uttam.<br></font></span></div>
<br></div></div>_______________________________________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/grass-user" rel="noreferrer" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br></blockquote></div><br></div>
</blockquote></div><br></div>