From f035e236a97bc739ab8b162d33e71fac4766586b Mon Sep 17 00:00:00 2001
From: Martin Fleischmann <martin@martinfleischmann.net>
Date: Thu, 6 Feb 2025 18:43:04 +0100
Subject: [PATCH] use min_scalar_type in rasterize (#97)

---
 xvec/zonal.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/xvec/zonal.py b/xvec/zonal.py
index 715ded1..15e1dc6 100644
--- a/xvec/zonal.py
+++ b/xvec/zonal.py
@@ -52,7 +52,7 @@ def _zonal_stats_rasterize(
 
     transform = acc._obj.rio.transform()
     length = len(geometry)
-    dtype = np.int16 if length < np.iinfo(np.int16).max else np.int32
+    dtype = np.min_scalar_type(length + 1)
 
     labels = features.rasterize(
         zip(geometry, range(length), strict=False),