Fastest way to randomly permute only a subset of indices (akin to non-black pixels) in a numpy array

Question

I have a 3D numpy array of shape, e.g.: (500,500,3), with only some [i,j,:] cells that are zero across all 3 channels (imagine an RGB image with all channels being zero = bl