code: def read_region_as_imagecoord(s, header):
      
      rp = RegionParser()
      ss = rp.parse(s)
      sss1 = rp.convert_attr(ss)
      sss2 = _check_wcs(sss1)
      sss3 = rp.sky_to_image(sss2, header)
  
      shape_list = rp.filter_shape(sss3)
      return ShapeList(shape_list)
docs: Read region as image coordinates.
  
      Parameters
      ----------
      s : str
          Region string
      header : `~astropy.io.fits.Header`
          FITS header
  
      Returns
      -------
      shapes : `~pyregion.ShapeList`
          List of `~pyregion.Shape`