Первичная обработка данных SRTM в ГИС SAGA - История изменений

Darsvid: /* Фильтрация */

2013-12-10T13:52:12Z

Фильтрация

← Предыдущая версия		Версия от 15:52, 10 декабря 2013
Строка 166:		Строка 166:
	# <tt>[http://www.webalice.it/alper78/saga_mod/saga_mod.html Grid – Filter (Perego 2009)]</tt><br/>		# <tt>[http://www.webalice.it/alper78/saga_mod/saga_mod.html Grid – Filter (Perego 2009)]</tt><br/>
	# <tt>[http://www.isprs.org/proceedings/XXXIII/congress/part3/935_XXXIII-part3.pdf Grid – Filter => DTM-Filter (slope-based)]</tt>, [http://ssrebelious.blogspot.com/2013/03/dsm-to-dem-conversion.html пример]<br/>		# <tt>[http://www.isprs.org/proceedings/XXXIII/congress/part3/935_XXXIII-part3.pdf Grid – Filter => DTM-Filter (slope-based)]</tt>, [http://ssrebelious.blogspot.com/2013/03/dsm-to-dem-conversion.html пример]<br/>
	# <tt>[http://switch.dl.sourceforge.net/project/saga-gis/SAGA%20-%20Documentation/GGA115/gga115_09.pdf Grid-Filter => Multi Direction Lee Filter]</tt></br>		# <tt>[http://switch.dl.sourceforge.net/project/saga-gis/SAGA%20-%20Documentation/GGA115/gga115_09.pdf Grid-Filter => Multi Direction Lee Filter]</tt><br/>
	# <tt>[http://www.cs.cf.ac.uk/meshfiltering/ Grid-Filter => Mesh Denoise]</tt>, [http://personalpages.manchester.ac.uk/staff/neil.mitchell/mdenoise/ пример]</ref>		# <tt>[http://www.cs.cf.ac.uk/meshfiltering/ Grid-Filter => Mesh Denoise]</tt>, [http://personalpages.manchester.ac.uk/staff/neil.mitchell/mdenoise/ пример]</ref>

Darsvid: /* Фильтрация */ добавила про Mesh Denoise

2013-12-10T13:50:59Z

Фильтрация: добавила про Mesh Denoise

← Предыдущая версия		Версия от 15:50, 10 декабря 2013
Строка 166:		Строка 166:
	# <tt>[http://www.webalice.it/alper78/saga_mod/saga_mod.html Grid – Filter (Perego 2009)]</tt><br/>		# <tt>[http://www.webalice.it/alper78/saga_mod/saga_mod.html Grid – Filter (Perego 2009)]</tt><br/>
	# <tt>[http://www.isprs.org/proceedings/XXXIII/congress/part3/935_XXXIII-part3.pdf Grid – Filter => DTM-Filter (slope-based)]</tt>, [http://ssrebelious.blogspot.com/2013/03/dsm-to-dem-conversion.html пример]<br/>		# <tt>[http://www.isprs.org/proceedings/XXXIII/congress/part3/935_XXXIII-part3.pdf Grid – Filter => DTM-Filter (slope-based)]</tt>, [http://ssrebelious.blogspot.com/2013/03/dsm-to-dem-conversion.html пример]<br/>
	# <tt>[http://switch.dl.sourceforge.net/project/saga-gis/SAGA%20-%20Documentation/GGA115/gga115_09.pdf Grid-Filter => Multi Direction Lee Filter]</tt></ref>		# <tt>[http://switch.dl.sourceforge.net/project/saga-gis/SAGA%20-%20Documentation/GGA115/gga115_09.pdf Grid-Filter => Multi Direction Lee Filter]</tt></br>
			# <tt>[http://www.cs.cf.ac.uk/meshfiltering/ Grid-Filter => Mesh Denoise]</tt>, [http://personalpages.manchester.ac.uk/staff/neil.mitchell/mdenoise/ пример]</ref>

	== Гидрологическая коррекция ==		== Гидрологическая коррекция ==

Vveco: /* Обрезка фрагмента */

2013-10-21T14:05:34Z

Обрезка фрагмента

