[GRASS5] GRASS and lat/long with sphere

Glynn Clements glynn.clements at virgin.net
Tue Oct 1 09:31:31 EDT 2002


Markus Neteler wrote:

> > The function get_a_e2_f() in get_ellipse.c requires that e2 (the
> > eccentricity-squared) is greater than zero (otherwise it returns zero,
> > indicating an error).
> > 
> > Changing it to allow e2 to be equal to zero doesn't result in any
> > problems which are immediately obvious. However, it may be that some
> > code can't cope with this (e.g. computing the reciprocal). In all
> > probability, the only way to find out is to make the change and wait
> > for bug reports.
> 
> May I ask you to apply the change? As far as I understand it should
> only users of "sphere" ellipsoid. So it is no severe change (I have
> already another tester in mind, hi Alessandro F.!).

Done.

The following programs all contain direct references to functions
whose names contain the words "ellipsoid" or "spheroid", and so are
prime candidates for testing.

etc/bin/cmd:
	d.geodesic
	d.where
	m.datum.shift
	m.gc2ll
	m.ll2gc
	m.ll2u
	m.proj2
	m.region.ll
	m.tiger.region
	m.u2ll
	r.buffer
	r.in.ll
	r.in.utm
	r.surf.idw
	v.apply.census
	v.in.tig.basic
	v.in.tig.lndmk
	v.out.sdts

etc/bin/inter:
	g.setproj
	m.proj

The functions in question are:

CC_get_spheroid            
CC_spheroid_name           
CC_u2ll_spheroid           
CC_u2ll_spheroid_parameters
G_ellipsoid_name           
G_get_ellipsoid_by_name    
G_get_ellipsoid_parameters 
G_get_spheroid_by_name     

-- 
Glynn Clements <glynn.clements at virgin.net>




More information about the grass-dev mailing list