code: def test_gcps_calculate_transform():
      src_gcps = [
          GroundControlPoint(row=0, col=0, x=156113, y=2818720, z=0),
          GroundControlPoint(row=0, col=800, x=338353, y=2785790, z=0),
          GroundControlPoint(
              row=800, col=800, x=297939, y=2618518, z=0
          ),
          GroundControlPoint(row=800, col=0, x=115698, y=2651448, z=0),
      ]
      _, width, height = calculate_default_transform(
          "epsg:3857", "epsg:4326", width=800, height=800, gcps=src_gcps
      )
      assert width == 1087
      assert height == 895
  
imports: import os
  import logging
  import numpy
  import pytest
  from numpy.testing import assert_almost_equal
  import rasterio
  from rasterio._err import CPLE_BaseError, CPLE_AppDefinedError
  from rasterio._warp import _calculate_default_transform
  from rasterio.control import GroundControlPoint
  from rasterio.crs import CRS
  from rasterio.errors import CRSError
  from rasterio.transform import from_bounds
  from rasterio.warp import calculate_default_transform, transform_bounds
  from tests.conftest import gdal_version
label_desc: Use image node_id 12 for calling a global function with example usage: test_gcps_calculate_transform() without return types