dtype problem with polygonize

[mdbartos]

dtype problem with polygonize

[mdbartos]

Hi David, thanks for letting me know. A check should probably be added for the datatype.

grid.catchment should probably also return the same datatype as the source flow direction file (unless the nodata and pour_value values don't fit into that datatype).

Note that if you call polygonize with data=None (the default), the input will default to grid.mask. If you previously used grid.clip_to('catch') to clip the mask to the catchment, this should give you an easy way of generating the outline of the catchment:
https://github.com/mdbartos/pysheds/blob/master/examples/shapefiles.ipynb

Looking at the polygonize function, I wonder if the input behavior should be changed to be more consistent with the other functions (e.g. accept a named parameter like 'catch' or 'dem').