[GRASS-user] r.out.gdal outputs incorrect projection file

Hello,

I have used r.out.gdal (GRASS 6.0.1) to import a number of geotiffs into
a GRASS location that is GCS (e.g. Projection 3 - Latitude, Longitude),
with the WGS84 datum. Both the location, and the geotiffs, are in that
coordinate system so that I can easily import without having to use the
-o argument. When I export the files as an AAIGrid, Uint16 ascii file, a
projection file is generated with the following information:

PROJCS["Sinusoidal(Sanson-Flamsteed)",GEOGCS["Unknown",DATUM["D_unknown",SPHEROID["Unknown",6378137,"inf"]],
PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Sinusoidal"],PARAMETER["longitude_of_center",0],
PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["meters",1]]

Why is it showing the projection as Sinusoidal?

Thanks, William

Bill Hudspeth wrote:

I have used r.out.gdal (GRASS 6.0.1) to import a number of geotiffs into
a GRASS location that is GCS (e.g. Projection 3 - Latitude, Longitude),
with the WGS84 datum. Both the location, and the geotiffs, are in that
coordinate system so that I can easily import without having to use the
-o argument. When I export the files as an AAIGrid, Uint16 ascii file, a
projection file is generated with the following information:

PROJCS["Sinusoidal(Sanson-Flamsteed)",GEOGCS["Unknown",DATUM["D_unknown",SPHEROID["Unknown",6378137,"inf"]],
PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Sinusoidal"],PARAMETER["longitude_of_center",0],
PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["meters",1]]

Why is it showing the projection as Sinusoidal?
Thanks, William

William,

(although some projection-expert knows this stuff better, two questions
to check again from the start...)

a. how did you create the location (exact command if possible)?

b. what do "g.region -p" and "g.proj -p" report for (=from within) the
location used to export (using r.out.gdal) the data?

# "b" should be somthing like in [*]

Nikos
---

[*]
# e.g. region settings
g.region -p

projection: 3 (Latitude-Longitude)
zone: 0
datum: wgs84
ellipsoid: wgs84
north: 36:27:38N
south: 35:52:34N
west: 27:30:58E
east: 28:14:55E
nsres: 0:00:01
ewres: 0:00:01
rows: 2104
cols: 2637
cells: 5548248

# e.g. projection settings
g.proj -p
-PROJ_INFO-------------------------------------------------
name : Lat/Lon
proj : ll
datum : wgs84
ellps : wgs84
no_defs : defined
-PROJ_UNITS------------------------------------------------
unit : degree
units : degrees
meters : 1.0

On Tue, Aug 11, 2009 at 1:09 AM, Nikos
Alexandris<nikos.alexandris@felis.uni-freiburg.de> wrote:

Bill Hudspeth wrote:

I have used r.out.gdal (GRASS 6.0.1) to import a number of geotiffs into

Note:
In 6.0.1, r.out.gdal is a script is based on 'gdal_translate'
(in 6.4.0 it is a C program).

...

PROJCS["Sinusoidal(Sanson-Flamsteed)",GEOGCS["Unknown",DATUM["D_unknown",SPHEROID["Unknown",6378137,"inf"]],
PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Sinusoidal"],PARAMETER["longitude_of_center",0],
PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["meters",1]]

Why is it showing the projection as Sinusoidal?

...

b. what do "g.region -p" and "g.proj -p" report for (=from within) the
location used to export (using r.out.gdal) the data?

Please also try
g.proj -w

which is the above WKT style.

Markus