#include<bits/stdc++.h>
#define int long long
#define endl '\n'
using namespace std;
const int mod=1e9+7;
void solve()
{
    int n,m,ans=1;
    cin>>n>>m;
    string s[n+1];
    vector<vector<int>> mp(n+1,vector<int>(m+1));
    for(int i=1;i<=n;i++)
    {
        cin>>s[i];
        s[i]=' '+s[i];
        for(int j=1;j<=m;j++) mp[i][j]=s[i][j]-'0';
    }
    // for(int i=1;i<=n;i++)
    // {
    //     for(int j=1;j<=m;j++) cout<<mp[i][j]<<' ';
    //     cout<<endl;
    // }
    vector<set<int>> line(m+1);
    for(int i=1;i<=m;i++)
    {
        for(int j=1;j<=n;j++)
        {
            if(mp[j][i]) line[i].insert(j),line[m-i+1].insert(j);
        }
    }
    // for(int i=1;i<=m;i++)
    // {
    //     for(auto x:line[i]) cerr<<x<<' ';
    //     cerr<<endl;
    // }
    vector<int> num(m+1);
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            num[j]+=mp[i][j];
            num[m-j+1]+=mp[i][j];
        }
    }
    // for(int i=1;i<=m;i++) cout<<num[i]<<' ';
    // cout<<endl;
    vector<int> fa(n+1);
    for(int i=1;i<=n;i++) fa[i]=i;
    function<int(int)> getfa=[&](int x)
    {
        if(x==fa[x]) return x;
        return fa[x]=getfa(fa[x]);
    };
    for(int i=1;i<=m;i++)
    {
        if(line[i].size()>2) {cout<<0<<endl;return;}
        if(line[i].size()==2)
        {
            if(m==1) ans=0;
            int x=*line[i].begin();
            int y=*(++line[i].begin());
            int temp=num[i]-mp[x][m-i+1],temp2=num[i]-mp[y][m-i+1];
            if(temp>2 || temp2>2) {cout<<0<<endl;return;}
            int fx=getfa(x),fy=getfa(y);
            if(fx==fy) continue;
            fa[fy]=fx;
        }
    }
    for(int i=1;i<=n;i++)
        if(getfa(i)==i) ans=(ans*2)%mod;
    cout<<ans<<endl;
}
signed main()
{
    cin.tie(0)->sync_with_stdio(0);
    int t=1;
    cin>>t;
    while(t--) solve();
    return 0;
}