# test zonal stats on points

import os, sys
import numpy as np
sys.path.append('./')
import zonepy as zp

gdb = os.path.join('examples', 'data', 'zone10km.gpkg')
ras = os.path.join('examples', 'data', 'zone10kmCellNum.tif')
lyrName = 'modelgrid' 
inputfile = os.path.join('examples', 'data', 'outputfilepoly.csv')
outTiff = os.path.join('examples', 'data', 'statRaster.tif')

zs = zp.ZoneClass(gdb, ras, lyrName=lyrName, fldname=None , 
                 buffDist=0, fact=1, 
                 outND=np.nan, nd_thresh=100, filenm=None)

zs.writeRaster(outTiff=outTiff, inputfile=inputfile)
