如何通过GeoDjango读取PostGIS中MultiLineStringM的三维数据？

问题背景 #

我通过线性参考（Linear Referencing）构建铁路线路，以MultiLineStringM格式存储在FME创建的PostGIS数据库中，每条线路对应特定编码，沿线标注了资产的起止公里数（km_from、km_to）。我只有读取权限，需要根据给定起止值计算LineString，但用GeoDjango访问时，发现原生模型不支持M值，尝试设置dim=3的MultiLineStringField将M值映射为Z维度，读取到的仍是二维MultiLineString。

可行解决方案 #

方案1：原生SQL查询提取M值 #

直接通过RawQuerySet或GeoDjango的数据库函数获取带M值的几何表示：

from django.contrib.gis.db.models.functions import AsEWKT
from .models import RailwayLine

# 用AsEWKT获取包含M维度的文本格式
lines = RailwayLine.objects.annotate(geom_ewkt=AsEWKT('geom')).all()
for line in lines:
    # EWKT格式示例：SRID=3128;MULTILINESTRINGM((x1 y1 m1, x2 y2 m2), ...)
    print(line.geom_ewkt)

# 原生SQL提取点的M值
raw_query = """
    SELECT id, ST_M((ST_DumpPoints(geom)).geom) AS m_value
    FROM your_railway_table
    WHERE line_code = %s
"""
results = RailwayLine.objects.raw(raw_query, ['LINE_001'])
for res in results:
    print(res.m_value)

方案2：自定义支持M维度的字段 #

扩展GeoDjango的字段类型，让它识别并保留M值：

from django.contrib.gis.db.models.fields import MultiLineStringField
from django.contrib.gis.geos import MultiLineString

class MultiLineStringMField(MultiLineStringField):
    def db_type(self, connection):
        return 'MULTILINESTRINGM'

    def from_db_value(self, value, expression, connection):
        if value is None:
            return value
        # 初始化时指定包含M维度
        return MultiLineString(value, dim=3, has_m=True)

# 模型中使用自定义字段
class RailwayLine(models.Model):
    line_code = models.CharField(max_length=50)
    geom = MultiLineStringMField(srid=3128, blank=True, null=True, spatial_index=True)

注意：需确保GEOS库版本支持M维度解析，且Django数据库适配器能正确处理该类型。

方案3：显式转换M值为Z维度 #

如果要继续用Z维度映射M值，需在查询时通过PostGIS函数强制转换：

from django.contrib.gis.db.models.functions import ST_Force3DZ
from .models import RailwayLine

# 将M值转换为Z值后读取
lines = RailwayLine.objects.annotate(geom_3dz=ST_Force3DZ('geom')).all()
for line in lines:
    for ls in line.geom_3dz:
        for point in ls:
            print(f"X: {point.x}, Y: {point.y}, Z(原M值): {point.z}")

内容的提问来源于stack exchange，提问作者Richard