net2 0.0.1__py3-none-any.whl

net2/net2.py

@@ -0,0 +1,210 @@
+#!/usr/bin/python3
+import numpy as np
+import geopandas as gpd
+import shapely
+import networkx as nx
+
+# Prepare network
+def prepare(G, ids_to_int=True):
+    # IDs to 'int'
+    if ids_to_int:
+        mapping = {i: int(i) for i in G.nodes()}
+        G = nx.relabel_nodes(G, mapping)
+    # nodes
+    for i in G.nodes:
+    # 'geometry' to 'shapely'
+        if 'geometry' in G.nodes[i]:
+            if isinstance(G.nodes[i]['geometry'], str):
+                G.nodes[i]['geometry'] = shapely.from_wkt((G.nodes[i]['geometry']))
+    # edges
+    for u,v in G.edges:
+        # 'geometry' to 'shapely'
+        if 'geometry' in G[u][v]:
+            if isinstance(G[u][v]['geometry'], str):
+                G[u][v]['geometry'] = shapely.from_wkt((G[u][v]['geometry']))
+        # 'length' to 'float'
+        if 'length' in G[u][v]:
+            G[u][v]['length'] = float(G[u][v]['length'])
+        # rename 'travel_time' to 'time'
+        if 'travel_time' in G[u][v]:
+            G[u][v]['time'] = G[u][v]['travel_time']
+            del G[u][v]['travel_time']
+        # 'time' to 'float'
+        if 'time' in G[u][v]:
+            G[u][v]['time'] = float(G[u][v]['time'])
+    return G
+
+# Calculate 'pos' dictionary with node 'x' and 'y'
+def pos(G):
+    return {i: [G.nodes[i]['geometry'].x, G.nodes[i]['geometry'].y] for i in G.nodes}
+
+# Network nodes to 'GeoDataFrame'
+def nodes_to_gdf(network, crs=None):
+    geom = []
+    node_id = []
+    for i in network.nodes:
+        geom.append(network.nodes[i]['geometry'])
+        node_id.append(i)
+    nodes = gpd.GeoDataFrame({'id':node_id, 'geometry':geom}, crs=crs)
+    return nodes
+
+# Network edges to 'GeoDataFrame'
+def edges_to_gdf(network, crs=None):
+    edges = nx.to_pandas_edgelist(network)
+    edges = gpd.GeoDataFrame(edges, crs=crs)
+    return edges
+
+# Nearest node to given point
+def nearest_node(G: nx.Graph, pnt: shapely.geometry.Point) -> tuple:
+    """
+    Find the nearest network node to the specified point
+    
+    :param G: `networkx` Network
+    :param pnt: `shapely` geometry of type 'Point'
+    """
+    min_distance = float('inf')
+    nearest_node = None
+    for i in G.nodes:
+        distance = pnt.distance(G.nodes[i]['geometry'])
+        if distance < min_distance:
+            min_distance = distance
+            nearest_node = i
+    return nearest_node, min_distance
+
+# Nearest edge to given point
+def nearest_edge(G: nx.Graph, pnt: shapely.geometry.Point) -> tuple:
+    """
+    Find the nearest network edge to the specified point
+    
+    :param G: `networkx` Network
+    :param pnt: `shapely` geometry of type 'Point'
+    """
+    min_distance = float('inf')
+    nearest_edge = None
+    for i in G.edges:
+        geom_edge = G.edges[i]['geometry']
+        distance = pnt.distance(geom_edge)
+        if distance < min_distance:
+            min_distance = distance
+            nearest_edge = i
+    return nearest_edge, min_distance
+
+# Split edge by point
+def split_edge(G, node_id, e, pnt_on_line, buffer_size):
+    pnt_on_line_b = pnt_on_line.buffer(buffer_size)
+    first_seg, buff_seg, last_seg = shapely.ops.split(G.edges[e]['geometry'], pnt_on_line_b).geoms
+    G.edges[e]['geometry'] = shapely.LineString(list(first_seg.coords) + list(pnt_on_line.coords) + list(last_seg.coords))
+    lines = shapely.ops.split(G.edges[e]['geometry'], pnt_on_line)
+    p = G.nodes[e[0]]['geometry']
+    line1 = filter(shapely.prepared.prep(p).intersects, lines.geoms)
+    line1 = list(line1)[0]
+    p = G.nodes[e[1]]['geometry']
+    line2 = filter(shapely.prepared.prep(p).intersects, lines.geoms)
+    line2 = list(line2)[0]
+    G.add_edge(e[0], node_id, geometry=line1, length=line1.length, time=G.edges[e]['time'] * (line1.length / G.edges[e]['geometry'].length))
+    G.add_edge(node_id, e[1], geometry=line2, length=line2.length, time=G.edges[e]['time'] * (line2.length / G.edges[e]['geometry'].length))
+    G.remove_edge(*e)
+    return G
+# Check if number is similar
+def is_same(a, b, treshold=0.001):
+    return abs(a - b) < treshold
+# Add node to network
+def add_node(network, pnt, buffer_size=1e-8):
+    G = network.copy()
+    # Detect nearest edge
+    e, dist = nearest_edge(G, pnt)
+    # Detect nearest node (from within the nearest edge)
+    node_id, dist = nearest_node(G.subgraph(e), pnt)
+    # If 'pnt' is on edisting node -> return that node
+    if dist == 0:
+        return G, node_id, dist
+    # Detect nearest point on the edge
+    pnt_on_line = shapely.ops.nearest_points(G.edges[e]['geometry'], pnt)
+    pnt_on_line = pnt_on_line[0]
+    # If 'pnt_on_line' is on edisting node -> return that node
+    if is_same(pnt_on_line.x, G.nodes[node_id]['geometry'].x) and is_same(pnt_on_line.y, G.nodes[node_id]['geometry'].y):
+        return G, node_id, dist
+    # Else - create new node
+    node_ids = [i for i in G.nodes if isinstance(i, (int, float))]
+    if(len(node_ids) > 0):
+        node_id = min(node_ids)-1
+    else:
+        node_id = -1
+    if node_id >= 0:
+        node_id = -1
+    G.add_node(node_id, geometry=pnt_on_line)
+    # Split edge
+    G = split_edge(G, node_id, e, pnt_on_line, buffer_size)
+    e = (e[1], e[0])
+    if e in G.edges and G.edges[e]['geometry'].intersects(pnt_on_line.buffer(buffer_size)):
+        G = split_edge(G, node_id, e, pnt_on_line, buffer_size)
+    return G, node_id, dist
+
+# Shortest path 1 (between existing nodes)
+def route1(network, node_start, node_end, weight):
+    try:
+        route = nx.shortest_path(network, node_start, node_end, weight)
+        weight_sum = nx.path_weight(network, route, weight=weight)
+        return {'route': route, 'weight': weight_sum}
+    except:
+        return {'route': np.nan, 'weight': np.nan}
+
+# Shortest path 2 (inserting new nodes)
+def route2(network, start, end, weight):
+    network, node_start, dist_start = add_node(network, start)
+    network, node_end, dist_end = add_node(network, end)
+    try:
+        route = nx.shortest_path(network, node_start, node_end, weight)
+        weight_sum = nx.path_weight(network, route, weight=weight)
+        return {'route': route, 'weight': weight_sum, 'dist_start': dist_start, 'dist_end': dist_end, 'network': network}
+    except:
+        return {'route': np.nan, 'weight': np.nan, 'dist_start': dist_start, 'dist_end': dist_end, 'network': network}
+
+# Shortest path 3 (inserting new nodes + walk to nodes + comparison w/ walking mode)
+def route3(network, start, end, time_weight, walking_speed=1.4): # m/s
+    import math
+    dist = math.sqrt(((start.x - end.x) ** 2) + ((start.y - end.y) ** 2))
+    time_walking = dist / walking_speed
+    try:
+        result = route2(network, start, end, time_weight)
+        time_driving_and_walking = result['dist_start']/walking_speed + result['weight'] + result['dist_end']/walking_speed
+    except:
+        return {'weight': np.nan}
+    if time_driving_and_walking <= time_walking:
+        return {'weight': time_driving_and_walking, 'mode': 'walking+driving'}
+    else:
+        return {'weight': time_walking, 'mode': 'walking'}
+
+# Create regular grid
+def create_grid(bounds, res, crs=None):
+    xmin, ymin, xmax, ymax = bounds
+    cols = list(np.arange(int(np.floor(xmin)), int(np.ceil(xmax+res)), res))
+    rows = list(np.arange(int(np.floor(ymin)), int(np.ceil(ymax+res)), res))
+    rows.reverse()
+    polygons = []
+    for x in cols:
+        for y in rows:
+            polygons.append(
+                shapely.Polygon([(x,y), (x+res, y), (x+res, y-res), (x, y-res)])
+            )
+    grid = gpd.GeoDataFrame({'geometry': polygons}, crs=crs)
+    sel = grid.intersects(shapely.box(*bounds))
+    grid = grid[sel]
+    return grid
+
+# Route to 'GeoDataFrame'
+def route_to_gdf(network, route, crs=None):
+    route_edges = nx.path_graph(route).edges
+    if len(route) == 1:
+        result = []
+        pnt = network.nodes[route[0]]['geometry']
+        line = shapely.LineString([pnt, pnt])
+        result = gpd.GeoDataFrame([{'from': route[0], 'to': route[0], 'geometry': line}], crs=crs)
+    if len(route) > 1:
+        result = []
+        for u,v in route_edges:
+            x = network.edges[u,v]['geometry']
+            result.append({'from': u, 'to': v, 'geometry': x})
+        result = gpd.GeoDataFrame(result, crs=crs)
+    return result
+

