Установить CRS для чтения файла с помощью растерио

Question

Я читаю изображение jpg и связанный с ним файл мира в Python с помощью Rasterio следующим образом:

import rasterio
with rasterio.open('/path/to/file.jpg') as src:
    print(src.width, src.height)
    print(src.crs)
    print(src.indexes)

Файл изображения и связанный с ним файл мира читаются правильно, однако CRS не определен ( Я думаю, это потому, что мировой файл не содержит CRS ).Вот вывод:

5000 5000
None
(1, 2, 3)

Как вручную установить CRS в Rasterio после чтения файла?

Answer 1

Не видя файл мира, я не знаю наверняка, если это правильно в деталях, но я использовал следующее, чтобы добавить transform и CRS в растр после чтения в файл с файлом мира:

from affine import Affine
import rasterio.crs

a, d, b, e, c, f = np.loadtxt(world_filename)    # order depends on convention
transform = Affine(a, b, c, d, e, f)
crs = rasterio.crs.CRS({"init": "epsg:4326"})    # or whatever CRS you know the image is in    
with rasterio.open('/path/to/file.jpg') as src:
    src.transform = transform
    src.crs = crs