Расчет площадей поверхностей с учетом рельефа средствами GRASS - История изменений

KolesovDmitry: /* Расчет площадей */

2012-03-19T05:18:17Z

Расчет площадей

← Предыдущая версия		Версия от 07:18, 19 марта 2012
Строка 115:		Строка 115:
	* Estimated Region Surface Area: итоговая расчетная площадь территории.		* Estimated Region Surface Area: итоговая расчетная площадь территории.

	[[Файл:~~Ambiguous_triangilation~~.png‎\|thumb\|На рисунке показаны два различных варианта построения треугольников при расчете площади поверхности (синими точками обозначены узлы решетки растра -- 9 ячеек)]]		[[Файл:Ambiguous_triangilation_1.png‎\|thumb\|На рисунке показаны два различных варианта построения треугольников при расчете площади поверхности (синими точками обозначены узлы решетки растра -- 9 ячеек)]]

	[[Файл:Ambiguous triangilation3D.png\|thumb\|На рисунке показаны два различных варианта построения треугольников при расчете площади поверхности (ячейки растра представлены вертикальными столбиками -- 4 ячейки)]]		[[Файл:Ambiguous triangilation3D.png\|thumb\|На рисунке показаны два различных варианта построения треугольников при расчете площади поверхности (ячейки растра представлены вертикальными столбиками -- 4 ячейки)]]

