title: ValueError: non-scalar numpy.ndarray cannot be used for fill - CuPy
body: <p>When I try to use the library <strong>CuPy</strong> with <strong>osgeo</strong>, I'm facing this error:</p>
  <pre><code>ValueError: non-scalar numpy.ndarray cannot be used for fill
  </code></pre>
  <p>I'm trying to fill this array:</p>
  <pre><code>im = cupy.zeros([ds.RasterYSize, ds.RasterXSize, ds.RasterCount], dtype=np.float32)
  
      for x in range(1, ds.RasterCount + 1):
          band = ds.GetRasterBand(x)
          im[:, :, x - 1] = band.ReadAsArray() #error here
  </code></pre>
  <p>Any fixes for this?</p>