← Предыдущая версия		Версия от 16:05, 21 октября 2013
Строка 53:		Строка 53:
	[[Файл:saga_dem_prep_05.png\|thumb\|700px\|center]]		[[Файл:saga_dem_prep_05.png\|thumb\|700px\|center]]

	После активации модуля в окне сообщений появится уведомление <span style="color:green">''Interactive module execution has been started''</span>. Для начала ввода координат инструментом [[Файл:saga_georef_action.png]] – <tt>Action</tt> щелкните в любой точке карты. В появившемся окне можно ввести координаты прямоугольника, охватывающего область интереса, – в нашем случае введем координаты углов листа топокарты в десятичных градусах.		После активации модуля в окне сообщений появится уведомление <span style="color:green">''Interactive module execution has been started''</span>. Для начала ввода координат инструментом [[Файл:saga_georef_action.png]] – <tt>Action</tt> щелкните в любой точке карты. В появившемся окне можно ввести координаты прямоугольника, охватывающего область интереса, – в нашем случае введем координаты углов листа топокарты в десятичных градусах.

	<span style="color:red">'''NB'''</span> По мере ввода значения координат автоматически корректируются программой в соответствии с разрешением (рядами и колонками) растра.		<span style="color:red">'''NB'''</span> По мере ввода значения координат автоматически корректируются программой в соответствии с разрешением (рядами и колонками) растра.
			Также нужный участок можно выделить, щелкнув и протянув по нужной области мышью.

	Нажав <tt>Okay</tt>, вы увидите, что на вкладке <tt>Data</tt> появился новый растр – обратите внимание, как отличаются его координаты, число рядов и колонок от исходного.		Нажав <tt>Okay</tt>, вы увидите, что на вкладке <tt>Data</tt> появился новый растр – обратите внимание, как отличаются его координаты, число рядов и колонок от исходного.

Александр Мурый в 11:59, 24 июня 2013

2013-06-24T11:59:27Z

← Предыдущая версия		Версия от 13:59, 24 июня 2013
Строка 1:		Строка 1:
	{{Статья\|~~Черновик~~}}		{{Статья\|Опубликована\|saga-srtm-preprocessing}}
	{{Аннотация\|Возможная последовательность шагов по подготовке данных SRTM к тематическому анализу в ГИС SAGA.}}		{{Аннотация\|Возможная последовательность шагов по подготовке данных SRTM к тематическому анализу в ГИС SAGA.}}

Александр Мурый в 11:16, 24 июня 2013

2013-06-24T11:16:29Z

Александр Мурый: /* Перепроецирование и ресамплинг */

2013-06-24T11:03:03Z

Перепроецирование и ресамплинг

Darsvid: /* Фильтрация */

2013-06-23T16:25:01Z

Фильтрация

← Предыдущая версия		Версия от 18:25, 23 июня 2013
Строка 121:		Строка 121:
	[[Файл:saga_dem_prep_18.png\|thumb\|700px\|center]]		[[Файл:saga_dem_prep_18.png\|thumb\|700px\|center]]

	Рисунок изолиний хорошо передает общие черты рельефа, даже в сравнении с топографической картой, но он содержит много мелких неровностей. ~~Эти неровности~~ как раз и ~~есть той~~ шумовой компонентой, речь о которой шла ранее. ~~Соответственно, прежде~~ чем проводить анализ ЦМВ, шум необходимо устранить с помощью фильтрации. Из-за сложности характера распределения полностью избавиться от шума не удастся, поэтому главная задача фильтрации – максимальная элиминация шумовой компоненты без утраты характерных черт рельефа местности.		Рисунок изолиний хорошо передает общие черты рельефа, даже в сравнении с топографической картой, но он содержит много мелких неровностей. Они как раз и являются шумовой компонентой, речь о которой шла ранее. Прежде чем проводить анализ ЦМВ, шум необходимо устранить с помощью фильтрации. Из-за сложности характера распределения полностью избавиться от шума не удастся, поэтому главная задача фильтрации – максимальная элиминация шумовой компоненты без утраты характерных черт рельефа местности.

	Воспользуемся простым однородным фильтром <tt>Grid – Filter => Simple Filter</tt>. Суть работы фильтра состоит в следующем: он получает новые значения ячеек растра в соответствии с некоторым математическим выражением, т.е. пересчитывает значения центральной ячейки на основе значений ее соседей. Результат фильтрации зависит от параметров скользящего окна, представленных в группе <tt>Options</tt> пунктами <tt>Search Mode</tt>, <tt>Filter</tt> и <tt>Radius</tt>. <tt>Search Mode</tt> и <tt>Radius</tt> совместно контролируют число соседних ячеек растра, которые будут учитываться при расчете нового значения центральной ячейки.		Воспользуемся простым однородным фильтром <tt>Grid – Filter => Simple Filter</tt>. Суть работы фильтра состоит в следующем: он получает новые значения ячеек растра в соответствии с некоторым математическим выражением, т.е. пересчитывает значения центральной ячейки на основе значений ее соседей. Результат фильтрации зависит от параметров скользящего окна, представленных в группе <tt>Options</tt> пунктами <tt>Search Mode</tt>, <tt>Filter</tt> и <tt>Radius</tt>. <tt>Search Mode</tt> и <tt>Radius</tt> совместно контролируют число соседних ячеек растра, которые будут учитываться при расчете нового значения центральной ячейки.

