zField[Ly1[mask].astype("int"), Lx1[mask].astype("int")] = 1.