building-eload 0.4.0__py3-none-any.whl

building_eload/__init__.py

@@ -0,0 +1,38 @@
+from pathlib import Path
+
+# Utilisation de Path(__file__).parents[1] pour naviguer vers le bon répertoire parent
+data_path = {
+    "data": Path(__file__).parents[1] / "data",
+    "bdtopo": Path(__file__).parents[1] / "data" / "bdtopo",
+    "elmas": Path(__file__).parents[1] / "data" / "elmas",
+    "simulation": Path(__file__).parents[1] / "data" / "simulation",
+    "activity_calendar": Path(__file__).parents[1] / "data" / "activity_calendar",
+    "validation": Path(__file__).parents[1] / "data" / "validation",
+    "climate": Path(__file__).parents[1] / "data" / "climate",
+    "representative_districts": Path(__file__).parents[1] / "data" / "representative_districts",
+}
+
+plot_path = Path(__file__).parents[1] / "plot"
+
+
+download_link = {
+    "bd_topo_11": "https://data.geopf.fr/telechargement/download/BDTOPO/BDTOPO_3-4_TOUSTHEMES_GPKG_LAMB93_R11_2024-09-15/BDTOPO_3-4_TOUSTHEMES_GPKG_LAMB93_R11_2024-09-15.7z",
+    "bd_topo_24": "https://data.geopf.fr/telechargement/download/BDTOPO/BDTOPO_3-4_TOUSTHEMES_GPKG_LAMB93_R24_2024-12-15/BDTOPO_3-4_TOUSTHEMES_GPKG_LAMB93_R24_2024-12-15.7z",
+    "bd_topo_27": "https://data.geopf.fr/telechargement/download/BDTOPO/BDTOPO_3-4_TOUSTHEMES_GPKG_LAMB93_R27_2025-03-15/BDTOPO_3-4_TOUSTHEMES_GPKG_LAMB93_R27_2025-03-15.7z",
+    "bd_topo_28": "https://data.geopf.fr/telechargement/download/BDTOPO/BDTOPO_3-4_TOUSTHEMES_GPKG_LAMB93_R28_2025-03-15/BDTOPO_3-4_TOUSTHEMES_GPKG_LAMB93_R28_2025-03-15.7z",
+    "bd_topo_32": "https://data.geopf.fr/telechargement/download/BDTOPO/BDTOPO_3-4_TOUSTHEMES_GPKG_LAMB93_R32_2025-03-15/BDTOPO_3-4_TOUSTHEMES_GPKG_LAMB93_R32_2025-03-15.7z",
+    "bd_topo_44": "https://data.geopf.fr/telechargement/download/BDTOPO/BDTOPO_3-4_TOUSTHEMES_GPKG_LAMB93_R44_2025-03-15/BDTOPO_3-4_TOUSTHEMES_GPKG_LAMB93_R44_2025-03-15.7z",
+    "bd_topo_52": "https://data.geopf.fr/telechargement/download/BDTOPO/BDTOPO_3-4_TOUSTHEMES_GPKG_LAMB93_R52_2025-03-15/BDTOPO_3-4_TOUSTHEMES_GPKG_LAMB93_R52_2025-03-15.7z",
+    "bd_topo_53": "https://data.geopf.fr/telechargement/download/BDTOPO/BDTOPO_3-4_TOUSTHEMES_GPKG_LAMB93_R53_2025-03-15/BDTOPO_3-4_TOUSTHEMES_GPKG_LAMB93_R53_2025-03-15.7z",
+    "bd_topo_75": "https://data.geopf.fr/telechargement/download/BDTOPO/BDTOPO_3-4_TOUSTHEMES_GPKG_LAMB93_R75_2025-03-15/BDTOPO_3-4_TOUSTHEMES_GPKG_LAMB93_R75_2025-03-15.7z",
+    "bd_topo_76": "https://data.geopf.fr/telechargement/download/BDTOPO/BDTOPO_3-4_TOUSTHEMES_GPKG_LAMB93_R76_2025-03-15/BDTOPO_3-4_TOUSTHEMES_GPKG_LAMB93_R76_2025-03-15.7z",
+    "bd_topo_84": "https://data.geopf.fr/telechargement/download/BDTOPO/BDTOPO_3-4_TOUSTHEMES_GPKG_LAMB93_R84_2025-03-15/BDTOPO_3-4_TOUSTHEMES_GPKG_LAMB93_R84_2025-03-15.7z",
+    "bd_topo_93": "https://data.geopf.fr/telechargement/download/BDTOPO/BDTOPO_3-4_TOUSTHEMES_GPKG_LAMB93_R93_2025-03-15/BDTOPO_3-4_TOUSTHEMES_GPKG_LAMB93_R93_2025-03-15.7z",
+    "bd_topo_94": "https://data.geopf.fr/telechargement/download/BDTOPO/BDTOPO_3-4_TOUSTHEMES_GPKG_LAMB93_R94_2025-03-15/BDTOPO_3-4_TOUSTHEMES_GPKG_LAMB93_R94_2025-03-15.7z",
+    "district_list": "https://drive.google.com/uc?export=download&id=1YbTrxEMXHyTfcblaDWeBTbZMilL4QZP7",
+    "elmas": "https://figshare.com/ndownloader/files/41895786",
+    "occupant_diaries": "https://drive.google.com/uc?export=download&id=1Q_QEOoAvtA0T_70EpA4T2JpDh-VCFh1B",
+    "ORE": "https://drive.google.com/uc?export=download&id=1IAmwRaE_HxVmkik49CStng0vlXMigjmg",
+    "ENEDIS_REGIONAL": "TODO",
+    "ENEDIS_NATIONAL": "TODO",
+}

