Django项目中如何从PostGIS获取栅格数据并计算指定位置的栅格值

现有代码错误点 #

• 数据查询逻辑错误：判断指定name的栅格存在后，调用了ndvimodel.objects.all()拉取全表数据，而非过滤对应name的数据集，会遍历到无关栅格数据。
• 坐标传入顺序错误：OGR/GEOS的Point构造遵循(经度, 纬度)的参数顺序，当前代码写反了经纬度位置，导致点位坐标完全偏移，无法匹配栅格位置。
• 缺少像素值查询逻辑：仅获取了rast字段对象，未调用栅格值查询方法，无法拿到点位对应的像素值。
• 缺少坐标系校验：如果传入经纬度的坐标系和栅格的坐标系不匹配，需要先做坐标转换，否则点位无法匹配到正确的栅格位置。

修正后代码 #

from django.contrib.gis.db.models.functions import ST_Value
from django.contrib.gis.geos import Point
import json

class dataviewView(views.View):
    def post(self, request):
        data = json.loads(request.body)
        name = data['name']
        lat = data['lat']
        lon = data['long']
        # srid根据你实际使用的坐标系填写，WGS84对应4326
        point = Point(lon, lat, srid=4326)
        
        # 若栅格坐标系和点位坐标系不一致，先做转换，示例为转成3857坐标系
        # point.transform(3857)
        
        # 直接过滤指定name的栅格，通过ST_Value注解拿到对应点位像素值
        target_raster = ndvimodel.objects.filter(name=name).annotate(
            pixel_value=ST_Value('rast', point)
        ).first()
        
        if target_raster:
            # 输出查询到的像素值
            print(f"对应点位像素值：{target_raster.pixel_value}")
            # 获取完整GDALRaster对象
            raster_data = target_raster.rast
            print(f"栅格对象类型：{type(raster_data)}")

额外校验项 #

• 确保ndvimodel模型中的rast字段是通过django.contrib.gis.db.models.RasterField定义的，且数据库已开启PostGIS栅格扩展。
• 确认数据库中存储的栅格数据坐标系和代码中配置的点位坐标系一致，或已做正确的坐标转换。
• 若栅格数据是分块存储的，需要确保查询逻辑覆盖所有可能的分块，避免点位落在未查询的分块中导致无返回值。

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