#include <bits/stdc++.h>
using namespace std;

#define F(a, b) fixed << setprecision(b) << a
typedef long long LL;
const LL MOD = 1e9 + 7;

struct DSU{

    vector<int> f, siz;

    DSU(){}
    DSU(int n){
        init(n);
    }

    void init(int n){
        f.resize(n + 1);
        iota(f.begin() + 1, f.end(), 1);
        siz.assign(n + 1, 1);
        return;
    }

    int find(int x){
        return x == f[x] ? f[x] : f[x] = find(f[x]);
    }

    bool same(int x, int y){
        return find(x) == find(y);
    }

    bool merge(int x, int y){
        x = find(x); y = find(y);
        if(x == y) return false;
        siz[x] += siz[y];
        f[y] = x;
        return true;
    }

    int size(int x){
        return siz[find(x)];
    }

};

LL qpow(LL base, LL b, LL p){
    LL ans = 1;
    while(b){
        if(b & 1) ans = ans * base % p;
        base = base * base % p;
        b >>= 1;
    }
    return ans;
}

void solve(){
    int r, c;
    cin >> r >> c;
    vector<vector<char>> a(r + 1, vector<char>(c + 1));
    for(int i = 1; i <= r; i++){
        for(int j = 1; j <= c; j++){
            cin >> a[i][j];
        }
    }

    DSU dsu(r + 1);
    vector<int> ag(r + 1);

    auto Ag = [&](int x){
        return ag[dsu.find(x)] = dsu.find(ag[dsu.find(x)]);
    };

    for(int j = 1; j <= c / 2; j++){
        int num = 0;
        int xr = 0, xc = 0, yr = 0, yc = 0;
        for(int i = 1; i <= r; i++){
            if(a[i][j] == '1'){
                num++;
                if(xr) yr = i;
                else xr = i;
                if(xc) yc = j;
                else xc = j;
            }
            if(a[i][c + 1 - j] == '1'){
                num++;
                if(yr) xr = i;
                else yr = i;
                if(yc) xc = c + 1 - j;
                else yc = c + 1 - j;
            }
        }
        if(num >= 3){
            cout << "0\n";
            return;
        }
        // cout << xr << " " << xc << endl;
        // cout << yr << " " << yc << endl;
        if(num == 2){
            if(xr != yr){
                if(xc != yc){
                    if(Ag(xr) && Ag(yr)){
                        if(Ag(xr) != dsu.find(yr) && Ag(yr) != dsu.find(xr)){
                            dsu.merge(xr, yr);
                            dsu.merge(Ag(xr), Ag(yr));
                        }else{
                            cout << "0\n";
                            return;
                        }
                    }else dsu.merge(xr, yr);
                }else{
                    if(dsu.same(xr, yr)){
                        cout << "0\n";
                        return;
                    }else{
                        if(Ag(yr)) dsu.merge(Ag(yr), xr);
                        else ag[dsu.find(yr)] = dsu.find(xr);
                        if(Ag(xr)) dsu.merge(Ag(xr), yr);
                        else ag[dsu.find(xr)] = dsu.find(yr);
                    }
                }
            }
        }
    }
    if(c & 1){
        int num = 0;
        for(int i = 1; i <= r; i++)
            if(a[i][c / 2 + 1] == '1') num++;
        if(num >= 2){
            cout << "0\n";
            return;
        }
    }

    vector<bool> vis1(r + 1), vis2(r + 1);
    int fnum = 0;
    for(int i = 1; i <= r; i++){
        if(!vis1[dsu.find(i)]){
            fnum++;
            vis1[dsu.find(i)] = true;
        }
        if(!vis2[dsu.find(i)]){
            if(Ag(i)){
                fnum--;
                vis2[Ag(i)] = true;
            }
            vis2[dsu.find(i)] = true;
        }
    }

    cout << qpow(2LL, fnum, MOD) << "\n";
    return;
}

int main(){
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    int T;
    cin >> T;
    while(T--){
        solve();
    }
    return 0;
}