I have a data sets with the following characteristics:
Numer of rows = 8384
Number of columns = 9102
XY corner coordinates :
Lower left: -4462500.000, -3999500.000
Upper left: -4462500.000, 4384500.000
Upper right: 4639500.000, 4384500.000
Lower right: 4639500.000, -3999500.000
Projection used: Lambert azimutal equal area
units = meters
pixel size = 1000 mts
Radius of sphere of influence = 6,370,997, mts
Longitud of origen = 100 00 00W
Latitud of origen = 45 00 00N
False easting = 0.0
False northing = 0.0
I can get it with “r.in.bin …” in “location 1”, but i need it in lat-long projection and when i try to get it in “location 2” using the following parameters into the “PROJ_INFO” file ;
type of projection= ll (lat-lon), datum= wgs84, gridnads= conus, ellips= wgs84.
With g.region n=90 s=-90 w=-180 e=180 res=0.008333 (30 secs) and using the “r.proj …” comand , i obtain my transformation with the following coordinates of the raster new image; n=40 s=-36 w=-49 e=51res=0.008333, but it’s wrong because the coordinates must be aproximately n=80 s=10 w=-170 e=-40.
How can i solve my problem?, which are the correct parameters in the “PROJ_INFO” file?, i try with ellips = sphere and can’t get results, where i must to definite lat, lon of origen and false northing and easting?
The data sets are int the urls:
http://lpdaac.usgs.gov/gtopo30/hydro/namerica.asp
and can be of the type : Download HYDRO1k DEM for North America as a tar file (28.2 MB).
I wait that anyone can help me in this.
Miguel
Do You Yahoo!?
Todo lo que quieres saber de Estados Unidos, América Latina y el resto del Mundo.
Visíta Yahoo! Noticias.