Create zonalstatsfun.py