Преобразование растерио и аффинность

Question

Я пытаюсь сделать базовую фильтрацию изображений. Ive включил дословный фрагмент из поваренной книги растерио (я удалил.astype() из медианного фильтра). Проблема в том, что мои входные и выходные растры должны иметь одинаковую степень, но не иметь. Преобразование и аффинное различны для ввода и вывода. Это ожидаемое поведение? Нужно ли что-то делать для аффинного преобразования и преобразования, чтобы вывод был таким же, как ввод?

Python 2.7.11 | Anaconda 4.0.0 (64-разрядная версия)| (по умолчанию, 16 февраля 2016 г., 09:58:36) [MSC v.1500 64 бит (AMD64)] на win32

rasterio == 0.36.0

import rasterio
from scipy.signal import medfilt

path = "map.tif"
output = "map2.tif"

with rasterio.open(path) as src:
    array = src.read()
    profile = src.profile

# apply a 5x5 median filter to each band
filtered = medfilt(array, (1, 5, 5))

# Write to tif, using the same profile as the source
with rasterio.open(output, 'w', **profile) as dst:
    dst.write(filtered)

    print profile
    print dst.profile

>>> {'count': 1, 'crs': CRS({'init': u'epsg:3857'}), 'interleave': 'band', 'dtype': 'float64', 'affine': Affine(100.0, 0.0, -13250000.0, 0.0, 100.0, 3980000.0), 'driver': u'GTiff', 'transform': (-13250000.0, 100.0, 0.0, 3980000.0, 0.0, 100.0), 'height': 1700, 'width': 1700, 'tiled': False, 'nodata': None}
>>> {'count': 1, 'crs': CRS({'init': u'epsg:3857'}), u'interleave': 'band', 'dtype': 'float64', 'affine': Affine(-13250000.0, 100.0, 0.0, 3980000.0, 0.0, 100.0), 'driver': u'GTiff', 'transform': (0.0, -13250000.0, 100.0, 100.0, 3980000.0, 0.0), 'height': 1700, 'width': 1700, u'tiled': False, 'nodata': None}

Answer 1

Растровые документы включают в себя историю аффинно-трансформационного использования, которое может оказаться полезным. Я имел обыкновение иметь пару строк, таких как следующие, чтобы справиться с этим:

out_profile = src.profile.copy()
out_affine = out_profile.pop("affine")
out_profile["transform"] = out_affine

# then, write the output raster

with rasterio.open(output, 'w', **out_profile) as dst:
    dst.write(filtered)

Я думаю, что это то, что здесь необходимо.