code: def test_identity_gcps():
      """Define an identity transform using GCPs"""
      # Tile: [53, 96, 8]
      src_crs = dst_crs = "EPSG:3857"
      width = height = 1000
      left, bottom, right, top = (
          -11740727.544603072,
          4852834.0517692715,
          -11584184.510675032,
          5009377.085697309,
      )
      # For comparison only, these are not used to calculate the transform.
      transform = from_bounds(left, bottom, right, top, width, height)
  
      # Define 4 ground control points at the corners of the image.
      gcps = [
          GroundControlPoint(row=0, col=0, x=left, y=top, z=0.0),
          GroundControlPoint(row=0, col=1000, x=right, y=top, z=0.0),
          GroundControlPoint(
              row=1000, col=1000, x=right, y=bottom, z=0.0
          ),
          GroundControlPoint(row=1000, col=0, x=left, y=bottom, z=0.0),
      ]
  
      # Compute an output transform.
      (
          res_transform,
          res_width,
          res_height,
      ) = _calculate_default_transform(
          src_crs, dst_crs, height=height, width=width, gcps=gcps
      )
  
      assert res_width == width
      assert res_height == height
      for res, exp in zip(res_transform, transform):
          assert round(res, 3) == round(exp, 3)
  
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 6 for calling a global function with example usage: test_identity_gcps() without return types