通过基于统计数据汇总5x5像素来缩小TIFF图像

Question

如何降低10m分辨率的TIFF映像，并创建一个50m的新图像，每个像素是第一个图像的统计信息？

初始TIFF映像是二进制分类图 - 意味着每个像素（10m）属于“水”（值= 0）或类“ ICE”（值= 1）的类。 我想创建一个新图像，其中每个像素是初始地图的5 x 5块中水百分比，这意味着新图像的每个像素都将具有50 m的分辨率，并且表示“水的比率或百分比” “以前地图的每5x5像素上的像素。您可以在此处看到示例：示例

这是一个映像示例（可以从Google驱动器下载）：

Answer 1

您的图像以一种相当奇怪的格式保存，使用32位浮子仅表示可以用一点点表示的两类数据，因此我用 imagemagick 使用：

magick YOURIMAGE.TIF -alpha off image.png

许多Python库会在您的实际TIFF上结结巴巴，因此也许要考虑使用另一种编写方式。

完成此操作后，代码可能看起来像这样：

#!/usr/bin/env python3

from PIL import Image
import numpy as np

# Set size of tiles for subsampling
tileX = tileY = 5

# Open image and convert to greyscale and thence to Numpy array
im = Image.open('image.png').convert('L')
na = np.array(im)

# Round height and width down to next lower multiple of tile sizes
h = (na.shape[0] // tileY) * tileY
w = (na.shape[1] // tileX) * tileX

# Create empty output array to fill
res = np.empty((h//tileY,w//tileX), np.uint8)

pxPerTile = tileX * tileY
for yoffset in range(0,h,tileY):
    for xoffset in range(0,w,tileX):
        # Count ice pixels in this 5x5 tile
        nonZero = np.count_nonzero(na[yoffset:yoffset+tileY, xoffset:xoffset+tileX])
        percent = int((100.0 * (pxPerTile - nonZero))/pxPerTile)
        res[yoffset//tileY, xoffset//tileX] = percent

# Make Numpy array back into PIL Image and save
Image.fromarray(res.astype(np.uint8)).save('result.png')

---

在反射上，您可以更快地进行操作，而更简单地使用

Answer 2

我在

#!/usr/bin/python3

import sys
import pyvips

image = pyvips.Image.new_from_file(sys.argv[1])

# label (0 == water, 1 == ice) is in the first band
label = image[0]

# average 5x5 areas
label = label.shrink(5, 5)

# turn into a percentage of water
water_percent = 100 * (1 - label)

# ... and save
water_percent.write_to_file(sys.argv[2])

版本

$ ./average.py ~/pics/meltPondClassifiedAndS12.tif x.png

一个 （相反）输出：