shapeFile = shaperead("shapeFiles/London.shp");
shapeFileR = shaperead("shapeFiles/LondonRoads.shp");
totalBox = getBoundingBox(shapeFile);
boxSize = 1e3;
avgLat = mean(totalBox(3:4));
avgLong = mean(totalBox(1,2));
longDist = CoorToDist(totalBox(1),totalBox(2),avgLat,avgLat);
latDist = CoorToDist(avgLong,avgLong,totalBox(3),totalBox(4));
kmToLat = 1/(latDist/1e3);
kmToLong = 1/(longDist/1e3);
boxScale = [kmToLong,kmToLat,kmToLong,kmToLat];
xSize = 1000;
ySize = 1000;
boundingBox = [3,3,1,1].*boxScale;
polygons = GenerateMap(boundingBox,shapeFile,1000,[xSize,ySize],totalBox);
polygonsR = GenerateMap(boundingBox,shapeFileR,200,[xSize,ySize],totalBox);
img = ones(xSize,ySize,3);
img = insertShape(img,'FilledPolygon',polygons,'Opacity',1);
img = insertShape(img,'Line',polygonsR,'LineWidth',3,'Opacity',1);
imshow(img);