ObjectNat 0.1.5__tar.gz → 0.2.1__tar.gz

{objectnat-0.1.5 → objectnat-0.2.1}/PKG-INFO

@@ -1,18 +1,18 @@
 Metadata-Version: 2.1
 Name: ObjectNat
-Version: 0.1.5
+Version: 0.2.1
 Summary: ObjectNat is an open-source library created for geospatial analysis created by IDU team
 License: BSD-3-Clause
-Author: Danila
+Author: DDonnyy
 Author-email: 63115678+DDonnyy@users.noreply.github.com
-Requires-Python: >=3.10,<4.0
+Requires-Python: >=3.10,<3.13
 Classifier: License :: OSI Approved :: BSD License
 Classifier: Programming Language :: Python :: 3
 Classifier: Programming Language :: Python :: 3.10
 Classifier: Programming Language :: Python :: 3.11
 Classifier: Programming Language :: Python :: 3.12
-Requires-Dist: dongraphio (>=0.3.13,<0.4.0)
 Requires-Dist: geopandas (>=0.14.3,<0.15.0)
+Requires-Dist: iduedu (>=0.1.4,<0.2.0)
 Requires-Dist: joblib (>=1.4.2,<2.0.0)
 Requires-Dist: networkit (>=11.0,<12.0)
 Requires-Dist: networkx (>=3.2.1,<4.0.0)
@@ -20,7 +20,7 @@ Requires-Dist: numpy (>=1.23.5,<2.0.0)
 Requires-Dist: pandarallel (>=1.6.5,<2.0.0)
 Requires-Dist: pandas (>=2.2.0,<3.0.0)
 Requires-Dist: population-restorator (>=0.2.3,<0.3.0)
-Requires-Dist: provisio (>=0.1.7,<0.2.0)
+Requires-Dist: pulp (>=2.8.0,<3.0.0)
 Requires-Dist: scikit-learn (>=1.4.0,<2.0.0)
 Requires-Dist: tqdm (>=4.66.2,<5.0.0)
 Description-Content-Type: text/markdown
