voxcity 0.4.4__py3-none-any.whl → 0.4.5__py3-none-any.whl

voxcity/geoprocessor/mesh.py

@@ -4,7 +4,7 @@ import matplotlib.colors as mcolors
 import matplotlib.cm as cm
 import matplotlib.pyplot as plt
 
-def create_voxel_mesh(voxel_array, class_id, meshsize=1.0, building_id_grid=None):
+def create_voxel_mesh(voxel_array, class_id, meshsize=1.0, building_id_grid=None, mesh_type=None):
     """
     Create a mesh from voxels preserving sharp edges, scaled by meshsize.
 
@@ -18,6 +18,11 @@ def create_voxel_mesh(voxel_array, class_id, meshsize=1.0, building_id_grid=None
         The real-world size of each voxel in meters, for x, y, and z.
     building_id_grid : np.ndarray (2D), optional
         2D grid of building IDs, shape (X, Y). Used when class_id=-3 (buildings).
+    mesh_type : str, optional
+        Type of mesh to create:
+        - None (default): create meshes for boundaries between different classes
+        - 'building_solar': only create meshes for boundaries between buildings (-3) 
+                           and void (0) or trees (-2)
 
     Returns
     -------
@@ -81,7 +86,7 @@ def create_voxel_mesh(voxel_array, class_id, meshsize=1.0, building_id_grid=None
             (x,   y-1, z)     # Bottom
         ]
 
-        # Only create faces where there's a transition between this class and "not this class"
+        # Only create faces where there's a transition based on mesh_type
         for face_idx, adj_coord in enumerate(adjacent_coords):
             try:
                 # If adj_coord is outside array bounds, it's a boundary => face is visible
@@ -89,7 +94,13 @@ def create_voxel_mesh(voxel_array, class_id, meshsize=1.0, building_id_grid=None
                     is_boundary = True
                 else:
                     adj_value = voxel_array[adj_coord]
-                    is_boundary = (adj_value == 0 or adj_value != class_id)
+                    
+                    if mesh_type == 'open_air' and class_id == -3:
+                        # For building_solar, only create faces at boundaries with void (0) or trees (-2)
+                        is_boundary = (adj_value == 0 or adj_value == -2)
+                    else:
+                        # Default behavior - create faces at any class change
+                        is_boundary = (adj_value == 0 or adj_value != class_id)
             except IndexError:
                 # Out of range => boundary
                 is_boundary = True

voxcity/simulator/view.py

@@ -1708,7 +1708,8 @@ def get_surface_view_factor(voxel_data, meshsize, **kwargs):
             voxel_data, 
             building_class_id, 
             meshsize,
-            building_id_grid=building_id_grid
+            building_id_grid=building_id_grid,
+            mesh_type='open_air'
         )
         if building_mesh is None or len(building_mesh.faces) == 0:
             print("No surfaces found in voxel data for the specified class.")

voxcity/utils/visualization.py

@@ -1270,7 +1270,7 @@ def visualize_numerical_grid_on_basemap(grid, rectangle_vertices, meshsize, valu
     plt.tight_layout()
     plt.show()
 
-def visualize_numerical_grid_gdf_on_basemap(gdf, value_name="value", cmap='viridis', vmin=None, vmax=None,
+def visualize_numerical_gdf_on_basemap(gdf, value_name="value", cmap='viridis', vmin=None, vmax=None,
                             alpha=0.6, figsize=(12, 8), basemap='CartoDB light',
                             show_edge=False, edge_color='black', edge_width=0.5):
     """Visualizes a GeoDataFrame with numerical values on a basemap.
