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 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]
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]
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