GDAL+MapInfo: разрезаем изображение прямоугольной сеткой - История изменений

Максим Дубинин: /* Создание контура изображения */

2012-10-06T07:36:14Z

Создание контура изображения

← Предыдущая версия		Версия от 09:36, 6 октября 2012
Строка 16:		Строка 16:

	Для начала создадим внешний контур изображения. Тут есть два способа — простой и сложный, выбор зависит от свойств конкретных данных и частностей задачи.		Для начала создадим внешний контур изображения. Тут есть два способа — простой и сложный, выбор зависит от свойств конкретных данных и частностей задачи.
	* '''Простой~~'''~~ способ ~~– простая~~ рамка изображения~~, которая~~ создаётся утилитой '''gdaltindex'''.		===Простой способ===
	Например, для всех растров в папках, начинающихся с «MRGN_» :		Простая рамка изображения создаётся утилитой '''gdaltindex'''.

			Например, для всех растров в папках, начинающихся с «MRGN_»:
	for /d %b in (MRGN_) do @for %a in (%b\.tif) do if exist "%~fa" gdaltindex -skip_different_projection MRGN_SP5_2011_index.shp "%~fa"		for /d %b in (MRGN_) do @for %a in (%b\.tif) do if exist "%~fa" gdaltindex -skip_different_projection MRGN_SP5_2011_index.shp "%~fa"
	для одной папки с растрами:		для одной папки с растрами:
	gdaltindex -skip_different_projection index.shp "*.tif"		gdaltindex -skip_different_projection index.shp "*.tif"

	* '''Сложный~~'''~~ способ~~, который учитывает~~, что ~~растр~~ может занимать только часть ~~изображения~~, а часть растра нас не интересует, т.к. ~~имеет~~ «черное» поле (NODATA=0).		===Сложный способ===
			Нужно учитывать, что значимая часть растра может занимать только часть матрицы, оставшаяся часть растра нас не интересует, т.к. представляет собой «черное» поле (NODATA).

	1. Cоздадим новый растр из предположения, что есть NODATA=0 и все остальное > 0. Используем утилиту '''gdal_calc''', чтобы получить ненулевые значения (в примере обработка целой папки с растрами в фрмате TIF):		1. Cоздадим новый растр из предположения, что есть NODATA=0 и все остальное > 0. Используем утилиту '''gdal_calc''', чтобы получить ненулевые значения (в примере обработка целой папки с растрами в фрмате TIF):

Максим Дубинин: /* Создание контура изображения */

2012-10-06T07:31:06Z

Создание контура изображения

← Предыдущая версия		Версия от 09:31, 6 октября 2012
Строка 19:		Строка 19:
	Например, для всех растров в папках, начинающихся с «MRGN_» :		Например, для всех растров в папках, начинающихся с «MRGN_» :
	for /d %b in (MRGN_) do @for %a in (%b\.tif) do if exist "%~fa" gdaltindex -skip_different_projection MRGN_SP5_2011_index.shp "%~fa"		for /d %b in (MRGN_) do @for %a in (%b\.tif) do if exist "%~fa" gdaltindex -skip_different_projection MRGN_SP5_2011_index.shp "%~fa"
	для одной папки с растрами ~~синтаксис имеет вид~~:		для одной папки с растрами:
	gdaltindex -skip_different_projection index.shp "*.tif"		gdaltindex -skip_different_projection index.shp "*.tif"

	* '''Сложный''' способ, который учитывает, что растр может занимать только часть изображения, а часть растра нас не интересует, т.к. имеет «черное» поле (NODATA=0).		* '''Сложный''' способ, который учитывает, что растр может занимать только часть изображения, а часть растра нас не интересует, т.к. имеет «черное» поле (NODATA=0).

	1. Cоздадим новый растр из предположения, что есть NODATA=0 и все остальное > 0. Используем утилиту '''gdal_calc''', чтобы ~~вычленить~~ ненулевые значения (в примере ~~снова~~ обработка целой папки с растрами TIF):		1. Cоздадим новый растр из предположения, что есть NODATA=0 и все остальное > 0. Используем утилиту '''gdal_calc''', чтобы получить ненулевые значения (в примере обработка целой папки с растрами в фрмате TIF):
	for /r %a in (.tif) do gdal_calc -A "%~fa" --outfile="%~dpna_c.tif" --calc="(A>0)127+(A==0)*0" --NoDataValue=0 --type=Byte –overwrite		for /r %a in (.tif) do gdal_calc -A "%~fa" --outfile="%~dpna_c.tif" --calc="(A>0)127+(A==0)*0" --NoDataValue=0 --type=Byte –overwrite

	2. Сжимаем полученные файлы с помощью опции COMPRESS=CCITTFAX4 с очевидным уменьшением разрядности до 1 бита через '''gdal_translate''':		2. Сжимаем полученные файлы с помощью опции COMPRESS=CCITTFAX4 с очевидным уменьшением разрядности до 1 бита через '''gdal_translate''':
	for /d %b in (MRGN_*) do @for %a in (%b\%~b_c.tif) do gdal_translate -of "GTiff" -a_nodata 127 -co "NBITS=1" -co "TFW=YES" -co "COMPRESS=CCITTFAX4" "%~fa" "%~dpna_msk.tif"		for /d %b in (MRGN_*) do @for %a in (%b\%~b_c.tif) do gdal_translate -of "GTiff" -a_nodata 127 -co "NBITS=1" -co "TFW=YES" -co "COMPRESS=CCITTFAX4" "%~fa" "%~dpna_msk.tif"