@@ -1298,7 +1298,7 @@ def visualize_numerical_grid_gdf_on_basemap(gdf, value_name="value", cmap='virid
     fig, ax = plt.subplots(figsize=figsize)
     
     # Plot the GeoDataFrame
-    gdf_web.plot(column='value',
+    gdf_web.plot(column=value_name,
                  ax=ax,
                  alpha=alpha,
                  cmap=cmap,
@@ -1325,7 +1325,7 @@ def visualize_numerical_grid_gdf_on_basemap(gdf, value_name="value", cmap='virid
     plt.tight_layout()
     plt.show()
 
-def visualize_point_grid_on_basemap(point_gdf, value_name='value', **kwargs):
+def visualize_point_gdf_on_basemap(point_gdf, value_name='value', **kwargs):
     """Visualizes a point GeoDataFrame on a basemap with colors based on values.
     
     Args:

{voxcity-0.4.4.dist-info → voxcity-0.4.5.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.2
 Name: voxcity
-Version: 0.4.4
+Version: 0.4.5
 Summary: voxcity is an easy and one-stop tool to output 3d city models for microclimate simulation by integrating multiple geospatial open-data
 Author-email: Kunihiko Fujiwara <kunihiko@nus.edu.sg>
 Maintainer-email: Kunihiko Fujiwara <kunihiko@nus.edu.sg>

{voxcity-0.4.4.dist-info → voxcity-0.4.5.dist-info}/RECORD

@@ -16,22 +16,22 @@ voxcity/exporter/obj.py,sha256=0RBFPMKGRH6uNmCLIwAoYFko1bOZKtTSwg7QnoPMud0,21593
 voxcity/geoprocessor/__init_.py,sha256=JzPVhhttxBWvaZ0IGX2w7OWL5bCo_TIvpHefWeNXruA,133
 voxcity/geoprocessor/draw.py,sha256=8Em2NvazFpYfFJUqG9LofNXaxdghKLL_rNuztmPwn8Q,13911
 voxcity/geoprocessor/grid.py,sha256=t-KAdYOk1iUeBuOMcoAb5e4hzJjMIfVItTdFP6FJs1A,44340
-voxcity/geoprocessor/mesh.py,sha256=z6e-BigN_K1CGpS9mopkf121pbCnBFJMu32XYfqpLI4,10506
+voxcity/geoprocessor/mesh.py,sha256=r3cRPLgpbhjwgESBemHWWJ5pEWl2KdkRhID6mdLhios,11171
 voxcity/geoprocessor/network.py,sha256=opb_kpUCAxDd1qtrWPStqR5reYZtVe96XxazNSen7Lk,18851
 voxcity/geoprocessor/polygon.py,sha256=8Vb2AbkpKYhq1kk2hQMc-gitmUo9pFIe910v4p1vP2g,37772
 voxcity/geoprocessor/utils.py,sha256=1BRHp-DDeOA8HG8jplY7Eo75G3oXkVGL6DGONL4BA8A,19815
 voxcity/simulator/__init_.py,sha256=APdkcdaovj0v_RPOaA4SBvFUKT2RM7Hxuuz3Sux4gCo,65
 voxcity/simulator/solar.py,sha256=amJ5GJtzdG9NE7oxrLj4EMRxGvDY9s2BQdDJ2paFfn0,56709
 voxcity/simulator/utils.py,sha256=sEYBB2-hLJxTiXQps1_-Fi7t1HN3-1OPOvBCWtgIisA,130
-voxcity/simulator/view.py,sha256=-CM-pPbfthg_meyMdcZs9Ar3CHaZhw2fRf_dM6ioN_Q,74875
+voxcity/simulator/view.py,sha256=YufbLuDXrLg1d1dedM6pVyiJ7uHsqY8F2sLLnIoJvB4,74910
 voxcity/utils/__init_.py,sha256=nLYrj2huBbDBNMqfchCwexGP8Tlt9O_XluVDG7MoFkw,98
 voxcity/utils/lc.py,sha256=RwPd-VY3POV3gTrBhM7TubgGb9MCd3nVah_G8iUEF7k,11562
 voxcity/utils/material.py,sha256=Vt3IID5Ft54HNJcEC4zi31BCPqi_687X3CSp7rXaRVY,5907
-voxcity/utils/visualization.py,sha256=UJ2wyFFbuI6glZJ8adyiVtNlZFUgleC278S-r2OO33g,85240
+voxcity/utils/visualization.py,sha256=SF8W7sqvBl3sZbB5noWCY9ic2D34Gq01VZYJ9NDNZ4Y,85237
 voxcity/utils/weather.py,sha256=CFPtoqRTajwMRswswDChwQ3BW1cGsnA3orgWHgz7Ehg,26304
-voxcity-0.4.4.dist-info/AUTHORS.rst,sha256=m82vkI5QokEGdcHof2OxK39lf81w1P58kG9ZNNAKS9U,175
-voxcity-0.4.4.dist-info/LICENSE,sha256=-hGliOFiwUrUSoZiB5WF90xXGqinKyqiDI2t6hrnam8,1087
-voxcity-0.4.4.dist-info/METADATA,sha256=73PRyC5f-BPeFk8FpmZvaUvUVYGVDqAHMK-_Rj_oQs0,25527
-voxcity-0.4.4.dist-info/WHEEL,sha256=52BFRY2Up02UkjOa29eZOS2VxUrpPORXg1pkohGGUS8,91
-voxcity-0.4.4.dist-info/top_level.txt,sha256=00b2U-LKfDllt6RL1R33MXie5MvxzUFye0NGD96t_8I,8
-voxcity-0.4.4.dist-info/RECORD,,
+voxcity-0.4.5.dist-info/AUTHORS.rst,sha256=m82vkI5QokEGdcHof2OxK39lf81w1P58kG9ZNNAKS9U,175
+voxcity-0.4.5.dist-info/LICENSE,sha256=-hGliOFiwUrUSoZiB5WF90xXGqinKyqiDI2t6hrnam8,1087
+voxcity-0.4.5.dist-info/METADATA,sha256=YDINX8e1EqSKFcCMTlPdFYj-wrWH7zIazT1fE89MkrQ,25527
+voxcity-0.4.5.dist-info/WHEEL,sha256=52BFRY2Up02UkjOa29eZOS2VxUrpPORXg1pkohGGUS8,91
+voxcity-0.4.5.dist-info/top_level.txt,sha256=00b2U-LKfDllt6RL1R33MXie5MvxzUFye0NGD96t_8I,8
+voxcity-0.4.5.dist-info/RECORD,,