@@ -28,59 +28,64 @@ Description-Content-Type: text/markdown
 # ObjectNat - Meta Library
 
 [![Code style: black](https://img.shields.io/badge/code%20style-black-000000.svg)](https://github.com/psf/black)
+[![PyPI version](https://img.shields.io/pypi/v/objectnat.svg)](https://pypi.org/project/objectnat/)
 
+- [РИДМИ (Russian)](README_ru.md)
 <p align="center">
-<img src="https://i.ibb.co/FWtHNQv/logo.png" alt="logo" width="400">
+<img src="https://github.com/user-attachments/assets/d3878cce-8eba-4f96-8458-9a798d436120" alt="logo" width="400">
 </p>
 
 #### **ObjectNat** is an open-source library created for geospatial analysis created by **IDU team**
 
 ## ObjectNat Components
 
-- [dongraphio](https://github.com/DDonnyy/dongraphio) : `dongraphio` provides graph functions
-- [provisio](https://github.com/DDonnyy/provisio) : `provisio` provides main provisio fuctions
+- [IduEdu](https://github.com/DDonnyy/IduEdu) : `IduEdu` provides graph functions
 - [population-restorator](https://github.com/kanootoko/population-restorator) : `restorator` provides city resettlement
 
 ## Features and how to use
 
-1. **[City graph from OSM](./examples/graph_generator.ipynb)** - Function to assemble a road, pedestrian, and public
-   transport graph from OpenStreetMap (OSM) and creating Intermodal graph.
+1. **[City graph from OSM (IduEdu)](https://github.com/DDonnyy/IduEdu/blob/main/examples/get_any_graph.ipynb)** - Functions to assemble a road, pedestrian,
+   and public transport graph from OpenStreetMap (OSM) and creating Intermodal graph. 
 
-   <img src="https://i.ibb.co/VpsPgL1/Graph-generator-v1.webp" alt="Graph-generator-v1" height="250">
+   <img src="https://github.com/user-attachments/assets/8dc98da9-8462-415e-8cc8-bdfca788e206" alt="IntermodalGraph" height="250">
 
 2. **[Adjacency matrix](./examples/calculate_adjacency_matrix.ipynb)** - Calculate adjacency matrix based on the provided
    graph and edge weight type (time or distance). The intermodal graph can be obtained using the previous example.
+
 3. **[Isochrones,transport accessibility](./examples/isochrone_generator.ipynb)** - Function for generating isochrones to
    analyze transportation accessibility from specified starting coordinates. Isochrones can be constructed based on
    pedestrian, automobile, or public transport graphs, or a combination thereof.
 
-   <img src="https://i.ibb.co/QM0tmZ2/isochrones-from-2-points.webp" alt="isochrones-from-2-points" height="250">
+   <img src="https://github.com/user-attachments/assets/37f308a5-db56-497d-b080-4edef3584fe5" alt="isochrones" height="250">
 
 4. **[Population restoration](./examples/restore_population.ipynb)** - Function for resettling population into the provided
    layer of residential buildings. This function distributes people among dwellings based on the total city population
    and the living area of each house.
 5. **[Service provision](./examples/calculate_provision.ipynb)** - Function for calculating the provision of residential
-   buildings and population with services. In case of missing data, this function utilizes previously described
-   functionality to retrieve the necessary information.
+   buildings and population with services. 
+
+   <img src="https://github.com/user-attachments/assets/5f2b3c55-9a02-4d70-80f4-503b77023eda" alt="ProvisionSchools" height="250">
 
-   <img src="https://i.ibb.co/CW7Xj5F/Burger-Provision5min.webp" alt="Burger-Provision5min" height="250">
-   
 6. **[Visibility analysis](./examples/visibility_analysis.ipynb)** - Function to get a quick estimate of visibility from a
    given point(s) to buildings within a given distance. Also, there is a visibility catchment area calculator for a
    large
    urban area. This function is designed to work with at least 1000 points spaced 10-20 meters apart for optimal
    results. Points can be generated using a road graph and random point distribution along edges.
 
-   <img src="https://i.ibb.co/LxcGTfN/visibility-from-point.webp" alt="visibility-from-point" height="250"> 
-
-   <img src="https://i.ibb.co/zNRzXc5/visibility-catchment-area.webp" alt="visibility-catchment-area" height="250">
+   <img src="https://github.com/user-attachments/assets/2927ac86-01e8-4b0e-9ea8-72ad81c13cf5" alt="visibility-from-point" height="250"> 
 
+   <img src="https://github.com/user-attachments/assets/b5b0d4b3-a02f-4ade-8772-475703cd6435" alt="visibility-catchment-area" height="250">
+   
 7. **[Point clusterization](./examples/point_clusterization.ipynb)** - Function to generate cluster polygons for given
    points based on a specified minimum distance and minimum points per cluster. Optionally, calculate the relative ratio
    between types of services within the clusters.
 
-   <img src="https://i.ibb.co/fF3c4YC/service-clusterization.webp" alt="service-clusterization" height="250">
-
+   <img src="https://github.com/user-attachments/assets/2a9ad722-87d2-4954-9612-5ac3765aa824" alt="service-clusterization" height="250">
+   
+8. **[Living buildings from OSM](./examples/download_buildings_from_osm.ipynb)** - This function downloads building geometries from OpenStreetMap (OSM) for a specified territory and assigns attributes to each building. Specifically, it determines whether a building is residential (`is_living` attribute) and estimates the approximate number of inhabitants (`approximate_pop` attribute).
+   
+   <img src="https://github.com/user-attachments/assets/d60dcd85-1a2e-4342-aae4-561aeda18858" alt="Living buildings" height="250">
+   
 ## Installation
 
 **ObjectNat** can be installed with ``pip``:
@@ -88,7 +93,16 @@ Description-Content-Type: text/markdown
 ```
 pip install ObjectNat
 ```
+### Configuration changes
 
+```python
+from objectnat import config
+
+config.set_timeout(10)  # Timeout for overpass queries
+config.change_logger_lvl('INFO')  # To mute all debug msgs
+config.set_enable_tqdm(False)  # To mute all tqdm's progress bars
+config.set_overpass_url('http://your.overpass-api.de/interpreter/URL')
+```
 ## Contacts
 
 - [NCCR](https://actcognitive.org/) - National

{objectnat-0.1.5 → objectnat-0.2.1}/README.md

@@ -1,59 +1,64 @@
 # ObjectNat - Meta Library
 
 [![Code style: black](https://img.shields.io/badge/code%20style-black-000000.svg)](https://github.com/psf/black)
+[![PyPI version](https://img.shields.io/pypi/v/objectnat.svg)](https://pypi.org/project/objectnat/)
 
+- [РИДМИ (Russian)](README_ru.md)
 <p align="center">
-<img src="https://i.ibb.co/FWtHNQv/logo.png" alt="logo" width="400">
+<img src="https://github.com/user-attachments/assets/d3878cce-8eba-4f96-8458-9a798d436120" alt="logo" width="400">
 </p>
 
 #### **ObjectNat** is an open-source library created for geospatial analysis created by **IDU team**
 
 ## ObjectNat Components
 
-- [dongraphio](https://github.com/DDonnyy/dongraphio) : `dongraphio` provides graph functions
-- [provisio](https://github.com/DDonnyy/provisio) : `provisio` provides main provisio fuctions
+- [IduEdu](https://github.com/DDonnyy/IduEdu) : `IduEdu` provides graph functions
 - [population-restorator](https://github.com/kanootoko/population-restorator) : `restorator` provides city resettlement
 
 ## Features and how to use
 
-1. **[City graph from OSM](./examples/graph_generator.ipynb)** - Function to assemble a road, pedestrian, and public
-   transport graph from OpenStreetMap (OSM) and creating Intermodal graph.
+1. **[City graph from OSM (IduEdu)](https://github.com/DDonnyy/IduEdu/blob/main/examples/get_any_graph.ipynb)** - Functions to assemble a road, pedestrian,
+   and public transport graph from OpenStreetMap (OSM) and creating Intermodal graph. 
 
-   <img src="https://i.ibb.co/VpsPgL1/Graph-generator-v1.webp" alt="Graph-generator-v1" height="250">
+   <img src="https://github.com/user-attachments/assets/8dc98da9-8462-415e-8cc8-bdfca788e206" alt="IntermodalGraph" height="250">
 
 2. **[Adjacency matrix](./examples/calculate_adjacency_matrix.ipynb)** - Calculate adjacency matrix based on the provided
    graph and edge weight type (time or distance). The intermodal graph can be obtained using the previous example.
+
 3. **[Isochrones,transport accessibility](./examples/isochrone_generator.ipynb)** - Function for generating isochrones to
    analyze transportation accessibility from specified starting coordinates. Isochrones can be constructed based on
    pedestrian, automobile, or public transport graphs, or a combination thereof.
 
-   <img src="https://i.ibb.co/QM0tmZ2/isochrones-from-2-points.webp" alt="isochrones-from-2-points" height="250">
+   <img src="https://github.com/user-attachments/assets/37f308a5-db56-497d-b080-4edef3584fe5" alt="isochrones" height="250">
 
 4. **[Population restoration](./examples/restore_population.ipynb)** - Function for resettling population into the provided
    layer of residential buildings. This function distributes people among dwellings based on the total city population
    and the living area of each house.
 5. **[Service provision](./examples/calculate_provision.ipynb)** - Function for calculating the provision of residential
-   buildings and population with services. In case of missing data, this function utilizes previously described
-   functionality to retrieve the necessary information.
+   buildings and population with services. 
+
+   <img src="https://github.com/user-attachments/assets/5f2b3c55-9a02-4d70-80f4-503b77023eda" alt="ProvisionSchools" height="250">
 
-   <img src="https://i.ibb.co/CW7Xj5F/Burger-Provision5min.webp" alt="Burger-Provision5min" height="250">
-   
 6. **[Visibility analysis](./examples/visibility_analysis.ipynb)** - Function to get a quick estimate of visibility from a
    given point(s) to buildings within a given distance. Also, there is a visibility catchment area calculator for a
    large
    urban area. This function is designed to work with at least 1000 points spaced 10-20 meters apart for optimal
    results. Points can be generated using a road graph and random point distribution along edges.
 
-   <img src="https://i.ibb.co/LxcGTfN/visibility-from-point.webp" alt="visibility-from-point" height="250"> 
-
-   <img src="https://i.ibb.co/zNRzXc5/visibility-catchment-area.webp" alt="visibility-catchment-area" height="250">
+   <img src="https://github.com/user-attachments/assets/2927ac86-01e8-4b0e-9ea8-72ad81c13cf5" alt="visibility-from-point" height="250"> 
 
+   <img src="https://github.com/user-attachments/assets/b5b0d4b3-a02f-4ade-8772-475703cd6435" alt="visibility-catchment-area" height="250">
+   
 7. **[Point clusterization](./examples/point_clusterization.ipynb)** - Function to generate cluster polygons for given
    points based on a specified minimum distance and minimum points per cluster. Optionally, calculate the relative ratio
    between types of services within the clusters.
 
-   <img src="https://i.ibb.co/fF3c4YC/service-clusterization.webp" alt="service-clusterization" height="250">
-
+   <img src="https://github.com/user-attachments/assets/2a9ad722-87d2-4954-9612-5ac3765aa824" alt="service-clusterization" height="250">
+   
+8. **[Living buildings from OSM](./examples/download_buildings_from_osm.ipynb)** - This function downloads building geometries from OpenStreetMap (OSM) for a specified territory and assigns attributes to each building. Specifically, it determines whether a building is residential (`is_living` attribute) and estimates the approximate number of inhabitants (`approximate_pop` attribute).
+   
+   <img src="https://github.com/user-attachments/assets/d60dcd85-1a2e-4342-aae4-561aeda18858" alt="Living buildings" height="250">
+   
 ## Installation
 
 **ObjectNat** can be installed with ``pip``:
@@ -61,7 +66,16 @@
 ```
 pip install ObjectNat
 ```
+### Configuration changes
 
+```python
+from objectnat import config
+
+config.set_timeout(10)  # Timeout for overpass queries
+config.change_logger_lvl('INFO')  # To mute all debug msgs
+config.set_enable_tqdm(False)  # To mute all tqdm's progress bars
+config.set_overpass_url('http://your.overpass-api.de/interpreter/URL')
+```
 ## Contacts
 
 - [NCCR](https://actcognitive.org/) - National

{objectnat-0.1.5 → objectnat-0.2.1}/pyproject.toml

@@ -1,28 +1,29 @@
 [tool.poetry]
 name = "ObjectNat"
-version = "0.1.5"
+version = "0.2.1"
 description = "ObjectNat is an open-source library created for geospatial analysis created by IDU team"
 license = "BSD-3-Clause"
-authors = ["Danila <63115678+DDonnyy@users.noreply.github.com>"]
+authors = ["DDonnyy <63115678+DDonnyy@users.noreply.github.com>"]
 readme = "README.md"
 
 packages = [{ include = "objectnat", from = "src" }]
 
 [tool.poetry.dependencies]
-python = "^3.10"
-geopandas = "^0.14.3"
-tqdm = "^4.66.2"
-networkit = "^11.0"
+python = ">=3.10,<3.13"
 numpy = "^1.23.5"
 pandas = "^2.2.0"
+geopandas = "^0.14.3"
+tqdm = "^4.66.2"
+pandarallel = "^1.6.5"
 networkx = "^3.2.1"
+networkit = "^11.0"
+pulp = "^2.8.0"
 population-restorator = "^0.2.3"
-dongraphio = "^0.3.13"
-provisio = "^0.1.7"
+iduedu = "^0.1.4"
 joblib = "^1.4.2"
-pandarallel = "^1.6.5"
 scikit-learn = "^1.4.0"
 
+
 [tool.poetry.group.dev.dependencies]
 black = "^24.2.0"
 pylint = "^3.0.3"
@@ -46,6 +47,11 @@ disable = [
     "duplicate-code",
     "missing-module-docstring",
     "missing-function-docstring",
+    "too-many-locals",
+    "too-many-branches",
+    "too-many-statements",
+    "too-many-arguments",
+    "cyclic-import"
 ]
 good-names = [
 
@@ -59,3 +65,4 @@ use_parentheses = true
 ensure_newline_before_comments = true
 line_length = 120
 split_on_trailing_comma = true
+skip = ["__init__.py"]

objectnat-0.2.1/src/objectnat/__init__.py

@@ -0,0 +1,13 @@
+"""
+ObjectNat
+========
+
+
+ObjectNat is an open-source library created for geospatial analysis created by IDU team.
+
+Homepage https://github.com/DDonnyy/ObjectNat.
+"""
+
+from ._config import config
+from ._api import *
+from ._version import VERSION as __version__

objectnat-0.1.5/src/objectnat/__init__.py → objectnat-0.2.1/src/objectnat/_api.py

@@ -1,16 +1,12 @@
-__version__ = "0.1.5"
+# pylint: disable=unused-import,wildcard-import,unused-wildcard-import
+from iduedu import *
 
-from dongraphio.enums import GraphType
-
-from .methods.adjacency_matrix import get_adjacency_matrix
 from .methods.balanced_buildings import get_balanced_buildings
 from .methods.cluster_points_in_polygons import get_clusters_polygon
 from .methods.coverage_zones import get_isochrone_zone_coverage, get_radius_zone_coverage
-from .methods.demands import get_demands
 from .methods.isochrones import get_accessibility_isochrones
 from .methods.living_buildings_osm import download_buildings
-from .methods.osm_graph import get_intermodal_graph_from_osm
-from .methods.provision import NoOsmIdException, NoWeightAdjacencyException, get_provision
+from .methods.provision.provision import clip_provision, get_service_provision, recalculate_links
 from .methods.visibility_analysis import (
     calculate_visibility_catchment_area,
     get_visibilities_from_points,

objectnat-0.2.1/src/objectnat/_config.py

@@ -0,0 +1,67 @@
+import sys
+from typing import Literal
+
+from iduedu import config as iduedu_config
+from loguru import logger
+
+
+class Config:
+    """
+    A configuration class to manage global settings for the application, such as Overpass API URL,
+    timeouts, and logging options.
+
+    Attributes
+    ----------
+    overpass_url : str
+        URL for accessing the Overpass API. Defaults to "http://lz4.overpass-api.de/api/interpreter".
+    timeout : int or None
+        Timeout in seconds for API requests. If None, no timeout is applied.
+    enable_tqdm_bar : bool
+        Enables or disables progress bars (via tqdm). Defaults to True.
+    logger : Logger
+        Logging instance to handle application logging.
+
+    Methods
+    -------
+    change_logger_lvl(lvl: Literal["TRACE", "DEBUG", "INFO", "WARN", "ERROR"])
+        Changes the logging level to the specified value.
+    set_overpass_url(url: str)
+        Sets a new Overpass API URL.
+    set_timeout(timeout: int)
+        Sets the timeout for API requests.
+    set_enable_tqdm(enable: bool)
+        Enables or disables progress bars in the application.
+    """
+
+    def __init__(
+        self,
+        overpass_url="http://lz4.overpass-api.de/api/interpreter",
+        timeout=None,
+        enable_tqdm_bar=True,
+    ):
+        self.overpass_url = overpass_url
+        self.timeout = timeout
+        self.enable_tqdm_bar = enable_tqdm_bar
+        self.logger = logger
+        self.iduedu_config = iduedu_config
+
+    def change_logger_lvl(self, lvl: Literal["TRACE", "DEBUG", "INFO", "WARN", "ERROR"]):
+        self.logger.remove()
+        self.logger.add(sys.stderr, level=lvl)
+        self.iduedu_config.change_logger_lvl(lvl)
+
+    def set_overpass_url(self, url: str):
+        self.overpass_url = url
+        self.iduedu_config.set_overpass_url(url)
+
+    def set_timeout(self, timeout: int):
+        self.timeout = timeout
+        self.iduedu_config.set_timeout(timeout)
+
+    def set_enable_tqdm(self, enable: bool):
+        self.enable_tqdm_bar = enable
+        self.iduedu_config.set_enable_tqdm(enable)
+
+
+config = Config()
+config.change_logger_lvl("INFO")

objectnat-0.2.1/src/objectnat/_version.py

@@ -0,0 +1 @@
+VERSION = "0.2.1"

{objectnat-0.1.5 → objectnat-0.2.1}/src/objectnat/methods/balanced_buildings.py

@@ -1,7 +1,10 @@
 import geopandas as gpd
 import population_restorator.balancer.houses as b_build
 import population_restorator.balancer.territories as b_terr
-from loguru import logger
+
+from objectnat import config
+
+logger = config.logger
 
 
 def get_balanced_buildings(

{objectnat-0.1.5 → objectnat-0.2.1}/src/objectnat/methods/cluster_points_in_polygons.py

@@ -2,9 +2,12 @@ from typing import Literal
 
 import geopandas as gpd
 import pandas as pd
-from loguru import logger
 from sklearn.cluster import DBSCAN, HDBSCAN
 
+from objectnat import config
+
+logger = config.logger
+
 
 def _get_cluster(services_select, min_dist, min_point, method):
     services_coords = pd.DataFrame(

{objectnat-0.1.5 → objectnat-0.2.1}/src/objectnat/methods/coverage_zones.py

@@ -2,8 +2,6 @@ from typing import Literal
 
 import geopandas as gpd
 import networkx as nx
-import pandas as pd
-from dongraphio import GraphType
 
 from .isochrones import get_accessibility_isochrones
 
@@ -48,7 +46,6 @@ def get_isochrone_zone_coverage(
     weight_type: Literal["time_min", "length_meter"],
     weight_value: int,
     city_graph: nx.Graph,
-    graph_type: list[GraphType],
 ) -> tuple[gpd.GeoDataFrame, gpd.GeoDataFrame | None, gpd.GeoDataFrame | None]:
     """
     Create isochrones for each service location based on travel time/distance.
@@ -56,15 +53,13 @@ def get_isochrone_zone_coverage(
     Parameters
     ----------
     services : gpd.GeoDataFrame
-        GeoDataFrame containing the service locations.
+        Layer containing the service locations.
     weight_type : str
         Type of weight used for calculating isochrones, either "time_min" or "length_meter".
     weight_value : int
         The value of the weight, representing time in minutes or distance in meters.
     city_graph : nx.Graph
         The graph representing the city's transportation network.
-    graph_type : list[GraphType]
-        List of graph types to be used for isochrone calculations.
 
     Returns
     -------
@@ -76,10 +71,10 @@ def get_isochrone_zone_coverage(
     >>> import networkx as nx
     >>> import geopandas as gpd
     >>> from shapely.geometry import Point
-    >>> from dongraphio import GraphType
+     >>> from iduedu import get_intermodal_graph
 
-    >>> # Create a sample city graph or download it from osm with get_intermodal_graph_from_osm()
-    >>> city_graph = nx.MultiDiGraph()
+    >>> # Create a sample city graph with get_intermodal_graph()
+    >>> graph = get_intermodal_graph(polygon=my_territory_polygon)
 
     >>> # Create a sample GeoDataFrame for services
     >>> services = gpd.read_file('services.geojson')
@@ -87,19 +82,9 @@ def get_isochrone_zone_coverage(
     >>> # Define parameters
     >>> weight_type = "time_min"
     >>> weight_value = 10
-    >>> graph_type = [GraphType.PUBLIC_TRANSPORT, GraphType.WALK]
 
     >>> # Get isochrone zone coverage
-    >>> isochrone_zones = get_isochrone_zone_coverage(services, weight_type, weight_value, city_graph, graph_type)
-    >>> isochrone_zones[0] # represent isochrones geodataframe
+    >>> isochrones, pt_stops, pt_routes = get_isochrone_zone_coverage(services, weight_type, weight_value, city_graph)
     """
-    assert services.crs == city_graph.graph["crs"], "CRS not match"
-    points = services.geometry.representative_point()
-    iso, routes, stops = get_accessibility_isochrones(
-        points, graph_type, weight_value, weight_type, city_graph, points.crs.to_epsg()
-    )
-    services_ = services.copy()
-    iso = gpd.GeoDataFrame(
-        pd.concat([iso.drop(columns=["point", "point_number"]), services_.drop(columns=["geometry"])], axis=1)
-    )
+    iso, routes, stops = get_accessibility_isochrones(services, weight_value, weight_type, city_graph)
     return iso, routes, stops

objectnat-0.2.1/src/objectnat/methods/isochrones.py

@@ -0,0 +1,140 @@
+from typing import Literal
+
+import geopandas as gpd
+import networkx as nx
+import numpy as np
+import pandas as pd
+from osmnx import graph_to_gdfs
+from pyproj import CRS
+from scipy.spatial import KDTree
+from shapely import Point
+from shapely.ops import unary_union
+
+from objectnat import config
+
+logger = config.logger
+
+
+def get_accessibility_isochrones(
+    points: gpd.GeoDataFrame,
+    weight_value: float,
+    weight_type: Literal["time_min", "length_meter"],
+    graph_nx: nx.Graph,
+) -> tuple[gpd.GeoDataFrame, gpd.GeoDataFrame | None, gpd.GeoDataFrame | None]:
+    """
+    Calculate accessibility isochrones from a gpd.GeoDataFrame based on the provided city graph.
+
+    Isochrones represent areas that can be reached from a given point within a specific time or distance,
+    using a graph that contains road and transport network data.
+
+    Parameters
+    ----------
+    points : gpd.GeoDataFrame
+        A GeoDataFrame containing the geometry from which accessibility isochrones should be calculated.
+        The CRS of this GeoDataFrame must match the CRS of the provided graph.
+    weight_value : float
+        The maximum distance or time threshold for calculating isochrones.
+    weight_type : Literal["time_min", "length_meter"]
+        The type of weight to use for distance calculations. Either time in minutes ("time_min") or distance
+        in meters ("length_meter").
+    graph_nx : nx.Graph
+        A NetworkX graph representing the city network.
+        The graph must contain the appropriate CRS and, for time-based isochrones, a speed attribute.
+
+    Returns
+    -------
+    tuple[gpd.GeoDataFrame, gpd.GeoDataFrame | None, gpd.GeoDataFrame | None]
+        A tuple containing:
+        - isochrones : GeoDataFrame with the calculated isochrone geometries.
+        - public transport stops (if applicable) : GeoDataFrame with public transport stops within the isochrone, or None if not applicable.
+        - public transport routes (if applicable) : GeoDataFrame with public transport routes within the isochrone, or None if not applicable.
+
+    Examples
+    --------
+    >>> from iduedu import get_intermodal_graph
+    >>> graph = get_intermodal_graph(polygon=my_territory_polygon)
+    >>> points = gpd.GeoDataFrame(geometry=[Point(30.33, 59.95)], crs=4326).to_crs(graph.graph['crs'])
+    >>> isochrones, pt_stops, pt_routes = get_accessibility_isochrones(points, weight_value=15, weight_type="time_min", graph_nx=my_graph)
+
+    """
+
+    assert points.crs == CRS.from_epsg(
+        graph_nx.graph["crs"]
+    ), f'CRS mismatch , points.crs = {points.crs.to_epsg()}, graph["crs"] = {graph_nx.graph["crs"]}'
+
+    nodes_with_data = list(graph_nx.nodes(data=True))
+    logger.info("Calculating isochrones distances...")
+    coordinates = np.array([(data["x"], data["y"]) for node, data in nodes_with_data])
+    tree = KDTree(coordinates)
+
+    target_coord = [(p.x, p.y) for p in points.representative_point()]
+    distances, indices = tree.query(target_coord)
+
+    nearest_nodes = [nodes_with_data[idx][0] for idx in indices]
+    del nodes_with_data
+    dist_nearest = pd.DataFrame(data=distances, index=nearest_nodes, columns=["dist"])
+    speed = 0
+    if graph_nx.graph["type"] in ["walk", "intermodal"] and weight_type == "time_min":
+        try:
+            speed = graph_nx.graph["walk_speed"]
+        except KeyError:
+            logger.warning("There is no walk_speed in graph, set to the default speed - 83.33 m/min")
+            speed = 83.33
+        dist_nearest = dist_nearest / speed
+    elif weight_type == "time_min":
+        speed = 20 * 1000 / 60
+        dist_nearest = dist_nearest / speed
+
+    if (dist_nearest > weight_value).all().all():
+        raise RuntimeError(
+            "The point(s) lie further from the graph than weight_value, it's impossible to "
+            "construct isochrones. Check the coordinates of the point(s)/their projection"
+        )
+
+    data = []
+    for source in nearest_nodes:
+        dist, path = nx.single_source_dijkstra(graph_nx, source, weight=weight_type, cutoff=weight_value)
+        for target_node, way in path.items():
+            source = way[0]
+            distance = dist.get(target_node, np.nan)
+            data.append((source, target_node, distance))
+    del dist
+    dist_matrix = pd.DataFrame(data, columns=["source", "destination", "distance"])
+    del data
+    dist_matrix = dist_matrix.pivot_table(index="source", columns="destination", values="distance", sort=False)
+
+    dist_matrix = dist_matrix.add(dist_nearest.dist, axis=0)
+    dist_matrix = dist_matrix.mask(dist_matrix >= weight_value, np.nan)
+    dist_matrix.dropna(how="all", inplace=True)
+
+    results = []
+    logger.info("Building isochrones geometry...")
+    for _, row in dist_matrix.iterrows():
+        geometry = []
+        for node_to, value in row.items():
+            if not pd.isna(value):
+                node = graph_nx.nodes[node_to]
+                point = Point(node["x"], node["y"])
+                geometry.append(
+                    point.buffer(round((weight_value - value) * speed * 0.8, 2))
+                    if weight_type == "time_min"
+                    else point.buffer(round((weight_value - value) * 0.8, 2))
+                )
+        geometry = unary_union(geometry)
+        results.append(geometry)
+
+    isochrones = gpd.GeoDataFrame(data=points, geometry=results, crs=graph_nx.graph["crs"])
+    isochrones["weight_type"] = weight_type
+    isochrones["weight_value"] = weight_value
+
+    isochrones_subgraph = graph_nx.subgraph(dist_matrix.columns)
+    nodes = pd.DataFrame.from_dict(dict(isochrones_subgraph.nodes(data=True)), orient="index")
+    if "desc" in nodes.columns and "stop" in nodes["desc"].unique():
+        pt_nodes = nodes[nodes["desc"] == "stop"]
+        nodes, edges = graph_to_gdfs(isochrones_subgraph.subgraph(pt_nodes.index))
+        nodes.reset_index(drop=True, inplace=True)
+        nodes = nodes[["desc", "route", "geometry"]]
+        edges.reset_index(drop=True, inplace=True)
+        edges = edges[["type", "route", "geometry"]]
+        return isochrones, nodes, edges
+    return isochrones, None, None