research

Author: shrddr

make_buildings.py

@@ -13,42 +13,37 @@
 def mil(fp):
     return math.floor(fp*1000000)
 
-if __name__ == "__main__": 
+if __name__ == "__main__":
 #    box = (27.4631,53.9359,27.4989,53.9566)
     box_minsk = (27.4013,53.8157,27.7827,53.9739)
-#    box_linetype = (27.643,53.843,27.647,53.846)
-#    box_hrodna = (23.7483,53.5909,23.9145,53.7544)
-    
-    irl_destroyed = []
-    irl_completed = []
+#    box_hrodna = (23.7483,53.5909,23.9145,53.7544)   
 
     q = loaders.Querier()
     ways = q.query_buildings(*box_minsk)
 
     IMZ = 18
-#    target = helpers.cleandir(f"buildings/yes")
-#    for wayid,nodes in ways:
-#        print(wayid)
-#        
-#        # fetch every tile the building touches
-#        for i,node in enumerate(nodes):
-#            fname = layers.maxar.gettile_wgs(node, IMZ, skipedge=True, edge=16)
-#            if fname is None:
-#                continue
-#            dst = str(target)
-#            shutil.copy(fname, dst)
+    target = helpers.cleandir(f"buildings/yes")
+    for wayid,nodes in ways:
+        print(wayid)
+        
+        # fetch every tile the building has a node at
+        for i,node in enumerate(nodes):
+            fname = layers.maxar.gettile_wgs(node, IMZ, skipedge=True, edge=16)
+            if fname is None:
+                continue
+            dst = str(target)
+            shutil.copy(fname, dst)
 
         # OR focus on the building and crop
 #        image = layers.maxar.tiles_way(nodes, IMZ, pad_pct=0.25, pad_px=48)
 #        dst = str(target / f"m{wayid}.jpg")
 #        cv2.imwrite(dst, image)
 
     mp = helpers.MercatorPainter(layers.maxar, *box_minsk, IMZ)
-    for wayid,nodes in ways:
+    for _, nodes in ways:
         mp.add_polyline_wgs(nodes, width=1)
-#    mp.show_fixedwindow(600, 600)
-#    mp.random_check(151202,84361)
-#    cv2.destroyAllWindows()
+    # TODO: add restricted zones where imagery is outdated
+    # use cv2.polyline(is_closed=True)
 
     target = helpers.cleandir(f"buildings/no")
 

readme.md

@@ -110,7 +110,7 @@ Task: classify [roof:shape](https://wiki.openstreetmap.org/wiki/Simple_3D_buildi
 
 ![Roofshapes](readme-images/roofshapes.png)
 
-Input: varies (see table)
+Input: varies (see table). As expected, the more the better!
 
 | Input size    | Frozen, 8 epoch train | Unfreeze, 8 epochs | 8 more epochs |
 | ------------- | --------------------- | ------------------ | ------------- |
@@ -127,13 +127,17 @@ Task: classify tiles with any type of building(s) vs. tiles with no buildings at
 
 ![Buildings: yes/no](readme-images/buildings_yes_no.png)
 
-Input: 5700 images of both categories, 20% validation
+Input: 5700 images of both categories, 20% validation.
 
-|                              | Frozen, 4 epoch train | Unfreeze, 6 epochs | 8 more epochs |
-| ---------------------------- | --------------------- | ------------------ | ------------- |
-| 256 no flip no rotation      | 3.0%                  | 3.0%               |               |
-| 256 no flip full 20 rotation | 2.3%                  | 2.7%               |               |
-| 256 no flip full 40 rotation |                       |                    |               |
-|                              |                       |                    |               |
-|                              |                       |                    |               |
-|                              |                       |                    |               |
+## First attempt
+
+For every node of every building, fetch the tile it belongs to. Don't use the tile if the node is closer than 16px to the edge (that would present only a tiny part of the building).
+
+|                                  | Frozen, 4 epoch train | Unfreeze, 6 epochs | 8 more epochs |
+| -------------------------------- | --------------------- | ------------------ | ------------- |
+| size=256 no flip, max_rotate=0   | 3.0%                  | 3.0%               |               |
+| size=256 no flip, max_rotate=20  | 2.3%                  | 2.7%               |               |
+| size=256, no flip, max_rotate=40 | 3.3%                  | 3.3%               |               |
+|                                  |                       |                    |               |
+|                                  |                       |                    |               |
+|                                  |                       |                    |               |