/[thuban]/branches/WIP-pyshapelib-bramz/test/test_proj.py
ViewVC logotype

Diff of /branches/WIP-pyshapelib-bramz/test/test_proj.py

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1818 by bh, Mon Oct 13 15:54:14 2003 UTC revision 2696 by bernhard, Sun Sep 17 23:19:21 2006 UTC
# Line 14  __version__ = "$Revision$" Line 14  __version__ = "$Revision$"
14  # $Id$  # $Id$
15    
16  import unittest  import unittest
17    import locale
18  import os  import os
19    
20    import localessupport
21  import xmlsupport  import xmlsupport
22  import support  import support
23  support.initthuban()  support.initthuban()
# Line 23  support.initthuban() Line 25  support.initthuban()
25  from Thuban import _  from Thuban import _
26  from Thuban.Model.proj import Projection, ProjFile, \  from Thuban.Model.proj import Projection, ProjFile, \
27       PROJ_UNITS_METERS, PROJ_UNITS_DEGREES       PROJ_UNITS_METERS, PROJ_UNITS_DEGREES
28    from Thuban.Model.messages import PROJECTION_ADDED, PROJECTION_REMOVED, \
29         PROJECTION_REPLACED
30  import Thuban.Model.resource as resource  import Thuban.Model.resource as resource
31    
32  from xmlsupport import sax_eventlist  from xmlsupport import sax_eventlist
# Line 54  class TestProjection(unittest.TestCase, Line 57  class TestProjection(unittest.TestCase,
57                                   (3620891.3077618643, 0.0,                                   (3620891.3077618643, 0.0,
58                                    3875381.8535437919, 252962.10480170773),                                    3875381.8535437919, 252962.10480170773),
59                                   epsilon = 1e-5)                                   epsilon = 1e-5)
60            self.assertFloatSeqEqual(proj.InverseBBox((3620891.3077618643, 0.0,
61                                                       3875381.8535437919,
62                                                       252962.10480170773)),
63                                     (-0.018341599754143501, 0.0,
64                                      2.017992992681688, 2.0377390677846736),
65                                     epsilon = 1e-5)
66    
67          # GetName()          # GetName()
68          self.assertEquals(proj.GetName(), _("Unknown"))          self.assertEquals(proj.GetName(), _("Unknown"))
# Line 79  class TestProjection(unittest.TestCase, Line 88  class TestProjection(unittest.TestCase,
88          self.assertEquals(proj.GetParameter("south"), "south")          self.assertEquals(proj.GetParameter("south"), "south")
89    
90      def test_get_projection_units_geo(self):      def test_get_projection_units_geo(self):
91          """Test Projection.GetProjectedUnits() for geographic projection"""          """Test Projection.GetProjectedUnits() for geographic projection.
92            Test for the alias 'longlat' as well.
93            """
94          proj = Projection(["proj=latlong", "to_meter=0.017453292519943295",          proj = Projection(["proj=latlong", "to_meter=0.017453292519943295",
95                             "ellps=clrk66"])                             "ellps=clrk66"])
96          self.assertEquals(proj.GetProjectedUnits(), PROJ_UNITS_DEGREES)          self.assertEquals(proj.GetProjectedUnits(), PROJ_UNITS_DEGREES)
97            proj = Projection(["proj=longlat", "to_meter=0.017453292519943295",
98                               "ellps=clrk66"])
99            self.assertEquals(proj.GetProjectedUnits(), PROJ_UNITS_DEGREES)
100    
101        def test_lc_numeric_robustness(self):
102            """Test if an LC_NUMERIC local with comma as decimal_point will work.
103    
104            Some versions of proj are not robust against this.
105            Starting with Python 2.4 there is a different behaviour when
106            calling C extensions and now they will see changes locales settings
107            which might tickle the bug in proj.
108            """
109            params = ["proj=latlong", "to_meter=0.01745", "ellps=clrk66"]
110    
111            oldlocale = localessupport.setdecimalcommalocale()
112            if oldlocale == None:
113                raise support.SkipTest(
114                        "No locale with comma as decimal_point found.")
115    
116            proj = Projection(params)
117            result1 =  proj.Forward(1,1)
118    
119            locale.setlocale(locale.LC_NUMERIC, "C")
120            proj = Projection(params)
121            result2= proj.Forward(1,1)
122    
123            locale.setlocale(locale.LC_NUMERIC, oldlocale)
124            self.assertFloatSeqEqual(result1, result2, epsilon = 1e-5 )
125    
126      def test_get_projection_units_normal(self):      def test_get_projection_units_normal(self):
127          """Test Projection.GetProjectedUnits() for normal projection"""          """Test Projection.GetProjectedUnits() for normal projection"""
# Line 114  class TestProjection(unittest.TestCase, Line 153  class TestProjection(unittest.TestCase,
153          self.assertEquals(proj.EPSGCode(), "42")          self.assertEquals(proj.EPSGCode(), "42")
154    
155    
 class ProjFileTest(unittest.TestCase, support.FileTestMixin):  
