Форумы GIS-Lab.info

Здравствуйте

Помогите решить задачу обрезки растра по векторному слою. Растр чернобелый(формат*.jpg, или *.bmp). векторный слой - таблица mapinfo.
нашел информацию про команды

gdal_translate с параметром [-projwin ulx uly lrx lry]. Я так понимаю, что экстент выходного растра прямоугольник? а как поступить, если нужно обрезать по произвольноиу контуру?

попробовал gdal_rasterize, выходной файл был записан без изменений.
кроме того windows ругнулась на gdal_rasterize.exe.
можно ли как то использовать в этой команде не shp-файлы, а tab?
не воспринимаются русские буквы в имени файлов? можно это как-нибудь обойти?

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

Да, охват - прямоугольник. Насчет кириллицы и tab - так переведите просто, чтобы без кириллицы и в shape. Вот еще два способа, попробуйте.

Способ 1 (возможно cutline поддерживает и tab, пробуйте)

gdal_translate -projwin ulx uly lrx lry rasterin.tif raster_clip.tif
gdalwarp -cutline shape_file raster_clip.tif raster_clip2.tif

Способ 2
http://trac.osgeo.org/gdal/wiki/FAQRast ... eGDAL1.8.0

Большое спасибо за ответ
попробовал
gdalwarp -cutline shape_file raster_clip.tif raster_clip2.tif
tab принимает в обработку
но выходной растр пустой. и возникает предупреждение: vector layer has a srs, source rastr dataset does not.
это изза несовпадения проекций? как узнать проекцию растра?

Посмотреть: Если автоматически не определилась проекция, то задайте вручную:

gdalinfo выдает Coordynate system is ''
это значит растр не имеет географической системы координат чтоли?

это значит что система координат не указана явным образом

пропишите систему координат, сделайте так, чтобы таб и растры были в одной системе координат (т.е. если они в разной, перепроецируйте либо растр либо вектор).

Вот тут http://asgs.narod.ru/programs.html лежит программа TiffMosaic. Вообще, как видно из названия, разрабатывалась для создания мозаики в MapInfo из ортофотопланов или космоснимков, но можно использовать и для обрезки растров, и как раз по вектору. Может и пригодится кому

Как определить охват векторного объекта/слоя, чтобы потом по нему обрезать?

Пример результата:

подскажите еще, как отсюда вытащить экстент? чтобы загнать их в gdal_translate -projwin ulx uly lrx lry

Добрый день! Объясните кто-нибудь популярно и доступно, пошагово на выполнить обрезку растра и какие в этом процессе "подводные камни".Help!

типо читайте выше

На самом деле подводные камни есть. при сложном векторе дает артефакты.
артефактов не дает, но никак не могу заставить заливать отрезаемое белым цветом (вернее заливать то заливает(-dstnodata "255 255 255" ), но вместе с полями прихватывает все пиксели, где хотя бы в 1 канале есть 255). По-моему это баг.

Denis Rykov писал(а):Не уверен, что правильно, но как вариант:
1. Открываем в QGIS вектор, по которому нужно произвести обрезку. В свойствах символики выставляем сплошную заливку. Для определенности будем считать, что цвет заливки 255 0 0 (красный).
2. File - Save As Image..., сохраняем в tiff.
3. Создаем новый проект.
4. Добавляем растр, который нужно обрезать и сверху накладываем растр, созданный по вектору.
5. Заходим в свойства растра, который был создан по вектору и задаем 100% прозрачность для цвета 255 0 0.
6. File - Save As Image.
Пример результата:

Все сделала, сработало! Интересно, что растр, созданный по вектору был в виде определенной фигуры, обрисованной по интересующей территории, а в итоге получился прямоугольный растр, в который попала вся обрисованная территория. В-принципе устраивает, но почему так получилось? В вашем примере обрезка была четко по заданному полю.

Озарян Юлия писал(а):

Denis Rykov писал(а):Не уверен, что правильно, но как вариант:
1. Открываем в QGIS вектор, по которому нужно произвести обрезку. В свойствах символики выставляем сплошную заливку. Для определенности будем считать, что цвет заливки 255 0 0 (красный).
2. File - Save As Image..., сохраняем в tiff.
3. Создаем новый проект.
4. Добавляем растр, который нужно обрезать и сверху накладываем растр, созданный по вектору.
5. Заходим в свойства растра, который был создан по вектору и задаем 100% прозрачность для цвета 255 0 0.
6. File - Save As Image.
Пример результата:

Все сделала, сработало! Интересно, что растр, созданный по вектору был в виде определенной фигуры, обрисованной по интересующей территории, а в итоге получился прямоугольный растр, в который попала вся обрисованная территория. В-принципе устраивает, но почему так получилось? В вашем примере обрезка была четко по заданному полю.

Еще, в качестве варианта, все элементарно делается в SAGA:
в инструментах вкладки: shapes\grid\spatial extent\ clip grid with rectangle