code: def compute_radius(wcs):
      
  
      ra, dec = wcs.wcs.crval
      img_center = SkyCoord(ra=ra * u.degree, dec=dec * u.degree)
      wcs_foot = wcs.calc_footprint()
      img_corners = SkyCoord(ra=wcs_foot[:, 0] * u.degree,
                             dec=wcs_foot[:, 1] * u.degree)
      radius = img_center.separation(img_corners).max().value
  
      return radius
docs: Compute the radius from the center to the furthest edge of the WCS.