Darsvid в 15:35, 18 июня 2013

2013-06-18T15:35:30Z

← Предыдущая версия		Версия от 17:35, 18 июня 2013
Строка 1:		Строка 1:
	{{Статья\|Черновик}}		{{Статья\|Черновик}}
	{{Аннотация\|~~Одна из возможных последовательностей~~ шагов по подготовке данных SRTM к тематическому анализу}}		{{Аннотация\|Возможная последовательность шагов по подготовке данных SRTM к тематическому анализу}}

	Глобальная цифровая модель высот Shuttle Radar Topography Mission (далее – ЦМВ SRTM), находящаяся в открытом доступе с 2003 года, – общедоступный набор геоданных, активно применяющийся в прикладных исследованиях различной направленности. Ее популярность обуславливается простотой получения, практически глобальным охватом и соответствием запросам среднемасштабного картографирования. По разным оценкам<ref>соотношение между точностью данных SRTM и различными картографическими масштабами детально рассмотрено в публикациях<br/>		Глобальная цифровая модель высот Shuttle Radar Topography Mission (далее – ЦМВ SRTM), находящаяся в открытом доступе с 2003 года, – общедоступный набор геоданных, активно применяющийся в прикладных исследованиях различной направленности. Ее популярность обуславливается простотой получения, практически глобальным охватом и соответствием запросам среднемасштабного картографирования. По разным оценкам<ref>соотношение между точностью данных SRTM и различными картографическими масштабами детально рассмотрено в публикациях<br/>

Darsvid в 15:35, 18 июня 2013

2013-06-18T15:35:02Z

← Предыдущая версия		Версия от 17:35, 18 июня 2013
Строка 1:		Строка 1:
	{{Статья\|Черновик}}		{{Статья\|Черновик}}
	{{Аннотация\|~~Последовательность~~ шагов по подготовке данных SRTM к анализу}}		{{Аннотация\|Одна из возможных последовательностей шагов по подготовке данных SRTM к тематическому анализу}}

	Глобальная цифровая модель высот Shuttle Radar Topography Mission (далее – ЦМВ SRTM), находящаяся в открытом доступе с 2003 года, – общедоступный набор геоданных, активно применяющийся в прикладных исследованиях различной направленности. Ее популярность обуславливается простотой получения, практически глобальным охватом и соответствием запросам среднемасштабного картографирования. По разным оценкам<ref>соотношение между точностью данных SRTM и различными картографическими масштабами детально рассмотрено в публикациях<br/>		Глобальная цифровая модель высот Shuttle Radar Topography Mission (далее – ЦМВ SRTM), находящаяся в открытом доступе с 2003 года, – общедоступный набор геоданных, активно применяющийся в прикладных исследованиях различной направленности. Ее популярность обуславливается простотой получения, практически глобальным охватом и соответствием запросам среднемасштабного картографирования. По разным оценкам<ref>соотношение между точностью данных SRTM и различными картографическими масштабами детально рассмотрено в публикациях<br/>

Darsvid в 15:33, 18 июня 2013

2013-06-18T15:33:15Z

