Update networks/legacy_roads/house.py

House.py

@@ -4,8 +4,6 @@ import numpy as np
 import math
 import matplotlib.pyplot as plt
 
-from pyxtension.streams import stream
-
 
 class House:
     def __init__(self, editor, coordinates_min, coordinates_max, direction, list_block):
@@ -51,38 +49,37 @@ class House:
 
         x_min += 1
         z_min += 1
-        x_max -= 2
+        x_max -= 1
-        z_max -= 2
+        z_max -= 1
-        if x_min + 1 > x_max :
+        if x_min + 1 > x_max - 1:
-            x = np.random.randint(x_max, x_min + 1)
+            x = np.random.randint(x_max - 1, x_min + 1)
         else:
-            x = np.random.randint(x_min + 1, x_max)
+            x = np.random.randint(x_min + 1, x_max - 1)
 
-        if z_min + 1 > z_max:
+        if z_min + 1 > z_max - 1:
-            z = np.random.randint(z_max, z_min + 1)
+            z = np.random.randint(z_max - 1, z_min + 1)
         else:
-            z = np.random.randint(z_min + 1, z_max)
+            z = np.random.randint(z_min + 1, z_max - 1)
 
         width = perimeter_width // 2
         depth = perimeter_depth // 2
         height = y_max - y_min
 
-        if x + width - 1 > x_max:
+        if x + width - 1 > x_max - 1:
-            x = x_max - width
+            x = x_max - width - 1
-        if z + depth - 1 > z_max:
+        if z + depth - 1 > z_max - 1:
-            z = z_max - depth
+            z = z_max - depth - 1
 
         x_plan3d = x - x_min
         z_plan3d = z - z_min
 
-        geometry.placeCuboid(self.editor, (x, y_min, z), (x + width - 1, y_min, z + depth - 1), self.floor)
+        for i in range(0, width - 1):
-        x_range = slice(x_plan3d, x_plan3d + width - 1)
+            for j in range(0, depth - 1):
-        z_range = slice(z_plan3d, z_plan3d + depth - 1)
+                self.editor.placeBlock((x + i, y_min, z + j), self.floor)
-        self.grid3d[x_range, 0, z_range] = True, 1
+                self.grid3d[x_plan3d + i, 0, z_plan3d + j] = True, 1
-        
         self.skeleton.append((x, z, width - 1, depth - 1, height))
-        print("Coordinates of the corners: ", (x, z), (x, z + depth - 1), 
+        print("Coordinates of the corners: ", (x, z), (x, z + depth - 1), (x + width - 1, z),
-              (x + width - 1, z), (x + width - 1, z + depth - 1))
+              (x + width - 1, z + depth - 1))
 
         x_min -= 1
         x_max -= 1
@@ -107,47 +104,54 @@ class House:
                 else:
                     new_width = np.random.randint(width - 2, 5)
 
-                if max(x_min+1, x-new_width) > min(x_max-new_width, x+width):
+                if max(x_min+1, x-new_width) > min(x_max-new_width-1, x+width):
-                    new_x = np.random.randint(min(x_max - new_width, x + width), 
+                    new_x = np.random.randint(
-                                              max(x_min + 1, x - new_width))
+                        min(x_max - new_width - 1, x + width), max(x_min + 1, x - new_width))
                 else:
-                    new_x = np.random.randint(max(x_min + 1, x - new_width), 
+                    new_x = np.random.randint(
-                                              min(x_max - new_width, x + width))
+                        max(x_min + 1, x - new_width), min(x_max - new_width - 1, x + width))
 
-                if max(z_min+1, z-new_depth) > min(z_max-new_depth, z+depth):
+                if max(z_min+1, z-new_depth) > min(z_max-new_depth-1, z+depth):
-                    new_z = np.random.randint(min(z_max - new_depth, z + depth), 
+                    new_z = np.random.randint(
-                                              max(z_min + 1, z - new_depth))
+                        min(z_max - new_depth - 1, z + depth), max(z_min + 1, z - new_depth))
                 else:
-                    new_z = np.random.randint(max(z_min + 1, z - new_depth), 
+                    new_z = np.random.randint(
-                                              min(z_max - new_depth, z + depth))
+                        max(z_min + 1, z - new_depth), min(z_max - new_depth - 1, z + depth))
 
                 new_x_plan3d = new_x - x_min - 1
                 new_z_plan3d = new_z - z_min + 1
 
                 adjacent_blocks = 0