156    
157      """Base class for the proj file tests"""  class TestProjFileSimple:
158    
159      def filename(self):      def test_init(self):
160          """Return the filename for the test"""          """Test ProjFile coinstructor"""
161          return self.temp_file_name(self.id() + ".proj")          proj_file = ProjFile("some_filename")
162            self.assertEquals(proj_file.GetFilename(), "some_filename")
163            self.assertEquals(len(proj_file.GetProjections()), 0)
164    
165        def test_set_filename(self):
166            """Test ProjFile.SetFilename()"""
167            proj_file = ProjFile("some_filename")
168            proj.SetFilename("other_name")
169            self.assertEquals(proj_file.GetFilename(), "other_name")
170    
 class TestProjFile(ProjFileTest, xmlsupport.ValidationTest):  
171    
172      """Test cases for reading and writing projection files.  class TestProjFile(unittest.TestCase, support.SubscriberMixin):
     """  
173    
174      def test(self):      """Test cases for ProjFile objects"""
         """Test ProjFile"""  
   
         proj0 = Projection(["proj=tmerc", "ellps=clrk66"])  
         proj1 = Projection(["proj=utm", "ellps=clrk66"])  
         proj2 = Projection(["proj=lcc", "ellps=clrk66",  
                             "lat_1=0", "lat_2=20"])  
   
         eq = self.assertEquals  
   
   
         #  
         # __init__()  
         # GetFilename()  
         # SetFilename()  
         #  
         for name in ["", "hello_world"]:  
             projFile = ProjFile(name)  
             eq(projFile.GetFilename(), name)  
   
             projFile.SetFilename("XXX")  
             projFile.SetFilename(name)  
             eq(projFile.GetFilename(), name)  
   
         # initial number of projections should be 0  
         eq(len(projFile.GetProjections()), 0)  
   
         #  
         # Add()  
         # Remove()  
         #  
         projFile.Add(proj0)  
         eq(len(projFile.GetProjections()), 1)  
         projFile.Remove(proj0)  
         eq(len(projFile.GetProjections()), 0)  
   
         # try to remove something that doesn't exist  
         self.assertRaises(ValueError, projFile.Remove, None)  
         self.assertRaises(ValueError, projFile.Remove, proj0)  
   
         projFile.Add(proj0)  
         projFile.Add(proj1)  
         projFile.Add(proj2)  
         eq(len(projFile.GetProjections()), 3)  
   
         # GetProjections() -- tests order  
         projs = projFile.GetProjections()  
         eq(projs[0], proj0)  
         eq(projs[1], proj1)  
         eq(projs[2], proj2)  
