使用rioxarray的rio.clip工具结合Shapefile裁剪NetCDF栅格并查找重叠区域的技术问询

我最近在处理栅格数据时遇到了一个需求：想用Shapefile配合rioxarray里的rio.clip工具裁剪NetCDF栅格，同时找出两者的重叠区域。

目前我先用一个简化的测试多边形做验证（实际场景里我会用州或者县一级的Shapefile），这个测试多边形的GeoJSON如下：

{"coordinates": [[[-103.623047,39.808536],[-95.932617,39.808536],[-95.932617,42.90816],[-103.623047,42.90816],[-103.623047,39.808536]]], "type": "Polygon"}

关于要处理的栅格，相关描述如下（原文内容截断）：它是通过读取...

备注：内容来源于stack exchange，提问作者deadpickle