In 2019 I implemented a Hapke photometric model using Python and Numba (there are many flavors) from the 2012 textbook by Hapke for a Space Weathering class project. I used Numba to gain massive speed improvements without using complications like cython.
Example output plot of Henyey-Greenstein asymmetry function
Python Hapke Model
