gdal读取postgis数据库中的栅格数据（python）

---

前言

gdal读取tif格式的栅格影像非常方便，它还支持读取postgis的栅格影像。不过，一开始我以为读取方式与tif会有较大差异，所以就上网搜索gdal如何读取postgis栅格。很不幸，内网几乎没有，外网没查过，一是不会魔法上网，二是硬看英文太累。功夫不负有心人，最后还是找到了利用gdal读取postgis栅格的JAVA实现（https://blog.youkuaiyun.com/qq_24309981/article/details/84184772）。由于不同语言下，gdal的函数（方法）名称几乎是一致的，所以这也算解决了我的问题。

---

一、向postgis导入栅格

读取之前的让你的数据库有栅格影像。因为gdal暂不支持向postgis导入栅格，因此得用其他方法------最简单的就是用postgis自带的程序raster2pgsql。具体使用在这里不赘述，本人学习的是这篇文章：https://blog.youkuaiyun.com/theonegis/article/details/55050943
这篇文章在内网被转发了好多次，我也不知道谁是原创了。

二、gdal读取postgis栅格

其实很简单，就是文件路径跟读取tif时设置不一样，其他都一样。

 \\tif的情形
 path=r'D:abc\def\xxx.tif'
 \\