Использование Python для преобразования шейп-файла в Geotiff - PullRequest
Новая
       126

Использование Python для преобразования шейп-файла в Geotiff

0 голосов
/ 05 марта 2019

У меня есть файл .shp, и я хотел бы преобразовать его в GEOTIFF.Мой файл формы состоит из большого многоугольника со множеством многоугольников внутри.Я использую следующий код, но выходной TIF состоит только из большого многоугольника.

from osgeo import ogr, gdal
import subprocess

InputVector = Shapefile
OutputImage = OutputfileName

gdalformat = 'GTiff'
datatype = gdal.GDT_Byte
burnVal = 1 
Shapefile = ogr.Open(InputVector)
Shapefile_layer = Shapefile.GetLayer()


Output = gdal.GetDriverByName(gdalformat).Create(OutputImage, RasterXSize, RasterYSize, 1, datatype, options=['COMPRESS=DEFLATE'])
Output.SetProjection(Projection)
Output.SetGeoTransform(GeoTransform) 

Band = Output.GetRasterBand(1)
Band.SetNoDataValue(0)
gdal.RasterizeLayer(Output, [1], Shapefile_layer, burn_values=[burnVal])

subprocess.call("gdaladdo --config COMPRESS_OVERVIEW DEFLATE "+OutputImage+" 2 4 8 16 32 64", shell=True)

Я не уверен, что я здесь не так делаю.

Спасибо

1 Ответ

0 голосов
/ 12 марта 2019

Моим первым предположением было бы убедиться, что у большого многоугольника есть «дыры» / NoData, где встречаются меньшие многоугольники.Я постараюсь быть более ясным с примером из инструмента стирания ArcMap.То, что вы хотите, чтобы ваш большой многоугольник выглядел так: с синим галочком (с наложением, очевидно, с стиранием), я думаю, что ваш большой полигон в настоящее время выглядит так же, как и с красным галочком, в результате чего ваш растр будет состоять только избольшой полигон.Мои извинения, если я полностью неправ, тогда мы найдем другое решение.

enter image description here

Если это так (ваш большой многоугольник похож насиняя галочка), сжигание всех полигонов до 1 также может создать путаницу.Вы можете создать новое поле в своем шейп-файле (например, UID) и присвоить ему уникальный числовой идентификатор (он должен быть числовым).Вы можете растеризовать на основе нового поля, например: 

gdal.RasterizeLayer(Output, [1], Shapefile_layer, options = ['ATTRIBUTE=UID'])
...