← Предыдущая версия		Версия от 17:33, 18 июня 2013
Строка 7:		Строка 7:
	# Karwel, A., Ewiak, I. [http://www.isprs.org/proceedings/XXXVII/congress/7_pdf/2_WG-VII-2/19.pdf Estimation of the accuracy of the SRTM terrain model on the area of Poland] // The International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences. – Vol. XXXVII, Part B7. – Beijing, 2008. – p. 169-172 <br />		# Karwel, A., Ewiak, I. [http://www.isprs.org/proceedings/XXXVII/congress/7_pdf/2_WG-VII-2/19.pdf Estimation of the accuracy of the SRTM terrain model on the area of Poland] // The International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences. – Vol. XXXVII, Part B7. – Beijing, 2008. – p. 169-172 <br />
	# Ozah, A.P., Kufoniyi, O. [http://www.isprs.org/proceedings/XXXVII/congress/7_pdf/7_WG-VII-7/09.pdf Accuracy assessment of contour interpolation from 1:50 000 topographical maps and SRTM data for 1:25 000 topographical mapping] // The International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences. – Vol. XXXVII, Part B7. – Beijing, 2008 – p. 1347-1353		# Ozah, A.P., Kufoniyi, O. [http://www.isprs.org/proceedings/XXXVII/congress/7_pdf/7_WG-VII-7/09.pdf Accuracy assessment of contour interpolation from 1:50 000 topographical maps and SRTM data for 1:25 000 topographical mapping] // The International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences. – Vol. XXXVII, Part B7. – Beijing, 2008 – p. 1347-1353
	</ref> детальность рельефа, представленного данными SRTM, колеблется в диапазоне масштабов от 50 000 до 100 000~~, удовлетворяя потребности среднемасштабного картографирования~~.		</ref> детальность рельефа, представленного данными SRTM, колеблется в диапазоне масштабов от 50 000 до 100 000.

	Четвертое поколение данных SRTM<ref>Jarvis A., Reuter H., Nelson A., Guevara E. Hole-filled seamless SRTM data V.4. International Centre for Tropical Agriculture (CIAT). – 2008. – http://srtm.csi.cgiar.org</ref> прошло несколько стадий обработки, позволивших повысить исходное качество. Основной целью этих улучшений было заполнение пробелов, характерных для территорий со сложным рельефом, поверхностей занятых водными объектами и прочих типов местностей, плохо поддающихся радарной съемке (например, пустынь). Для этого применялись несколько интерполяционных алгоритмов, а в роли вспомогательных источников использовались локальные и национальные ЦМР более высокого разрешения.		Четвертое поколение данных SRTM<ref>Jarvis A., Reuter H., Nelson A., Guevara E. Hole-filled seamless SRTM data V.4. International Centre for Tropical Agriculture (CIAT). – 2008. – http://srtm.csi.cgiar.org</ref> прошло несколько стадий обработки, позволивших повысить исходное качество. Основной целью этих улучшений было заполнение пробелов, характерных для территорий со сложным рельефом, поверхностей занятых водными объектами и прочих типов местностей, плохо поддающихся радарной съемке (например, пустынь). Для этого применялись несколько интерполяционных алгоритмов, а в роли вспомогательных источников использовались локальные и национальные ЦМР более высокого разрешения.