net2-0.0.1.dist-info/METADATA

@@ -0,0 +1,35 @@
+Metadata-Version: 2.4
+Name: net2
+Version: 0.0.1
+Summary: Package for working with spatial networks
+Project-URL: Homepage, https://github.com/michaeldorman/net2/
+Project-URL: Issues, https://github.com/michaeldorman/net2/issues
+Author-email: Michael Dorman <dorman@post.bgu.ac.il>
+License: Copyright (c) 2026 The Python Packaging Authority
+        
+        Permission is hereby granted, free of charge, to any person obtaining a copy
+        of this software and associated documentation files (the "Software"), to deal
+        in the Software without restriction, including without limitation the rights
+        to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+        copies of the Software, and to permit persons to whom the Software is
+        furnished to do so, subject to the following conditions:
+        
+        The above copyright notice and this permission notice shall be included in all
+        copies or substantial portions of the Software.
+        
+        THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+        IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+        FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+        AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+        LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+        OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+        SOFTWARE.
+License-File: LICENSE
+Classifier: Operating System :: OS Independent
+Classifier: Programming Language :: Python :: 3
+Requires-Python: >=3.12
+Description-Content-Type: text/markdown
+
+# `net2`
+
+`net2` is a Python package for working with spatial networks.

net2-0.0.1.dist-info/RECORD

@@ -0,0 +1,6 @@
+net2/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+net2/net2.py,sha256=V3l5FgLYT87N_HdBdIVnjI_dKmoCb9SQ1oQTOSlLTwY,8169
+net2-0.0.1.dist-info/METADATA,sha256=sthiv9FcJWmiFs67My1CCKwQ_Nnf0tjkV5kSDqVIuec,1761
+net2-0.0.1.dist-info/WHEEL,sha256=WLgqFyCfm_KASv4WHyYy0P3pM_m7J5L9k2skdKLirC8,87
+net2-0.0.1.dist-info/licenses/LICENSE,sha256=CUUUEnItJwG4KAPWBman2XvVZgg0TPMn6FwA9GlVxek,1074
+net2-0.0.1.dist-info/RECORD,,

net2-0.0.1.dist-info/WHEEL

@@ -0,0 +1,4 @@
+Wheel-Version: 1.0
+Generator: hatchling 1.28.0
+Root-Is-Purelib: true
+Tag: py3-none-any

net2-0.0.1.dist-info/licenses/LICENSE

@@ -0,0 +1,19 @@
+Copyright (c) 2026 The Python Packaging Authority
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.