Step 1 安装 Python3

# 下载 Python3 安装包
# 解压缩 .tgz 文件
tar -xvf Python-3.9.6.tgz
# 编译安装 Python3
./configure
make
make install

Step 2 安装 conda

# 下载安装包 https://docs.conda.io/en/latest/miniconda.html
# 赋予执行权限
./iniconda3-py39_4.9.2-Linux-x86_64.sh

Step 3 安装 gdal

开源的空间数据处理程序

conda install -c conda-forge gdal

Step 4 安装 rasterio

MapBox 在 gdal 基础上开发的栅格工具

sudo add-apt-repository ppa:ubuntugis/ppa
sudo apt-get update
sudo apt-get install python-numpy gdal-bin libgdal-dev
pip install rasterio

Step 5 安装 rio-rgbify

MapBox 发布的将 dem 栅格编码为 rgb 栅格的 rasterio 插件

pip install rio-rgbify

Step6 数据预处理

  1. 确定坐标系
    GeoTiff 的坐标系必须是 WGS84 Web 墨卡托 (EPSG:3857)。 通常 Geo Tiff 格式的灰度栅格图片看起来是黑白的。

如果不清楚坐标系的话，可以使用 rasterio 提供的命令行来获取坐标系：

rio info --indent 2 xxxx.tif
  1. 转换坐标系，清理负数值
  • 如果坐标系不是 EPSG:3857，就需要进行坐标转换；
  • 用 负数值表示无数据 ，而 Terrain-RGB 无法表示负值，所以需要进一步处理；
# 同时进行坐标系转换和清楚表示无数据的负值
gdalwarp -t_srs EPSG:3857 -dstnodata None -co TILED=YES -co COMPRESS=DEFLATE -co BIGTIFF=IF_NEEDED xxxx.tif xxxx_n.tif
  1. 将转换后的文件转换为 Terrain-RGB 格式
    Terrain-RGB 用 3 个 byte 通过 rgb 三通道来表示高程， 比原来的灰度 tiff 要小很多
  • 将 灰度数据 转换为 RGB 数据
    高度计算公式：
height = -10000 + ((R * 256 * 256 + G * 256 + B) * 0.1)

因此设置 ribify 的参数 base value 的参数为 -10000 ， interval 为 0.1 ，继续输入以下命令：

rio rgbify -b -10000 -i 0.1 xxxx_n.tif xxxx_n_rgb.tif

Step7 通过 GeoServer 发布切片服务

  1. 创建工作区
  2. 创建数据存储
    数据源选择 GeoTIFF
  3. 发布图层
  4. 创建图层组

Step 8 查看结果