Задачи на сфере: прямая геодезическая задача - История изменений

Игорь Белов: /* Решение прямой задачи средствами PROJ.4 */

2020-05-10T05:46:52Z

Решение прямой задачи средствами PROJ.4

← Предыдущая версия		Версия от 07:46, 10 мая 2020
Строка 251:		Строка 251:
	</syntaxhighlight>		</syntaxhighlight>

	== Решение прямой задачи средствами PROJ.4 ==		== Решение прямой задачи средствами PROJ ==

	В пакет PROJ.4 входит программа '''geod''', предназначенная для решения прямых и обратных геодезических задач на сфере. Так выглядит команда обработки файла '''dir.dat''':		В пакет PROJ входит программа '''geod''', предназначенная для решения прямых и обратных геодезических задач на сфере. Так выглядит команда обработки файла '''dir.dat''':

	<syntaxhighlight lang="bash">		<syntaxhighlight lang="bash">
	$ geod +a=6371000 -f "%f" +units=km dir.dat		$ geod +a=6371000 +b=6371000 -f "%f" +units=km dir.dat
	</syntaxhighlight>		</syntaxhighlight>

ErnieBoyd: /* Альтернативные методы */

2014-06-21T06:16:24Z

Альтернативные методы

← Предыдущая версия		Версия от 08:16, 21 июня 2014
Строка 275:		Строка 275:
	К наиболее надёжным относится следующий способ:		К наиболее надёжным относится следующий способ:

	<math>\begin{~~array}{rcl~~}		<math>\begin{align}
	\xi & = & \cos \sigma \cos \varphi_1 - \sin \sigma \sin \varphi_1 \cos \alpha_1 \\		\xi & = \cos \sigma \cos \varphi_1 - \sin \sigma \sin \varphi_1 \cos \alpha_1 \\
	\eta & = & \sin \sigma \sin \alpha_1 \\		\eta & = \sin \sigma \sin \alpha_1 \\
	\operatorname{tg\,} (\lambda_2 - \lambda_1) & = & \dfrac{\eta}{\xi} \\		\operatorname{tg\,} (\lambda_2 - \lambda_1) & = \dfrac{\eta}{\xi} \\
	\operatorname{tg\,} \varphi_2 & = & \dfrac{\sin \varphi_1 \cos \sigma + \cos \varphi_1 \sin \sigma \cos \alpha_1}{\xi \cos (\lambda_2 - \lambda_1) + \eta \sin (\lambda_2 - \lambda_1)}		\operatorname{tg\,} \varphi_2 & = \dfrac{\sin \varphi_1 \cos \sigma + \cos \varphi_1 \sin \sigma \cos \alpha_1}{\xi \cos (\lambda_2 - \lambda_1) + \eta \sin (\lambda_2 - \lambda_1)}
	\end{~~array~~}</math>		\end{align}</math>

	В сферической тригонометрии углы и стороны должны быть в диапазоне [0, 180°]. Алгоритмизация формул требует анализа и обработки случаев, когда входные величины не попадают в эти рамки.		В сферической тригонометрии углы и стороны должны быть в диапазоне [0, 180°]. Алгоритмизация формул требует анализа и обработки случаев, когда входные величины не попадают в эти рамки.

ErnieBoyd: /* Преобразование декартовых координат в сферические */

2014-06-21T06:15:14Z

Преобразование декартовых координат в сферические

← Предыдущая версия		Версия от 08:15, 21 июня 2014
Строка 148:		Строка 148:
	=== Преобразование декартовых координат в сферические ===		=== Преобразование декартовых координат в сферические ===

	: <math>\begin{~~array}{rcl~~}		: <math>\begin{align}
	\lambda & = & \operatorname{arctg\,} \dfrac{y}{x} \\		\lambda & = \operatorname{arctg\,} \dfrac{y}{x} \\
	\varphi & = & \operatorname{arctg\,} \dfrac{z}{\sqrt{x^2 + y^2}}		\varphi & = \operatorname{arctg\,} \dfrac{z}{\sqrt{x^2 + y^2}}
	\end{~~array~~}</math>		\end{align}</math>

	В данном случае в роли сферических координат ''φ'', ''λ'' окажутся ''φ''₂, ''λ''₂.		В данном случае в роли сферических координат ''φ'', ''λ'' окажутся ''φ''₂, ''λ''₂.

ErnieBoyd: /* Вращение вокруг оси */

2014-06-21T06:14:30Z

Вращение вокруг оси

← Предыдущая версия		Версия от 08:14, 21 июня 2014
Строка 110:		Строка 110:
	Представим оператор вращения вокруг оси ''X'' на угол ''θ'' в следующем виде:		Представим оператор вращения вокруг оси ''X'' на угол ''θ'' в следующем виде:

	: <math>\begin{~~array}{rcl~~}		: <math>\begin{align}
	x' & = & x \\		x' & = x \\
	y' & = & y \cos \theta + z \sin \theta \\		y' & = y \cos \theta + z \sin \theta \\
	z' & = & - y \sin \theta + z \cos \theta		z' & = - y \sin \theta + z \cos \theta
	\end{~~array~~}</math>		\end{align}</math>

	Операторы вращения вокруг осей ''Y'' и ''Z'' получаются перестановкой символов.		Операторы вращения вокруг осей ''Y'' и ''Z'' получаются перестановкой символов.

ErnieBoyd: /* Преобразование сферических координат в декартовы */

2014-06-21T06:13:42Z

Преобразование сферических координат в декартовы

