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)
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
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)
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 <spareeth@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
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)
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.
On Tue, Sep 22, 2015 at 11:59 PM, sajid pareeth <spareeth@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
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)
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.
On Tue, Sep 22, 2015 at 11:59 PM, sajid pareeth <spareeth@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
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)
However I noticed that the projected boundary values of N, S, W and E in GRASS are not same as the values obtained from projecting the same image in ENVI image processing software and in gdalwarp.
Projection of the scene in Lambert Azimuthal Equal Area is similar (but not exactly equal) in GRASS and gdal, but in ENVI it is totally different. I just did the comparison to see if things are similar from different softwares in case they have to be ported to other platforms.
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.
On Tue, Sep 22, 2015 at 11:59 PM, sajid pareeth <spareeth@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
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)
However I noticed that the projected boundary values of N, S, W and E in GRASS are not same as the values obtained from projecting the same image in ENVI image processing software and in gdalwarp.
Projection of the scene in Lambert Azimuthal Equal Area is similar (but not exactly equal) in GRASS and gdal, but in ENVI it is totally different. I just did the comparison to see if things are similar from different softwares in case they have to be ported to other platforms.
Can you print the projection information in both GRASS (g.proj -p or -w) and ENVI? Are they really the same?