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