Can we optimize geotop calculations by using .data instead of .values?

[dbrakenhoff]

We call .values quite a bit in the code, which forces the dask array to read data and return a numpy array for a particular variable. If we use the .data attribute we can delay computation until the very end, and use dask support for chunking etc?

Something worth checking out for large datasets.