Максим Дубинин: /* Постановка задачи */

2012-10-06T07:14:03Z

Постановка задачи

← Предыдущая версия		Версия от 09:14, 6 октября 2012
Строка 1:		Строка 1:
	== Постановка задачи ==		== Постановка задачи ==
	Допустим, что есть большое изображение и есть векторная сетка (полигоны), по которой бы хотелось это изображение разрезать. Сетка должна ~~совпадать с изображением по проекции (или системе~~ координат). Может и не совпадать, но это потребует дополнительных действий.		Допустим, что есть большое изображение и есть векторная сетка (полигоны), по которой бы хотелось это изображение разрезать. Сетка должна иметь такую же систему координат, как и изображение. Может и не совпадать, но это потребует дополнительных действий.

Максим Дубинин: /* Создание контура изображения */

2012-10-06T07:09:01Z

Создание контура изображения

← Предыдущая версия		Версия от 09:09, 6 октября 2012
Строка 22:		Строка 22:
	gdaltindex -skip_different_projection index.shp "*.tif"		gdaltindex -skip_different_projection index.shp "*.tif"

	* '''Сложный''' способ, который учитывает, что растр может занимать только часть изображения, а часть растра нас не интересует, т.к. имеет «черное» поле (NODATA=0). ~~В этом случае, я делаю еще «пару танцев с бубном»:~~		* '''Сложный''' способ, который учитывает, что растр может занимать только часть изображения, а часть растра нас не интересует, т.к. имеет «черное» поле (NODATA=0).

	1. Cоздадим новый растр из предположения, что есть NODATA=0 и все остальное > 0. Используем утилиту '''gdal_calc''', чтобы вычленить ненулевые значения (в примере снова обработка целой папки с растрами TIF):		1. Cоздадим новый растр из предположения, что есть NODATA=0 и все остальное > 0. Используем утилиту '''gdal_calc''', чтобы вычленить ненулевые значения (в примере снова обработка целой папки с растрами TIF):
	for /r %a in (.tif) do gdal_calc -A "%~fa" --outfile="%~dpna_c.tif" --calc="(A>0)127+(A==0)*0" --NoDataValue=0 --type=Byte –overwrite		for /r %a in (.tif) do gdal_calc -A "%~fa" --outfile="%~dpna_c.tif" --calc="(A>0)127+(A==0)*0" --NoDataValue=0 --type=Byte –overwrite

Boris: /* Создание контура изображения */

2012-10-05T23:01:47Z

Создание контура изображения

← Предыдущая версия		Версия от 01:01, 6 октября 2012
Строка 68:		Строка 68:

	Повторим, что нам необходимо было получить:		Повторим, что нам необходимо было получить:

	- 4 координаты 2-х углов сетки, которые лежат только в пределах обрезаемого (скорее разбираемого) растра и имеют ту же СК, что и обрезаемый растр.		- 4 координаты 2-х углов сетки, которые лежат только в пределах обрезаемого (скорее разбираемого) растра и имеют ту же СК, что и обрезаемый растр.
	Если это у нас есть, то вперед - интерационо разберем растр на части.		Если это у нас есть, то вперед - интерационо разберем растр на части.

