I am performing an operation in which I have to calculate sum of all pixel intensity, if the pixel value is >0. Currently, I am iterating through every pixel, which is slower. Can anyone please suggest a faster way to improve my code?

```
h, w = image.shape[:]
total_intensity = 0
for x in range(h):
    for y in range(w):
        if (image[x][y] > 0):
            total_intensity += image[x][y]

```