Примеры использования инструментов GDAL - История изменений

Максим Дубинин: /* Создание растров уменьшенного разрешения (т.н. quicklook, preview) */

2023-11-04T13:42:50Z

Создание растров уменьшенного разрешения (т.н. quicklook, preview)

← Предыдущая версия		Версия от 15:42, 4 ноября 2023
Строка 286:		Строка 286:

	<pre>gdalwarp -tr 8000 8000 -r average input.tif output.tif</pre>		<pre>gdalwarp -tr 8000 8000 -r average input.tif output.tif</pre>

			или, если на выходе нужен растр определенного размера:

			<pre>gdalwarp -ts 8000 8000 -r average input.tif output.tif</pre>

	==Ссылки по теме==		==Ссылки по теме==
	*[http://www.gdal.org/index.html Библиотеки GDAL/OGR]		*[http://www.gdal.org/index.html Библиотеки GDAL/OGR]
	*[http://gis-lab.info/qa/ogr2ogr-examples.html Примеры использования ogr2ogr]		*[http://gis-lab.info/qa/ogr2ogr-examples.html Примеры использования ogr2ogr]

Александр Мурый: /* Работа с цветовой таблицей */

2017-08-09T09:43:46Z

Работа с цветовой таблицей

← Предыдущая версия		Версия от 11:43, 9 августа 2017
Строка 188:		Строка 188:

	==Работа с цветовой таблицей==		==Работа с цветовой таблицей==
	Если на входе есть растр в ~~Grayscale~~ а на выходе нужно получить изображение с индексированными цветами (палитрой) то можно воспользоваться подходом через VRT.		Если на входе есть растр с цветами в шкале серого (grayscale), а на выходе нужно получить изображение с индексированными цветами (палитрой),о можно воспользоваться подходом через VRT.

			=== С использованием VRT ===

	1. Создаем файл VRT:		1. Создаем файл VRT:
Строка 210:		Строка 212:
	</pre>		</pre>

	Каждая Entry - это цвет кодированный RGB соответствующий по порядку значениям от 0 до 255. Если цветов меньше 255, то остальные будут добавлены автоматически с RGB = 0,0,0 (черный).		Каждая Entry - это цвет, кодированный в RGB и соответствующий по порядку значениям от 0 до 255. Если цветов меньше 255, то остальные будут добавлены автоматически с RGB = 0,0,0 (черный).

	После того как таблица отредактирована, нужно сохранить ее обратно в растр:		После того как таблица отредактирована, нужно сохранить ее обратно в растр:

Максим Дубинин: /* Примеры конвертации */

2016-02-22T10:38:14Z

Примеры конвертации

← Предыдущая версия		Версия от 12:38, 22 февраля 2016
Строка 82:		Строка 82:
	*EHdr (rw+): ESRI .hdr Labelled		*EHdr (rw+): ESRI .hdr Labelled

	==~~Примеры конвертации~~==		==Конвертация==

	Извлечь три канала с номерами 1, 2, 3 в новый файл из исходного с перекомбинацией, в котором каналов может быть больше.		Извлечь три канала с номерами 1, 2, 3 в новый файл из исходного с перекомбинацией, в котором каналов может быть больше.
Строка 102:		Строка 102:
	Конвертация 16 битного одноканального растра в 8 битный:		Конвертация 16 битного одноканального растра в 8 битный:

	<pre>gdal_translate -scale -ot Byte ~~input~~.tif output.tif</pre>		<pre>gdal_translate -scale -ot Byte input_16bit.tif output.tif</pre>

	Пакетная конвертация всех JPG в TIF (Windows):		Пакетная конвертация всех JPG в TIF (Windows):

Максим Дубинин: /* Конвертация с перепроецированием */

2016-02-06T05:45:06Z

Конвертация с перепроецированием

← Предыдущая версия		Версия от 07:45, 6 февраля 2016
Строка 152:		Строка 152:
	gdalwarp позволяет не только конвертировать данные из одного формата в другой, но и одновременно произвести перепроецирование данных из одной системы координат в другую. Для этого используются параметры:		gdalwarp позволяет не только конвертировать данные из одного формата в другой, но и одновременно произвести перепроецирование данных из одной системы координат в другую. Для этого используются параметры:

	*-a_srs используется для указания системы координат для данных		*-a_srs используется для указания системы координат (СК) для данных
	*-s_srs используется для перезаписи информации о системе координат		*-s_srs используется для перезаписи информации о системе координат
	*-t_srs перепроецирования данных в требуемую систему координат		*-t_srs перепроецирования данных в требуемую систему координат

	~~Например~~, ~~перепроецировать из проекции растра в проекцию Альберса~~ можно так:		В самом простейшем случае это делается следующим образом (ключ -tr указывает разрешение целевого растра):

			<pre>gdalwarp.exe -tr 0.0083 0.0083 -t_srs "EPSG:4326" in.tif out.tif</pre>

			Если EPSG-кода у СК конечного растра нет, то указать целевую СК можно так:

	<pre>gdalwarp.exe -t_srs "+proj=aea +lat_1=52 +lat_2=64 +lat_0=0 +lon_0=45 +x_0=8500000 +y_0=0 +ellps=krass +units=m +towgs84=28,-130,-95,0,0,0,0 +no_defs" in.tif out.tif</pre>		<pre>gdalwarp.exe -t_srs "+proj=aea +lat_1=52 +lat_2=64 +lat_0=0 +lon_0=45 +x_0=8500000 +y_0=0 +ellps=krass +units=m +towgs84=28,-130,-95,0,0,0,0 +no_defs" in.tif out.tif</pre>

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

2016-01-30T04:35:55Z

Переназначение системы координат

← Предыдущая версия		Версия от 06:35, 30 января 2016
Строка 139:		Строка 139:

	<pre>gdal_edit -a_srs "EPSG:4326" input.tif</pre>		<pre>gdal_edit -a_srs "EPSG:4326" input.tif</pre>

			Если код неизвестен, можно указать описание системы координат в формате Proj.4:

			<pre>gdal_edit -a_srs "+proj=sinu +lon_0=0 +x_0=0 +y_0=0 +a=6371007.181 +b=6371007.181 +units=m +no_defs" input.tif</pre>

			Узнать строку описания системы координат можно следующим образом:

			<pre>gdalsrsinfo template.tif</pre>

	===Конвертация с перепроецированием===		===Конвертация с перепроецированием===

Antonio.2: remove duplicate gdal_calc

2016-01-20T15:08:53Z

remove duplicate gdal_calc

← Предыдущая версия		Версия от 17:08, 20 января 2016
Строка 14:		Строка 14:
	Утилиты GDAL предназначены для конвертации растровых данных из одного формата в другой и выполнения над ними различных операций. Установить утилиты GDAL для Windows можно с [http://gis-lab.info/qa/qgis-osgeo4w.html помощью OSGeo4W]. Удобный визуальный интерфейс для утилит имеется в [http://gis-lab.info/qa/gdaltools.html GDALTools для QGIS]. В комплект GDAL входят следующие утилиты:		Утилиты GDAL предназначены для конвертации растровых данных из одного формата в другой и выполнения над ними различных операций. Установить утилиты GDAL для Windows можно с [http://gis-lab.info/qa/qgis-osgeo4w.html помощью OSGeo4W]. Удобный визуальный интерфейс для утилит имеется в [http://gis-lab.info/qa/gdaltools.html GDALTools для QGIS]. В комплект GDAL входят следующие утилиты:

	*[http://www.gdal.org/gdal_calc.html gdal_calc] - растровый калькулятор, арифметические операции с растрами;		*[http://www.gdal.org/gdal_calc.html gdal_calc.py] - растровый калькулятор, арифметические операции с растрами;
	*[http://www.gdal.org/gdal-config.html gdal-config] - получить опции необходимые для создания ПО использующего GDAL;		*[http://www.gdal.org/gdal-config.html gdal-config] - получить опции необходимые для создания ПО использующего GDAL;
	*[http://www.gdal.org/gdal_contour.html gdal_contour] - получение изолиний по цифровым моделям рельефа (ЦМР);		*[http://www.gdal.org/gdal_contour.html gdal_contour] - получение изолиний по цифровым моделям рельефа (ЦМР);
	*[http://www.gdal.org/gdal_calc.html gdal_calc] - растровый калькулятор, арифметические операции с растрами;
	*[http://www.gdal.org/gdal_edit.html gdal_edit.py ] - редактировать информацию о растре;		*[http://www.gdal.org/gdal_edit.html gdal_edit.py ] - редактировать информацию о растре;
	*[http://www.gdal.org/gdal_fillnodata.html gdal_fillnodata.py] - заполнение областей имеющих значение NODATA;		*[http://www.gdal.org/gdal_fillnodata.html gdal_fillnodata.py] - заполнение областей имеющих значение NODATA;

Максим Дубинин: /* Стэки, композиты, мозаики */

2016-01-04T22:53:57Z

Стэки, композиты, мозаики

← Предыдущая версия		Версия от 00:53, 5 января 2016
Строка 110:		Строка 110:
	==Стэки, композиты, мозаики==		==Стэки, композиты, мозаики==

	Создание композитного изображения из серии отдельных растров, каждый из которых в своем файле TIF. Разрешение выходного файла устанавливается по первому их растров. Таким образом, если первый канал 15 м, а остальные 30 м, то последние будут пересчитаны на 15 м. Чтобы указать, что ~~производится помещение каждого растра~~ в ~~свой~~ слой, а не мозаицирование, используется ключ -separate:		Создание композитного изображения из серии отдельных растров, каждый из которых в своем файле TIF. Разрешение выходного файла устанавливается по первому их растров. Таким образом, если первый канал 15 м, а остальные 30 м, то последние будут пересчитаны на 15 м. Чтобы указать, что каждый исходный растр должен попасть в отдельный слой (layer-stacking), а не мозаицирование, используется ключ -separate:

	<pre>gdal_merge.py -o output.tif band1.tif band2.tif band3.tif band4.tif band5.tif -separate</pre>		<pre>gdal_merge.py -o output.tif band1.tif band2.tif band3.tif band4.tif band5.tif -separate</pre>

Максим Дубинин в 13:31, 8 декабря 2015

2015-12-08T13:31:35Z

← Предыдущая версия		Версия от 15:31, 8 декабря 2015
Строка 25:		Строка 25:
	*[http://www.gdal.org/gdal_rasterize.html gdal_rasterize] - растеризация векторных данных;		*[http://www.gdal.org/gdal_rasterize.html gdal_rasterize] - растеризация векторных данных;
	*[http://www.gdal.org/gdal_retile.html gdal_retile.py] - создать новый набор тайлов и/или перестроить пирамидные слои;		*[http://www.gdal.org/gdal_retile.html gdal_retile.py] - создать новый набор тайлов и/или перестроить пирамидные слои;
			*[http://www.gdal.org/gdal_grid.html gdal_grid] - создание ЦМР из векторных данных;
	*[http://www.gdal.org/gdal_sieve.html gdal_sieve.py] - фильтрация осколочных объектов растра;		*[http://www.gdal.org/gdal_sieve.html gdal_sieve.py] - фильтрация осколочных объектов растра;
	*[http://www.gdal.org/gdal_translate.html gdal_translate] - конвертация растров из формата в формат;		*[http://www.gdal.org/gdal_translate.html gdal_translate] - конвертация растров из формата в формат;

Максим Дубинин: /* HDF4 */

2015-08-19T23:26:51Z

HDF4

← Предыдущая версия		Версия от 01:26, 20 августа 2015
Строка 235:		Строка 235:
	==HDF4==		==HDF4==

	В HDF4 распространяется множество данных дистанционного зондирования, например MODIS и ASTER~~. Использовать -geoloc для перепроецирования не нужно~~.		В HDF4 распространяется множество данных дистанционного зондирования, например MODIS и ASTER.

	Импорт данных ASTER L1A:		Импорт данных ASTER L1A:

	<pre>gdalwarp -overwrite -of GTiff HDF4_EOS:EOS_SWATH:"110601_081441.hdf":VNIR_Band1:ImageData b1.tif </pre>		<pre>gdalwarp -overwrite -of GTiff HDF4_EOS:EOS_SWATH:"110601_081441.hdf":VNIR_Band1:ImageData b1.tif </pre>

			Импорт данных MODIS Level 3,4:

			<pre>gdal_translate HDF4_EOS:EOS_GRID:"MCD12Q1.A2001001.h00v08.051.2014287161513.hdf":MOD12Q1:Land_Cover_Type_1 output.tif</pre>

			Если одна из частей названия SDS (массива данных внутри HDF) имеет пробелы, ее нужно взять в кавычки. Использовать -geoloc для перепроецирования не нужно.

	==Расчеты==		==Расчеты==

Максим Дубинин: /* Расчеты */

2015-08-17T19:52:49Z

Расчеты

← Предыдущая версия		Версия от 21:52, 17 августа 2015
Строка 251:		Строка 251:

	<pre>gdal_calc.bat -A input.tif --outfile=output.tif --calc="1*(logical_and(A>=249,A<255)) " --NoDataValue=0'</pre>		<pre>gdal_calc.bat -A input.tif --outfile=output.tif --calc="1*(logical_and(A>=249,A<255)) " --NoDataValue=0'</pre>

			Следует обратить внимание, что часть ELSE выражения по умолчанию всегда будет сбрасывать остальные значения в 0. При этом, часть --NoDataValue=0 только говорит к какому значению нужно "приклеить" ярлык NODATA. Таким образом, если в выражении выше сказать, например, --NoDataValue=255, то, несмотря на то, что NODATA установится на 255, значения не удовлетворяющие условию в calc вовсе не станут равны 255, а останутся равны 0. Что бы установить остальные значения в число отличное от 0, нужно сделать это явным образом, например:

			<pre>gdal_calc --overwrite -A input.tif --outfile=input.tif --calc="1(A==0) + 255(A>0)" --NoDataValue=255</pre>

			Присвоить пикселям со значением 0 значение 1, а остальным присвоить 255, установив на это значение тег NODATA.

	==Создание растров уменьшенного разрешения (т.н. quicklook, preview)==		==Создание растров уменьшенного разрешения (т.н. quicklook, preview)==