[Proj] Proj may give erroneous results depending on locale

Peter S Galbraith psg at debian.org
Tue Jul 20 10:55:42 PDT 2004


Hi, a Debian user (Radoslaw Grzanka <radekg2 at poczta.onet.pl>) reported
the following bug when his locale is set to pl_PL :

: Hello,
:   I have encountered bug which resulted in proj giving different results
: depending on locale set.
:   As you can see in the signature my locale is pl_PL and when used with
: gettext (or gtk+) gives wrong coordinates after transformation.
:   I'm attaching simple testcase that shows the problem.
:   The bug is also occures in testing vesion of package (4.4.7) and was
: really hard to trace.
: 
: Best Regards,
:   Radoslaw Grzanka

When compiling the following C++ code and linking with either 4.4.7 or
4.4.8 versions of proj, the bug appears when running the executable (on
Debian systems).  I'm unfortunately not an expert in proj so cannot
track this down.  I don't think this is Debian-specific (but I could be
wrong).

--------cut-------------------------
//// Simple testcase
#include <locale.h>
#include <string>
#include <iostream>
#include <projects.h>

const std::string projparams("+x_0=115073.415747 +y_0=-6288530.367751"
        " +proj=lcc +lat_1=37.000000 +lat_2=65.000000 +lon_0=8.199763"
        " +ellps=WGS84 +to_meter=63.506519 +no_defs");

projUV proj_input;              // Input Values
projUV proj_output;             // Output Values
projPJ pj;

int main()
{
//////////////////////////////////
//// TEST 1 - without locales set.
        pj = pj_init_plus(projparams.c_str());

        proj_input.v = 53.5736;         // Lattitude
        proj_input.u = 7.42731;         // Longitude
        proj_input.v *= DEG_TO_RAD;     // To Rads
        proj_input.u *= DEG_TO_RAD;     // To Rads

        proj_output = pj_fwd(proj_input, pj);

        std::cout << "Lattitude: " << proj_output.v
                << " Longitude: " << proj_output.u << std::endl;

        pj_free(pj);
//// RESULT - Lattitude: 3197.17 Longitude: 1029.99 - GOOD.
///////////////////////////////////
////     TEST 2 - with locales set!
        setlocale (LC_ALL, "");                 // SETTING LOCALES!

        pj = pj_init_plus(projparams.c_str());

        proj_input.v = 53.5736;         // Lattitude
        proj_input.u = 7.42731;         // Longitude
        proj_input.v *= DEG_TO_RAD;     // To Rads
        proj_input.u *= DEG_TO_RAD;     // To Rads

        proj_output = pj_fwd(proj_input, pj);

        std::cout << "Lattitude: " << proj_output.v
                << " Longitude: " << proj_output.u << std::endl;

        pj_free(pj);
//// RESULT - Lattitude: 3220.56 Longitude: 1242.12 - BAD.
}
--------cut-------------------------

I had to add pl_PL as a possible locale on my system and generate the
proper locale files, then set pl_PL as my locale, then I could reproduce
the bug.

Thanks,
-- 
Peter S. Galbraith, Debian Developer          <psg at debian.org>
                                 http://people.debian.org/~psg
GPG key 1024/D2A913A1 - 97CE 866F F579 96EE  6E68 8170 35FF 799E



More information about the Proj mailing list