Hi Uttam,

if I understand your problem correctly, you can use the -p flag of r.proj
to print input map's bounds in the current projection. With this
information you can set the region boundaries appropriately using g.region.

Afterwards you can 'import' your raster into your laea-projection location
using the r.proj as mentioned before.

best,
Johannes

On Wed, Sep 23, 2015 at 8:21 PM, Uttam Sinha <uttamsinh...@gmail.com> wrote:

>
>
> Hi,
>
> I tried the way you told. I created a new location in laea without setting
> the default boundary. Then from this mapset I used
>
> r.proj location*=input_location* mapset=*input_mapset* input=*input_raster
> name* output=
>
> *output_raster*
> I get an error:
>
>
> ERROR: Input raster map is outside current region
>
> Any suggestions.
>
> Uttam.
>
>
> On Tue, Sep 22, 2015 at 11:59 PM, sajid pareeth <spare...@gmail.com>
> wrote:
>
>> Hi Uttam
>>
>> If I understand you correctly you want to reproject a landsat image from
>> sinusoidal to laea projection!!
>> You dont have to manually set the bounding box while reprojecting a
>> raster, as GRASS will take care of it.
>>
>> First you create a new location and mapset with desired output projection
>> , which in this case is Lamberts Azimuthal Equal Area.
>> Then get into this mapset, use the following r.proj command.
>> r.proj location*=input_location* mapset=*input_mapset* input=*input_raster
>> name* output=
>> *output_raster*
>>
>> More details here: https://grass.osgeo.org/grass70/manuals/r.proj.html
>>
>> If you then need to set your working area to a particular bounding box,
>> use g.region
>>
>> g.region raster=*output_raster*
>>
>> regards
>>
>> Sajid
>>
>> On Tue, Sep 22, 2015 at 11:07 PM, Uttam Sinha <uttamsinh...@gmail.com>
>> wrote:
>>
>>>
>>>
>>> Hi All,
>>>
>>> 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.
>>>
>>> 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.
>>>
>>> 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?
>>>
>>> 2.) How do I reproject the images from Sinusoidal to Lamberts without
>>> georeferencing the image at source location (in sinusoidal projection)
>>>
>>> Appreciate any reply.
>>>
>>> Uttam.
>>>
>>> _______________________________________________
>>> grass-user mailing list
>>> grass-user@lists.osgeo.org
>>> http://lists.osgeo.org/mailman/listinfo/grass-user
>>>
>>
>>
>
> _______________________________________________
> grass-user mailing list
> grass-user@lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-user
>
_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Reply via email to