code: def test_rasterize_existing_output(tmpdir, runner, basic_feature):
      """
      Create a rasterized output, then rasterize additional pixels into it.
      The final result should include rasterized pixels from both
      """
  
      truth = np.zeros(DEFAULT_SHAPE)
      truth[2:4, 2:4] = 1
      truth[4:6, 4:6] = 1
  
      output = str(tmpdir.join("test.tif"))
      result = runner.invoke(
          main_group,
          [
              "rasterize",
              output,
              "--dimensions",
              DEFAULT_SHAPE[0],
              DEFAULT_SHAPE[1],
              "--bounds",
              bbox(0, 10, 10, 0),
          ],
          input=json.dumps(basic_feature),
          catch_exceptions=False,
      )
  
      assert result.exit_code == 0
      assert os.path.exists(output)
  
      coords = np.array(basic_feature["geometry"]["coordinates"]) + 2
      basic_feature["geometry"]["coordinates"] = coords.tolist()
  
      result = runner.invoke(
          main_group,
          [
              "rasterize",
              "--overwrite",
              "-o",
              output,
              "--dimensions",
              DEFAULT_SHAPE[0],
              DEFAULT_SHAPE[1],
          ],
          input=json.dumps(basic_feature),
      )
  
      assert result.exit_code == 0
  
      with rasterio.open(output) as out:
          assert np.array_equal(truth, out.read(1, masked=False))
  
imports: import json
  import os
  import numpy
  import rasterio
  from rasterio.rio.main import main_group
label_desc: Use image node_id 10 for calling a global function with example usage: test_rasterize_existing_output(tmpdir, runner, basic_feature) without return types