Ортокоррекция данных OrbView-3 с помощью GDAL - История изменений

Bishop: Исправление про артефакты

2012-06-17T16:48:22Z

Исправление про артефакты

← Предыдущая версия		Версия от 18:48, 17 июня 2012
Строка 292:		Строка 292:
	На изображении видно, что при выполнении ортокоррекции без учета геоида отмечается сдвиг трека в сторону от дороги. Результат выполнения ортокоррекции с учетом геоида полностью идентичен результатам работы ПО ENVI EX.		На изображении видно, что при выполнении ортокоррекции без учета геоида отмечается сдвиг трека в сторону от дороги. Результат выполнения ортокоррекции с учетом геоида полностью идентичен результатам работы ПО ENVI EX.

	Однако, на изображении, полученном в результате ортокорреции при помощи GDAL, имеются артефакты. Причем, ортокорреция, выполненная в программе wxGIS, тоже основанной на библиотеке GDAL, не приводит к появлению артефактов. Данное необъяснимое поведение библиотеки (версия GDAL 1.9.0, released 2011/12/29) требует отдельного изучения.		Однако, на изображении, полученном в результате ортокорреции при помощи GDAL, имеются артефакты. Причем, ортокорреция, выполненная в программе wxGIS, тоже основанной на библиотеке GDAL, не приводит к появлению артефактов. <strike>Данное необъяснимое поведение библиотеки (версия GDAL 1.9.0, released 2011/12/29) требует отдельного изучения</strike>. Для устранения артефактов в командую строку необходимо добавить ключ -et 0.0. Пример команды:

			> gdalwarp -dstnodata 0 -srcnodata 0 -overwrite -t_srs epsg:32635 -wo "INIT_DEST=NO_DATA" -et 0.0 -rpc -to "RPC_DEM=D:\temp\rb\DEM_merged.tif" -to "RPC_HEIGHT=22.0157"
			D:\temp\rb\3v050909p0000897861a520004700712m_001631680.org.tif D:\temp\rb\3v050909p0000897861a520004700712m_001631680.rec2.tif

			Данный ключ отключает излишнюю оптимизацию при интерполяции растра.

	[[Файл:Gdal ortho paradox.png\|740px]]		[[Файл:Gdal ortho paradox.png\|740px]]

Максим Дубинин в 08:57, 27 февраля 2012

2012-02-27T08:57:28Z

← Предыдущая версия		Версия от 10:57, 27 февраля 2012
Строка 1:		Строка 1:
	{{Статья\|~~Черновик~~}}		{{Статья\|Опубликована\|orbview3-ortho-gdal}}
			{{Аннотация\|Как произвести ортокоррекцию данных OrbView-3 с помощью GDAL}}

	== Введение ==		== Введение ==

Максим Дубинин: переименовал Ортокоррекция космических снимков с КА OrbView-3 с помощью GDAL в Ортокоррекция данных OrbView-3 с помощью GDAL

2012-02-27T08:51:32Z

переименовал Ортокоррекция космических снимков с КА OrbView-3 с помощью GDAL в Ортокоррекция данных OrbView-3 с помощью GDAL

← Предыдущая версия	Версия от 10:51, 27 февраля 2012
(нет различий)

Максим Дубинин: переименовал Ортокоррекция космических снимков с КА OrbView-3 в Ортокоррекция космических снимков с КА OrbView-3 с помощью GDAL

2012-02-27T08:47:58Z

переименовал Ортокоррекция космических снимков с КА OrbView-3 в Ортокоррекция космических снимков с КА OrbView-3 с помощью GDAL

← Предыдущая версия	Версия от 10:47, 27 февраля 2012
(нет различий)

Voltron: опечатки

2012-02-26T19:15:34Z

опечатки