-                for i in range(new_x_plan3d, new_x_plan3d + new_width) :
+                for i in range(new_x_plan3d, new_x_plan3d + new_width):
-                    for j in range(new_z_plan3d, new_z_plan3d + new_depth) :
+                    for j in range(new_z_plan3d, new_z_plan3d + new_depth):
-                        if self.grid3d[     i - 1, 0, j    ]['bool'] and self.grid3d[i - 1, 0, j    ]['int'] == 1 or \
+                        if self.grid3d[i - 1, 0, j]['bool'] and self.grid3d[i - 1, 0, j]['int'] == 1 or \
-                                self.grid3d[i + 1, 0, j    ]['bool'] and self.grid3d[i + 1, 0, j    ]['int'] == 1 or \
+                                self.grid3d[i + 1, 0, j]['bool'] and self.grid3d[i + 1, 0, j]['int'] == 1 or \
                                 self.grid3d[i, 0, j - 1]['bool'] and self.grid3d[i, 0, j - 1]['int'] == 1 or \
                                 self.grid3d[i, 0, j + 1]['bool'] and self.grid3d[i, 0, j + 1]['int'] == 1:
                             adjacent_blocks += 1
 
-                new_x_range = slice(new_x_plan3d, new_x_plan3d + new_width)
+                if adjacent_blocks < 3:
-                new_z_range = slice(new_z_plan3d, new_z_plan3d + new_depth)
-                if adjacent_blocks < 3 or np.any(self.grid3d[new_x_range, 0, new_z_range]['bool']) :
                     continue
 
-                geometry.placeCuboid(self.editor, (new_x, y_min, new_z), (new_x+new_width-1, y_min, new_z+new_depth-1), self.floor)
+                if not np.any(
-                new_x_plan3d += 1
+                        self.grid3d[new_x_plan3d:new_x_plan3d + new_width, 0, new_z_plan3d:new_z_plan3d + new_depth][
-                new_z_plan3d -= 1 
+                            'bool']):
-                new_x_range = slice(new_x_plan3d, new_x_plan3d + new_width)
+                    new_x_plan3d = new_x - x_min
-                new_z_range = slice(new_z_plan3d, new_z_plan3d + new_depth) 
+                    new_z_plan3d = new_z - z_min
-                self.grid3d[new_x_range, 0, new_z_range] = True, 2
+                    for i in range(0, new_width):
+                        for j in range(0, new_depth):
+                            self.grid3d[new_x_plan3d + i, 0,
+                                        new_z_plan3d + j] = True, 2
 
-                self.skeleton.append((new_x, new_z, new_width, new_depth, height))
+                            if i == 0 or i == new_width - 1 or j == 0 or j == new_depth - 1:
+                                continue
+                            else:
+                                self.editor.placeBlock(
+                                    (new_x + i, y_min, new_z + j), self.floor)
+
+                    self.skeleton.append(
+                        (new_x, new_z, new_width, new_depth, height))
                     break
-                    
             else:
                 print("Failed to place rectangle after 100000 attempts.")
 
@@ -1245,7 +1249,6 @@ class House:
             self.placeStairs()
 
 
-
 if __name__ == "__main__":
     editor = Editor(buffering=True)
     buildArea = editor.getBuildArea()

main.py

@@ -15,8 +15,7 @@ from networks.geometry.Point3D import Point3D
 from networks.geometry.Point2D import Point2D
 from networks.geometry.Circle import Circle
 
-from PIL import Image as img
+from PIL import Image
-from PIL.Image import Image
 from utils.JsonReader import JsonReader
 from utils.YamlReader import YamlReader
 from buildings.Building import Building
@@ -45,7 +44,6 @@ def main():
     time_remove_tree = time.time() - start_time
     print(f"[TIME] Remove tree {time_remove_tree}")
 
-    """
     start_time = time.time()
     smooth_terrain('./world_maker/data/heightmap.png',
                    './world_maker/data/heightmap_smooth.png', './world_maker/data/smooth_sobel_watermap.png')
@@ -119,17 +117,17 @@ def main():
                       entranceDirection[random.randint(0, 3)], blocks)
         house.build()
     time_houses = time.time() - start_time
-    print(f"[TIME] Houses {time_houses}") """
+    print(f"[TIME] Houses {time_houses}")
 
     print("[GDMC] Done!\n\n")
 
     print(f"[TIME] Total {time.time() - start_time_all}")
     print(f"[TIME] World_maker {time_world_maker}")
     print(f"[TIME] Remove tree {time_remove_tree}")
