postgis存储遥感影像

### 使用 PostGIS 存储和管理遥感影像数据

#### 安装必要的扩展和支持库
为了能够有效地存储和处理遥感影像，在 PostgreSQL 数据库中除了安装基本的 PostGIS 扩展外，还需要额外支持栅格数据类型的模块。这可以通过如下命令完成：

yum install -y postgis24_10  # 安装PostGIS及其依赖项[^1]

接着，在目标数据库内激活 `postgis` 及其栅格功能：

CREATE EXTENSION postgis;
CREATE EXTENSION postgis_raster;  -- 启用栅格数据支持[^3]

#### 导入遥感影像文件至数据库
对于导入操作，通常会利用 GDAL 工具集中的 `raster2pgsql` 命令行工具来实现。该工具允许用户将常见的地理空间图像格式转换成适合存入 PostGIS 的形式。

假设有一个 GeoTIFF 文件名为 `landsat.tif` 需要上传到名为 `public.rasters` 表里，则可按下面的方式执行：

raster2pgsql -s 4326 -I -C -M landsat.tif public.rasters | psql -d your_database_name

这里 `-s 4326` 参数指定了输入坐标系 EPSG 编码；而其他选项则用于优化性能以及设置索引等特性。

#### 查询与分析栅格数据
一旦完成了上述准备工作之后，就可以通过 SQL 语句来进行各种查询和计算工作了。例如获取特定区域内的平均反射率值（假定波段编号为1），可以编写如下脚本：

SELECT rid,
       AVG((ST_SummaryStats(rast, 1)).mean) OVER () as avg_reflectance
FROM public.rasters
WHERE ST_Intersects(ST_MakeEnvelope(minx, miny, maxx, maxy, 4326), rast);

此代码片段展示了如何基于给定边界框筛选符合条件的数据记录，并对其内部像素强度求均值得出统计结果。