building_eload/core/building_loader.py

@@ -0,0 +1,127 @@
+# building_loader.py
+
+import logging
+import warnings
+from pathlib import Path
+from typing import Optional
+
+import buildingdata
+import geopandas as gpd
+from shapely import geometry
+
+from building_eload import data_path
+
+logger = logging.getLogger(__name__)
+
+
+
+class BuildingLoader:
+    """Load buildings from a district and filter them based on intersection with the district geometry."""
+
+    INTERSECTION_THRESHOLD = 0.5
+    BUILDING_LAYER_NAME = "batiment"
+    GIS_ENGINE = "pyogrio"
+    BUILDING_FILE_PREFIX = "BDTOPO_R"
+    BUILDING_FILE_SUFFIX = ".gpkg"
+    DEFAULT_DISTRICT_DATA_FILE = None
+    DEFAULT_BUILDING_FOOTPRINT_FOLDER = Path(data_path["bdtopo"])
+
+    def __init__(self, district_id: str, district_data_file: str = None, building_footprint_folder: str = None):
+
+        self.district_id = str(district_id)
+        self.district_data_file = (
+            Path(district_data_file) if district_data_file is not None else self.DEFAULT_DISTRICT_DATA_FILE
+        )
+        self.building_footprint_folder = Path(
+            building_footprint_folder or self.DEFAULT_BUILDING_FOOTPRINT_FOLDER
+        )
+        logger.debug(f"{self.logger_header} District data file: {self.district_data_file}.")
+        logger.debug(f"{self.logger_header} Building footprint folder: {self.building_footprint_folder}.")
+
+    @property
+    def logger_header(self):
+        if self.district_id is None:
+            return "[BuildingLoader]"
+        return f"[District {self.district_id}]"
+
+    # -------------------------------------------------------
+    # Main function
+    # -------------------------------------------------------
+
+    def load_buildings(self) -> Optional[gpd.GeoDataFrame]:
+        """Load buildings for a specific district."""
+        logger.info(f"{self.logger_header} Loading buildings.")
+
+        district_data = self._load_district(self.district_id)
+        if district_data.empty:
+            logger.warning(f"{self.logger_header} District {self.district_id} not found in district file.")
+            return None
+        district_geometry = self._get_district_geometry(district_data)
+        buildings = self._load_buildings_in_district_bounds(district_geometry, district_data)
+        filtered_buildings = self._filter_buildings_on_bounds(buildings, district_geometry)
+        filtered_buildings["district"] = self.district_id
+
+        
+        logger.info(f"{self.logger_header} Loaded {len(filtered_buildings)} buildings.")
+
+        return filtered_buildings
+
+    # -------------------------------------------------------
+    # Utilis
+    # -------------------------------------------------------
+
+    def _load_district(self, district_id: str) -> gpd.GeoDataFrame:
+        district_id = str(district_id)
+        if self.district_data_file is not None:
+            districts = gpd.read_parquet(self.district_data_file)
+        else:
+            districts = buildingdata.get_districts()
+
+        district_codes = districts["district"].astype(str).str.zfill(len(district_id))
+        return districts[district_codes == district_id]
+
+    @staticmethod
+    def _get_district_geometry(district_data: gpd.GeoDataFrame) -> geometry.Polygon:
+        return district_data["geometry"].iloc[0]
+
+    def _load_buildings_in_district_bounds(
+        self, district_geometry: geometry.Polygon, district_data: gpd.GeoDataFrame
+    ) -> gpd.GeoDataFrame:
+        bounds = district_geometry.convex_hull.bounds
+        building_path = self._get_building_path(district_data)
+
+        if not building_path.exists():
+            logger.debug(
+                f"{self.logger_header} Local BDTOPO file {building_path} not found; "
+                "loading buildings through buildingdata."
+            )
+            return buildingdata.get_bdtopo(self.district_id)
+
+        with warnings.catch_warnings():
+            warnings.filterwarnings("ignore")
+            logger.debug(f"{self.logger_header} Loading buildings from {building_path}.")
+            return gpd.read_file(
+                building_path,
+                layer=self.BUILDING_LAYER_NAME,
+                engine=self.GIS_ENGINE,
+                bbox=bounds,
+            )
+
+    def _filter_buildings_on_bounds(
+        self, buildings: gpd.GeoDataFrame, district_geometry: geometry.Polygon
+    ) -> gpd.GeoDataFrame:
+        buildings["intersection_share"] = (
+            buildings["geometry"].intersection(district_geometry).area / buildings["geometry"].area
+        )
+        filtered_buildings = buildings.loc[
+            buildings.intersection_share >= self.INTERSECTION_THRESHOLD
+        ].copy()
+        filtered_buildings.drop(columns=["intersection_share"], inplace=True)
+        return filtered_buildings
+
+    def _get_building_path(self, district_data: gpd.GeoDataFrame) -> Path:
+        region_id = int(district_data["region"].iloc[0])
+        return (
+            self.building_footprint_folder
+            / f"{self.BUILDING_FILE_PREFIX}{region_id}{self.BUILDING_FILE_SUFFIX}"
+        )