How to compute mean for non zero elements of a 3d volume along one of the dimensions

Question

I have a 3d volume HxWxD and I want to compute the mean only for values that are non zero along the D dimension to obtain a 2d HxW mea