Boris: /* Создание контура изображения */

2012-10-05T23:01:11Z

Создание контура изображения

← Предыдущая версия		Версия от 01:01, 6 октября 2012
Строка 67:		Строка 67:
	Все «танцы с бубном» завершены.		Все «танцы с бубном» завершены.

	Повторим, что ~~собственно~~ нам ~~обязательно надо~~ получить:		Повторим, что нам необходимо было получить:
	* 4 координаты углов сетки;		- 4 координаты 2-х углов сетки, которые лежат только в пределах обрезаемого (скорее разбираемого) растра и имеют ту же СК, что и обрезаемый растр.
	* эти координаты должны лежать только в пределах обрезаемого (скорее разбираемого) растра;		Если это у нас есть, то вперед - интерационо разберем растр на части.
	* иметь ту же СК, что и обрезаемый растр.

	== Обрезка растра ==		== Обрезка растра ==

Boris: /* Обрезка растра */

2012-10-05T22:57:27Z

Обрезка растра

← Предыдущая версия		Версия от 00:57, 6 октября 2012
Строка 75:		Строка 75:
	Само использование GDAL, а именно, '''gdal_translate''', гораздо короче:		Само использование GDAL, а именно, '''gdal_translate''', гораздо короче:
	# Создаем набор команд для вычисления контуров в виде BAT-файла.		# Создаем набор команд для вычисления контуров в виде BAT-файла.
	# Команда выполняется в каталоге, где лежит файл Q1.MID (так я назвал экспорт координат углов), файлы выполнения сохраняются в каталог root\1\, исходный файл находится по адресу '''"f:\test\test.tif"''':		# Команда выполняется в каталоге, где лежит файл '''Q1.MID''' (так я назвал экспорт координат углов), файлы выполнения сохраняются в каталог '''"root\1\"''', исходный файл находится по адресу '''"f:\test\test.tif"''':
	for /f "tokens=1-5 delims=," %a in (Q1.MID) do @echo @if not exist "root\1\test_%~a.tif" gdal_translate -of GTIFF -projwin %b %c %d %e -a_nodata 0 -co "TFW=YES" -co "compress=LZW" "f:\test\test.tif" "root\1\test_%~a.tif">>root\1\do_%~a.bat		for /f "tokens=1-5 delims=," %a in (Q1.MID) do @echo @if not exist "root\1\test_%~a.tif" gdal_translate -of GTIFF -projwin %b %c %d %e -a_nodata 0 -co "TFW=YES" -co "compress=LZW" "f:\test\test.tif" "root\1\test_%~a.tif">>root\1\do_%~a.bat
	# Можно сразу же запускать команды на выполнение в цикле for командной строки DOS, но создание набора BAT-файлов позволяет воспользоваться параллельным вычислением фрагментов или вернуться к прерванному вычислению, если понадобиться.		# Можно сразу же запускать команды на выполнение в цикле for командной строки DOS, но создание набора BAT-файлов позволяет воспользоваться параллельным вычислением фрагментов или вернуться к прерванному вычислению, если понадобиться.

Boris: /* Обрезка растра */

2012-10-05T22:43:32Z

Обрезка растра

