code: def test_write_crs_transform(tmpdir):
      name = str(tmpdir.join("test_write_crs_transform.tif"))
      a = np.ones((100, 100), dtype=rasterio.ubyte) * 127
      transform = affine.Affine(
          300.0379266750948,
          0.0,
          101985.0,
          0.0,
          -300.041782729805,
          2826915.0,
      )
  
      with rasterio.open(
          name,
          "w",
          driver="GTiff",
          width=100,
          height=100,
          count=1,
          crs={
              "units": "m",
              "no_defs": True,
              "datum": "WGS84",
              "proj": "utm",
              "zone": 18,
          },
          transform=transform,
          dtype=rasterio.ubyte,
      ) as s:
          s.write(a, indexes=1)
      assert s.crs.to_epsg() == 32618
      info = subprocess.check_output(["gdalinfo", name]).decode("utf-8")
      # make sure that pixel size is nearly the same as transform
      # (precision varies slightly by platform)
      assert re.search(
          r"Pixel Size = \(300.03792\d+,-300.04178\d+\)", info
      )
  
imports: import logging
  import re
  from pathlib import Path
  import subprocess
  import affine
  import numpy
  import pytest
  from .conftest import requires_gdal35, gdal_version
  import rasterio
  from rasterio.drivers import blacklist
  from rasterio.enums import MaskFlags, Resampling
  from rasterio.env import Env, GDALVersion
  from rasterio.errors import RasterioIOError
label_desc: Use image node_id 13 for calling a global function with example usage: test_write_crs_transform(tmpdir) without return types