← Предыдущая версия		Версия от 21:15, 26 февраля 2012
Строка 123:		Строка 123:
	Для нашего случая система координат будет '''WGS 84 / UTM zone 35N''' или '''epsg:32635'''.		Для нашего случая система координат будет '''WGS 84 / UTM zone 35N''' или '''epsg:32635'''.

	~~Для получения~~ необходимых данных о рельефе рассмотрим ~~способ~~ с данными '''SRTM'''. Для получения файла рельефа необходимо:		Получение необходимых данных о рельефе рассмотрим на примере работы с данными '''SRTM'''. Для получения файла рельефа необходимо:
	<ol>		<ol>
	<li> На странице [http://srtm.csi.cgiar.org/SELECTION/inputCoord.asp SRTM Data Selection Options] выписываем номера рядов и столбцов с территорией съемки.</li>		<li> На странице [http://srtm.csi.cgiar.org/SELECTION/inputCoord.asp SRTM Data Selection Options] выписываем номера рядов и столбцов с территорией съемки.</li>
Строка 134:		Строка 134:
	Команда для операционной системы UNIX:		Команда для операционной системы UNIX:
	gdalbuildvrt srtm.vrt srtm*tif		gdalbuildvrt srtm.vrt srtm*tif
	Команда для операционной системы ~~UNIX~~:		Команда для операционной системы Windows:
	gdalbuildvrt.exe srtm.vrt srtm*tif		gdalbuildvrt.exe srtm.vrt srtm*tif
	</ol>		</ol>

Максим Дубинин в 18:41, 26 февраля 2012

2012-02-26T18:41:48Z

Bishop в 09:51, 26 февраля 2012

2012-02-26T09:51:22Z

← Предыдущая версия		Версия от 11:51, 26 февраля 2012
Строка 6:		Строка 6:

	Сами космические снимки с КА OrbView-3 свободно распространяются и могут быть получены через [[Каталог данных OrbView-3\|каталог данных OrbView-3]].		Сами космические снимки с КА OrbView-3 свободно распространяются и могут быть получены через [[Каталог данных OrbView-3\|каталог данных OrbView-3]].
	Однако, в исходном виде снимки не имеют географической привязки - в формат распространения TIFF не внедрены специальные теги для описания привязки, т.е. формат не представляет из себя классического GeoTIFF. В комплекте поставки имеются все необходимые метаданные для осуществления привязки и ортокоррекции снимков. Как это осуществить будет рассказано далее.		Однако, в исходном виде, снимки не имеют географической привязки - в формат распространения TIFF не внедрены специальные теги для описания привязки, т.е. формат не представляет из себя классического GeoTIFF. В комплекте поставки имеются все необходимые метаданные для осуществления привязки и ортокоррекции снимков. Как это осуществить будет рассказано далее.

	== Программное обеспечение и данные ==		== Программное обеспечение и данные ==
Строка 279:		Строка 279:
	D:\temp\rb\3v050909p0000897861a520004700712m_001631680.org.tif D:\temp\rb\3v050909p0000897861a520004700712m_001631680.rec.tif		D:\temp\rb\3v050909p0000897861a520004700712m_001631680.org.tif D:\temp\rb\3v050909p0000897861a520004700712m_001631680.rec.tif

	Модель RPC использует высотные отметки в системе координат WGS84, а большинство данных о рельефе (DEM) используют относительные отметки геоида. Для устранения этого расхождения необходимо определить превышения геоида над эллипсоидом WGS84 и использовать это значение в команде ортокоррекции. Для определения превышения, возьмем координаты средней точки снимка (27d 2'55.34"E, 55d 1'14.26"N) и загрузим на этот онлайн [http://www.unavco.org/community_science/science-support/geoid/geoid.html ресурс] или на этот онлайн [http://geographiclib.sourceforge.net/cgi-bin/GeoidEval ресурс] (~~Такой~~ расчет может быть выполнен как с помощью онлайн ресурсов, так и с использованием программного обеспечения, например, proj4). В результате получим '''22.0157 м''' для модели EGM96.		Модель RPC использует высотные отметки в системе координат WGS84, а большинство данных о рельефе (DEM) используют относительные отметки геоида. Для устранения этого расхождения необходимо определить превышения геоида над эллипсоидом WGS84 и использовать это значение в команде ортокоррекции. Для определения превышения, возьмем координаты средней точки снимка (27d 2'55.34"E, 55d 1'14.26"N) и загрузим на этот онлайн [http://www.unavco.org/community_science/science-support/geoid/geoid.html ресурс] или на этот онлайн [http://geographiclib.sourceforge.net/cgi-bin/GeoidEval ресурс] (такой расчет может быть выполнен как с помощью онлайн ресурсов, так и с использованием программного обеспечения, например, proj4). В результате получим '''22.0157 м''' для модели EGM96.

	Для ортокоррекции с учетом геоида выполним следующую команду.		Для ортокоррекции с учетом геоида выполним следующую команду.
Строка 295:		Строка 295:
	На изображении видно, что при выполнении ортокоррекции без учета геоида отмечается сдвиг трека в сторону от дороги. Результат выполнения ортокоррекции с учетом геоида полностью идентичен результатам работы ПО ENVI EX.		На изображении видно, что при выполнении ортокоррекции без учета геоида отмечается сдвиг трека в сторону от дороги. Результат выполнения ортокоррекции с учетом геоида полностью идентичен результатам работы ПО ENVI EX.

	Однако, на изображении, полученном в результате ортокорреции при помощи GDAL имеются артефакты. Причем, ортокорреция, выполненная в программе wxGIS, тоже основанной на библиотеке GDAL, не приводит к появлению артефактов. Данное необъяснимое поведение библиотеки (версия GDAL 1.9.0, released 2011/12/29) требует отдельного изучения.		Однако, на изображении, полученном в результате ортокорреции при помощи GDAL, имеются артефакты. Причем, ортокорреция, выполненная в программе wxGIS, тоже основанной на библиотеке GDAL, не приводит к появлению артефактов. Данное необъяснимое поведение библиотеки (версия GDAL 1.9.0, released 2011/12/29) требует отдельного изучения.

	[[Файл:Gdal ortho paradox.png\|740px]]		[[Файл:Gdal ortho paradox.png\|740px]]

Bishop: Окончательно

2012-02-26T00:00:30Z

Окончательно

Внесённые изменения

Bishop: переписал первых 2 главы

2012-02-25T19:41:59Z

переписал первых 2 главы

@@ Строка 1: / Строка 1: @@
 {{Статья|Черновик}}
-=== Краткое руководство ===
+== Введение ==
-Три простых шага для получения компактного восьмибитного растра из исходного:
+В данной статье рассматривается способ ортокоррекции космических снимков с КА OrbView-3 при помощи библиотеки [http://gdal.org GDAL].
-. Скачиваем нужные участки рельефа SRTM в GeoTIFF, для орторектификации (распрямления) по ним снимков. Для этого:
+Сами космические снимки с КА OrbView-3 свободно распространяются и могут быть получены через [[Каталог данных OrbView-3|каталог данных OrbView-3]].
-* на странице http://srtm.csi.cgiar.org/SELECTION/inputCoord.asp смотрим на картинку и выписываем куда-нибудь номера рядов и столбцов с интересующей нас местностью;
+Однако, в исходном виде снимки не имеют географической привязки - в формат распространения TIFF не внедрены специальные теги для описания привязки, т.е. формат не представляет из себя классического GeoTIFF. В комплекте поставки имеются все необходимые метаданные для осуществления привязки и ортокоррекции снимков. Как это осуществить будет рассказано далее.
-. Для удобства объединяем все файлы SRTM в формате GeoTIFF в единый виртуальный растр:
+== Программное обеспечение и данные ==
    gdalbuildvrt srtm.vrt srtm*tif
 . Для получения компактного восьмибитного геотиффа из исходного сырого делаем следующее:

Bishop: переименовал Орторектификация OrbView-3 в Ортокоррекция космических снимков с КА OrbView-3

2012-02-25T17:52:26Z

переименовал Орторектификация OrbView-3 в Ортокоррекция космических снимков с КА OrbView-3

← Предыдущая версия

Версия от 19:52, 25 февраля 2012

← Предыдущая версия		Версия от 20:41, 26 февраля 2012
Строка 5:		Строка 5:
	В данной статье рассматривается способ ортокоррекции космических снимков с КА OrbView-3 при помощи библиотеки [http://gdal.org GDAL].		В данной статье рассматривается способ ортокоррекции космических снимков с КА OrbView-3 при помощи библиотеки [http://gdal.org GDAL].

	~~Сами космические~~ снимки с КА OrbView-3 свободно распространяются и могут быть получены через [[Каталог данных OrbView-3\|каталог данных OrbView-3]].		Космические снимки с КА OrbView-3 свободно распространяются и могут быть получены через [[Каталог данных OrbView-3\|каталог данных OrbView-3]].
	Однако, в исходном виде, снимки не имеют географической привязки - в формат распространения TIFF не внедрены специальные теги для описания привязки, т.е. формат не представляет из себя классического GeoTIFF. В комплекте поставки имеются все необходимые метаданные для осуществления привязки и ортокоррекции снимков. Как это осуществить будет рассказано далее.		Однако, в исходном виде, снимки не имеют географической привязки - в формат распространения (TIFF) не внедрены специальные теги для описания привязки, т.е. формат не представляет из себя классического GeoTIFF. Однако, в комплекте поставки имеются все необходимые метаданные для осуществления привязки и ортокоррекции снимков. Как это осуществить будет рассказано далее.

	== Программное обеспечение и данные ==		== Программное обеспечение и данные ==

	Для выполнения ортокоррекции космического снимка ~~необходимы следующие файлы~~:		Для выполнения ортокоррекции космического снимка необходимо:
	# Свободно распространяемое программное обеспечение [http://gdal.org/gdalwarp.html gdalwarp], [http://gdal.org/gdalbuildvrt.html gdalbuildvrt], [http://gdal.org/gdal_translate.html gdal_translate] из состава библиотеки GDAL.		# Свободно распространяемое программное обеспечение [http://gdal.org/gdalwarp.html gdalwarp], [http://gdal.org/gdalbuildvrt.html gdalbuildvrt], [http://gdal.org/gdal_translate.html gdal_translate] из состава библиотеки GDAL. Для получения программного обеспечения можно перейти по [http://trac.osgeo.org/gdal/wiki/DownloadingGdalBinaries ссылке], где перечислены варианты загрузки.
	~~# Космический снимок с КА OrbView-3~~		# Космический снимок с КА OrbView-3. Снимки на необходимую территорию ищем в [[Каталог данных OrbView-3\|каталоге]].
	~~# Рельеф на территорию съемки в формате [[DEM]] (Digital elevation model)~~		# Рельеф на территорию съемки в формате в виде ЦМР (Цифровой модели рельефа)

	Для получения программного обеспечения можно перейти по [http://trac.osgeo.org/gdal/wiki/DownloadingGdalBinaries ссылке], где перечислены варианты загрузки.

	Снимки на необходимую территорию ищем в [[Каталог данных OrbView-3\|каталоге]].

	В качестве рельефа могут быть использованы: [http://gis-lab.info/qa/srtm.html SRTM], [http://gis-lab.info/qa/aster-gdem.html ASTER GDEM] v.1 и v.2 и др.		В качестве рельефа могут быть использованы: [http://gis-lab.info/qa/srtm.html SRTM], [http://gis-lab.info/qa/aster-gdem.html ASTER GDEM] v.1 и v.2 и др.
	Рельеф загружаем ~~от сюда~~:		Рельеф загружаем отсюда:
	* [http://gis-lab.info/data/srtm-tif SRTM]		* [http://gis-lab.info/data/srtm-tif SRTM]
	* [http://www.gdem.aster.ersdac.or.jp/index.jsp ASTER GDEM]		* [http://www.gdem.aster.ersdac.or.jp/index.jsp ASTER GDEM]

	Для демонстрации операций по ~~ортокоррекци~~ можно воспользоваться подготовленными наборами данных (снимки и рельеф) по [http://gis-lab.info/data/orbview-3/samples/kursk/ Курской области] и [http://gis-lab.info/data/orbview-3/samples/rb/ Республике Беларусь].		Для демонстрации операций по ортокоррекции можно воспользоваться подготовленными наборами данных (снимки и рельеф) по районам [http://gis-lab.info/data/orbview-3/samples/kursk/ Курской области] и [http://gis-lab.info/data/orbview-3/samples/rb/ Республике Беларусь].

	Набор данных по Курской области включает в себя:		Набор данных по Курской области включает в себя: