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

[Nikos Alexandris]

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