模拟星系

from astropy.modeling.models import Gaussian2D
import numpy as np
import matplotlib.pyplot as plt
x,y = np.mgrid[0:200,0:200]
model=Gaussian2D(1,100,100,20,50,np.pi/4)
galaxy_model=model(x,y)
plt.imshow(galaxy_model,cmap='jet')

　　

from：https://www.youtube.com/watch?v=djkx4nT59EM