Тайлы из Landsat 8 - История изменений

Zverik: /* Паншарпенинг */ опечатка

2015-08-28T18:27:04Z

Паншарпенинг: опечатка

← Предыдущая версия		Версия от 20:27, 28 августа 2015
Строка 33:		Строка 33:
	-lum ${BASE}_B3.TIF 0.5 -lum ${BASE}_B4.TIF 0.5		-lum ${BASE}_B3.TIF 0.5 -lum ${BASE}_B4.TIF 0.5

	Что тут происходит? Первые три параметра передают панхроматический канал (восьмой), значение «нет данных» и имя выходного файла. Следующая тройка с <code>-~~ref~~</code> — каналы для красного, зелёного и синего цветов. В Landsat 8 им соответствуют 4, 3 и 2, но можно выбрать другие, из [http://gis-lab.info/qa/landsat-bandcomb.html невидимых областей спектра].		Что тут происходит? Первые три параметра передают панхроматический канал (восьмой), значение «нет данных» и имя выходного файла. Следующая тройка с <code>-rgb</code> — каналы для красного, зелёного и синего цветов. В Landsat 8 им соответствуют 4, 3 и 2, но можно выбрать другие, из [http://gis-lab.info/qa/landsat-bandcomb.html невидимых областей спектра].

	Наконец, параметры <code>-lum</code> корректируют наложение цвета на панхроматический канал. [https://github.com/gina-alaska/dans-gdal-scripts/wiki/Gdal_landsat_pansharp Описание алгоритма] слишком сложное, мы взяли коэффициенты из [https://gist.github.com/dwtkns/d18828165160384d0fe1 этих функций] Дерека Уоткинса.		Наконец, параметры <code>-lum</code> корректируют наложение цвета на панхроматический канал. [https://github.com/gina-alaska/dans-gdal-scripts/wiki/Gdal_landsat_pansharp Описание алгоритма] слишком сложное, мы взяли коэффициенты из [https://gist.github.com/dwtkns/d18828165160384d0fe1 этих функций] Дерека Уоткинса.

Zverik: /* Изображение */ предупреждения

2015-02-23T11:32:31Z

Изображение: предупреждения

← Предыдущая версия		Версия от 13:32, 23 февраля 2015
Строка 64:		Строка 64:

	convert result.tif -quality 85 result.jpg		convert result.tif -quality 85 result.jpg

			Программа выдаст несколько предупреждений об «unknown field with tag». Это нормально: она не умеет читать поля TIFF с геопривязкой, и пропускает их.

	== Тайлы ==		== Тайлы ==

Zverik: /* Коррекция цвета */ предупреждение

2015-02-23T11:29:25Z

Коррекция цвета: предупреждение

← Предыдущая версия		Версия от 13:29, 23 февраля 2015
Строка 38:		Строка 38:

	== Коррекция цвета ==		== Коррекция цвета ==

			Снимок после паншарпенинга получится либо весь в белой дымке, либо невозможно тёмный. Для его использования обязательна коррекция цвета, хотя бы методом автоматического растягивания гистограммы.

	=== Автоматически ===		=== Автоматически ===

	~~Белую дымку удалить непросто.~~ В GDAL нет команды для обработки цвета, подобной вкладке в QGIS, но тот же аляскинском наборе содержит скрипт, выдающий приемлемый результат:		В GDAL нет команды для обработки цвета, подобной вкладке в QGIS, но тот же аляскинском наборе содержит скрипт, выдающий приемлемый результат:

	gdal_contrast_stretch -ndv 0 -percentile-range 0.05 0.95 result.tif result2.tif		gdal_contrast_stretch -ndv 0 -percentile-range 0.05 0.95 result.tif result2.tif

Максим Дубинин в 12:24, 13 февраля 2015

2015-02-13T12:24:53Z

← Предыдущая версия		Версия от 14:24, 13 февраля 2015
Строка 1:		Строка 1:
	{{Статья\|~~Опубликовано~~\|landsat-tiles}}		{{Статья\|Опубликована\|landsat-tiles}}

	{{Аннотация\|Упрощённая инструкция по получению и обработке снимков Landsat 8 и подготовке из них тайлов}}		{{Аннотация\|Упрощённая инструкция по получению и обработке снимков Landsat 8 и подготовке из них тайлов}}

Максим Дубинин в 12:18, 13 февраля 2015

2015-02-13T12:18:26Z

← Предыдущая версия		Версия от 14:18, 13 февраля 2015
Строка 1:		Строка 1:
	{{Статья\|~~Черновик~~}}		{{Статья\|Опубликовано\|landsat-tiles}}
	{{Аннотация\|Упрощённая инструкция по получению и обработке снимков Landsat 8 и подготовке из них тайлов}}		{{Аннотация\|Упрощённая инструкция по получению и обработке снимков Landsat 8 и подготовке из них тайлов}}

Zverik: window

2015-02-02T07:09:18Z

window

← Предыдущая версия		Версия от 09:09, 2 февраля 2015
Строка 3:		Строка 3:

	Главное свойство данных ДЗЗ Landsat 8 — свободная лицензия, которая позволяет обрисовывать снимки для, в частности, OpenStreetMap. Эта статья не вдаётся в теорию, но содержит шаги для поиска и обработки снимков Landsat 8 и их превращения в тайлы для публикации или обрисовки в JOSM.		Главное свойство данных ДЗЗ Landsat 8 — свободная лицензия, которая позволяет обрисовывать снимки для, в частности, OpenStreetMap. Эта статья не вдаётся в теорию, но содержит шаги для поиска и обработки снимков Landsat 8 и их превращения в тайлы для публикации или обрисовки в JOSM.

			Некоторые шаги будет непросто выполнить в ОС Windows. Для пользователей этой системы советуем [http://gis-lab.info/qa/qgis-osgeo4w.html установить OSGeo4W] и делать [http://courses.neteler.org/processing-landsat-8-data-in-grass-gis-7-rgb-composites-and-pan-sharpening/ паншарпенинг и цветовую коррекцию в GRASS].

	== Получение снимка ==		== Получение снимка ==

Zverik: /* Коррекция цвета */ привязка

2015-02-01T21:56:48Z

Коррекция цвета: привязка

← Предыдущая версия		Версия от 23:56, 1 февраля 2015
Строка 35:		Строка 35:

	== Коррекция цвета ==		== Коррекция цвета ==

			=== Автоматически ===

	Белую дымку удалить непросто. В GDAL нет команды для обработки цвета, подобной вкладке в QGIS, но тот же аляскинском наборе содержит скрипт, выдающий приемлемый результат:		Белую дымку удалить непросто. В GDAL нет команды для обработки цвета, подобной вкладке в QGIS, но тот же аляскинском наборе содержит скрипт, выдающий приемлемый результат:
Строка 41:		Строка 43:

	Доли отбрасываемых цветов можно покрутить на ±0.03: чем ближе к 0-100, тем темнее получается картинка, но тем больше деталей сохраняется.		Доли отбрасываемых цветов можно покрутить на ±0.03: чем ближе к 0-100, тем темнее получается картинка, но тем больше деталей сохраняется.

			=== В графическом редакторе ===

			Раз tif — обычная картинка с доп. полями, можно открыть его в [http://imagej.nih.gov/ij/ ImageJ] или ином редакторе и подкрутить контраст вручную, акцентируя нужные полосы. Чтобы сохранить геопривязку, понадобятся программы из пакета ''libgeotiff'' (иногда ''geotiff-bin''). Экспорт геопривязки в файл:

			listgeo result.tif > result.geo

			Восстановление геопривязки для отредактированного файла (сохранять нужно, понятно, тоже в формате TIFF):

			geotifcp -g result.geo edited.tif edited_geo.tif

	== Изображение ==		== Изображение ==

Zverik: /* Паншарпенинг */ переделал на переменную

2015-02-01T21:20:36Z

Паншарпенинг: переделал на переменную

← Предыдущая версия		Версия от 23:20, 1 февраля 2015
Строка 23:		Строка 23:
	-rgb B4.TIF -rgb B3.TIF -rgb B2.TIF -lum B3.TIF 0.5 -lum B4.TIF 0.5		-rgb B4.TIF -rgb B3.TIF -rgb B2.TIF -lum B3.TIF 0.5 -lum B4.TIF 0.5

	Здесь опущены длинные идентификаторы снимков, вида «LC81850172014281LGN00». ~~С ними будет длиннее и менее понятно~~:		Здесь опущены длинные идентификаторы снимков, вида «LC81850172014281LGN00». Для удобства пригодится такой скрипт:

	gdal_landsat_pansharp -pan ~~LC81850172014281LGN00_B8~~.TIF -ndv 0 -o result.tif \		BASE=LC81850172014281LGN00
	-rgb ~~LC81850172014281LGN00_B4~~.TIF -rgb ~~LC81850172014281LGN00_B3~~.TIF -rgb ~~LC81850172014281LGN00_B2~~.TIF \		gdal_landsat_pansharp -pan ${BASE}_B8.TIF -ndv 0 -o result.tif \
	-lum ~~LC81850172014281LGN00_B3~~.TIF 0.5 -lum ~~LC81850172014281LGN00_B4~~.TIF 0.5		-rgb ${BASE}_B4.TIF -rgb ${BASE}_B3.TIF -rgb ${BASE}_B2.TIF \
			-lum ${BASE}_B3.TIF 0.5 -lum ${BASE}_B4.TIF 0.5

	Что тут происходит? Первые три параметра передают панхроматический канал (восьмой), значение «нет данных» и имя выходного файла. Следующая тройка с <code>-ref</code> — каналы для красного, зелёного и синего цветов. В Landsat 8 им соответствуют 4, 3 и 2, но можно выбрать другие, из [http://gis-lab.info/qa/landsat-bandcomb.html невидимых областей спектра].		Что тут происходит? Первые три параметра передают панхроматический канал (восьмой), значение «нет данных» и имя выходного файла. Следующая тройка с <code>-ref</code> — каналы для красного, зелёного и синего цветов. В Landsat 8 им соответствуют 4, 3 и 2, но можно выбрать другие, из [http://gis-lab.info/qa/landsat-bandcomb.html невидимых областей спектра].

Zverik: /* Паншарпенинг */ лишний пробел

2015-02-01T21:16:53Z

Паншарпенинг: лишний пробел

← Предыдущая версия		Версия от 23:16, 1 февраля 2015
Строка 20:		Строка 20:
	У нас есть 3 канала с разрешением 30 м/пк и один — 15 м/пк. Паншарпенинг — это объединение их всех в один цветной снимок максимального разрешения. Раньше его делали длинным цепочками команд GDAL или ImageMagick, но нам повезло. Умные люди из GINA (университет Аляски) написали [https://github.com/gina-alaska/dans-gdal-scripts/ скрипт], и теперь достаточно одной команды. Набор аляскинских скриптов есть в большинстве дистрибутивов, поищите пакет ''dans-gdal-scripts''. Команда такая:		У нас есть 3 канала с разрешением 30 м/пк и один — 15 м/пк. Паншарпенинг — это объединение их всех в один цветной снимок максимального разрешения. Раньше его делали длинным цепочками команд GDAL или ImageMagick, но нам повезло. Умные люди из GINA (университет Аляски) написали [https://github.com/gina-alaska/dans-gdal-scripts/ скрипт], и теперь достаточно одной команды. Набор аляскинских скриптов есть в большинстве дистрибутивов, поищите пакет ''dans-gdal-scripts''. Команда такая:

	gdal_landsat_pansharp -pan B8.TIF -ndv 0 -o result.tif \		gdal_landsat_pansharp -pan B8.TIF -ndv 0 -o result.tif \
	-rgb B4.TIF -rgb B3.TIF -rgb B2.TIF -lum B3.TIF 0.5 -lum B4.TIF 0.5		-rgb B4.TIF -rgb B3.TIF -rgb B2.TIF -lum B3.TIF 0.5 -lum B4.TIF 0.5

Zverik: /* Получение снимка */ tar

2015-02-01T21:14:14Z

Получение снимка: tar

← Предыдущая версия		Версия от 23:14, 1 февраля 2015
Строка 10:		Строка 10:
	http://libra.developmentseed.org		http://libra.developmentseed.org

	Напоминаем, что каждый снимок весит 500-800 мегабайт в архиве (2 гигабайта вне его), и для обработки потребуются ещё гигабайта четыре.		Напоминаем, что каждый снимок весит 500-800 мегабайт в архиве (2 гигабайта вне его), и для обработки потребуются ещё гигабайта четыре. Под линуксом файл распаковывается такой командой:

			tar -xjvf LC*.tar.bz

	В архиве вы получите [http://geektimes.ru/post/183416/ 11 каналов] и слой QA, на котором отмечены точки без информации о земной поверхности (например, облака). Мы будем работать с каналами 2, 3, 4 (синий, зелёный, красный) и 8 (ч/б с 15-метровым разрешением).		В архиве вы получите [http://geektimes.ru/post/183416/ 11 каналов] и слой QA, на котором отмечены точки без информации о земной поверхности (например, облака). Мы будем работать с каналами 2, 3, 4 (синий, зелёный, красный) и 8 (ч/б с 15-метровым разрешением).