← Предыдущая версия		Версия от 13:16, 24 июня 2013
Строка 1:		Строка 1:
	{{Статья\|Черновик}}		{{Статья\|Черновик}}
	{{Аннотация\|Возможная последовательность шагов по подготовке данных SRTM к тематическому анализу}}		{{Аннотация\|Возможная последовательность шагов по подготовке данных SRTM к тематическому анализу в ГИС SAGA.}}

	Глобальная цифровая модель высот Shuttle Radar Topography Mission (далее – ЦМВ SRTM), находящаяся в открытом доступе с 2003 года, – общедоступный набор геоданных, активно применяющийся в прикладных исследованиях различной направленности. Ее популярность обуславливается простотой получения, практически глобальным охватом и соответствием запросам среднемасштабного картографирования. По разным оценкам<ref>~~соотношение~~ между точностью данных SRTM и различными картографическими масштабами детально рассмотрено в публикациях<br/>		Глобальная цифровая модель высот Shuttle Radar Topography Mission (далее – ЦМВ SRTM), находящаяся в открытом доступе с 2003 года, — общедоступный набор геоданных, активно применяющийся в прикладных исследованиях различной направленности. Ее популярность обуславливается простотой получения, практически глобальным охватом и соответствием запросам среднемасштабного картографирования. По разным оценкам, <ref>Соотношение между точностью данных SRTM и различными картографическими масштабами детально рассмотрено в публикациях:<br/>
	# Карионов Ю.И. [http://www.geoprofi.ru/default.aspx?mode=binary&id=1168 Оценка точности матрицы высот SRTM] // Геопрофи. – 2010, №10, c.48-51<br/>		* Карионов Ю.И. [http://www.geoprofi.ru/default.aspx?mode=binary&id=1168 Оценка точности матрицы высот SRTM] // Геопрофи. – 2010, №10, c.48-51<br/>
	# Jarvis, A. [http://srtm.csi.cgiar.org/PDF/Jarvis4.pdf Practical use of SRTM data in the tropics: Comparisons with digital elevation models generated from cartographic data] / A. Jarvis, J. Rubiano, A. Nelson, A. Farrow and M. Mulligan. – Cali, CO: Centro Internacional de Agricultura Tropical (CIAT), 2004.– 32 p. (Working document no. 198) <br />		* Jarvis, A. [http://srtm.csi.cgiar.org/PDF/Jarvis4.pdf Practical use of SRTM data in the tropics: Comparisons with digital elevation models generated from cartographic data] / A. Jarvis, J. Rubiano, A. Nelson, A. Farrow and M. Mulligan. – Cali, CO: Centro Internacional de Agricultura Tropical (CIAT), 2004.– 32 p. (Working document no. 198) <br />
	# Karwel, A., Ewiak, I. [http://www.isprs.org/proceedings/XXXVII/congress/7_pdf/2_WG-VII-2/19.pdf Estimation of the accuracy of the SRTM terrain model on the area of Poland] // The International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences. – Vol. XXXVII, Part B7. – Beijing, 2008. – p. 169-172 <br />		* Karwel, A., Ewiak, I. [http://www.isprs.org/proceedings/XXXVII/congress/7_pdf/2_WG-VII-2/19.pdf Estimation of the accuracy of the SRTM terrain model on the area of Poland] // The International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences. – Vol. XXXVII, Part B7. – Beijing, 2008. – p. 169-172 <br />
	# Ozah, A.P., Kufoniyi, O. [http://www.isprs.org/proceedings/XXXVII/congress/7_pdf/7_WG-VII-7/09.pdf Accuracy assessment of contour interpolation from 1:50 000 topographical maps and SRTM data for 1:25 000 topographical mapping] // The International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences. – Vol. XXXVII, Part B7. – Beijing, 2008 – p. 1347-1353		* Ozah, A.P., Kufoniyi, O. [http://www.isprs.org/proceedings/XXXVII/congress/7_pdf/7_WG-VII-7/09.pdf Accuracy assessment of contour interpolation from 1:50 000 topographical maps and SRTM data for 1:25 000 topographical mapping] // The International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences. – Vol. XXXVII, Part B7. – Beijing, 2008 – p. 1347-1353
	</ref> детальность рельефа, представленного данными SRTM, ~~колеблется~~ в ~~диапазоне~~ масштабов от 50 ~~000 до 100~~ 000.		</ref> детальность рельефа, представленного данными SRTM, в целом соответствует таковой на топографических картах масштабов 1:100 000 — 1:50 000.

	Четвертое поколение данных SRTM<ref>Jarvis A., Reuter H., Nelson A., Guevara E. Hole-filled seamless SRTM data V.4. International Centre for Tropical Agriculture (CIAT). – 2008. – http://srtm.csi.cgiar.org</ref> прошло несколько стадий обработки, позволивших повысить исходное качество. Основной целью этих улучшений было заполнение пробелов, характерных для территорий со сложным рельефом, поверхностей занятых водными объектами и прочих типов местностей, плохо поддающихся радарной съемке (например, пустынь). Для этого применялись несколько интерполяционных алгоритмов, а в роли вспомогательных источников использовались локальные и национальные ЦМР более высокого разрешения.		Четвертое поколение данных SRTM<ref>Jarvis A., Reuter H., Nelson A., Guevara E. Hole-filled seamless SRTM data V.4. International Centre for Tropical Agriculture (CIAT). – 2008. – http://srtm.csi.cgiar.org</ref> прошло несколько стадий обработки, позволивших повысить исходное качество. Основной целью этих улучшений было заполнение пробелов, характерных для территорий со сложным рельефом, поверхностей занятых водными объектами и прочих типов местностей, плохо поддающихся радарной съемке (например, пустынь). Для этого применялись несколько интерполяционных алгоритмов, а в роли вспомогательных источников использовались локальные и национальные ЦМР более высокого разрешения.
Строка 162:		Строка 162:
	Поскольку в данном случае нам удалось получить удовлетворительный для учебных целей результат с первого раза, дальнейшая фильтрация проводиться не будет. В более сложных случаях допускается многоразовое использование фильтраций, когда в качестве исходного на каждом последующем этапе используется слой, полученный на предыдущем.		Поскольку в данном случае нам удалось получить удовлетворительный для учебных целей результат с первого раза, дальнейшая фильтрация проводиться не будет. В более сложных случаях допускается многоразовое использование фильтраций, когда в качестве исходного на каждом последующем этапе используется слой, полученный на предыдущем.

	Кроме того, мы использовали самый простой вариант фильтра, но библиотеки <tt>SAGA Grid – Filter</tt> и <tt>Grid – Filter (Perego 2009)</tt> содержат 20 различных модулей фильтров, которые могут использоваться как по отдельности, так и совместно в зависимости от особенностей распределения шумовой компоненты, типов артефактов, характера рельефа и т.д. <ref> ~~подробнее~~ о работе некоторых фильтров<br/>		Кроме того, мы использовали самый простой вариант фильтра, но библиотеки <tt>SAGA Grid – Filter</tt> и <tt>Grid – Filter (Perego 2009)</tt> содержат 20 различных модулей фильтров, которые могут использоваться как по отдельности, так и совместно в зависимости от особенностей распределения шумовой компоненты, типов артефактов, характера рельефа и т.д. <ref>Подробнее о работе некоторых фильтров:<br/>
	# <tt>[http://www.webalice.it/alper78/saga_mod/saga_mod.html Grid – Filter (Perego 2009)]</tt><br/>		# <tt>[http://www.webalice.it/alper78/saga_mod/saga_mod.html Grid – Filter (Perego 2009)]</tt><br/>
	# <tt>[http://www.isprs.org/proceedings/XXXIII/congress/part3/935_XXXIII-part3.pdf Grid – Filter => DTM-Filter (slope-based)]</tt>, [http://ssrebelious.blogspot.com/2013/03/dsm-to-dem-conversion.html пример]<br/>		# <tt>[http://www.isprs.org/proceedings/XXXIII/congress/part3/935_XXXIII-part3.pdf Grid – Filter => DTM-Filter (slope-based)]</tt>, [http://ssrebelious.blogspot.com/2013/03/dsm-to-dem-conversion.html пример]<br/>
Строка 170:		Строка 170:

	Модули гидрологической коррекции объединены библиотекой <tt>Terrain Analysis – Preprocessing</tt>, которая содержит различные инструменты для заполнения локальных впадин (sinks) для облегчения дальнейшего гидрологического анализа.		Модули гидрологической коррекции объединены библиотекой <tt>Terrain Analysis – Preprocessing</tt>, которая содержит различные инструменты для заполнения локальных впадин (sinks) для облегчения дальнейшего гидрологического анализа.
	Воспользуемся модулем <tt>Terrain Analysis – Preprocessing => Fill Sinks (Planchon/ Darboux, 2001)</tt><ref>~~ознакомиться~~ с использованным алгоритмом гидрологической коррекции можно в оригинальной публикации Planchon O., Darboux F. [https://www.rocq.inria.fr/estime/DYNAS/PDF/planchon01b.pdf A fast, simple and versatile algorithm to fill the depressions of digital elevation models] // Catena, 2002, № 46(2-3), р. 159-176</ref>. Принцип его действия состоит в следующем: вместо постепенного заполнения локальных понижений, вся поверхность сначала «заливается» слоем воды, а после ее излишек удаляется, оставляя после себя заполненные понижения. Гибкости алгоритму дает возможность заполнять впадины как до горизонтальной поверхности, так и с сохранением незначительного уклона (например, 0,01°). Первый вариант удобен в том случае, если нужно оценить объем впадины, второй – для оконтуривания дренажной сети. Если в дальнейшем планируется проводить гидрологический анализ ЦМР, то следует выбрать второй вариант.		Воспользуемся модулем <tt>Terrain Analysis – Preprocessing => Fill Sinks (Planchon/ Darboux, 2001)</tt><ref>Ознакомиться с использованным алгоритмом гидрологической коррекции можно в оригинальной публикации Planchon O., Darboux F. [https://www.rocq.inria.fr/estime/DYNAS/PDF/planchon01b.pdf A fast, simple and versatile algorithm to fill the depressions of digital elevation models] // Catena, 2002, № 46(2-3), р. 159-176</ref>. Принцип его действия состоит в следующем: вместо постепенного заполнения локальных понижений, вся поверхность сначала «заливается» слоем воды, а после ее излишек удаляется, оставляя после себя заполненные понижения. Гибкости алгоритму дает возможность заполнять впадины как до горизонтальной поверхности, так и с сохранением незначительного уклона (например, 0,01°). Первый вариант удобен в том случае, если нужно оценить объем впадины, второй – для оконтуривания дренажной сети. Если в дальнейшем планируется проводить гидрологический анализ ЦМР, то следует выбрать второй вариант.

	[[Файл:saga_dem_prep_23.png\|thumb\|700px\|center]]		[[Файл:saga_dem_prep_23.png\|thumb\|700px\|center]]

← Предыдущая версия		Версия от 13:03, 24 июня 2013
Строка 75:		Строка 75:
	[[Файл:saga_dem_prep_09.png\|600px\|center]]		[[Файл:saga_dem_prep_09.png\|600px\|center]]

	== Перепроецирование и ~~ресамплинг~~ ==		== Перепроецирование и ресэмплинг ==

	Данные SRTM распространяются в географической СК на основе эллипсоида WGS 84, поэтому для дальнейшего анализа их необходимо перевести в ~~спрецированную~~ СК. Для этого воспользуемся уже знакомым по привязке топокарт модулем <tt>Projection – Proj.4 => Proj.4 (Dialog, Grid)</tt>. В его диалоговом окне сначала введем параметры исходной проекции <tt>Source Projection Parameters</tt>: выберем географическую СК и датум WGS 84, оставив без изменений прочее.		Данные SRTM распространяются в географической СК на основе эллипсоида WGS-84, поэтому для дальнейшего анализа их необходимо перевести в спроецированную СК. Для этого воспользуемся уже знакомым по привязке топокарт модулем <tt>Projection – Proj.4 => Proj.4 (Dialog, Grid)</tt>. В его диалоговом окне сначала введем параметры исходной проекции <tt>Source Projection Parameters</tt>: выберем географическую СК и датум WGS-84, оставив без изменений прочее.

	[[Файл:saga_dem_prep_10.png\|thumb\|500px\|center]]		[[Файл:saga_dem_prep_10.png\|thumb\|500px\|center]]
Строка 89:		Строка 89:
	[[Файл:saga_dem_prep_12.png\|thumb\|500px\|center]]		[[Файл:saga_dem_prep_12.png\|thumb\|500px\|center]]

	После нажатия <tt>Okay</tt>, мы вернемся в основной диалог, где в качестве метода ~~передискретиации~~ значений выберем билинейную интерполяцию. Данный способ хорошо подходит для континуальных данных (~~каковыми~~ и ~~есть ЦМВ~~), поскольку определяет новое значение ячейки на основе средневзвешенного расстояния от центров 4-х ближайших исходных ячеек, что в свою очередь приводит к незначительному сглаживанию данных.		После нажатия <tt>Okay</tt>, мы вернемся в основной диалог, где в качестве метода передискретации значений выберем билинейную интерполяцию. Данный способ хорошо подходит для континуальных данных (в каковым относится и поле высот), поскольку определяет новое значение ячейки на основе средневзвешенного расстояния от центров 4-х ближайших исходных ячеек, что, в свою очередь, приводит к незначительному сглаживанию данных.

	[[Файл:saga_dem_prep_13.png\|thumb\|500px\|center]]		[[Файл:saga_dem_prep_13.png\|thumb\|500px\|center]]

	В результате перед вам появится окно в котором нужно будет ввести номер зоны UTM, а после – окно с параметрами растра в новой СК (охват, пространственное разрешение, соответствующее количество ~~рядов~~ и ~~колонок~~). При этом по умолчанию предлагается размер ячейки ≈90 м, как это и заявлено для данных ~~ЦМВ~~ SRTM. Но для дальнейшего анализа такое пространственное разрешение не очень удобно, поэтому воспользуемся диалогом для его изменения.		В результате перед вам появится окно, в котором нужно будет ввести номер зоны UTM, а после — окно с параметрами растра в новой СК (охват, пространственное разрешение, соответствующее количество строк и столбцов матрицы). При этом по умолчанию предлагается размер ячейки ≈90 м, как это и заявлено для данных SRTM. Но для дальнейшего анализа такое пространственное разрешение не очень удобно, поэтому воспользуемся диалогом для его изменения.

	Один из простых способов определиться с размером ячейки рассмотрен в статье [http://champs.cecs.ucf.edu/Library/Journal_Articles/pdfs/Hengl_Finding%20the%20right%20pixel%20size.pdf Hengl, 2006]. Согласно предложенному правилу, размер ячейки должен быть равен 0,5 мм аналоговой карты в масштабе исследования. Т.е., если в качестве рабочего масштаба мы выберем масштаб топокарты 1:100 000 (кроме того, именно этому масштабу согласно большинству выводов соответствуют данные SRTM), размер ячейки растра составит 50 м. Обратите внимание, что при вводе числа автоматически пересчитываются и другие значения.		Один из простых способов определиться с размером ячейки рассмотрен в статье [http://champs.cecs.ucf.edu/Library/Journal_Articles/pdfs/Hengl_Finding%20the%20right%20pixel%20size.pdf Hengl, 2006]. Согласно предложенному правилу, размер ячейки должен быть равен 0,5 мм аналоговой карты в масштабе исследования. Т.е., если в качестве рабочего масштаба мы выберем масштаб топокарты 1:100 000 (кроме того, именно этому масштабу согласно большинству выводов соответствуют данные SRTM), размер ячейки растра составит 50 м. Обратите внимание, что при вводе числа автоматически пересчитываются и другие значения.
Строка 105:		Строка 105:
	После сообщения <span style="color:green">''Module execution succeeded''</span> на вкладке <tt>Data</tt> появится новый элемент, который через контекстное меню слоя <tt>Save As…</tt> следует сохранить в рабочую папку проекта под именем ''srtm_m-37-121_utm.sgrd''. Исходный растр ''srtm_m-37-121_gcs'' теперь можно закрыть через контекстное меню слоя <tt>Close</tt>.		После сообщения <span style="color:green">''Module execution succeeded''</span> на вкладке <tt>Data</tt> появится новый элемент, который через контекстное меню слоя <tt>Save As…</tt> следует сохранить в рабочую папку проекта под именем ''srtm_m-37-121_utm.sgrd''. Исходный растр ''srtm_m-37-121_gcs'' теперь можно закрыть через контекстное меню слоя <tt>Close</tt>.

	Для проверки результатов, загрузите рабочий лист топографической карты М-37-121 в проекции UTM и двойным щелчком откройте его в новом окне, в качестве параметра цветового отображения в свойствах объекта установите <tt>Type: RGB</tt>. В ~~эту~~ же ~~карту~~ откройте слой ''srtm_m-37-121_utm''. Если все сделано верно, то благодаря тому, что лист топокарты и растр ЦМВ имеют общую проекцию и отвечают одной и той же территории, слои наложатся. На вкладке свойств объекта (справа) в разделе <tt>Display</tt> установите значение <tt>Transparency [%]: 50</tt>, нажмите <tt>Apply</tt> – это сделает слой ЦМВ наполовину прозрачным и вы сможете лучше оценить взаимное соответствие топокарты и данных SRTM.		Для проверки результатов, загрузите рабочий лист топографической карты М-37-121 в проекции UTM и двойным щелчком откройте его в новом окне, в качестве параметра цветового отображения в свойствах объекта установите <tt>Type: RGB</tt>. В этой же карте откройте слой ''srtm_m-37-121_utm''. Если все сделано верно, то благодаря тому, что лист топокарты и растр ЦМВ имеют общую проекцию и отвечают одной и той же территории, слои наложатся. На вкладке свойств объекта (справа) в разделе <tt>Display</tt> установите значение <tt>Transparency [%]: 50</tt>, нажмите <tt>Apply</tt> – это сделает слой ЦМВ наполовину прозрачным и вы сможете лучше оценить взаимное соответствие топокарты и данных SRTM.

	[[Файл:saga_dem_prep_16.png\|thumb\|600px\|center]]		[[Файл:saga_dem_prep_16.png\|thumb\|600px\|center]]