<div dir="ltr">Hi Uttam, <div><br></div><div>you may want to first do v.in.region in your input location, this creates a vector with the limits of your input location (set g.region to your raster of interest first), then v.proj the output of v.in.region in your target location, set region there with the imported vector and just then try to reproject the raster image. </div><div><br></div><div>so, steps would be:</div><div><br></div><div># in input_location</div><div>g.region raster=your_raster</div><div><br></div><div>v.in.region output=vector_region_of_interest</div><div><br></div><div># in target location </div><div>v.proj location=input_location mapset=input_mapset input=vector_region_of_interest output=projected_boundaries</div><div><br></div><div>g.region vector=projected_boundaries</div><div><br></div><div>r.proj location=input_location mapset=input_mapset input=input_raster output=output_raster</div><div><br></div><div>You may also wanna check r.import if you use grass7.1 (<a href="https://grass.osgeo.org/grass71/manuals/r.import.html">https://grass.osgeo.org/grass71/manuals/r.import.html</a>) or as an add-on in grass 7.0 (<a href="https://grass.osgeo.org/grass70/manuals/addons/r.import.html">https://grass.osgeo.org/grass70/manuals/addons/r.import.html</a>). </div><div>You can install it with g.extension r.import</div><div><br></div><div>HTH, </div><div>Vero </div></div><div class="gmail_extra"><br><div class="gmail_quote">2015-09-23 15:21 GMT-03:00 Uttam Sinha <span dir="ltr"><<a href="mailto:uttamsinha09@gmail.com" target="_blank">uttamsinha09@gmail.com</a>></span>:<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><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<span class=""><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></span>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.<span class="HOEnZb"><font color="#888888"><br><br></font></span></div><span class="HOEnZb"><font color="#888888">Uttam.<br><div><div><br></div></div></font></span></div><div class="HOEnZb"><div class="h5"><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>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><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>
</div></div><br>_______________________________________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org">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>