传入波段号1时ST_Contour报错“波段号需在1到0之间”的问题

按以下步骤排查解决：

• 第一步：确认栅格波段数
  先执行SQL验证你的栅格确实包含1个波段：

SELECT ST_NumBands(my_raster_column) FROM my_schema.my_table WHERE my_other_column IS NOT NULL;

如果返回值不是1，说明你的栅格数据本身存在问题，需先修复数据。

• 第二步：简化调用，使用默认波段号
  ST_Contour的bandnum参数默认值就是1，可省略该参数，只传入栅格和等高线间隔：

SELECT ST_Contour(my_raster_column, 5) FROM my_schema.my_table WHERE my_other_column IS NOT NULL;

如果此方式能正常运行，说明之前的参数传递存在位置混淆问题。

• 第三步：用关键字参数明确指定参数
  若必须指定波段号和偏移量，使用关键字参数避免位置错误（PostgreSQL支持通过参数名绑定值）：

SELECT ST_Contour(
    rast := my_raster_column,
    bandnum := 1,
    interval := 5.0,
    offset := 0.0
) FROM my_schema.my_table WHERE my_other_column IS NOT NULL;

• 第四步：调整Python代码的参数传递方式
  使用GeoAlchemy2时，同样用关键字参数调用func.ST_Contour，确保参数映射正确：

from sqlalchemy import select
from geoalchemy2 import func
# 假设你的模型类是MyTable，对应my_schema.my_table
from your_models import MyTable

async def fetch_contours():
    async with your_async_session() as session:
        stmt = select(
            func.ST_Contour(
                rast=MyTable.my_raster_column,
                bandnum=1,
                interval=5.0,
                offset=0.0
            )
        ).where(MyTable.my_other_column.isnot(None))
        result = await session.execute(stmt)
        # 后续处理结果逻辑

• 第五步：检查PostGIS版本
  若以上方法均无效，可能是你使用的PostGIS版本存在bug，建议升级到3.x系列的稳定版本（PostGIS 2.4之后才引入ST_Contour，旧版本可能存在参数逻辑问题）。

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