Astropy:确定一个点(RA，DEC)是否在给定正方形区域角点坐标的正方形区域内

from astropy.coordinates import SkyCoord
import astropy.units as u

# 定义正方形区域的角点坐标
corner1 = SkyCoord(ra=ra1*u.deg, dec=dec1*u.deg)
corner2 = SkyCoord(ra=ra2*u.deg, dec=dec2*u.deg)
corner3 = SkyCoord(ra=ra3*u.deg, dec=dec3*u.deg)
corner4 = SkyCoord(ra=ra4*u.deg, dec=dec4*u.deg)

# 定义待检查的点的坐标
point = SkyCoord(ra=ra*u.deg, dec=dec*u.deg)

# 检查点是否在正方形区域内
is_inside = corner1.contains(point) and corner2.contains(point) and corner3.contains(point) and corner4.contains(point)

if is_inside:
    print("点在正方形区域内")
else:
    print("点不在正方形区域内")