1 |
# Copyright (c) 2002, 2003 by Intevation GmbH |
# Copyright (c) 2002, 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. |
15 |
# $Id$ |
# $Id$ |
16 |
|
|
17 |
import unittest |
import unittest |
18 |
|
import locale |
19 |
import os |
import os |
20 |
|
|
21 |
|
import localessupport |
22 |
import xmlsupport |
import xmlsupport |
23 |
import support |
import support |
24 |
support.initthuban() |
support.initthuban() |
58 |
(3620891.3077618643, 0.0, |
(3620891.3077618643, 0.0, |
59 |
3875381.8535437919, 252962.10480170773), |
3875381.8535437919, 252962.10480170773), |
60 |
epsilon = 1e-5) |
epsilon = 1e-5) |
61 |
|
self.assertFloatSeqEqual(proj.InverseBBox((3620891.3077618643, 0.0, |
62 |
|
3875381.8535437919, |
63 |
|
252962.10480170773)), |
64 |
|
(-0.018341599754143501, 0.0, |
65 |
|
2.017992992681688, 2.0377390677846736), |
66 |
|
epsilon = 1e-5) |
67 |
|
|
68 |
# GetName() |
# GetName() |
69 |
self.assertEquals(proj.GetName(), _("Unknown")) |
self.assertEquals(proj.GetName(), _("Unknown")) |
89 |
self.assertEquals(proj.GetParameter("south"), "south") |
self.assertEquals(proj.GetParameter("south"), "south") |
90 |
|
|
91 |
def test_get_projection_units_geo(self): |
def test_get_projection_units_geo(self): |
92 |
"""Test Projection.GetProjectedUnits() for geographic projection""" |
"""Test Projection.GetProjectedUnits() for geographic projection. |
93 |
|
Test for the alias 'longlat' as well. |
94 |
|
""" |
95 |
proj = Projection(["proj=latlong", "to_meter=0.017453292519943295", |
proj = Projection(["proj=latlong", "to_meter=0.017453292519943295", |
96 |
"ellps=clrk66"]) |
"ellps=clrk66"]) |
97 |
self.assertEquals(proj.GetProjectedUnits(), PROJ_UNITS_DEGREES) |
self.assertEquals(proj.GetProjectedUnits(), PROJ_UNITS_DEGREES) |
98 |
|
proj = Projection(["proj=longlat", "to_meter=0.017453292519943295", |
99 |
|
"ellps=clrk66"]) |
100 |
|
self.assertEquals(proj.GetProjectedUnits(), PROJ_UNITS_DEGREES) |
101 |
|
|
102 |
|
def test_lc_numeric_robustness(self): |
103 |
|
"""Test if an LC_NUMERIC local with comma as decimal_point will work. |
104 |
|
|
105 |
|
Some versions of proj are not robust against this. |
106 |
|
Starting with Python 2.4 there is a different behaviour when |
107 |
|
calling C extensions and now they will see changes locales settings |
108 |
|
which might tickle the bug in proj. |
109 |
|
""" |
110 |
|
params = ["proj=latlong", "to_meter=0.01745", "ellps=clrk66"] |
111 |
|
|
112 |
|
oldlocale = localessupport.setdecimalcommalocale() |
113 |
|
if oldlocale == None: |
114 |
|
raise support.SkipTest( |
115 |
|
"No locale with comma as decimal_point found.") |
116 |
|
|
117 |
|
proj = Projection(params) |
118 |
|
#print proj.work_around_broken_proj |
119 |
|
result1 = proj.Forward(1.2,3.2) |
120 |
|
|
121 |
|
locale.setlocale(locale.LC_NUMERIC, "C") |
122 |
|
proj = Projection(params) |
123 |
|
result2= proj.Forward(1.2,3.2) |
124 |
|
|
125 |
|
locale.setlocale(locale.LC_NUMERIC, oldlocale) |
126 |
|
self.assertFloatSeqEqual(result1, result2, epsilon = 1e-5 ) |
127 |
|
|
128 |
def test_get_projection_units_normal(self): |
def test_get_projection_units_normal(self): |
129 |
"""Test Projection.GetProjectedUnits() for normal projection""" |
"""Test Projection.GetProjectedUnits() for normal projection""" |
173 |
|
|
174 |
class TestProjFile(unittest.TestCase, support.SubscriberMixin): |
class TestProjFile(unittest.TestCase, support.SubscriberMixin): |
175 |
|
|
176 |
"""Test cases for reading and writing projection files. |
"""Test cases for ProjFile objects""" |
|
""" |
|
177 |
|
|
178 |
def setUp(self): |
def setUp(self): |
179 |
self.clear_messages() |
self.clear_messages() |
247 |
|
|
248 |
class ProjFileReadTests(ProjFileTest): |
class ProjFileReadTests(ProjFileTest): |
249 |
|
|
250 |
"""Test read ProjFile objects from files""" |
"""Test read ProjFile objects from files |
251 |
|
|
252 |
|
The files only cover error handling and the system projection file. |
253 |
|
""" |
254 |
|
|
255 |
def test_read_non_existing_file(self): |
def test_read_non_existing_file(self): |
256 |
"""Test read_proj_file with non-existing file""" |
"""Test read_proj_file with non-existing file""" |
264 |
As currently written this only works on unix-like systems and |
As currently written this only works on unix-like systems and |
265 |
not e.g. on MS Windows. |
not e.g. on MS Windows. |
266 |
""" |
""" |
267 |
|
if os.name != "posix": |
268 |
|
raise support.SkipTest("Test only works on posix systems") |
269 |
filename = self.filename() |
filename = self.filename() |
270 |
file = open(filename, "w") |
file = open(filename, "w") |
271 |
file.close() |
file.close() |
281 |
self.assertRaises(SAXParseException, resource.read_proj_file, filename) |
self.assertRaises(SAXParseException, resource.read_proj_file, filename) |
282 |
|
|
283 |
def test_get_system_proj_file(self): |
def test_get_system_proj_file(self): |
284 |
"""Test resource.get_system_proj_file() |
"""Test resource.get_system_proj_file(DEFAULT_PROJ_FILE) |
285 |
|
|
286 |
This is primarily to test whether the system proj file contains |
This is primarily to test whether the system proj file contains |
287 |
invalid projection paramers and whether the proj file is not |
invalid projection paramers and whether the proj file is not |
288 |
empty |
empty |
289 |
""" |
""" |
290 |
projfile, warnings = resource.get_system_proj_file() |
projfile, warnings\ |
291 |
|
= resource.get_system_proj_file(resource.DEFAULT_PROJ_FILE) |
292 |
self.assertEquals(warnings, []) |
self.assertEquals(warnings, []) |
293 |
self.assert_(len(projfile.GetProjections()) > 0) |
self.assert_(len(projfile.GetProjections()) > 0) |
294 |
|
|
295 |
|
# see whether it got cached and we get the same projfile object |
296 |
|
# when we read the file again |
297 |
|
projfile2, warnings \ |
298 |
|
= resource.get_system_proj_file(resource.DEFAULT_PROJ_FILE) |
299 |
|
self.assert_(projfile is projfile2) |
300 |
|
|
301 |
|
|
302 |
class WriteProjFileTests(ProjFileTest): |
class WriteProjFileTests(ProjFileTest): |
303 |
|
|
366 |
self.doTestWrite(pf, file_contents) |
self.doTestWrite(pf, file_contents) |
367 |
|
|
368 |
|
|
369 |
class TestLoadingProjFile(support.FileLoadTestCase): |
class ProjFileLoadTestCase(support.FileLoadTestCase): |
370 |
|
|
371 |
|
"""Base class for the test cases that read specific test files""" |
372 |
|
|
373 |
file_extension = ".proj" |
file_extension = ".proj" |
374 |
|
|
375 |
|
def tearDown(self): |
376 |
|
"""Clear the cache explicitly""" |
377 |
|
resource.clear_proj_file_cache() |
378 |
|
|
379 |
|
|
380 |
|
class TestLoadingProjFile(ProjFileLoadTestCase): |
381 |
|
|
382 |
file_contents = '''\ |
file_contents = '''\ |
383 |
<?xml version="1.0" encoding="UTF-8"?> |
<?xml version="1.0" encoding="UTF-8"?> |
384 |
<!DOCTYPE projectionlist SYSTEM "projfile.dtd"> |
<!DOCTYPE projectionlist SYSTEM "projfile.dtd"> |
430 |
"proj=tmerc", "units=m", |
"proj=tmerc", "units=m", |
431 |
"x_0=400000.000", "y_0=0.000"]) |
"x_0=400000.000", "y_0=0.000"]) |
432 |
|
|
433 |
|
def test_caching(self): |
434 |
|
# test whether the projfile cache works |
435 |
|
projfile, warnings = resource.read_proj_file(self.filename()) |
436 |
|
projfile2, warnings = resource.read_proj_file(self.filename()) |
437 |
|
|
438 |
|
# Both projfiles should be the same object |
439 |
|
self.assert_(projfile2 is projfile) |
440 |
|
|
441 |
class TestLoadingProjFileWithEmptyProjectionlist(support.FileLoadTestCase): |
# If we clear the cache we should get a new one. |
442 |
|
resource.clear_proj_file_cache() |
443 |
|
projfile3, warnings = resource.read_proj_file(self.filename()) |
444 |
|
self.assert_(projfile3 is not projfile) |
445 |
|
|
446 |
|
|
447 |
|
class TestLoadingProjFileWithEmptyProjectionlist(ProjFileLoadTestCase): |
448 |
|
|
|
file_extension = ".proj" |
|
449 |
file_contents = '''\ |
file_contents = '''\ |
450 |
<?xml version="1.0" encoding="UTF-8"?> |
<?xml version="1.0" encoding="UTF-8"?> |
451 |
<!DOCTYPE projectionlist SYSTEM "projfile.dtd"> |
<!DOCTYPE projectionlist SYSTEM "projfile.dtd"> |
462 |
self.assertEquals(len(projfile.GetProjections()), 0) |
self.assertEquals(len(projfile.GetProjections()), 0) |
463 |
|
|
464 |
|
|
465 |
class TestProjFileWithInvalidParameters(unittest.TestCase, |
class TestProjFileWithInvalidParameters(ProjFileLoadTestCase): |
|
support.FileLoadTestCase): |
|
466 |
|
|
|
file_extension = ".proj" |
|
467 |
file_contents = '''\ |
file_contents = '''\ |
468 |
<?xml version="1.0" encoding="UTF-8"?> |
<?xml version="1.0" encoding="UTF-8"?> |
469 |
<!DOCTYPE projectionlist SYSTEM "projfile.dtd"> |
<!DOCTYPE projectionlist SYSTEM "projfile.dtd"> |