← Предыдущая версия		Версия от 08:13, 21 июня 2014
Строка 73:		Строка 73:
	=== Преобразование сферических координат в декартовы ===		=== Преобразование сферических координат в декартовы ===

	: <math>\begin{~~array}{rcl~~}		: <math>\begin{align}
	x & = & \cos \varphi \cos \lambda \\		x & = \cos \varphi \cos \lambda \\
	y & = & \cos \varphi \sin \lambda \\		y & = \cos \varphi \sin \lambda \\
	z & = & \sin \varphi		z & = \sin \varphi
	\end{~~array~~}</math>		\end{align}</math>

	В данном случае в качестве сферических координат ''φ'', ''λ'' подставим углы (90° − ''σ''), (180° − ''α''₁).		В данном случае в качестве сферических координат ''φ'', ''λ'' подставим углы (90° − ''σ''), (180° − ''α''₁).

ErnieBoyd: /* Преобразование декартовых координат в сферические */

2014-03-23T08:16:33Z

Преобразование декартовых координат в сферические

← Предыдущая версия		Версия от 10:16, 23 марта 2014
Строка 174:		Строка 174:
	double p;		double p;

	p = ~~sqrt~~(x[0] ~~* x[0] + x[1] *~~ x[1]);		p = hypot(x[0], x[1]);
	y[1] = atan2(x[1], x[0]);		y[1] = atan2(x[1], x[0]);
	y[0] = atan2(x[2], p);		y[0] = atan2(x[2], p);

	return ~~sqrt~~(p ~~* p + x[2] *~~ x[2]);		return hypot(p, x[2]);
	}		}
	</syntaxhighlight>		</syntaxhighlight>

ErnieBoyd: /* Ссылки */

2014-03-19T03:50:39Z

Ссылки

← Предыдущая версия		Версия от 05:50, 19 марта 2014
Строка 285:		Строка 285:

	== Ссылки ==		== Ссылки ==
	* [[Вычисление расстояния и начального азимута между двумя точками на сфере]]		* [http://gis-lab.info/qa/great-circles.html Вычисление расстояния и начального азимута между двумя точками на сфере]
	* [[Задачи на сфере: обратная геодезическая задача]]		* [http://gis-lab.info/qa/sphere-geodesic-invert-problem.html Задачи на сфере: обратная геодезическая задача]
	* [[Задачи на сфере: угловая засечка]]		* [http://gis-lab.info/qa/sphere-geodesic-angular-resection.html Задачи на сфере: угловая засечка]
	* [[Задачи на сфере: линейная засечка]]		* [http://gis-lab.info/qa/sphere-geodesic-linear-resection.html Задачи на сфере: линейная засечка]
	* [http://www.pm298.ru/sferich.php Краткий справочник по сферической тригонометрии]		* [http://www.pm298.ru/sferich.php Краткий справочник по сферической тригонометрии]
	* [http://trac.osgeo.org/proj/wiki/man_geod man_geod – PROJ.4]		* [http://trac.osgeo.org/proj/wiki/man_geod man_geod – PROJ.4]
	* [http://en.wikipedia.org/wiki/Earth_radius Earth radius]		* [http://en.wikipedia.org/wiki/Earth_radius Earth radius]
	* [http://gis-lab.info/docs/books/sphere-trigonometry/sphere-trigonometry.rar Степанов Н. Н. Сферическая тригонометрия]		* [http://gis-lab.info/docs/books/sphere-trigonometry/sphere-trigonometry.rar Степанов Н. Н. Сферическая тригонометрия]

ErnieBoyd в 17:41, 18 марта 2014

2014-03-18T17:41:55Z

← Предыдущая версия		Версия от 19:41, 18 марта 2014
Строка 1:		Строка 1:
	{{Статья\|Опубликована\|~~http://gis-lab.info/qa/~~sphere-geodesic-direct-problem}}		{{Статья\|Опубликована\|sphere-geodesic-direct-problem}}

	{{Аннотация\|Прямая геодезическая задача — это нахождение положения точки по координатам исходного пункта и значениям начального направления и расстояния.}}		{{Аннотация\|Прямая геодезическая задача — это нахождение положения точки по координатам исходного пункта и значениям начального направления и расстояния.}}

ErnieBoyd в 17:41, 18 марта 2014

2014-03-18T17:41:25Z

← Предыдущая версия		Версия от 19:41, 18 марта 2014
Строка 1:		Строка 1:
	{{Статья\|~~Черновик~~}}		{{Статья\|Опубликована\|http://gis-lab.info/qa/sphere-geodesic-direct-problem}}

	{{Аннотация\|Прямая геодезическая задача — это нахождение положения точки по координатам исходного пункта и значениям начального направления и расстояния.}}		{{Аннотация\|Прямая геодезическая задача — это нахождение положения точки по координатам исходного пункта и значениям начального направления и расстояния.}}

ErnieBoyd: /* Решение прямой задачи средствами PROJ.4 */

2014-03-18T06:17:02Z

Решение прямой задачи средствами PROJ.4

← Предыдущая версия		Версия от 08:17, 18 марта 2014
Строка 259:		Строка 259:
	</syntaxhighlight>		</syntaxhighlight>

	Параметр <tt>+a</tt> определяет радиус сферы, <tt>-f</tt> — формат вывода угловых величин, <tt>+units</tt> — единица измерения расстояний. В ~~результате~~ получим идентичный результат:		Параметр <tt>+a</tt> определяет радиус сферы, <tt>-f</tt> — формат вывода угловых величин, <tt>+units</tt> — единица измерения расстояний. В итоге получим идентичный результат:

	<pre>		<pre>