1. Pascal / Говнокод #10822

    +90

    1. 01
    2. 02
    3. 03
    4. 04
    5. 05
    6. 06
    7. 07
    8. 08
    9. 09
    10. 10
    11. 11
    12. 12
    13. 13
    14. 14
    15. 15
    16. 16
    17. 17
    18. 18
    19. 19
    20. 20
    21. 21
    22. 22
    23. 23
    24. 24
    25. 25
    26. 26
    27. 27
    28. 28
    29. 29
    30. 30
    31. 31
    32. 32
    33. 33
    34. 34
    35. 35
    36. 36
    37. 37
    38. 38
    39. 39
    40. 40
    41. 41
    42. 42
    43. 43
    44. 44
    45. 45
    46. 46
    47. 47
    48. 48
    49. 49
    50. 50
    51. 51
    52. 52
    53. 53
    54. 54
    55. 55
    56. 56
    57. 57
    58. 58
    59. 59
    60. 60
    61. 61
    62. 62
    63. 63
    64. 64
    65. 65
    66. 66
    67. 67
    function Vincenty(Lat1, Lon1, Lat2, Lon2: Extended): Extended;
const  // Параметры эллипсоида:
a = 6378245.0;
f = 1 / 298.3;
b = (1 - f) * a;
EPS = 0.5E-30;
var
APARAM, BPARAM, CPARAM, OMEGA, TanU1, TanU2,
Lambda, LambdaPrev, SinL, CosL, USQR, U1, U2,
SinU1, CosU1, SinU2, CosU2, SinSQSigma, CosSigma,
TanSigma, Sigma, SinAlpha, Cos2SigmaM, DSigma : Extended;
begin
lon1 := lon1 * (PI / 180); 
lat1 := lat1 * (PI / 180);
lon2 := lon2 * (PI / 180); 
lat2 := lat2 * (PI / 180); //Пересчет значений координат в радианы

TanU1 := (1 - f) * Tan(lat1); 
TanU2 := (1 - f) * Tan(lat2);
U1 := ArcTan(TanU1);  
U2 := ArcTan(TanU2);
SinCos(U1, SinU1, CosU1); 
SinCos(U2, SinU2, CosU2);
OMEGA := lon2 - lon1; 
lambda := OMEGA;

repeat //Начало цикла итерации

LambdaPrev:= lambda;
SinCos(lambda, SinL, CosL);
SinSQSigma := (CosU2 * SinL * CosU2 * SinL) +
  (CosU1 * SinU2 - SinU1 * CosU2 * CosL) *
  (CosU1 * SinU2 - SinU1 * CosU2 * CosL);

CosSigma := SinU1 * SinU2 + CosU1 * CosU2 * CosL;
TanSigma:= Sqrt(SinSQSigma) / CosSigma;

if TanSigma > 0  then Sigma := ArcTan(TanSigma)
  else Sigma := ArcTan(TanSigma) + Pi;

if SinSQSigma = 0  then SinAlpha := 0
   else SinAlpha := CosU1 * CosU2 * SinL / Sqrt(SinSQSigma);

if (Cos(ArcSin(SinAlpha)) * Cos(ArcSin(SinAlpha))) = 0  then Cos2SigmaM := 0
  else Cos2SigmaM:= CosSigma - (2 * SinU1 * SinU2 / (Cos(ArcSin(SinAlpha)) * Cos(ArcSin(SinAlpha))));

CPARAM:= (f / 16) * Cos(ArcSin(SinAlpha)) * Cos(ArcSin(SinAlpha)) *
  (4 + f * (4 - 3 * Cos(ArcSin(SinAlpha)) * Cos(ArcSin(SinAlpha))));

lambda := OMEGA + (1 - CPARAM) * f * SinAlpha * (ArcCos(CosSigma) +
  CPARAM * Sin(ArcCos(CosSigma)) * (Cos2SigmaM + CPARAM * CosSigma *
  (-1 + 2 * Cos2SigmaM * Cos2SigmaM)));

