1 |
# Copyright (c) 2001, 2003 by Intevation GmbH |
# Copyright (c) 2001, 2003, 2006 by Intevation GmbH |
2 |
# Authors: |
# Authors: |
3 |
# Bernhard Herzog <[email protected]> |
# Bernhard Herzog <[email protected]> |
4 |
|
# Bernhard Reiter <[email protected]> |
5 |
# |
# |
6 |
# This program is free software under the GPL (>=v2) |
# This program is free software under the GPL (>=v2) |
7 |
# Read the file COPYING coming with Thuban for details. |
# Read the file COPYING coming with Thuban for details. |
|
|
|
8 |
__version__ = "$Revision$" |
__version__ = "$Revision$" |
9 |
|
|
10 |
from types import StringTypes |
from types import StringTypes |
11 |
|
import locale |
12 |
|
|
13 |
from Thuban import _ |
from Thuban import _ |
14 |
from Thuban.Lib.connector import Publisher |
from Thuban.Lib.connector import Publisher |
22 |
PROJ_UNITS_METERS = 1 |
PROJ_UNITS_METERS = 1 |
23 |
PROJ_UNITS_DEGREES = 2 |
PROJ_UNITS_DEGREES = 2 |
24 |
|
|
25 |
|
def _do_we_have_to_work_around_broken_proj(): |
26 |
|
""" If we have a problematic locale, check if proj results are good. """ |
27 |
|
if locale.localeconv()['decimal_point'] != '.': |
28 |
|
params = ["proj=latlong", "to_meter=0.01745", "ellps=clrk66"] |
29 |
|
proj = BaseProjection(params) |
30 |
|
result1 = proj.Forward(1,1) |
31 |
|
|
32 |
|
savedlocale = locale.getlocale(locale.LC_NUMERIC) |
33 |
|
locale.setlocale(locale.LC_NUMERIC, "C") |
34 |
|
|
35 |
|
proj = BaseProjection(params) |
36 |
|
result2 = proj.Forward(1,1) |
37 |
|
|
38 |
|
locale.setlocale(locale.LC_NUMERIC, savedlocale) |
39 |
|
if result1 != result2: |
40 |
|
return True |
41 |
|
return False |
42 |
|
|
43 |
class Projection(BaseProjection): |
class Projection(BaseProjection): |
44 |
|
"""A proj4 projection object that remembers the parameters. |
45 |
|
|
46 |
"""A proj4 projection object that remembers the parameters""" |
The proj library is not robust against decimal_point != '.' locales. |
47 |
|
Since python 2.4 calls C extensions with the set locale, it can create |
48 |
|
a problem. It seems that calling |
49 |
|
self.assuregoodlocale() |
50 |
|
self.assureinitlocale() |
51 |
|
before BaseProjection.__init__() is enough to work around this. |
52 |
|
|
53 |
|
We assuming that the locale stays the same after a projection |
54 |
|
has been initialised |
55 |
|
and thus we can return to it in self.assureinitlocale(). |
56 |
|
""" |
57 |
|
|
58 |
def __init__(self, params, name = None, epsg = None): |
def __init__(self, params, name = None, epsg = None): |
59 |
"""Initialize the Projection |
"""Initialize the Projection |
62 |
|
|
63 |
params -- a list of 'parameter=value' strings |
params -- a list of 'parameter=value' strings |
64 |
|
|
65 |
name -- (optional) The name of the projectin. If None or omitted |
name -- (optional) The name of the projection. If None or omitted |
66 |
it defaults to 'Unknown' in the local language. |
it defaults to 'Unknown' in the local language. |
67 |
|
|
68 |
epsg -- (optional) The EPSG code as a string. |
epsg -- (optional) The EPSG code as a string. |
69 |
""" |
""" |
70 |
|
self.initlocale = locale.getlocale(locale.LC_NUMERIC) |
71 |
|
self.work_around_broken_proj = _do_we_have_to_work_around_broken_proj() |
72 |
|
|
73 |
|
self.assuregoodlocale() |
74 |
BaseProjection.__init__(self, params) |
BaseProjection.__init__(self, params) |
75 |
|
self.assureinitlocale() |
76 |
|
|
77 |
if name is None: |
if name is None: |
78 |
self.name = _("Unknown") |
self.name = _("Unknown") |
82 |
self.epsg = epsg |
self.epsg = epsg |
83 |
self.params = params |
self.params = params |
84 |
|
|
85 |
|
def assuregoodlocale(self): |
86 |
|
if self.work_around_broken_proj: |
87 |
|
locale.setlocale(locale.LC_NUMERIC, "C") |
88 |
|
|
89 |
|
def assureinitlocale(self): |
90 |
|
if self.work_around_broken_proj: |
91 |
|
locale.setlocale(locale.LC_NUMERIC, self.initlocale) |
92 |
|
|
93 |
def _transform_bbox(self, trafo, bbox): |
def _transform_bbox(self, trafo, bbox): |
94 |
# This is not really the correct way to determine the bbox of a |
# This is not really the correct way to determine the bbox of a |
95 |
# projected bbox, but for now it works well enough |
# projected bbox, but for now it works well enough |