невозможно обрезать растр (ETrF.etrf.tif) в python? - PullRequest
Новая
       28

невозможно обрезать растр (ETrF.etrf.tif) в python?

0 голосов
/ 16 июня 2019

У меня есть 50 растров, сохраненных как LC81520412017282LGN00_ETrF.etrf.tif, и шейп-файл многоугольника, сохраненный как .shp.

, как при попытке обрезать это изображение с помощью Clip (Управление данными) в arcgis, оно прекрасно работает.но для полной обработки всего изображения я написал скрипт, как показано ниже 

import os
import arcpy
arcpy.env.workspace = "E:/Project"
for root, dirs, files in os.walk('E:/Project/a/r'):
    for file in files:
        if os.path.splitext(file)[1] == '.tif':
           print (file)
           arcpy.Clip_management(file,"#",
           "E:/Project/be/test.gdb/%file%"
           ,"be/canal.shp","#","Clipping 
           Geometry", "NO_MAINTAIN_EXTENT")

. Приведенный выше скрипт выдает ошибку «неподдерживаемый файл или файл не существует».

1 Ответ

0 голосов
/ 22 июня 2019

управление клипами только для того, чтобы обрезать шейп-файл, чтобы "вырезать" растровый файл, используя шейп-файл, который вы должны использовать 

ExtractByMask (in_raster, in_mask_data)

вот пример 

import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
outExtractByMask = ExtractByMask("elevation", "mask.shp")
outExtractByMask.save("C:/sapyexamples/output/maskextract")
...