until Abs(lambda - LambdaPrev) < EPS; // Конец цикла итерации

USQR:= Cos(ArcSin(SinAlpha)) * Cos(ArcSin(SinAlpha)) *(a * a - b * b) / (b * b);
APARAM := 1 + (USQR / 16384) * (4096 + USQR * (-768 + USQR * (320 - 175 * USQR)));

BPARAM := (USQR / 1024) * (256 + USQR * (-128 + USQR * (74 - 47 * USQR)));
DSigma := BPARAM * SQRT(SinSQSigma) * (Cos2SigmaM + BPARAM / 4 *
(CosSigma * (-1 + 2 * Cos2SigmaM * Cos2SigmaM) - BPARAM / 6 * Cos2SigmaM *
(-3 + 4 * SinSQSigma) * (-3 + 4 * Cos2SigmaM * Cos2SigmaM)));

Result := b * APARAM * (Sigma - DSigma);
end; 

{ ©Drkb::04255 }

    Алгоритм расчёта километража между двумя точками на земной поверхности методом Винсенти, найден в drkb3.0. Там же весь этот ГК уместился в одной строчке:

    distance=sqrt(pow((lon1 - lon2)*111*COS(lat2/57.295781), 2) + pow((lat1) - lat)*111, 2));

    , чудноо... :)

    Запостил: Alegun, 09 Июня 2012

    Комментарии (17) RSS

    • ava guest 09.06.2012 07:35 # +3

      В одной строчке, кстати то ещё говнищще, да к тому же дающе неправильный результат.
      Ответить

    • ava roman-kashitsyn 09.06.2012 07:53 # +6

      > //Начало цикла итерации
      циклы высшего порядка
      Ответить

    • ava Pedofil 09.06.2012 12:18 # −12

      показать все, что скрытоПЕНИС!!!!!!!!!!!!!
      Ответить

    • ava TarasB 09.06.2012 17:22 # 0

      Я не вкурил, найти растояние между точками - это одна теорема косинусов, что за пиздец тут наворотили?
      Ответить

      • ava bormand 09.06.2012 17:47 # −2

        Тут же не сфера...
        Ответить

        • ava TarasB 11.06.2012 13:26 # +2

          А что тут? Неужели эту 1% эллисоидальность кто-то учитывает? Тогда в натуре без интеграла не решить.
          Ответить

          • ava 3.14159265 11.06.2012 14:52 # +8

            2 строчка же.
            >const // Параметры эллипсоида:

            Ну подумаешь самолет пролетев 1000км ошибется на 1% и приземлится за пару километров от аэродрома.
            Что тут такого?
            Ответить

            • ava TarasB 11.06.2012 15:21 # +1

              Да ничего, просто есть факторы и поважнее эллипсоидности, и их хрен учтёшь.
              Ответить

              • ava geust 11.06.2012 15:30 # +3

                > просто есть факторы и поважнее
                Какие? Лично Сталин что-ли?
                http://tinyurl.com/crs54oa
                Ответить

                • ava bormand 11.06.2012 15:34 # +2

                  Рельеф местности например.
                  Ответить

                  • ava geust 11.06.2012 15:39 # +2

                    http://instruct.uwo.ca/earth-sci/505/acadimages/geoid.jpg
                    Ответить

                    • ava TarasB 11.06.2012 15:41 # +1

                      Вот-вот, и поправки на геоид, сдаётся мне, окажутся не менее важными, чем поправки на эллипсность.
                      Ответить

      • ava Lure Of Chaos 09.06.2012 18:25 # +4

        о, если бы земля покоилась на трех слонах, плывущих на большой черепахе в безграничном океане вселенной...
        Ответить

        • ava TarasB 11.06.2012 15:22 # 0

          http://ru.wikipedia.org/wiki/Теоремы_косинусов_(сферическая_геометрия )
          Ответить

    Добавить комментарий