<div dir="ltr"><div><br><br></div>I have been trying for a long time. Something wrong is going on. I appreciate if you can correct me.<br><div><div class="gmail_quote"><br><div dir="ltr"><div><div><div><div><div><div></div></div><br></div>I have all my source images in sinusoidal projection in a location called "Landsat_sinu". <br><br></div>My target locatoin is "Landsat_laea" I tried using the default option in GRASS GUI as follows:<br><br></div>1.) <b>Select coordinate system parameters from a list</b><br><br></div>2.) Projection code: <b>laea</b> <br><div><div><br></div><div>3.) Datum with associated ellipsoid: <br></div><div>Central parallel: <b>30</b><br></div><div>Central Meridian: <b>-96</b><br></div><div>False easting: <b>0</b><br></div><div>False northing: <b>0</b><br><br></div><div>4.) Datum code: <b>nad83</b><br></div><div>5.) Select from the list of datum transformation: <b>1 used in whole nad83 region</b><br><br></div><span class=""><div>+proj=laea<br></div><div>+lat_0=30<br></div><div>+lon_0=-96<br></div><div>+x_0=0<br></div><div>+y_0=0<br></div></span><div>+no_defs<br></div><div>+a=6378137<br></div><div>+rf=298.257222101<br></div><div>+towgs84=0.000,0.000,0.000<br></div><div>to_meter=1<br></div><div>+datum=nad83<br><br></div><div>Finish.<br><br></div><div>Is everything ok with this for the below (required) target definition?<br></div><span class=""><div><br></div><div><b>*PROJ.4 : '+proj=laea +lat_0=30 +lon_0=-96 +x_0=0 +y_0=0 +datum=NAD83<br>
+units=m +no_defs '*</b><br><br></div></span><div>Looks like at some point I am making a small mistake. I must get this as the projected coordinates.<br><br><b>
<span style="color:rgb(0,0,255)">north=1064601.38274499<br>
south=830176.771464028<br>
east=-1960125.15283348<br>
west=-2265106.68066176</span></b><br><br><br></div><div>Any help is appreciated.<br></div><div><br></div><div>Uttam.<br></div><div><br></div></div></div><div class=""><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Oct 14, 2015 at 4:52 PM, Moritz Lennert <span dir="ltr"><<a href="mailto:mlennert@club.worldonline.be" target="_blank">mlennert@club.worldonline.be</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">[Please keep conversations on the list]<br>
<br>
On 14/10/15 23:19, Uttam Kumar wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
_<br>
<br>
_<br>
_*This information is for the original data (Source):*_<br>
<br>
*PROJ.4 : '+proj=sinu +lon_0=0 +x_0=0 +y_0=0 +a=6371007.181<br>
+b=6371007.181 +units=m +no_defs '*<span><br>
<br>
<br>
Coordinate System is:<br>
PROJCS["Sinusoidal",<br>
     GEOGCS["GCS_Unknown",<br>
         DATUM["Unknown",<br>
             SPHEROID["S_Unknown",6371007.181,0.0]],<br>
         PRIMEM["Greenwich",0.0],<br>
         UNIT["Degree",0.0174532925199433]],<br>
     PROJECTION["Sinusoidal"],<br>
     PARAMETER["False_Easting",0.0],<br>
     PARAMETER["False_Northing",0.0],<br>
     PARAMETER["longitude_of_center",0.0],<br>
     UNIT["Meter",1]]<br>
Origin = (-10801805.197699999436736,4130102.079069999977946)<br>
Pixel Size = (30.000000000000000,-30.000000000000000)<br>
<br></span>
*<br>
_This information is for the destination data (Target):_*<br>
<br>
*PROJ.4 : '+proj=laea +lat_0=30 +lon_0=-96 +x_0=0 +y_0=0 +datum=NAD83<br>
+units=m +no_defs '*<br>
</blockquote><span>
<br>
In your first mail you said:<br>
<br>
Lambert Azimuthal Equal Area, Datum: WGS84<br>
<br></span>
but here the datum is NAD83.<br>
<br>
Using your projection info, I get almost exactly the same results in GRASS as you report from ENVI:<br>
<br>
On 08/10/15 21:36, Uttam Kumar wrote:> Reprojected in ENVI (Lambert Azimuthal Equal Area, Datum: WGS84)<span><br>
> -------------------------------------------------------------------------<br>
> Corner Coordinates:<br>
><br>
> Upper Left  (-2265106.681W, 1064601.383N)<br>
> Lower Left  (-2265106.681W,  830181.383N)<br>
> Upper Right (-1960126.681W, 1064601.383N)<br>
> Lower Right (-1960126.681W,  830181.383N)<br>
<br>
<br></span>
GRASS 7.0.1 (test_LAEA):~ > <a href="http://v.info" rel="noreferrer" target="_blank">v.info</a> -g Image_Boundary_box<br>
north=1064601.38274499<br>
south=830176.771464028<br>
east=-1960125.15283348<br>
west=-2265106.68066176<br>
<br>
<br>
Don't have the time now to explore the exact differences, but either it's a question of definition of the bounding box, or it possibly be a question of which specific NAD83 transformation parameters were used.<span><font color="#888888"><br>
<br>
Moritz<br>
</font></span></blockquote></div><br></div>
</div></div></div><br></div></div>