← Предыдущая версия		Версия от 00:43, 6 октября 2012
Строка 75:		Строка 75:
	Само использование GDAL, а именно, '''gdal_translate''', гораздо короче:		Само использование GDAL, а именно, '''gdal_translate''', гораздо короче:
	# Создаем набор команд для вычисления контуров в виде BAT-файла.		# Создаем набор команд для вычисления контуров в виде BAT-файла.
	# Команда выполняется в каталоге, где лежит файл Q1.MID (так я назвал экспорт координат углов), файлы выполнения сохраняются в каталог root\1\, ~~в нем же должна быть воссоздана структура исходных каталогов~~:		# Команда выполняется в каталоге, где лежит файл Q1.MID (так я назвал экспорт координат углов), файлы выполнения сохраняются в каталог root\1\, исходный файл находится по адресу '''"f:\test\test.tif"''':
	for /f "tokens=1-5 delims=," %a in (Q1.MID) do @echo @if not exist "root\1\test_%~a.tif" gdal_translate -of GTIFF -projwin %b %c %d %e -a_nodata 0 -co "TFW=YES" -co "compress=LZW" "f:\test\test.tif" "root\1\test_%~a.tif">>root\1\do_%~a.bat		for /f "tokens=1-5 delims=," %a in (Q1.MID) do @echo @if not exist "root\1\test_%~a.tif" gdal_translate -of GTIFF -projwin %b %c %d %e -a_nodata 0 -co "TFW=YES" -co "compress=LZW" "f:\test\test.tif" "root\1\test_%~a.tif">>root\1\do_%~a.bat
	# Можно сразу же запускать команды на выполнение в цикле for командной строки DOS, но создание набора BAT-файлов позволяет воспользоваться параллельным вычислением фрагментов или вернуться к прерванному вычислению, если понадобиться.		# Можно сразу же запускать команды на выполнение в цикле for командной строки DOS, но создание набора BAT-файлов позволяет воспользоваться параллельным вычислением фрагментов или вернуться к прерванному вычислению, если понадобиться.

Boris: /* Создание контура изображения */

2012-10-05T22:40:39Z

Создание контура изображения

← Предыдущая версия		Версия от 00:40, 6 октября 2012
Строка 44:		Строка 44:

	Теперь откроем контуры планшетов в ГИС. Проделаем с ними следующие превращения (большинство избыточны):		Теперь откроем контуры планшетов в ГИС. Проделаем с ними следующие превращения (большинство избыточны):
	* Увеличим размер каждого планшета на некоторую величину (не обязательно, но если потом придется выводить растр с перекрытием, это даст запас пикселей по краям). Я увеличиваю на 5% ~~от размера планшета~~, ~~пересчитанных~~ в метры.		* Увеличим размер каждого планшета на некоторую величину (не обязательно, но если потом придется выводить растр с перекрытием, это даст запас пикселей по краям). Я увеличиваю размера планшета приблизительно на 5%, для простоты создания буфера я перед созданием буфера пересчитываю их в метры.
	* Переводим полученный вектор в проекцию (или систему координат) растра.		* Переводим полученный вектор в проекцию (или систему координат) растра.
	* Новый растр открываем и находим для каждого контура MBR (minimum bounding rectangle). В MapInfo это стандартная функция, которую можно применить к объектам слоя, затирая их (это команда на MapBasic, которая выполняется в окне MapBasic внутри MapInfo):		* Новый растр открываем и находим для каждого контура MBR (minimum bounding rectangle). В MapInfo это стандартная функция, которую можно применить к объектам слоя, затирая их (это команда на MapBasic, которая выполняется в окне MapBasic внутри MapInfo):

Boris: /* Небольшое отступление */

2012-10-05T22:38:12Z

Небольшое отступление

← Предыдущая версия		Версия от 00:38, 6 октября 2012
Строка 10:		Строка 10:
	=== Небольшое отступление ===		=== Небольшое отступление ===
	Самый известный вариант сетки — это сетка планшетов СК-42 (и аналогичных). Она стандартна, удобно создана для разных масштабов, и имеется у всех, кому нужна. С ней есть только один небольшой минус — она не прямоугольна. А смысл в том, чтобы резать изображение, не нарушая его пиксельной структуры — без пересчетов и искажений. Утилита '''gdalwarp''' позволяет обрезать изображение векторной сеткой, записанной в файл. Это не наш случай, т.к. он производит пересчет пикселей под векторный контур.		Самый известный вариант сетки — это сетка планшетов СК-42 (и аналогичных). Она стандартна, удобно создана для разных масштабов, и имеется у всех, кому нужна. С ней есть только один небольшой минус — она не прямоугольна. А смысл в том, чтобы резать изображение, не нарушая его пиксельной структуры — без пересчетов и искажений. Утилита '''gdalwarp''' позволяет обрезать изображение векторной сеткой, записанной в файл. Это не наш случай, т.к. он производит пересчет пикселей под векторный контур.
			Идея нашей обработки - извлечение прямоугольных участков растра, заданных двумя вершинами. Для этого у '''gdal_translate''' есть два механизма - задать координаты растра и задать координаты в проекции растра, а программа сама пересчитает эти координаты в координаты растра (и выведет их нам при работе).

	== Создание контура изображения ==		== Создание контура изображения ==