def setUp(self):
self.tempdir = TempDir()
self.outimage = os.path.join(self.tempdir.dir, 'test.imi')
self.mi = MapImage(self.outimage)
self.mi.open('w')
self.m = self.mi.createMap()
self.m.scale = 1e-3
def testCreateLayerSansBBox(self):
m = self.m
layer = Layer(m, name="The Layer", filename="layer", layertype=LayerTypePolyline)
m.addLayer(layer)
layer.open("w")
line = CellElementPolyline(((16.185, 58.5912), (16.186, 58.5915)), objtype=2)
layer.addCellElement(line)
self.mi.close()
## Read back
mi = MapImage(self.outimage)
mi.open('r')
m = mi.maps[0]
class TestBBox(myTestCase):
"""Test to create map without a priori bounding box"""
def setUp(self):
self.tempdir = TempDir()
self.outimage = os.path.join(self.tempdir.dir, 'test.imi')
self.mi = MapImage(self.outimage)
self.mi.open('w')
self.m = self.mi.createMap()
self.m.scale = 1e-3
def tearDown(self):
self.m.close()
self.mi.close()
def testCreateLayerSansBBox(self):
m = self.m
layer = Layer(m, name="The Layer", filename="layer", layertype=LayerTypePolyline)
m.addLayer(layer)
layer.open("w")
line = CellElementPolyline(((16.185, 58.5912), (16.186, 58.5915)), objtype=2)
layer.addCellElement(line)
self.mi.close()
## Read back
mi = MapImage(self.outimage)
mi.open('r')
m = mi.maps[0]
def testCreateMapImage(self):
imagefile = os.path.join(self.testdatadir.dir, 'test.imi')
mi = MapImage(imagefile)
mi.open('w')
m = mi.createMap()
mi.close()
img = Image(imagefile)
self.assertTrue('bmp4bit.ics' in img.listdir())
self.assertTrue('add_maps.cfg' in img.listdir())
self.assertTrue('00map.ini' in img.listdir())
def testAddPoint(self):
outimage = os.path.join(self.tempdir.dir, 'test.imi')
mi = MapImage(outimage)
mi.open('w')
m = mi.createMap()
m.scale = 1e-3
m.bbox = ((-3.0, -4.0), (2.0, 1.0))
m.open("w")
layer = Layer(m, name="The Layer", filename="layer", layertype=LayerTypePoint, nlevels = 1)
m.addLayer(layer)
layer.open("w")
point = CellElementPoint((-1.5006, -3.5006))
layer.addCellElement(point)
cell = layer.getCell(point.cellnum)
self.assertEqual(point.cellnum, 4)
point.discretizeGeometry(layer.getCell(point.cellnum))
self.assertAlmostEqual(point.x, -1.501)
self.assertAlmostEqual(point.y, -3.501)
mi.close()
## Read back
mi = MapImage(outimage)
mi.open('r')
m = mi.maps[0]
import sys
from magellan.mapdir import MapDirectory, Image
from magellan.Map import Map, MapImage
from magellan.Layer import Layer, LayerTypePolyline, DetailMapLayerStyle
from magellan.CellElement import CellElementPolyline
from magellan.SearchGroup import Feature,FeatureNormal,FeatureStreet, GroupNormal
import osgeo.ogr as ogr
from magellan.rsttable import toRSTtable
layerinfo = True
dbinfo = False
filename = sys.argv[1]
mi = MapImage(filename)
mi.open('r')
m = mi.maps[0]
m.open('r')
print "Group names: ", m.getGroupNames()
print ""
## Print layer element
print 80 * "*"
print "Layers:"
if layerinfo:
for layer in m.layers:
layer.open('r')
print "Layer: ", layer.name
import sys
from magellan.mapdir import MapDirectory, Image
from magellan.Map import Map, MapImage
from magellan.Layer import Layer, LayerTypePolyline, DetailMapLayerStyle
from magellan.CellElement import CellElementPolyline
from magellan.SearchGroup import Feature, FeatureNormal, FeatureStreet, GroupNormal
import osgeo.ogr as ogr
from magellan.rsttable import toRSTtable
layerinfo = True
dbinfo = False
filename = sys.argv[1]
mi = MapImage(filename)
mi.open('r')
m = mi.maps[0]
m.open('r')
print "Group names: ", m.getGroupNames()
print ""
## Print layer element
print 80 * "*"
print "Layers:"
if layerinfo:
for layer in m.layers:
layer.open('r')
print "Layer: ", layer.name
