Python 3, 135 127 127123 bytes
Answering my own question since there is no python solution yet. Expects a transposed matrix. Strongly inspired by Lynn's rotation method. I'm sure there is room for improvement!
def f(x):
x+=[[3 for y in x]]x+=[[3]*len(x)]
for _ in x*4*len(x):x=[map(lambda a,b:a|a*2&b,y,(0,*y))for y in zip(*x)][::-1]
return 3in x[0]