175    
176          projFile.Remove(proj2)      def setUp(self):
177          projFile.Remove(proj1)          self.clear_messages()
178            self.proj0 = Projection(["proj=tmerc", "ellps=clrk66"])
179            self.proj1 = Projection(["proj=utm", "ellps=clrk66"])
180            self.proj2 = Projection(["proj=lcc", "ellps=clrk66",
181                                     "lat_1=0", "lat_2=20"])
182            self.proj_file = ProjFile("some_filename")
183            for msg in [PROJECTION_ADDED, PROJECTION_REMOVED, PROJECTION_REPLACED]:
184                self.proj_file.Subscribe(msg, self.subscribe_with_params, msg)
185    
186        def tearDown(self):
187            self.clear_messages()
188            self.proj_file.Destroy()
189    
190        def test_add_remove(self):
191            """Test ProjFile.Add() and ProjFile.Remove()"""
192            self.proj_file.Add(self.proj0)
193            self.proj_file.Add(self.proj1)
194            self.assertEquals(self.proj_file.GetProjections(),
195                              [self.proj0, self.proj1])
196            self.check_messages([(self.proj0, PROJECTION_ADDED),
197                                 (self.proj1, PROJECTION_ADDED)])
198            self.clear_messages()
199    
200            self.proj_file.Remove(self.proj0)
201            self.assertEquals(self.proj_file.GetProjections(), [self.proj1])
202            self.check_messages([(self.proj0, PROJECTION_REMOVED)])
203    
204        def test_remove_non_existing(self):
205            """Test ProjFile.Remove(<proj not in projfile>)"""
206            self.assertRaises(ValueError, self.proj_file.Remove, self.proj0)
207            # Nothing happened, so no messages should have been sent
208            self.check_messages([])
209    
210        def test_replace(self):
211            """Test ProjFile.Replace()"""
212            self.proj_file.Add(self.proj0)
213            self.proj_file.Add(self.proj1)
214            self.clear_messages()
215    
216          # Replace()          # Replace()
217          projFile.Replace(proj0, proj1)          self.proj_file.Replace(self.proj0, self.proj2)
218          projs = projFile.GetProjections()          self.assertEquals(self.proj_file.GetProjections(),
219          eq(projs[0], proj1)                            [self.proj2, self.proj1])
220            self.check_messages([(self.proj0, self.proj2, PROJECTION_REPLACED)])
221          # replace a non-existent projection  
222          self.assertRaises(ValueError, projFile.Replace, None, proj2)      def test_replace_non_existing(self):
223          self.assertRaises(ValueError, projFile.Replace, proj0, proj2)          """Test ProjFile.Replace(<proj not in projfile>, <some proj>)"""
224            self.proj_file.Add(self.proj0)
225            self.proj_file.Add(self.proj1)
226            self.clear_messages()
227            self.assertRaises(ValueError,
228                              self.proj_file.Replace, self.proj2, self.proj0)
229            # All projections should still be there
230            self.assertEquals(self.proj_file.GetProjections(),
231                              [self.proj0, self.proj1])
232            # Nothing happened, so no messages should have been sent
233            self.check_messages([])
234    
235    
236    class ProjFileTest(unittest.TestCase, support.FileTestMixin,
237                       xmlsupport.ValidationTest):
238    
239        """Base class for the proj file tests that read or write files"""
240    
241        def filename(self):
242            """Return the filename for the test"""
243            return self.temp_file_name(self.id() + ".proj")
244    
245    
246    class ProjFileReadTests(ProjFileTest):
247    
248        """Test read ProjFile objects from files
249    
250        The files only cover error handling and the system projection file.
251        """
252    
253      def test_read_non_existing_file(self):      def test_read_non_existing_file(self):
254          """Test read_proj_file with non-existing file"""          """Test read_proj_file with non-existing file"""
# Line 203  class TestProjFile(ProjFileTest, xmlsupp Line 262  class TestProjFile(ProjFileTest, xmlsupp
262          As currently written this only works on unix-like systems and          As currently written this only works on unix-like systems and
263          not e.g. on MS Windows.          not e.g. on MS Windows.
264          """          """
265            if os.name != "posix":
266                raise support.SkipTest("Test only works on posix systems")
267          filename = self.filename()          filename = self.filename()
268          file = open(filename, "w")          file = open(filename, "w")
269          file.close()          file.close()
# Line 218  class TestProjFile(ProjFileTest, xmlsupp Line 279  class TestProjFile(ProjFileTest, xmlsupp
279          self.assertRaises(SAXParseException, resource.read_proj_file, filename)          self.assertRaises(SAXParseException, resource.read_proj_file, filename)
280    
281      def test_get_system_proj_file(self):      def test_get_system_proj_file(self):
282          """Test resource.get_system_proj_file()          """Test resource.get_system_proj_file(DEFAULT_PROJ_FILE)
283    
284          This is primarily to test whether the system proj file contains          This is primarily to test whether the system proj file contains
285          invalid projection paramers and whether the proj file is not          invalid projection paramers and whether the proj file is not
286          empty          empty
287          """          """
288          projfile, warnings = resource.get_system_proj_file()          projfile, warnings\
289                      = resource.get_system_proj_file(resource.DEFAULT_PROJ_FILE)
290          self.assertEquals(warnings, [])          self.assertEquals(warnings, [])
291          self.assert_(len(projfile.GetProjections()) > 0)          self.assert_(len(projfile.GetProjections()) > 0)
292    
293            # see whether it got cached and we get the same projfile object
294            # when we read the file again
295            projfile2, warnings \
296                       = resource.get_system_proj_file(resource.DEFAULT_PROJ_FILE)
297            self.assert_(projfile is projfile2)
298    
299    
300  class WriteProjFileTests(ProjFileTest, xmlsupport.ValidationTest):  class WriteProjFileTests(ProjFileTest):
301    
302      """Test cases for writing proj files"""      """Test cases for writing proj files"""
303    
# Line 296  class WriteProjFileTests(ProjFileTest, x Line 364  class WriteProjFileTests(ProjFileTest, x
364          self.doTestWrite(pf, file_contents)          self.doTestWrite(pf, file_contents)
365    
366    
367  class TestLoadingProjFile(support.FileLoadTestCase):  class ProjFileLoadTestCase(support.FileLoadTestCase):
368    
369        """Base class for the test cases that read specific test files"""
370    
371      file_extension = ".proj"      file_extension = ".proj"
372    
373        def tearDown(self):
374            """Clear the cache explicitly"""
375            resource.clear_proj_file_cache()
376    
377    
378    class TestLoadingProjFile(ProjFileLoadTestCase):
379    
380      file_contents = '''\      file_contents = '''\
381  <?xml version="1.0" encoding="UTF-8"?>  <?xml version="1.0" encoding="UTF-8"?>
382  <!DOCTYPE projectionlist SYSTEM "projfile.dtd">  <!DOCTYPE projectionlist SYSTEM "projfile.dtd">
# Line 350  class TestLoadingProjFile(support.FileLo Line 428  class TestLoadingProjFile(support.FileLo
428                                 "proj=tmerc", "units=m",                                 "proj=tmerc", "units=m",
429                                 "x_0=400000.000", "y_0=0.000"])                                 "x_0=400000.000", "y_0=0.000"])
430    
431        def test_caching(self):
432            # test whether the projfile cache works
433            projfile, warnings = resource.read_proj_file(self.filename())
434            projfile2, warnings = resource.read_proj_file(self.filename())
435    
436  class TestLoadingProjFileWithEmptyProjectionlist(support.FileLoadTestCase):          # Both projfiles should be the same object
437            self.assert_(projfile2 is projfile)
438    
439            # If we clear the cache we should get a new one.
440            resource.clear_proj_file_cache()
441            projfile3, warnings = resource.read_proj_file(self.filename())
442            self.assert_(projfile3 is not projfile)
443    
444    
445    class TestLoadingProjFileWithEmptyProjectionlist(ProjFileLoadTestCase):
446    
     file_extension = ".proj"  
