refactoring : terraforming.py

main.py

@@ -15,7 +15,8 @@ from networks.geometry.Point3D import Point3D
 from networks.geometry.Point2D import Point2D
 from networks.geometry.Circle import Circle
 
-from PIL import Image
+from PIL import Image as img
+from PIL.Image import Image
 from utils.JsonReader import JsonReader
 from utils.YamlReader import YamlReader
 from buildings.Building import Building
@@ -44,6 +45,7 @@ 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')
@@ -117,17 +119,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):
@@ -162,7 +164,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 = Image.open(
+            heightmap_smooth = img.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(host="http://gdmc.ale-pri.com:9000", buffering=  True) 
+    editor = Editor(buffering=  True) 
     
     if y1==y2:
         

world_maker/Skeleton.py

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

world_maker/World.py

@@ -29,8 +29,10 @@ class World:
         editor = Editor(buffering=True)
         buildArea = editor.getBuildArea()
 
-        self.coordinates_min = [min(buildArea.begin[i], buildArea.last[i]) for i in range(3)]
+        self.coordinates_min = [
-        self.coordinates_max = [max(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.length_x = self.coordinates_max[0] - self.coordinates_min[0] + 1
         self.length_y = self.coordinates_max[1] - self.coordinates_min[1] + 1
@@ -62,7 +64,8 @@ class World:
         Add block or list of block to the volume.
         """
 
-        self.volume[volumeCoordinates[0]][volumeCoordinates[1]][volumeCoordinates[2]] = None
+        self.volume[volumeCoordinates[0]][volumeCoordinates[1]
+                                          ][volumeCoordinates[2]] = None
 
     def getBlockFromCoordinates(self, coordinates):
         """
@@ -85,9 +88,11 @@ 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 = (Block.coordinates[0] + i, Block.coordinates[1] + j, Block.coordinates[2] + k)
+                        coordinates = (
+                            Block.coordinates[0] + i, Block.coordinates[1] + j, Block.coordinates[2] + k)
                         if self.isInVolume(coordinates):
-                            Block.addNeighbors([self.getBlockFromCoordinates(coordinates)])
+                            Block.addNeighbors(
+                                [self.getBlockFromCoordinates(coordinates)])
 
     def setVolume(self):
         """
@@ -99,7 +104,8 @@ 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([Block((x, y, z), editor.getBlock((x, y, z)).id)])
+                    self.addBlocks(
+                        [Block((x, y, z), editor.getBlock((x, y, z)).id)])
 
     def getData(self):
         """
@@ -120,8 +126,10 @@ class World:
 
         slice = editor.loadWorldSlice(buildRect)
 
-        heightmapData = list(np.array(slice.heightmaps["MOTION_BLOCKING_NO_LEAVES"], dtype=np.uint8))
+        heightmapData = list(
-        treesmapData = list(np.array(slice.heightmaps["MOTION_BLOCKING"], dtype=np.uint8))
+            np.array(slice.heightmaps["MOTION_BLOCKING_NO_LEAVES"], dtype=np.uint16))
+        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]):
@@ -148,7 +156,8 @@ class World:
 
                                     # print('getData for tree', xzStart[0] + x + i, k, xzStart[1] + z + j)
 
-                                    blockNeighbor = slice.getBlock((x + i, k, z + j))
+                                    blockNeighbor = slice.getBlock(
+                                        (x + i, k, z + j))
                                     if blockNeighbor.id not in lookup.TREES:
                                         height += k
                                         number += 1
@@ -162,7 +171,9 @@ class World:
                 else:
                     watermap.putpixel((x, z), 0)
 
-                self.addBlocks([Block((xzStart[0] + x, 100, xzStart[1] + z), block)])  # y set to 100 for 2D
+                # y set to 100 for 2D
+                self.addBlocks(
+                    [Block((xzStart[0] + x, 100, xzStart[1] + z), block)])
 
         return heightmap, watermap, treesmap
 
@@ -213,7 +224,8 @@ class World:
 
         slice = editor.loadWorldSlice(buildRect)
 
-        heightmapData = list(np.array(slice.heightmaps["MOTION_BLOCKING_NO_LEAVES"], dtype=np.uint8))
+        heightmapData = list(
+            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,7 +2,8 @@ from typing import Union
 
 import numpy as np
 from gdpc      import Editor, Block, geometry, lookup
-from PIL import Image
+from PIL       import Image as img
+from PIL.Image import Image
 from skimage   import morphology
 
 from world_maker.data_analysis import handle_import_image
@@ -16,38 +17,25 @@ 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]), max(
+    distance = (max(build_rectangle.end[0], build_rectangle.begin[0]) - min(build_rectangle.end[0], build_rectangle.begin[0]),
-        build_rectangle.end[1], build_rectangle.begin[1]) - min(build_rectangle.end[1], build_rectangle.begin[1]))
+                max(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 = Image.new("L", distance, 0)
+    removed_treesmap = img.new("L", distance, 0)
 
-    removed = []
+    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 and treesmap.getpixel((x, z)) > 0 and (x, z) not in removed:
+        if mask.getpixel((x, z)) != 0 and treesmap.getpixel((x, z)) > 0:
 
             tree_area = morphology.flood(treesmap, (z, x), tolerance=1)
-                blend = Image.blend(Image.fromarray(tree_area).convert(
+            removed_treesmap = img.fromarray(np.where(tree_area, treesmap, 0).astype(np.uint8))
-                    '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(
+            geometry.placeLine(editor, (start[0] + x, y+1, start[1] + z), (start[0] + x, y_top, start[1] + z), Block('air'))
-                    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.")
@@ -62,39 +50,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]), max(
+    distance = (max(build_rectangle.end[0], build_rectangle.begin[0]) - min(build_rectangle.end[0], build_rectangle.begin[0]),
-        build_rectangle.end[1], build_rectangle.begin[1]) - min(build_rectangle.end[1], build_rectangle.begin[1]))
+                max(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 = Image.new("RGB", distance, 0)
+    smooth_terrain_delta = img.new("RGB", distance, 0)
 
     slice = editor.loadWorldSlice(build_rectangle)
     smoothable_blocks = lookup.OVERWORLD_SOILS | lookup.OVERWORLD_STONES | lookup.SNOWS
 
-    for x in range(0, distance[0]):
+    for (x, z) in [(x, z) for x in range(distance[0]) for z in range(distance[1])] :
-        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(
+            geometry.placeLine(editor, (start[0] + x, y, start[1] + z), (start[0] + x, y_smooth, start[1] + z), Block('air'))
-                                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)
-                            editor.placeBlock(
-                                (start[0] + x, y_smooth, start[1] + z), block)
 
         else:
-                            geometry.placeLine(
+            geometry.placeLine(editor, (start[0] + x, y, start[1] + z), (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)
 
     smooth_terrain_delta.save('./world_maker/data/smooth_terrain_delta.png')
     print("[Smooth terrain] Done.")