ArcGIS光栅金字塔的Rasterio版本

Question

在ArcGIS中，当您加载一个栅格时，您可以选择构建金字塔。就像这样：

​

​

我发现可以使用以下Python代码和Rasterio库复制这些文件：

import rasterio
from rasterio.warp import Resampling

with rasterio.Env(TIFF_USE_OVR=True):

    with rasterio.open('new_new.tif', "w", **meta_data) as dest:
        dest.write(raster_data)
        overviews = [2,4,8,10,14,16]
        dest.build_overviews(overviews, Resampling.nearest)

但是，如果可能的话，我想要指导的是，我上面的代码中有多少“概述”将复制ArcGIS运行的金字塔工具？如果可能的话，我想尝试创建一个相同的输出。

Answer 1

您可以使用文档中的代码检查原始金字塔级别。

src = rasterio.open(path, 'r')
[src.overviews(i) for i in src.indexes]
src.tags(ns='rio_overview').get('resampling')

Arcmap使用视图级别(2,4,8,16)。据我所知，不可能使用rasterio编写带有压缩的概述文件，这样您就无法模仿ESRI通过压缩创建的ovr。