Изменить интервал сетки по изогнутой оси атропии ⇐ Python

[Anonymous]

Код: Выделить всё

import astropy.units as u
import numpy as np
from astropy.io import fits
from astropy.wcs import WCS
import matplotlib.pyplot as plt
from astropy.visualization import simple_norm
from astropy.coordinates import SkyCoord, Galactic

file = "hi4pi-hvc-vlsr-gal-car.fits"
hdul = fits.open(file)
hdr = hdul[0].header
data = hdul[0].data
wcs = WCS(hdr)
norm = simple_norm(data, "linear")

# Plot image on axes ax
fig, ax = plt.subplots(figsize=(10, 6), subplot_kw={"projection": wcs})
cmap = plt.get_cmap("RdBu", 20)
im = ax.imshow(data, cmap=cmap)
ax.set_xlabel("Galactic Longitude")
ax.set_ylabel("Galactic Latitude")

overlay = ax.get_coords_overlay("icrs")
major_ticks = np.arange(0, 360, 30)
#ax.set_xticks(major_ticks)
overlay.grid(color="grey", ls="dotted")
plt.colorbar(im, orientation="horizontal")

plt.savefig("area.png", dpi=300, bbox_inches="tight")
plt.show()

Я хочу установить сетку RA на np.arange(0,360,30), интервал в 30 градусов вместо интервала в 90 градусов.
и декабрь меткой справа лучше быть Dec, но я понятия не имею, как ее изменить.

файл соответствия находится здесь: Westmeier 2018

Код: Выделить всё

WCS Keywords

Number of WCS axes: 2
CTYPE : 'GLON' 'GLAT'
CRVAL : np.float64(180.000000024) np.float64(0.0)
CRPIX : np.float64(2162.0) np.float64(1080.0)
PC1_1 PC1_2  : np.float64(1.0) np.float64(0.0)
PC2_1 PC2_2  : np.float64(0.0) np.float64(1.0)
CDELT : np.float64(-0.0833333333) np.float64(0.0833333333)
NAXIS : 4323  2144

Я не знаю, как установить интервалы сетки в overlay.grid. Я проверил документацию, но там, похоже, нет никаких функций.

Подробнее здесь: https://stackoverflow.com/questions/790 ... urved-axis