mozzherinmetodichkaofcartography61

Приложение 5

Rem Программа расчета величин искажений в картографических проекциях

Rem Блок определения переменных

Dim latitude, longitude, delta_latitude, delta_longitude As Integer

Dim small_l1, small_l2, teta As Single

Dim scaling, large_L1, large_L2 As Long

Dim m, n, p, a, b, k, omega, alpha0, ro, beta As Single

Rem Блок ввода исходных данных

Input “Введите широту точки в градусах: ”, latitude Input “Введите долготу точки в градусах: “, longitude

Input “Введите частоту градусной сетки по параллели в градусах: “, delta_latitude Input “Введите частоту градусной сетки по меридиану в градусах: “, delta_longitude Input “Введите расстояние по касательной к меридиану в сантиметрах: “, small_l1 Input “Введите расстояние по касательной к параллели в сантиметрах: “, small_l2 Input “Введите угол между касательными к меридиану и параллели в градусах: “, teta Input “Введите знаменатель главного масштаба: “, scaling

Input “Введите длину дуги меридиана в метрах: “, large_L1 Input “Введите длину дуги параллели в метрах: “, large_L2

Rem Блок расчета величин искажений

m = CInt(small_l1 * scaling / large_L1) / 100 n = CInt(small_l2 * scaling / large_L2) / 100

p = CInt(m * n * Sin(teta * 3.14159 / 180) * 100) / 100

a = CInt((Sqr(m ^ 2 + 2 * p + n ^ 2) + Sqr(m ^ 2 - 2 * p + n ^ 2)) / 2 * 100) / 100 b = CInt((Sqr(m ^ 2 + 2 * p + n ^ 2) - Sqr(m ^ 2 - 2 * p + n ^ 2)) / 2 * 100) / 100 k = CInt(a / b * 100) / 100

omega = (a - b) / (a + b)

omega = CInt(2 * Atn(omega / Sqr(omega ^ 2 + 1)) * 180 / 3.14159 * 10) / 10 If m = b Then

alpha0 = 90 Else

alpha0 = CInt(Atn(b / a * Sqr((a ^ 2 - m ^ 2) / (m ^ 2 - b ^ 2))) * 180 / 3.14159 * 10) / 10 End If

ro = CInt(Sqr((p - 1) ^ 2 + (k - 1) ^ 2) * 100) / 100

beta = CInt(Atn((k - 1) / (p - 1)) * 180 / 3.14159 * 10) / 10

Rem Блок вывода результатов на экран

Print “Величина искажения длин по меридиану, m: “; using “##.##”; m Print “Величина искажения длин по параллели, n: “; using “##.##”; n Print “Величина искажения площади, p: “; using “##.##”; p

Print “Наибольшая величина искажения длин, a: “; using “##.##”; a Print “Наименьшая величина искажения длин, b: “; using “##.##”; b Print “Величина искажения форм, k: “; using “##.##”; k

Print “Величина искажения угла, omega: “; using “##.#”; omega

Print “Азимут направления наибольшего искажения длин, alpha0: “; using “##.#”; alpha0 Print “Вектор ро, ro: “; using “##.##”; ro

Print “Направление вектора ро, beta: “; using “##.#”; beta

End

91

Приложение 6

92

Карта 3

Нормальная цилиндрическая равноугольная проекция Меркатора на секущем по 30 параллелям цилиндре

95

Приложение 8

Rem Программа расчета координат промежуточных точек ортодромий и локсодромий

Rem Блок определения переменных

Dim phiA, phiB, lambdaA, lambdaB, A As Single

Dim alphaort, alphalox As Single

Dim delta_latitude, delta_longitude, lambda, phi As Integer

Rem Блок ввода исходных данных

Input “Введите широту начальной точки в градусах: ”, phiA Input “Введите долготу начальной точки в градусах: “, lambdaA Input “Введите широту конечной точки в градусах: ”, phiB Input “Введите долготу конечной точки в градусах: “, lambdaB

Input “Введите частоту градусной сетки по параллели в градусах: “, delta_latitude Input “Введите частоту градусной сетки по меридиану в градусах: “, delta_longitude

Rem Блок расчета координат промежуточных точек ортодромии

A = (Atn(1 / (Tan(phiB * 3.14159 / 180) / (Tan(phiA * 3.14159 / 180) * sin ((lambdaB - lambdaA)

*3.14159 / 180) - 1 / Tan((lambdaB - lambdaA) * 3.14159 / 180)))) * 180 / 3.14159 Input “Введите первый меридиан, с которого следует начать расчет ортодромии в

градусах: “, lambda

Do While lambda < lambdaB

phi = Atn((Tan(phiA * 3.14159 / 180) * Sin((A – lambdaA + lambda) * 3.14159 / 180)) / Sin(A * 3.14159 / 180)) * 180 / 3.14159

Print “Долгота промежуточной точки ортодромии: “; using “##”; lambda; Print “Широта промежуточной точки ортодромии: “; using “##.#”; phi; lambda = lambda + delta_longitude

Loop

alphaort = Sin(phiA * 3.14159 / 180) / Tan(A * 3.14159 / 180) * 180 / 3.14159 Print “Азимут линии ортодромии в начальной точке “; using “##.#”; alphaort

Rem Блок расчета координат промежуточных точек локсодромии

alphalox = Atn(3.14159 * (lambdaB - lambdaA) / (180 * (Log(45 + phiB / 2) - Log(45 + phiA / 2)))) * 180 / 3.14159

Input “Введите первую параллель, с которой следует начать расчет локсодромии в градусах: “, phi

Do While phi < phiB

lambda = lambdaA + Tan(alphalox * 3.14159 / 180) * (180 * (Log(45 + phi / 2) - Log(45 + phiA / 2)) / 3.14159)

Print “Широта промежуточной точки локсодромии: “; using “##”; phi; Print “Долгота промежуточной точки локсодромии: “; using “##.#”; lambda; lambda = lambda + delta_latitude

Loop

Print “Азимут линии локсодромии “; using “##.#”; alphalox

End

99