sij/mapnik
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

convert to degrees if projection is geographic

Browse source
This commit is contained in:
Artem Pavlenko 2008-02-20 12:44:34 +00:00
parent 24c605e807
commit f3a3ce04db
1 changed files with 16 additions and 6 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

22
src/projection.cpp
View file

@ -82,6 +82,11 @@ namespace mapnik {
p = pj_fwd(p,proj_);
x = p.u;
y = p.v;
if (pj_is_latlong(proj_))
{
x *=RAD_TO_DEG;
y *=RAD_TO_DEG;
}
}
void projection::inverse(double & x,double & y) const
@ -89,12 +94,17 @@ namespace mapnik {
#ifdef MAPNIK_THREADSAFE
mutex::scoped_lock lock(mutex_);
#endif
projUV p;
p.u = x;
p.v = y;
p = pj_inv(p,proj_);
x = RAD_TO_DEG * p.u;
y = RAD_TO_DEG * p.v;
if (pj_is_latlong(proj_))
{
x *=DEG_TO_RAD;
y *=DEG_TO_RAD;
}
projUV p;
p.u = x;
p.v = y;
p = pj_inv(p,proj_);
x = RAD_TO_DEG * p.u;
y = RAD_TO_DEG * p.v;
}
projection::~projection()