在Python中，如何在相同的区域切片FITS图像？

我有一个FITS图像，我想在Python中将其切成相等的垂直区域，以便分别获得每个区域中的计数数，并最终从中获得1D径向轮廓。我该怎么做？这是fits文件的信息：

from astropy.io import fits

hdu_list = fits.open('file_name')
hdu_list.info()


No.    Name      Ver    Type      Cards   Dimensions   Format
  0  PRIMARY       1 PrimaryHDU      27   (922, 731)   float32  

到目前为止，我获得的是将图像绘制如下：

image_data = fits.getdata('file_name', ext=0)

from matplotlib.colors import LogNorm

plt.imshow(image_data, norm=LogNorm())

plt.hist2d(x, y, bins=1, norm=LogNorm())

plt.colorbar()
plt.grid()

plt.xlim(197,682)
plt.ylim(115,550)

plt.xlabel('x')
plt.ylabel('y')