如何在numpy中对矩阵进行“短路”布尔或“ if”运算？

Question

我正在使用以下公式计算云覆盖率：

map = np.logical_or(
         np.logical_and(B02 > 0.175, (B03 - B04) / (B03 + B04) > 0),
            B02 > 0.39
     )

其中所有B02，B03和B04均为2000x2000双矩阵。有时我在B03 + B04中遇到零，并且numpy警告我有关错误的值。

如何以任何方式排除它们？