code: def test_expandim_in_groupby_aggregate_multiple_funcs():
      # https://github.com/geopandas/geopandas/pull/2296#issuecomment-1021966443
      # There are two calls to _constructor_expanddim here
      # SeriesGroupBy._aggregate_multiple_funcs() and
      # SeriesGroupBy._wrap_series_output() len(output) > 1
  
      s = GeoSeries.from_xy([0, 1, 2], [0, 1, 3])
  
      def union(s):
          return s.union_all()
  
      def total_area(s):
          return s.area.sum()
  
      grouped = s.groupby([0, 1, 0])
      agg = grouped.agg([total_area, union])
      assert_obj_no_active_geo_col(agg, GeoDataFrame, geo_colname=None)
      result = grouped.agg([union, total_area])
      assert_obj_no_active_geo_col(
          result, GeoDataFrame, geo_colname=None
      )
      assert_object(grouped.agg([total_area, total_area]), pd.DataFrame)
      assert_object(grouped.agg([total_area]), pd.DataFrame)
  
imports: import pandas
  import pyproj
  import pytest
  from shapely.geometry import Point
  import numpy
  from geopandas import GeoDataFrame, GeoSeries
  from geopandas.testing import assert_geodataframe_equal
label_desc: Use image node_id 17 for calling a global function with example usage: test_expandim_in_groupby_aggregate_multiple_funcs() and returns: s, s