code: def multi_point(point_gdf):
      """Create a multi-point GeoDataFrame."""
      multi_point = point_gdf.unary_union
      out_df = gpd.GeoDataFrame(
          geometry=gpd.GeoSeries(
              [
                  multi_point,
                  Point(2, 5),
                  Point(-11, -14),
                  Point(-10, -12),
              ]
          ),
          crs="epsg:4326",
      )
      out_df = out_df.rename(columns={0: "geometry"}).set_geometry(
          "geometry"
      )
      out_df["attr"] = ["tree", "another tree", "shrub", "berries"]
      return out_df
  
file_path: earthpy/earthpy/tests/test_clip.py
imports: import pytest
  import numpy
  from shapely.geometry import Polygon, Point, LineString
  import geopandas
  import earthpy.clip
label_desc: Use image node_id 5 for calling a global function with example usage: multi_point(point_gdf) and returns: out_df