Question

我有一个名为geodf的geopandas数据框。该数据框中每行的几何图形是多面。例如，如果我们通过地理数据框考虑几何列中的第一行：

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)

然后我通过我们的Crate 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'])

我在上面的查询中通过box

传递多边形.format(box)

我想要做的是在上面的查询中传递bound而不是box（请注意，上面的查询适用于box）。但是，当我尝试传递bound时，我收到以下错误：

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

我无法诊断上面的错误。我想知道为什么bounds不起作用以及我怎么能使它工作？我们不想使用边界框，因为我们的多边形边界中没有包含大量多余的区域

Answer 1

您可以将multipolygon表示为类似于geoJSON格式的dict 然后使用参数替换。我对形状不熟悉，所以我无法评论你如何获得dict表示，但是使用手工制作的字典，它看起来如下：

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())

通过Crate DB表

1 个答案: