EPSG WKT в растерио - PullRequest
Новая
       56

EPSG WKT в растерио

0 голосов
/ 05 мая 2018

Добрый вечер всем. У меня вопрос по поводу вариантов деформации в растерио. У меня есть два набора растровых данных из Албании. Первый растр - продукт отсечения, а второй - границы страны. Во-первых, у меня есть следующая проекция. введите описание изображения здесь

Для второго растра у меня есть следующая информация: введите описание изображения здесь

Моя цель - перепроектировать растр # 1 с пространственной информацией растра 2. Я попытался использовать растерио, и результат показывает следующую пространственную привязку: Krassovsky_1942_Transverse_Mercator, а датум неизвестен. Этот результат был получен после использования следующего кода: 

import numpy as np
import rasterio
from rasterio.warp import calculate_default_transform, reproject, Resampling


#-----------------------------------------------------------------------------  #
#                   OPENING IMAGE WITH THE CORRECT SPATIAL REFERENCE
#-----------------------------------------------------------------------------#
with rasterio.open(r"C:\Users\Roger\Documents\git\ArcpyRecipes\Test_rasterio\test_albania\raster\AL_020m_nat_buffer100m.tif") as proj:
dst_crs = proj.crs
print dst_crs
transform_proj, width_proj, height_proj = calculate

#-----------------------------------------------------------------------------#
#                   OPENING IMAGE WITH TO BE REPROJECTED
#-----------------------------------------------------------------------------#

with rasterio.open(r"C:\Users\Roger\Documents\git\ArcpyRecipes\Test_rasterio\test_albania\results\umprojected\ALB_adm0.tif") as src:
transform, width, height = calculate_default_transform(
    src.crs, dst_crs, src.width, src.height, *src.bounds)
kwargs = src.meta.copy()
kwargs.update({
    "crs": dst_crs,
    "transform": transform,
    "width": width,
    "height": height
})

#---------------------------------------------------------------------------
#                   CREATING THE OUTPUT IMAGE   
#---------------------------------------------------------------------------

with rasterio.open(r"C:\Users\Roger\Documents\git\ArcpyRecipes\Test_rasterio\test_albania\results\projected/Albita_projected_2.tif", "w", **kwargs) as dst:
    for i in range(1, src.count + 1):
        reproject(
            source = rasterio.band(src, i),
            destination=rasterio.band(dst, i),
            src_transform = src.transform,
            src_crs=src.crs,
            dst_transform = transform,
            dst_crs = dst_crs,
            resampling = Resampling.nearest)

Я прочитал, что код EPGS можно получить в ArcMap с помощью arcpy.Describe, но я хочу использовать модуль rasterio для выполнения этой задачи. Я думал, что опция (dataset.crs) после открытия изображения с растерио даст мне CRS для проецирования, но это не так. Может ли кто-нибудь помочь мне с этим?

Большое спасибо

...