通过板条箱DB表通过多元而不是边界框

Question

我有一个称为 geodf的地理数据框。该数据框架中每一行的几何形状是多角形。例如，如果我们通过GeodataFrame考虑几何列中的第一行：

bound = geodf['geometry'].iloc[0]
print(bound)

看起来像这样：

MULTIPOLYGON (((-86.46228799973933 34.31021100007911, -86.46447100007661 34.31018399970606, -86.46447800010341 34.31197299998977, -86.4623920000716 34.31198799958079

多重子本身很大，所以我尝试为其创建一个盒子

box = shapely.geometry.box(*geodf.bounds)

然后我通过板条箱DB服务器通过以下查询传递了box：

query = """
                 SELECT geohash, grid,shape 
                 FROM geo 
                 WHERE layer = 'pop' 
                AND MATCH(shape, '{}') USING within;
          """.format(box)
geoCursor.execute(query)
results = geoCursor.fetchall()
District = gpd.GeoDataFrame(results, columns = ['geohash', 'grid','shape'])

我在上面查询中通过 .format(box)

通过多边形box

我想做的是通过上述查询中的bound而不是box(请注意，上面的查询适用于box(。但是，当我尝试通过bound时，我会收到以下错误：

ProgrammingError: SQLActionException[UnhandledServerException: java.lang.IllegalArgumentException: Cannot convert Map "{type=MultiPolygon, coordinates=[[[[D@2b59d486}" to shape]

我无法诊断上述错误。我想知道为什么bounds不起作用，以及如何使它起作用？我们不想使用边界框，因为我们的多角边界中没有很多多余的区域

Answer 1

您可以代表多重子，因为dict类似于geojson格式然后使用参数替代。我不熟悉Shapely，所以我无法评论您如何获得命令表示，但是有了手动制作的词典，它看起来如下：

multipolygon = {
    "type": "MultiPolygon",
    "coordinates": [
        [
            [
                [ 40.0, 40.0 ],
                [ 20.0, 45.0 ],
                [ 45.0, 30.0 ],
                [ 40.0, 40.0 ]
            ]
        ],
        [
            [
                [ 20.0, 35.0 ],
                [ 10.0, 30.0 ],
                [ 10.0, 10.0 ],
                [ 30.0, 5.0 ],
                [ 45.0, 20.0 ],
                [ 20.0, 35.0 ]
            ],
            [
                [ 30.0, 20.0 ],
                [ 20.0, 15.0 ],
                [ 20.0, 25.0 ],
                [ 30.0, 20.0 ]
            ]
        ]
    ]
}
c.execute('select * from t where match(x, ?) using within', (multipolygon, ))
print(c.fetchall())