[GRASS-user] gps import: definition of WGS84

On 5/23/08 10:38 AM, "grass-user-request@lists.osgeo.org"
<grass-user-request@lists.osgeo.org> wrote:

Date: Fri, 23 May 2008 19:12:39 +0200
From: Moritz Lennert <mlennert@club.worldonline.be>
Subject: Re: [GRASS-user] gps import: definition of WGS84
To: Nikos Alexandris <nikos.alexandris@felis.uni-freiburg.de>
Cc: grass-user@lists.osgeo.org
Message-ID: <4836FB07.6060402@club.worldonline.be>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed

On 23/05/08 18:45, Nikos Alexandris wrote:

On Fri, 2008-05-23 at 18:39 +0200, Moritz Lennert wrote:

Hello,

I'm trying to figure out why we are seeing different results when
importing GPS data into different software (GRASS is doing great, it's
the others that are having the trouble). It must be linked to the
definition of the projection the GPS data comes in.

v.in.gpsbabel defines the default projection as

IN_PROJ="+proj=longlat +towgs84=0.000,0.000,0.000"

Now my question: is this equivalent to EPSG 4326, which is defined as

<4326> +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs <>

?

Dear Moritz,

I think the first definition is not complete. If the "+datum=WGS84"
parameter was at least present it should be ok since a datum includes
(or should include) all required information.

The "problem" is that it's the first that works.

To be more precise:

I have data in a GPS.

I have a raster topographic map of the area where the GPS data was
taken. This map is in ESPG 31370.

I create my EPSG 31370 location and use, import my raster and use
v.in.gpsbabel to import the gps data, not setting the proj= parameter.

The gps points and tracks are perfectly positioned compared to the topo map.

However, when I try to use QGIS (or gvSIG) with on-the-fly projection,

My bet is that this is the problem. We've had regular requests to implement
on-the-fly projection in GRASS but have not done it so far. It is
convenient, but tends to have errors that don't show up until you actually
try to do something that requires precise rectification and overlays. This
means that, unless you have something to compare it with, you can have
errors and not know it.

Michael

defining the projection of my tracks/points as EPSG 4326 (and I've tried
a few others), or as a custom projection defined as

+proj=longlat +towgs84=0.000,0.000,0.000 +ellps=WGS84

my topo map projection as 31370, and my project/view projection as
31370, I always see my GPS data several tens of meters off...

Moritz

__________________________________________
Michael Barton, Professor of Anthropology
Director of Graduate Studies
School of Human Evolution & Social Change
Center for Social Dynamics and Complexity
Arizona State University

phone: 480-965-6213
fax: 480-965-7671
www: http://www.public.asu.edu/~cmbarton

Hi all,

In settings -> options -> projection, try to set "prompt for projection" and
select the desired projection (define a custum one, if you prefer) every
time you import points (or whatever) using gpsbabel. Probably you set to use
a projection (different from the one defined with gpsbabel) as standard
projection for every "import" in the current project. Actually I'm not sure
about it, but just try, it could work :slight_smile:

Marco