ssb-sgis 0.2.7__tar.gz → 0.2.8__tar.gz

{ssb_sgis-0.2.7 → ssb_sgis-0.2.8}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: ssb-sgis
-Version: 0.2.7
+Version: 0.2.8
 Summary: GIS functions used at Statistics Norway.
 Home-page: https://github.com/statisticsnorway/ssb-sgis
 License: MIT

{ssb_sgis-0.2.7 → ssb_sgis-0.2.8}/pyproject.toml

@@ -1,6 +1,6 @@
 [tool.poetry]
 name = "ssb-sgis"
-version = "0.2.7"
+version = "0.2.8"
 description = "GIS functions used at Statistics Norway."
 authors = ["Statistics Norway <ort@ssb.no>"]
 license = "MIT"

{ssb_sgis-0.2.7 → ssb_sgis-0.2.8}/src/sgis/__init__.py

@@ -3,7 +3,9 @@ from .geopandas_tools.buffer_dissolve_explode import (
     buff,
     buffdiss,
     buffdissexp,
+    buffdissexp_by_cluster,
     dissexp,
+    dissexp_by_cluster,
 )
 from .geopandas_tools.general import (
     bounds_to_points,
@@ -23,6 +25,7 @@ from .geopandas_tools.general import (
 from .geopandas_tools.geometry_types import (
     get_geom_type,
     is_single_geom_type,
+    make_all_singlepart,
     to_single_geom_type,
 )
 from .geopandas_tools.neighbors import (

{ssb_sgis-0.2.7 → ssb_sgis-0.2.8}/src/sgis/geopandas_tools/buffer_dissolve_explode.py

@@ -17,6 +17,7 @@ for the following:
 from geopandas import GeoDataFrame, GeoSeries
 
 from .geometry_types import make_all_singlepart
+from .polygon_operations import get_polygon_clusters
 
 
 def _decide_ignore_index(kwargs: dict) -> tuple[dict, bool]:
@@ -190,6 +191,31 @@ def dissexp(
     )
 
 
+def dissexp_by_cluster(gdf: GeoDataFrame) -> GeoDataFrame:
+    return (
+        gdf.explode(ignore_index=True)
+        .pipe(get_polygon_clusters, cluster_col="cluster")
+        .pipe(dissexp, by="cluster")
+        .reset_index(drop=True)
+    )
+
+
+def buffdissexp_by_cluster(
+    gdf: GeoDataFrame,
+    distance: int | float,
+    *,
+    resolution: int = 50,
+    copy: bool = True,
+) -> GeoDataFrame:
+    return (
+        buff(gdf, distance, resolution=resolution, copy=copy)
+        .explode(ignore_index=True)
+        .pipe(get_polygon_clusters, cluster_col="cluster")
+        .pipe(dissexp, by="cluster")
+        .reset_index(drop=True)
+    )
+
+
 def buff(
     gdf: GeoDataFrame | GeoSeries,
     distance: int | float,

{ssb_sgis-0.2.7 → ssb_sgis-0.2.8}/src/sgis/geopandas_tools/general.py

@@ -2,6 +2,7 @@ import warnings
 
 import geopandas as gpd
 import numpy as np
+import pandas as pd
 from geopandas import GeoDataFrame, GeoSeries
 from geopandas.array import GeometryDtype
 from shapely import (
@@ -519,11 +520,17 @@ def to_lines(*gdfs: GeoDataFrame, copy: bool = True) -> GeoDataFrame:
     if len(lines) == 1:
         return lines[0]
 
-    unioned = lines[0].overlay(lines[1], how="union", keep_geom_type=True)
+    if len(lines[0]) and len(lines[1]):
+        unioned = lines[0].overlay(lines[1], how="union", keep_geom_type=True)
+    else:
+        unioned = pd.concat([lines[0], lines[1]], ignore_index=True)
 
     if len(lines) > 2:
         for line_gdf in lines[2:]:
-            unioned = unioned.overlay(line_gdf, how="union", keep_geom_type=True)
+            if len(line_gdf):
+                unioned = unioned.overlay(line_gdf, how="union", keep_geom_type=True)
+            else:
+                unioned = pd.concat([unioned, line_gdf], ignore_index=True)
 
     return make_all_singlepart(unioned, ignore_index=True)
 

{ssb_sgis-0.2.7 → ssb_sgis-0.2.8}/src/sgis/geopandas_tools/polygon_operations.py

@@ -23,6 +23,147 @@ from .neighbors import get_neighbor_indices
 from .overlay import clean_overlay
 
 
+def get_polygon_clusters(
+    *gdfs: GeoDataFrame | GeoSeries,
+    cluster_col: str = "cluster",
+    allow_multipart: bool = False,
+) -> GeoDataFrame | tuple[GeoDataFrame]:
+    """Find which polygons overlap without dissolving.
+
+    Devides polygons into clusters in a fast and precice manner by using spatial join
+    and networkx to find the connected components, i.e. overlapping geometries.
+    If multiple GeoDataFrames are given, the clusters will be based on all
+    combined.
+
+    This can be used instead of dissolve+explode, or before dissolving by the cluster
+    column. This has been tested to be a lot faster if there are many
+    non-overlapping polygons, but somewhat slower than dissolve+explode if most
+    polygons overlap.
+
+    Args:
+        gdfs: One or more GeoDataFrames of polygons.
+        cluster_col: Name of the resulting cluster column.
+        allow_multipart: Whether to allow mutipart geometries in the gdfs.
+            Defaults to False to avoid confusing results.
+
+    Returns:
+        One or more GeoDataFrames (same amount as was given) with a new cluster column.
+
+    Examples
+    --------
+
+    Create geometries with three clusters of overlapping polygons.
+
+    >>> import sgis as sg
+    >>> gdf = sg.to_gdf([(0, 0), (1, 1), (0, 1), (4, 4), (4, 3), (7, 7)])
+    >>> buffered = sg.buff(gdf, 1)
+    >>> gdf
+                                                geometry
+    0  POLYGON ((1.00000 0.00000, 0.99951 -0.03141, 0...
+    1  POLYGON ((2.00000 1.00000, 1.99951 0.96859, 1....
+    2  POLYGON ((1.00000 1.00000, 0.99951 0.96859, 0....
+    3  POLYGON ((5.00000 4.00000, 4.99951 3.96859, 4....
+    4  POLYGON ((5.00000 3.00000, 4.99951 2.96859, 4....
+    5  POLYGON ((8.00000 7.00000, 7.99951 6.96859, 7....
+
+    Add a cluster column to the GeoDataFrame:
+
+    >>> gdf = sg.get_polygon_clusters(gdf, cluster_col="cluster")
+    >>> gdf
+       cluster                                           geometry
+    0        0  POLYGON ((1.00000 0.00000, 0.99951 -0.03141, 0...
+    1        0  POLYGON ((2.00000 1.00000, 1.99951 0.96859, 1....
+    2        0  POLYGON ((1.00000 1.00000, 0.99951 0.96859, 0....
+    3        1  POLYGON ((5.00000 4.00000, 4.99951 3.96859, 4....
+    4        1  POLYGON ((5.00000 3.00000, 4.99951 2.96859, 4....
+    5        2  POLYGON ((8.00000 7.00000, 7.99951 6.96859, 7....
+
+    If multiple GeoDataFrames are given, all are returned with common
+    cluster values.
+
+    >>> gdf2 = sg.to_gdf([(0, 0), (7, 7)])
+    >>> gdf, gdf2 = sg.get_polygon_clusters(gdf, gdf2, cluster_col="cluster")
+    >>> gdf2
+       cluster                 geometry
+    0        0  POINT (0.00000 0.00000)
+    1        2  POINT (7.00000 7.00000)
+    >>> gdf
+       cluster                                           geometry
+    0        0  POLYGON ((1.00000 0.00000, 0.99951 -0.03141, 0...
+    1        0  POLYGON ((2.00000 1.00000, 1.99951 0.96859, 1....
+    2        0  POLYGON ((1.00000 1.00000, 0.99951 0.96859, 0....
+    3        1  POLYGON ((5.00000 4.00000, 4.99951 3.96859, 4....
+    4        1  POLYGON ((5.00000 3.00000, 4.99951 2.96859, 4....
+    5        2  POLYGON ((8.00000 7.00000, 7.99951 6.96859, 7....
+
+    Dissolving 'by' the cluster column will make the dissolve much
+    faster if there are a lot of non-overlapping polygons.
+
+    >>> dissolved = gdf.dissolve(by="cluster", as_index=False)
+    >>> dissolved
+       cluster                                           geometry
+    0        0  POLYGON ((0.99951 -0.03141, 0.99803 -0.06279, ...
+    1        1  POLYGON ((4.99951 2.96859, 4.99803 2.93721, 4....
+    2        2  POLYGON ((8.00000 7.00000, 7.99951 6.96859, 7....
+    """
+    if isinstance(gdfs[-1], str):
+        *gdfs, cluster_col = gdfs
+
+    concated = pd.DataFrame()
+    orig_indices = ()
+    for i, gdf in enumerate(gdfs):
+        if isinstance(gdf, GeoSeries):
+            gdf = gdf.to_frame()
+
+        if not isinstance(gdf, GeoDataFrame):
+            raise TypeError("'gdfs' should be GeoDataFrames or GeoSeries.")
+
+        if not allow_multipart and len(gdf) != len(gdf.explode(index_parts=False)):
+            raise ValueError(
+                "All geometries should be exploded to singlepart "
+                "in order to get correct polygon clusters. "
+                "To allow multipart geometries, set allow_multipart=True"
+            )
+
+        orig_indices = orig_indices + (gdf.index,)
+
+        gdf["i__"] = i
+
+        concated = pd.concat([concated, gdf], ignore_index=True)
+
+    neighbors = get_neighbor_indices(concated, concated)
+
+    edges = [(source, target) for source, target in neighbors.items()]
+
+    graph = nx.Graph()
+    graph.add_edges_from(edges)
+
+    component_mapper = {
+        j: i
+        for i, component in enumerate(nx.connected_components(graph))
+        for j in component
+    }
+
+    concated[cluster_col] = component_mapper
+
+    concated = _push_geom_col(concated)
+
+    n_gdfs = concated["i__"].unique()
+
+    if len(n_gdfs) == 1:
+        concated.index = orig_indices[0]
+        return concated.drop(["i__"], axis=1)
+
+    unconcated = ()
+    for i in n_gdfs:
+        gdf = concated[concated["i__"] == i]
+        gdf.index = orig_indices[i]
+        gdf = gdf.drop(["i__"], axis=1)
+        unconcated = unconcated + (gdf,)
+
+    return unconcated
+
+
 def eliminate_by_longest(
     gdf: GeoDataFrame,
     to_eliminate: GeoDataFrame,
@@ -215,146 +356,6 @@ def _eliminate_by_area(
     return eliminated
 
 
-def get_polygon_clusters(
-    *gdfs: GeoDataFrame | GeoSeries,
-    cluster_col: str = "cluster",
-    allow_multipart: bool = False,
-) -> GeoDataFrame | tuple[GeoDataFrame]:
-    """Find which polygons overlap without dissolving.
-
-    Devides polygons into clusters in a fast and precice manner by using spatial join
-    and networkx to find the connected components, i.e. overlapping geometries.
-    If multiple GeoDataFrames are given, the clusters will be based on all
-    combined.
-
-    This can be used instead of dissolve+explode, or before dissolving by the cluster
-    column. This has been tested to be a lot faster if there are many
-    non-overlapping polygons, but somewhat slower than dissolve+explode if most
-    polygons overlap.
-
-    Args:
-        gdfs: One or more GeoDataFrames of polygons.
-        cluster_col: Name of the resulting cluster column.
-        allow_multipart: Whether to allow mutipart geometries in the gdfs.
-            Defaults to False to avoid confusing results.
-
-    Returns:
-        One or more GeoDataFrames (same amount as was given) with a new cluster column.
-
-    Examples
-    --------
-
-    Create geometries with three clusters of overlapping polygons.
-
-    >>> import sgis as sg
-    >>> gdf = sg.to_gdf([(0, 0), (1, 1), (0, 1), (4, 4), (4, 3), (7, 7)])
-    >>> buffered = sg.buff(gdf, 1)
-    >>> gdf
-                                                geometry
-    0  POLYGON ((1.00000 0.00000, 0.99951 -0.03141, 0...
-    1  POLYGON ((2.00000 1.00000, 1.99951 0.96859, 1....
-    2  POLYGON ((1.00000 1.00000, 0.99951 0.96859, 0....
-    3  POLYGON ((5.00000 4.00000, 4.99951 3.96859, 4....
-    4  POLYGON ((5.00000 3.00000, 4.99951 2.96859, 4....
-    5  POLYGON ((8.00000 7.00000, 7.99951 6.96859, 7....
-
-    Add a cluster column to the GeoDataFrame:
-
-    >>> gdf = sg.get_polygon_clusters(gdf, cluster_col="cluster")
-    >>> gdf
-       cluster                                           geometry
-    0        0  POLYGON ((1.00000 0.00000, 0.99951 -0.03141, 0...
-    1        0  POLYGON ((2.00000 1.00000, 1.99951 0.96859, 1....
-    2        0  POLYGON ((1.00000 1.00000, 0.99951 0.96859, 0....
-    3        1  POLYGON ((5.00000 4.00000, 4.99951 3.96859, 4....
-    4        1  POLYGON ((5.00000 3.00000, 4.99951 2.96859, 4....
-    5        2  POLYGON ((8.00000 7.00000, 7.99951 6.96859, 7....
-
-    If multiple GeoDataFrames are given, all are returned with common
-    cluster values.
-
-    >>> gdf2 = sg.to_gdf([(0, 0), (7, 7)])
-    >>> gdf, gdf2 = sg.get_polygon_clusters(gdf, gdf2, cluster_col="cluster")
-    >>> gdf2
-       cluster                 geometry
-    0        0  POINT (0.00000 0.00000)
-    1        2  POINT (7.00000 7.00000)
-    >>> gdf
-       cluster                                           geometry
-    0        0  POLYGON ((1.00000 0.00000, 0.99951 -0.03141, 0...
-    1        0  POLYGON ((2.00000 1.00000, 1.99951 0.96859, 1....
-    2        0  POLYGON ((1.00000 1.00000, 0.99951 0.96859, 0....
-    3        1  POLYGON ((5.00000 4.00000, 4.99951 3.96859, 4....
-    4        1  POLYGON ((5.00000 3.00000, 4.99951 2.96859, 4....
-    5        2  POLYGON ((8.00000 7.00000, 7.99951 6.96859, 7....
-
-    Dissolving 'by' the cluster column will make the dissolve much
-    faster if there are a lot of non-overlapping polygons.
-
-    >>> dissolved = gdf.dissolve(by="cluster", as_index=False)
-    >>> dissolved
-       cluster                                           geometry
-    0        0  POLYGON ((0.99951 -0.03141, 0.99803 -0.06279, ...
-    1        1  POLYGON ((4.99951 2.96859, 4.99803 2.93721, 4....
-    2        2  POLYGON ((8.00000 7.00000, 7.99951 6.96859, 7....
-    """
-    if isinstance(gdfs[-1], str):
-        *gdfs, cluster_col = gdfs
-
-    concated = pd.DataFrame()
-    orig_indices = ()
-    for i, gdf in enumerate(gdfs):
-        if isinstance(gdf, GeoSeries):
-            gdf = gdf.to_frame()
-
-        if not isinstance(gdf, GeoDataFrame):
-            raise TypeError("'gdfs' should be one or more GeoDataFrames or GeoSeries.")
-
-        if not allow_multipart and len(gdf) != len(gdf.explode(index_parts=False)):
-            raise ValueError(
-                "All geometries should be exploded to singlepart "
-                "in order to get correct polygon clusters. "
-                "To allow multipart geometries, set allow_multipart=True"
-            )
-
-        orig_indices = orig_indices + (gdf.index,)
-        gdf["i__"] = i
-
-        concated = pd.concat([concated, gdf], ignore_index=True)
-
-    neighbors = get_neighbor_indices(concated, concated)
-
-    edges = [(source, target) for source, target in neighbors.items()]
-
-    graph = nx.Graph()
-    graph.add_edges_from(edges)
-
-    component_mapper = {
-        j: i
-        for i, component in enumerate(nx.connected_components(graph))
-        for j in component
-    }
-
-    concated[cluster_col] = component_mapper
-
-    concated = _push_geom_col(concated)
-
-    n_gdfs = concated["i__"].unique()
-
-    if len(n_gdfs) == 1:
-        concated.index = orig_indices[0]
-        return concated.drop(["i__"], axis=1)
-
-    unconcated = ()
-    for i in n_gdfs:
-        gdf = concated[concated["i__"] == i]
-        gdf.index = orig_indices[i]
-        gdf = gdf.drop(["i__"], axis=1)
-        unconcated = unconcated + (gdf,)
-
-    return unconcated
-
-
 def get_overlapping_polygons(
     gdf: GeoDataFrame | GeoSeries, ignore_index: bool = False
 ) -> GeoDataFrame | GeoSeries:

{ssb_sgis-0.2.7 → ssb_sgis-0.2.8}/src/sgis/geopandas_tools/to_geodataframe.py

@@ -3,6 +3,7 @@ from collections.abc import Iterator, Sized
 
 import geopandas as gpd
 import pandas as pd
+import shapely
 from geopandas import GeoDataFrame, GeoSeries
 from pandas.api.types import is_array_like, is_dict_like, is_list_like
 from shapely import Geometry, box, wkb, wkt
@@ -128,7 +129,7 @@ def to_gdf(
         raise TypeError("'to_gdf' doesn't accept GeoDataFrames as input type.")
 
     if isinstance(geom, GeoSeries):
-        geom_col = "geometry" if not geometry else geometry
+        geom_col = geometry if geometry else "geometry"
         return _geoseries_to_gdf(geom, geom_col, crs, **kwargs)
 
     geom_col = _find_geometry_column(geom, geometry)
@@ -144,8 +145,26 @@ def to_gdf(
         return GeoDataFrame({geom_col: geom}, geometry=geom_col, crs=crs, **kwargs)
 
     if not is_dict_like(geom):
-        geom = GeoSeries(_make_shapely_geoms(geom), index=index)
-        return GeoDataFrame({geom_col: geom}, geometry=geom_col, crs=crs, **kwargs)
+        if not hasattr(geom, "__iter__") and hasattr(geom, "__dict__"):
+            if all(attr in geom.__dict__ for attr in ["minx", "miny", "maxx", "maxy"]):
+                geom = GeoSeries(
+                    shapely.box(*(geom.minx, geom.miny, geom.maxx, geom.maxy)),
+                    index=index,
+                )
+                return GeoDataFrame(
+                    {geom_col: geom}, geometry=geom_col, crs=crs, **kwargs
+                )
+        if hasattr(geom, "__iter__") and all(isinstance(g, dict) for g in geom):
+            crs = crs if crs else _get_crs(geom)
+            geom = pd.concat(GeoSeries(_from_json(g)) for g in geom)
+            if index is not None:
+                geom.index = index
+            else:
+                geom = geom.reset_index(drop=True)
+            return GeoDataFrame({geom_col: geom}, geometry=geom_col, crs=crs, **kwargs)
+        else:
+            geom = GeoSeries(_make_shapely_geoms(geom), index=index)
+            return GeoDataFrame({geom_col: geom}, geometry=geom_col, crs=crs, **kwargs)
 
     # now we have dict, Series or DataFrame
 
@@ -177,10 +196,61 @@ def to_gdf(
     if geometry and geom_col not in geom or isinstance(geom, pd.DataFrame):
         raise ValueError("Cannot find geometry column(s)", geometry)
 
+    # geojson, __geo_interface__
+    if (
+        isinstance(geom, dict)
+        and sum(key in geom for key in ["type", "coordinates", "features"]) >= 2
+    ):
+        if "geometry" in geom:
+            geometry = "geometry"
+
+        crs = crs if crs else _get_crs(geom)
+        print(crs)
+        geom = GeoSeries(_from_json(geom), index=index)
+        return GeoDataFrame({geom_col: geom}, geometry=geom_col, crs=crs, **kwargs)
+
     geoseries = _series_like_to_geoseries(geom, index=index)
     return GeoDataFrame(geometry=geoseries, crs=crs, **kwargs)
 
 
+def _get_crs(geom):
+    if not is_dict_like(geom) and is_dict_like(geom[0]):
+        crss = list({_get_crs(g) for g in geom})
+        if len(crss) == 1:
+            return crss[0]
+        return None
+    if "properties" in geom:
+        return _get_crs(geom["properties"])
+    if "crs" in geom:
+        geom = geom["crs"]
+        while is_dict_like(geom):
+            if "properties" in geom:
+                geom = geom["properties"]
+            elif "name" in geom:
+                geom = geom["name"]
+            else:
+                return None
+        return geom
+    return None
+
+
+def _from_json(geom: dict):
+    if not isinstance(geom, dict) and isinstance(geom[0], dict):
+        return [_from_json(g) for g in geom]
+    if "geometry" in geom:
+        return _from_json(geom["geometry"])
+    if "features" in geom:
+        return _from_json(geom["features"])
+    coords = geom["coordinates"]
+    constructor = eval("shapely.geometry." + geom.get("type", Point))
+    try:
+        return constructor(coords)
+    except TypeError:
+        while len(coords) == 1:
+            coords = coords[0]
+        return constructor(coords)
+
+
 def _series_like_to_geoseries(geom, index):
     if index is None:
         index = geom.keys()
@@ -252,6 +322,8 @@ def _is_one_geometry(geom) -> bool:
 def _make_shapely_geoms(geom):
     if _is_one_geometry(geom):
         return _make_one_shapely_geom(geom)
+    if isinstance(geom, dict) and "coordinates" in geom:
+        return _from_json(geom)
     return (_make_one_shapely_geom(g) for g in geom)
 
 

{ssb_sgis-0.2.7 → ssb_sgis-0.2.8}/src/sgis/helpers.py

@@ -58,9 +58,9 @@ def get_name(var: object, n: int = 5) -> str | None:
     frame = inspect.currentframe().f_back.f_back
 
     for _ in range(n):
-        locals_ = frame.f_locals
-
-        names = [var_name for var_name, var_val in locals_.items() if var_val is var]
+        names = [
+            var_name for var_name, var_val in frame.f_locals.items() if var_val is var
+        ]
         if names and len(names) == 1:
             return names[0]
 

{ssb_sgis-0.2.7 → ssb_sgis-0.2.8}/src/sgis/maps/explore.py

@@ -4,6 +4,7 @@ This module holds the Explore class, which is the basis for the explore, samplem
 clipmap functions from the 'maps' module.
 """
 import warnings
+from numbers import Number
 from statistics import mean
 
 import branca as bc
@@ -15,10 +16,11 @@ from folium import plugins
 from geopandas import GeoDataFrame
 from IPython.display import display
 from jinja2 import Template
+from shapely import Geometry
 from shapely.geometry import LineString
 
 from ..geopandas_tools.general import clean_geoms, make_all_singlepart
-from ..geopandas_tools.geometry_types import get_geom_type
+from ..geopandas_tools.geometry_types import get_geom_type, to_single_geom_type
 from ..geopandas_tools.to_geodataframe import to_gdf
 from ..helpers import unit_is_degrees
 from .httpserver import run_html_server
@@ -110,6 +112,17 @@ class Explore(Map):
 
         super().__init__(*gdfs, column=column, **kwargs)
 
+        # remove columns not renerable by leaflet (list columns etc.)
+        new_gdfs = []
+        for gdf in self.gdfs:
+            cols_to_keep = [
+                col
+                for col in gdf.columns
+                if isinstance(gdf[col].iloc[0], (Number, str, Geometry))
+            ]
+            new_gdfs.append(gdf[cols_to_keep])
+        self._gdfs = new_gdfs
+
         self.popup = popup
         self.max_zoom = max_zoom
         self.smooth_factor = smooth_factor
@@ -266,6 +279,7 @@ class Explore(Map):
 
         new_gdfs = []
         for gdf in self._gdfs:
+            print(gdf)
             if get_geom_type(gdf) == "mixed" and not unit_is_degrees(gdf):
                 gdf[gdf._geometry_column_name] = gdf.buffer(0.01)
             gdf = make_all_singlepart(gdf)
@@ -615,7 +629,26 @@ class Explore(Map):
                     "supported as marker values"
                 )
 
-        gdf = clean_geoms(gdf)
+        gdf = clean_geoms(gdf).pipe(make_all_singlepart)
+        if get_geom_type(gdf) == "mixed":
+            if gdf.geom_type.str.lower().str.contains("polygon").any():
+                warnings.warn(
+                    "GeoJsonTooltip is not configured to render for GeoJson "
+                    "GeometryCollection geometries. Keeping only polygons."
+                )
+                gdf = to_single_geom_type(gdf, geom_type="polygon")
+            elif gdf.geom_type.str.lower().str.contains("line").any():
+                warnings.warn(
+                    "GeoJsonTooltip is not configured to render for GeoJson "
+                    "GeometryCollection geometries. Keeping only lines."
+                )
+                gdf = to_single_geom_type(gdf, geom_type="line")
+            else:
+                warnings.warn(
+                    "GeoJsonTooltip is not configured to render for GeoJson "
+                    "GeometryCollection geometries. Keeping only points."
+                )
+                gdf = to_single_geom_type(gdf, geom_type="point")
 
         # prepare tooltip and popup
         if isinstance(gdf, GeoDataFrame):