447      file_contents = '''\      file_contents = '''\
448  <?xml version="1.0" encoding="UTF-8"?>  <?xml version="1.0" encoding="UTF-8"?>
449  <!DOCTYPE projectionlist SYSTEM "projfile.dtd">  <!DOCTYPE projectionlist SYSTEM "projfile.dtd">
# Line 370  class TestLoadingProjFileWithEmptyProjec Line 460  class TestLoadingProjFileWithEmptyProjec
460          self.assertEquals(len(projfile.GetProjections()), 0)          self.assertEquals(len(projfile.GetProjections()), 0)
461    
462    
463  class TestProjFileWithInvalidParameters(unittest.TestCase,  class TestProjFileWithInvalidParameters(ProjFileLoadTestCase):
                                         support.FileLoadTestCase):  
464    
     file_extension = ".proj"  
465      file_contents = '''\      file_contents = '''\
466  <?xml version="1.0" encoding="UTF-8"?>  <?xml version="1.0" encoding="UTF-8"?>
467  <!DOCTYPE projectionlist SYSTEM "projfile.dtd">  <!DOCTYPE projectionlist SYSTEM "projfile.dtd">
# Line 413  class TestProjFileWithInvalidParameters( Line 501  class TestProjFileWithInvalidParameters(
501    
502    
503  if __name__ == "__main__":  if __name__ == "__main__":
504      unittest.main()      support.run_tests()
Legend:
Removed from v.1818  
changed lines
  Added in v.2696
[email protected]
 ViewVC Help
Powered by ViewVC 1.1.26