Максим Дубинин: переименовал Расчет площадей поверхностей с учетом рельефа средствами GRASS ГИС в [[Расчет площадей поверхностей с учетом рельефа средс...

2012-03-19T04:54:34Z

переименовал Расчет площадей поверхностей с учетом рельефа средствами GRASS ГИС в [[Расчет площадей поверхностей с учетом рельефа средс...

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

Версия от 06:54, 19 марта 2012

Максим Дубинин в 04:46, 19 марта 2012

2012-03-19T04:46:09Z

← Предыдущая версия		Версия от 06:46, 19 марта 2012
Строка 1:		Строка 1:
	{{Статья\|~~Черновик~~}}		{{Статья\|Опубликована\|3darea-grass}}
	{{Аннотация\|Время от времени возникает задача рассчитать площадь поверхности с учетом ее наклона, т.е. того, что эта поверхность находится в трехмерном пространстве. Эта статья рассматривает пример того, как эту задачу можно решить в GRASS GIS.}}		{{Аннотация\|Время от времени возникает задача рассчитать площадь поверхности с учетом ее наклона, т.е. того, что эта поверхность находится в трехмерном пространстве. Эта статья рассматривает пример того, как эту задачу можно решить в GRASS GIS.}}

KolesovDmitry: /* Расчет площадей */

2012-03-18T10:36:14Z

Расчет площадей

← Предыдущая версия		Версия от 12:36, 18 марта 2012
Строка 115:		Строка 115:
	* Estimated Region Surface Area: итоговая расчетная площадь территории.		* Estimated Region Surface Area: итоговая расчетная площадь территории.

	[[Файл:Ambiguous_triangilation.png‎\|thumb\|На рисунке показаны два различных варианта построения треугольников при расчете площади поверхности (синими точками обозначены узлы решетки растра)]]		[[Файл:Ambiguous_triangilation.png‎\|thumb\|На рисунке показаны два различных варианта построения треугольников при расчете площади поверхности (синими точками обозначены узлы решетки растра -- 9 ячеек)]]

			[[Файл:Ambiguous triangilation3D.png\|thumb\|На рисунке показаны два различных варианта построения треугольников при расчете площади поверхности (ячейки растра представлены вертикальными столбиками -- 4 ячейки)]]


	Таким образом, получаем, что величина 9.027246e+11 -- искомая площадь территории с учетом рельефа, 1.440478e+12 -- общая площадь без учета рельефа (если необходимо произвести сравнение площадей с учетом и без учета рельефа, то из общей площади 1.440478e+12 сначала нужно вычесть число 5.439092e+11 -- площадь пустых ячеек, которые в расчетах не участвовали).		Таким образом, получаем, что величина 9.027246e+11 -- искомая площадь территории с учетом рельефа, 1.440478e+12 -- общая площадь без учета рельефа (если необходимо произвести сравнение площадей с учетом и без учета рельефа, то из общей площади 1.440478e+12 сначала нужно вычесть число 5.439092e+11 -- площадь пустых ячеек, которые в расчетах не участвовали).

Максим Дубинин в 07:54, 18 марта 2012

2012-03-18T07:54:57Z

← Предыдущая версия		Версия от 09:54, 18 марта 2012
Строка 1:		Строка 1:
	{{Статья\|Черновик}}		{{Статья\|Черновик}}
	{{Аннотация\|Время от времени возникает задача рассчитать площадь поверхности с учетом ее наклона, т.е. того, что эта поверхность ~~погружена~~ в ~~трехмерное пространство~~. ~~Данная~~ статья рассматривает пример того, как эту задачу можно ~~реализовать~~ в GRASS GIS.}}		{{Аннотация\|Время от времени возникает задача рассчитать площадь поверхности с учетом ее наклона, т.е. того, что эта поверхность находится в трехмерном пространстве. Эта статья рассматривает пример того, как эту задачу можно решить в GRASS GIS.}}

			Для примера возьмем набор geosample (http://gis-lab.info/qa/geosample.html). Данный набор содержит растр высот (слой relief) и векторную карту охраняемых природных территорий (слой oopt). Требуется посчитать площади каждой из охраняемых территорий как с учетом рельефа, так и без этого.

	Для определенности возьмем набор geosample (http://gis-lab.info/qa/geosample.html). Данный набор содержит растр высот relief и векторную карту охраняемых природных территорий oopt.		== Перепроецирование ==
	~~Требуется посчитать площади каждой из охраняемых территорий как с учетом рельефа, так и без этого.~~

	~~== Перепроецирование ==~~		Поскольку набор geosample содержит данные в географической системе координат (широта/долгота), а измерения требуется произвести в метрической системе координат, то их предварительно нужно перепроецировать, например в равноплощадную проекцию Альберса. Для этого создадим еще один набор geosample_sib_aea в новой проекции со следующими параметрами:
	Поскольку набор geosample содержит данные в системе широта/долгота, то их предварительно нужно перепроецировать, например в равноплощадную проекцию Альберса. Для этого создадим еще один набор geosample_sib_aea в ~~требуемой~~ проекции со следующими параметрами:

	<syntaxhighlight lang="bash">		<syntaxhighlight lang="bash">
Строка 28:		Строка 27:
	</syntaxhighlight>		</syntaxhighlight>

	Посмотрим поближе на исходные данные:		Посмотрим поближе на исходные данные (слой relief):
	<syntaxhighlight lang="bash">		<syntaxhighlight lang="bash">
	> r.info relief		> r.info relief
Строка 60:		Строка 59:
	</syntaxhighlight>		</syntaxhighlight>

	Аналогично по oopt:		Аналогично по слою oopt:
	<syntaxhighlight lang="bash">		<syntaxhighlight lang="bash">
	> v.info oopt		> v.info oopt
Строка 96:		Строка 95:

	== Расчет площадей ==		== Расчет площадей ==
	GRASS содержит модуль r.surf.area (http://grass.osgeo.org/gdp/html_grass64/r.surf.area.html), который производит требуемые расчеты. ~~Расчитаем площади~~ с использованием данного модуля:		GRASS содержит модуль r.surf.area (http://grass.osgeo.org/gdp/html_grass64/r.surf.area.html), который производит требуемые расчеты. Рассчитаем общую площадь всей территории охвата слоя relief с использованием данного модуля:
	<syntaxhighlight lang="bash">		<syntaxhighlight lang="bash">
	> r.surf.area relief		> r.surf.area relief
Строка 109:		Строка 108:
	Done.		Done.
	</syntaxhighlight>		</syntaxhighlight>
	Остановимся подробнее на ~~результате~~:		Остановимся подробнее на результатах расчета:
	* Null value area ignored in calculation: площадь ячеек растра, которые не ~~учавствовали~~ в расчетах (были замаскированы).		* Null value area ignored in calculation: площадь ячеек растра, которые не участвовали в расчетах (были замаскированы).
	* Plan area used in calculation: площадь всего растра (площадь проекции растра на горизонтальную плоскость).		* Plan area used in calculation: площадь всего растра как площадь проекции растра на горизонтальную плоскость.
	* Surface Area Calculation(low, high, avg): площадь растра с учетом рельефа ~~рассчитывается~~ путем построения треугольников, вершины которых находятся в центре ячеек растра. Поскольку для соседних ячеек можно построить несколько вариантов триангуляции (см. пример на рисунке), то и результаты расчета площадей могут различаться. Поэтому алгоритм сохраняет минимальное, максимальное и среднее значение расчетной площади для каждой пары треугольников.		* Surface Area Calculation(low, high, avg): площадь растра с учетом рельефа. Рассчитывается путем построения треугольников, вершины которых находятся в центре ячеек растра. Поскольку для соседних ячеек можно построить несколько вариантов триангуляции (см. пример на рисунке), то и результаты расчета площадей могут различаться. Поэтому алгоритм сохраняет минимальное, максимальное и среднее ((минимальная + максимальная)/2) значение расчетной площади для каждой пары треугольников.
	* Current Region plan area: площадь текущей области (без учета рельефа)		* Current Region plan area: площадь текущей области (без учета рельефа).
	* Estimated Region Surface Area: итоговая расчетная площадь территории.		* Estimated Region Surface Area: итоговая расчетная площадь территории.

Строка 121:		Строка 120:

	Нужно заметить, что результаты расчетов модуля зависят от текущего разрешения (что вполне естественно, т.к. результат и должен зависеть от точности исходных данных).		Нужно заметить, что результаты расчетов модуля зависят от текущего разрешения (что вполне естественно, т.к. результат и должен зависеть от точности исходных данных).

	Посмотрим на текущее разрешение:		Посмотрим на текущее разрешение:
	<syntaxhighlight lang="bash">		<syntaxhighlight lang="bash">
Строка 265:		Строка 265:
	Done.		Done.
	</syntaxhighlight>		</syntaxhighlight>


	Таким же образом можно получить информацию и по остальным участкам.		Таким же образом можно получить информацию и по остальным участкам.
Строка 278:		Строка 277:
	done		done
	</syntaxhighlight>		</syntaxhighlight>


	Начальные строки файла area.txt представлены ниже:		Начальные строки файла area.txt представлены ниже:

KolesovDmitry: /* Расчет площадей */

2012-03-03T08:28:52Z

Расчет площадей

← Предыдущая версия		Версия от 10:28, 3 марта 2012
Строка 112:		Строка 112:
	* Null value area ignored in calculation: площадь ячеек растра, которые не учавствовали в расчетах (были замаскированы).		* Null value area ignored in calculation: площадь ячеек растра, которые не учавствовали в расчетах (были замаскированы).
	* Plan area used in calculation: площадь всего растра (площадь проекции растра на горизонтальную плоскость).		* Plan area used in calculation: площадь всего растра (площадь проекции растра на горизонтальную плоскость).
	* Surface Area Calculation(low, high, avg): площадь растра с учетом рельефа рассчитывается путем построения треугольников, вершины которых находятся в центре ячеек растра. Поскольку для соседних ячеек можно построить несколько вариантов триангуляции, то и результаты расчета площадей могут различаться. Поэтому алгоритм сохраняет минимальное, максимальное и среднее значение расчетной площади для каждой пары треугольников.		* Surface Area Calculation(low, high, avg): площадь растра с учетом рельефа рассчитывается путем построения треугольников, вершины которых находятся в центре ячеек растра. Поскольку для соседних ячеек можно построить несколько вариантов триангуляции (см. пример на рисунке), то и результаты расчета площадей могут различаться. Поэтому алгоритм сохраняет минимальное, максимальное и среднее значение расчетной площади для каждой пары треугольников.
	* Current Region plan area: площадь текущей области (без учета рельефа)		* Current Region plan area: площадь текущей области (без учета рельефа)
	* Estimated Region Surface Area: итоговая расчетная площадь территории.		* Estimated Region Surface Area: итоговая расчетная площадь территории.

KolesovDmitry: /* Расчет площадей */

2012-03-03T08:28:17Z

Расчет площадей

← Предыдущая версия		Версия от 10:28, 3 марта 2012
Строка 115:		Строка 115:
	* Current Region plan area: площадь текущей области (без учета рельефа)		* Current Region plan area: площадь текущей области (без учета рельефа)
	* Estimated Region Surface Area: итоговая расчетная площадь территории.		* Estimated Region Surface Area: итоговая расчетная площадь территории.

			[[Файл:Ambiguous_triangilation.png‎\|thumb\|На рисунке показаны два различных варианта построения треугольников при расчете площади поверхности (синими точками обозначены узлы решетки растра)]]

	Таким образом, получаем, что величина 9.027246e+11 -- искомая площадь территории с учетом рельефа, 1.440478e+12 -- общая площадь без учета рельефа (если необходимо произвести сравнение площадей с учетом и без учета рельефа, то из общей площади 1.440478e+12 сначала нужно вычесть число 5.439092e+11 -- площадь пустых ячеек, которые в расчетах не участвовали).		Таким образом, получаем, что величина 9.027246e+11 -- искомая площадь территории с учетом рельефа, 1.440478e+12 -- общая площадь без учета рельефа (если необходимо произвести сравнение площадей с учетом и без учета рельефа, то из общей площади 1.440478e+12 сначала нужно вычесть число 5.439092e+11 -- площадь пустых ячеек, которые в расчетах не участвовали).

KolesovDmitry: /* Расчет площадей */

2012-03-01T14:57:26Z

Расчет площадей

@@ Строка 96: / Строка 96: @@
 == Расчет площадей ==
-GRASS содержит модуль r.surf.area (http://grass.osgeo.org/gdp/html_grass64/r.surf.area.html), который производит требуемые расчеты.
+GRASS содержит модуль r.surf.area (http://grass.osgeo.org/gdp/html_grass64/r.surf.area.html), который производит требуемые расчеты. Расчитаем площади с использованием данного модуля:
 Нужно заметить, что результаты расчетов модуля зависят от текущего разрешения (что вполне естественно, т.к. результат и должен зависеть от точности исходных данных).
 Посмотрим на текущее разрешение:
@@ Строка 115: / Строка 136: @@
 cells:      3314597
 </syntaxhighlight>
 Огрубим разрешение:

Aleksandr Dezhin: + шаблон аннотации

2012-01-26T18:47:45Z

+ шаблон аннотации

← Предыдущая версия		Версия от 20:47, 26 января 2012
Строка 1:		Строка 1:
	{{Статья\|Черновик}}		{{Статья\|Черновик}}
	Время от времени возникает задача рассчитать площадь поверхности с учетом ее наклона, т.е. того, что эта поверхность погружена в трехмерное пространство. Данная статья рассматривает пример того, как эту задачу можно реализовать в GRASS GIS.		{{Аннотация\|Время от времени возникает задача рассчитать площадь поверхности с учетом ее наклона, т.е. того, что эта поверхность погружена в трехмерное пространство. Данная статья рассматривает пример того, как эту задачу можно реализовать в GRASS GIS.}}

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

2012-01-26T18:41:40Z

← Предыдущая версия		Версия от 20:41, 26 января 2012
Строка 1:		Строка 1:
			{{Статья\|Черновик}}
	Время от времени возникает задача рассчитать площадь поверхности с учетом ее наклона, т.е. того, что эта поверхность погружена в трехмерное пространство. Данная статья рассматривает пример того, как эту задачу можно реализовать в GRASS GIS.		Время от времени возникает задача рассчитать площадь поверхности с учетом ее наклона, т.е. того, что эта поверхность погружена в трехмерное пространство. Данная статья рассматривает пример того, как эту задачу можно реализовать в GRASS GIS.