-    #print(f"[TIME] Smooth terrain {time_smooth_terrain}")
+    print(f"[TIME] Smooth terrain {time_smooth_terrain}")
-    #print(f"[TIME] Roads {time_roads}")
+    print(f"[TIME] Roads {time_roads}")
-    #print(f"[TIME] Buildings {time_buildings}")
+    print(f"[TIME] Buildings {time_buildings}")
-    #print(f"[TIME] Houses {time_houses}")
+    print(f"[TIME] Houses {time_houses}")
 
 
 def get_height_building_from_center(center, position, length_world):
@@ -164,7 +162,7 @@ def set_roads(skeleton: Skeleton, origin):
         for j in range(len(skeleton.lines[i])):
             xyz = transpose_form_heightmap('./world_maker/data/heightmap.png',
                                            skeleton.coordinates[skeleton.lines[i][j]], origin)
-            heightmap_smooth = img.open(
+            heightmap_smooth = Image.open(
                 './world_maker/data/road_heightmap.png')
             skeleton.lines[i][j] = [xyz[0], heightmap_smooth.getpixel(
                 (skeleton.coordinates[skeleton.lines[i][j]][0], skeleton.coordinates[skeleton.lines[i][j]][-1])), xyz[2]]

networks/legacy_roads/house.py

@@ -48,7 +48,7 @@ def mur_sol(co1,co2,block):
         
         tailleZ=co2[2]-co1[2]
         midtailleZ=(tailleZ//2)+z1
-    editor = Editor(buffering=  True) 
+    editor = Editor(host="http://gdmc.ale-pri.com:9000", buffering=  True) 
     
     if y1==y2:
         

world_maker/Skeleton.py

@@ -4,8 +4,7 @@ from typing import List, Union
 
 import numpy as np
 from gdpc import Editor
-from PIL.Image import Image
+from PIL import Image, ImageDraw
-from PIL import Image as img, ImageDraw
 from skan.csr import skeleton_to_csgraph
 from skimage.morphology import skeletonize
 from networks.geometry.Point3D import Point3D
@@ -13,7 +12,7 @@ from networks.geometry.Point3D import Point3D
 
 def handle_import_image(image: Union[str, Image]) -> Image:
     if isinstance(image, str):
-        return img.open(image)
+        return Image.open(image)
     return image
 
 
@@ -204,7 +203,7 @@ class Skeleton:
         # xzDistance = (max(buildRect.end[0], buildRect.begin[0]) - min(buildRect.end[0], buildRect.begin[0]),
         #              max(buildRect.end[1], buildRect.begin[1]) - min(buildRect.end[1], buildRect.begin[1]))
 
-        heightmap = img.open(
+        heightmap = Image.open(
             "./world_maker/data/heightmap.png").convert('RGB')
         # roadsArea = Image.new("L", xzDistance, 0)
         # width, height = heightmap.size
@@ -265,9 +264,9 @@ class Skeleton:
 
     def road_area(self, name: str, radius: int = 10) -> Image:
         print("[Skeleton] Start mapping the road area...")
-        heightmap = img.open("./world_maker/data/heightmap.png")
+        heightmap = Image.open("./world_maker/data/heightmap.png")
         width, height = heightmap.size
-        road_area_map = img.new("L", (width, height), 0)
+        road_area_map = Image.new("L", (width, height), 0)
         road_area_map_draw = ImageDraw.Draw(road_area_map)
 
         # Lines

world_maker/World.py

@@ -29,10 +29,8 @@ class World:
         editor = Editor(buffering=True)
         buildArea = editor.getBuildArea()
 
-        self.coordinates_min = [
+        self.coordinates_min = [min(buildArea.begin[i], buildArea.last[i]) for i in range(3)]
-            min(buildArea.begin[i], buildArea.last[i]) for i in range(3)]
+        self.coordinates_max = [max(buildArea.begin[i], buildArea.last[i]) for i in range(3)]
-        self.coordinates_max = [
-            max(buildArea.begin[i], buildArea.last[i]) for i in range(3)]
 
         self.length_x = self.coordinates_max[0] - self.coordinates_min[0] + 1
         self.length_y = self.coordinates_max[1] - self.coordinates_min[1] + 1
@@ -64,8 +62,7 @@ class World:
         Add block or list of block to the volume.
         """
 
-        self.volume[volumeCoordinates[0]][volumeCoordinates[1]
+        self.volume[volumeCoordinates[0]][volumeCoordinates[1]][volumeCoordinates[2]] = None
-                                          ][volumeCoordinates[2]] = None
 
     def getBlockFromCoordinates(self, coordinates):
         """
@@ -88,11 +85,9 @@ class World:
             for j in range(-1, 2):
                 for k in range(-1, 2):
                     if not (i == 0 and j == 0 and k == 0):
-                        coordinates = (
+                        coordinates = (Block.coordinates[0] + i, Block.coordinates[1] + j, Block.coordinates[2] + k)
-                            Block.coordinates[0] + i, Block.coordinates[1] + j, Block.coordinates[2] + k)
                         if self.isInVolume(coordinates):
-                            Block.addNeighbors(
+                            Block.addNeighbors([self.getBlockFromCoordinates(coordinates)])
-                                [self.getBlockFromCoordinates(coordinates)])
 
     def setVolume(self):
         """
@@ -104,8 +99,7 @@ class World:
         for x in range(self.coordinates_min[0], self.coordinates_max[0] + 1):
             for y in range(self.coordinates_min[1], self.coordinates_max[1] + 1):
                 for z in range(self.coordinates_min[2], self.coordinates_max[2] + 1):
-                    self.addBlocks(
+                    self.addBlocks([Block((x, y, z), editor.getBlock((x, y, z)).id)])
-                        [Block((x, y, z), editor.getBlock((x, y, z)).id)])
 
     def getData(self):
         """
@@ -126,10 +120,8 @@ class World:
 
         slice = editor.loadWorldSlice(buildRect)
 
-        heightmapData = list(
+        heightmapData = list(np.array(slice.heightmaps["MOTION_BLOCKING_NO_LEAVES"], dtype=np.uint8))
-            np.array(slice.heightmaps["MOTION_BLOCKING_NO_LEAVES"], dtype=np.uint16))
+        treesmapData = list(np.array(slice.heightmaps["MOTION_BLOCKING"], dtype=np.uint8))
-        treesmapData = list(
-            np.array(slice.heightmaps["MOTION_BLOCKING"], dtype=np.uint16))
 
         for x in range(0, xzDistance[0]):
             for z in range(0, xzDistance[1]):
@@ -156,8 +148,7 @@ class World:
 
                                     # print('getData for tree', xzStart[0] + x + i, k, xzStart[1] + z + j)
 
-                                    blockNeighbor = slice.getBlock(
+                                    blockNeighbor = slice.getBlock((x + i, k, z + j))
-                                        (x + i, k, z + j))
                                     if blockNeighbor.id not in lookup.TREES:
                                         height += k
                                         number += 1
@@ -171,9 +162,7 @@ class World:
                 else:
                     watermap.putpixel((x, z), 0)
 
-                # y set to 100 for 2D
+                self.addBlocks([Block((xzStart[0] + x, 100, xzStart[1] + z), block)])  # y set to 100 for 2D
-                self.addBlocks(
-                    [Block((xzStart[0] + x, 100, xzStart[1] + z), block)])
 
         return heightmap, watermap, treesmap
 
@@ -224,8 +213,7 @@ class World:
 
         slice = editor.loadWorldSlice(buildRect)
 
-        heightmapData = list(
+        heightmapData = list(np.array(slice.heightmaps["MOTION_BLOCKING_NO_LEAVES"], dtype=np.uint8))
-            np.array(slice.heightmaps["MOTION_BLOCKING_NO_LEAVES"], dtype=np.uint8))
 
         for x in range(0, xzDistance[0]):
             for z in range(0, xzDistance[1]):

world_maker/terraforming.py

@@ -2,8 +2,7 @@ from typing import Union
 
 import numpy as np
 from gdpc import Editor, Block, geometry, lookup
-from PIL       import Image as img
+from PIL import Image
-from PIL.Image import Image
 from skimage import morphology
 
 from world_maker.data_analysis import handle_import_image
@@ -17,25 +16,38 @@ def remove_trees(heightmap: Union[str, Image], treesmap: Union[str, Image], mask
 
     start = build_rectangle.begin
 
-    distance = (max(build_rectangle.end[0], build_rectangle.begin[0]) - min(build_rectangle.end[0], build_rectangle.begin[0]),
+    distance = (max(build_rectangle.end[0], build_rectangle.begin[0]) - min(build_rectangle.end[0], build_rectangle.begin[0]), max(
-                max(build_rectangle.end[1], build_rectangle.begin[1]) - min(build_rectangle.end[1], build_rectangle.begin[1]))
+        build_rectangle.end[1], build_rectangle.begin[1]) - min(build_rectangle.end[1], build_rectangle.begin[1]))
 
     heightmap = handle_import_image(heightmap).convert('L')
     treesmap = handle_import_image(treesmap).convert('L')
     mask = handle_import_image(mask)
 
-    removed_treesmap = img.new("L", distance, 0)
+    removed_treesmap = Image.new("L", distance, 0)
 
-    for (x, z) in [(x, z) for x in range(distance[0]) for z in range(distance[1])] :
+    removed = []
+    for x in range(0, distance[0]):
+        for z in range(0, distance[1]):
 
-        if mask.getpixel((x, z)) != 0 and treesmap.getpixel((x, z)) > 0:
+            if mask.getpixel((x, z)) != 0 and treesmap.getpixel((x, z)) > 0 and (x, z) not in removed:
 
                 tree_area = morphology.flood(treesmap, (z, x), tolerance=1)
-            removed_treesmap = img.fromarray(np.where(tree_area, treesmap, 0).astype(np.uint8))
+                blend = Image.blend(Image.fromarray(tree_area).convert(
+                    'L'), removed_treesmap.convert('L'), 0.5)
 
+                array = np.array(blend.convert('L'))
+                bool_array = array > 1
+                removed_treesmap = Image.fromarray(bool_array)
+
+                removed.append((x, z))
+
+    for x in range(0, distance[0]):
+        for z in range(0, distance[1]):
+            if removed_treesmap.getpixel((x, z)) != 0:
                 y = heightmap.getpixel((x, z))
                 y_top = removed_treesmap.getpixel((x, z))
-            geometry.placeLine(editor, (start[0] + x, y+1, start[1] + z), (start[0] + x, y_top, start[1] + z), Block('air'))
+                geometry.placeLine(
+                    editor, (start[0] + x, y+1, start[1] + z), (start[0] + x, y_top, start[1] + z), Block('air'))
 
     removed_treesmap.save('./world_maker/data/removed_treesmap.png')
     print("[Remove tree] Done.")
@@ -50,39 +62,39 @@ def smooth_terrain(heightmap: Union[str, Image], heightmap_smooth: Union[str, Im
 
     start = build_rectangle.begin
 
-    distance = (max(build_rectangle.end[0], build_rectangle.begin[0]) - min(build_rectangle.end[0], build_rectangle.begin[0]),
+    distance = (max(build_rectangle.end[0], build_rectangle.begin[0]) - min(build_rectangle.end[0], build_rectangle.begin[0]), max(
-                max(build_rectangle.end[1], build_rectangle.begin[1]) - min(build_rectangle.end[1], build_rectangle.begin[1]))
+        build_rectangle.end[1], build_rectangle.begin[1]) - min(build_rectangle.end[1], build_rectangle.begin[1]))
 
     heightmap = handle_import_image(heightmap).convert('L')
     heightmap_smooth = handle_import_image(heightmap_smooth).convert('L')
     mask = handle_import_image(mask).convert('L')
 
-    smooth_terrain_delta = img.new("RGB", distance, 0)
+    smooth_terrain_delta = Image.new("RGB", distance, 0)
 
     slice = editor.loadWorldSlice(build_rectangle)
     smoothable_blocks = lookup.OVERWORLD_SOILS | lookup.OVERWORLD_STONES | lookup.SNOWS
 
-    for (x, z) in [(x, z) for x in range(distance[0]) for z in range(distance[1])] :
+    for x in range(0, distance[0]):
+        for z in range(0, distance[1]):
 
-        if mask.getpixel((x, z)) == 0:
+            if mask.getpixel((x, z)) != 0:
-            continue
                 y = heightmap.getpixel((x, z))
                 y_smooth = heightmap_smooth.getpixel((x, z))
                 delta = y - y_smooth
                 smooth_terrain_delta.putpixel((x, z), delta)
 
-        if delta == 0:
+                if delta != 0:
-            continue
                     block = slice.getBlock((x, y, z))
-
+                    if block.id in smoothable_blocks:
-        if block.id not in smoothable_blocks:
-            continue
                         if delta > 0:
-            geometry.placeLine(editor, (start[0] + x, y, start[1] + z), (start[0] + x, y_smooth, start[1] + z), Block('air'))
+                            geometry.placeLine(
-            editor.placeBlock((start[0] + x, y_smooth, start[1] + z), block)
+                                editor, (start[0] + x, y, start[1] + z), (start[0] + x, y_smooth, start[1] + z), Block('air'))
+                            editor.placeBlock(
+                                (start[0] + x, y_smooth, start[1] + z), block)
 
                         else:
-            geometry.placeLine(editor, (start[0] + x, y, start[1] + z), (start[0] + x, y_smooth, start[1] + z), block)
+                            geometry.placeLine(
+                                editor, (start[0] + x, y, start[1] + z), (start[0] + x, y_smooth, start[1] + z), block)
 
     smooth_terrain_delta.save('./world_maker/data/smooth_terrain_delta.png')
     print("[Smooth terrain] Done.")