zField[dy[mask].astype("int"), dx[mask].astype("int")] = 1.