Пересчет координат в проекциях Гаусса-Крюгера

Добрый день.
Подскажите пожалуйста как при переходе из плоских координат в географические учитывается Начальная широта, когда она не равна 0

 

Добрый вечер, после всех сообщений про СРУ Ваш вопрос меня обнадежил - есть в этой жизни .
Если начальная широта не равна О, то она считается после 6 лет в НИИГАиК или еще много где это узнать.
Что-то конкретно надо пересчитать, выложите. Пересчитают.

 

Поконкретнее бы о каких СК речь...
Думаю что как точка начала отсчёта этих прямоугольных плоских координат на север. Осевой меридиан, соответственно, на восток (иногда -500 км для положительных значений координат).

 

yeti, http://geodesist.ru/forum/threads/archaeosys.20703/page-9#post-543630

 

спасибо за предложение пересчитать, но я и сам могу пересчитать, мне нужен алгоритм . Я правильно понимаю в в ссылке выше рассказывается про перевод географических в плоские координаты? А на оборот. Вычисляем β затем В_х потом B, а где добавить В₀.

 

Попутно задам еще один вопрос
координаты:

это Колумбия, ближе примера не нашел их датум:

a, м	1 / f	ΔX, м	ΔY, м	ΔZ, м	L0	B0	m	false X, м	false Y, м
6378388.0	297	307.0	304.0	-318.0	-74.0809	4.599047	1.0000	1000000	1000000

перевожу в их градусы
1)если использовать преобразования через пространственные прямоугольные координаты, вместо недостающих параметров используя нули, то результат аналогичен мапинфо и геокалькулятору.
2)если же переводить сразу из геодезических в геодезические (по Молоденскому), то появляется ошибка которая при переходе к местным координатам выливается в 30 метров по х

вот часть того самого алгоритма, толи точность этого расчета такая, толи я что то не то делаю

Код:

da = a_sk - a_wgs
a_sr = (a_wgs + a_sk) / 2									
a_sr = a_sk
de2 = e2_sk - e2_wgs
e2_sr = (e2_wgs + e2_sk) / 2
 
		  N = a_sr * (1 - e2_sr * Sin(B) ^ 2) ^ -0.5
		  M = a_sr * (1 - e2_sr) * (1 - e2_sr * Sin(B) ^ 2) ^ -1.5
 
dB = ro / (M + H) * (N / a_sr * e2_sr * Sin(B) * Cos(B) * da _
+ (N ^ 2 / a_sr ^ 2 + 1) * N * Sin(B) * Cos(B) * de2 / 2 _
- (dx_sk * Cos(L) + dy_sk * Sin(L)) * Sin(B) + dz_sk * Cos(B))
 
dL = ro / ((N + H) * Cos(B)) * (-dx_sk * Sin(L) + dy_sk * Cos(L))
 
dH = -a_sr / N * da + N * Sin(B) ^ 2 * de2 / 2 _
+ (dx_sk * Cos(L) + dy_sk * Sin(L)) * Cos(B) + dz_sk * Sin(B)
 
 
		dB = (dB / 3600) * Pi / 180
		dL = (dL / 3600) * Pi / 180
 
B = B - dB
L = L - dL
H = H - dH

Здесь кстати эта начальная широта присутствует. когда переходишь от градусов к метрам я ее учел. посчитал длину дуги до точки и до начальной широты и потом взял их разность. А вот обратный пересчет не понял.

 

http://www.remotesensing.org/geotiff/proj_list/transverse_mercator.html

 

может, что называется, на пальцах объясните?!? ото я не могу...
это как наш В_х? типа В_х=В_х+В0, но результат не тот

п.с. про 3 и 7 параметров тоже интересно