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