这基本上是一个关于在 Django 代码中运行自定义 PostGIS 函数的问题。
我必须完成的任务是LineStrings根据它们与给定几何图形的 Frechet 距离进行过滤。Frechet 距离应该使用ST_FrechetDistancePostGIS 提供的函数来计算。

在另一个基于 SQLAlchemy 的项目中，我使用以下功能完成完全相同的任务（它正在运行）：

from geoalchemy2 import Geography, Geometry
from sqlalchemy import func, cast

def get_matched_segments(wkt: str, freche_threshold: float = 0.002):
    matched_segments = db_session.query(RoadElement).filter(
        func.ST_Dwithin(
            RoadElement.geom,
            cast(wkt, Geography),
            10
        )
    ).filter(
        (func.ST_FrechetDistance(
            cast(RoadElement.geom, Geometry),
            cast(wkt, Geometry),
            0.1
        ) < freche_threshold) |
        # Frechet Distance is sensitive to geometry direction
        (func.ST_FrechetDistance(
            cast(RoadElement.geom, Geometry),
            func.ST_Reverse(cast(wkt, Geometry)),
            0.1
        ) < freche_threshold)
    )
    return matched_segments

正如我所说，上面的函数正在运行，我想在 Django 中重新实现它。我必须添加额外的几何 SRS 转换，因为在基于 SQLite 的项目中，LineStrings 最初位于 EPSG:4326 中，而在 Django 中它们位于 EPSG:3857 中。以下是我想出的办法：

from django.db.models import Func, Value, Q, QuerySet, F
from django.contrib.gis.geos import GEOSGeometry


class HighwayOnlyMotor(models.Model):
    geom = LineStringField(srid=3857)

def get_matched_segments(wkt: str, freche_threshold: float = 0.002) -> QuerySet:
    linestring = GEOSGeometry(wkt, srid=4326)
    transform_ls = linestring.transform(3857, clone=True)
    linestring.reverse()
    frechet_annotation = HighwayOnlyMotor.objects.filter(
        geom__dwithin=(transform_ls, D(m=20))  
    ).annotate(
        fre_forward=Func(
            Func(F('geom'), Value(4326), function='ST_Transform'),
            Value(wkt),
            Value(0.1),
            function='ST_FrechetDistance'
        ),
        fre_backward=Func(
            Func(F('geom'), Value(4326), function='ST_Transform'),
            Value(linestring.wkt),
            Value(0.1),
            function='ST_FrechetDistance'
        )
    )
    matched_segments = frechet_annotation.filter(
        Q(fre_forward__lte=freche_threshold) |
        Q(fre_backward__lte=freche_threshold)
    )
    return matched_segments

它不起作用，因为frechet_annotationQuerySet 引发了异常：

django.db.utils.ProgrammingError: cannot cast type double precision to bytea
LINE 1: ...548 55.717805109,36.825235998 55.717761246)', 0.1)::bytea AS...
                                                             ^

似乎我错误地定义了“ST_FrechetDistance”计算。我该如何修复它？

更新

检查了 Django 编写的 SQL。它总体上是正确的，但尝试将结果转换为 会FrecheDistance破坏bytea它ST_FrechetDistance(...)::bytea。当我手动运行未进行转换的查询时bytea，SQL 可以正常工作。那么问题是如何避免将结果转换为bytea？