ism3d.utils.misc.makepb

ism3d.utils.misc.makepb(header, phasecenter=None, antsize=<Quantity 12. m>)

make a 2D Gaussian image approximated to the ALMA (or VLA?) primary beam

https://help.almascience.org/index.php?/Knowledgebase/Article/View/90/0/90

note: this is just a apprximate solution, assuming the the pointing is towards the reference pixel in the header

and use the first channel as the reference frequency.