Форумы GIS-Lab.info

KotAlex писал(а):желательно что-нибудь подходящее для Delphi / C++ Builder / MS Visual Studio.

gamm писал(а): версия 2015 MS Visual Studio нормально работает с Питоном (в том числе бесплатная). Язык, конечно мутноватый (хотя и не сильно мутнее Матлаба/R), но, поскольку нужны только вызовы подпрограмм, продраться можно, да и примеров много

Игорь Белов писал(а): Там есть всё, чего можно желать.

Ariki писал(а): ...

 geod -le

geod +ellps=sphere ...

geod +a=6371000 ...

Ariki писал(а):Не совсем. У меня выводит такие параметры для sphere:
sphere a=6370997.0 b=6370997.0 Normal Sphere (r=6370997)
Лучше задать +a и +b явно.

geod +a=6370997.0 ...

geod +ellps=sphere ...

  MERIT a=6378137.0      rf=298.257       MERIT 1983
    SGS85 a=6378136.0      rf=298.257       Soviet Geodetic System 85
    GRS80 a=6378137.0      rf=298.257222101 GRS 1980(IUGG, 1980)
    IAU76 a=6378140.0      rf=298.257       IAU 1976
     airy a=6377563.396    b=6356256.910    Airy 1830
   APL4.9 a=6378137.0.     rf=298.25        Appl. Physics. 1965
    NWL9D a=6378145.0.     rf=298.25        Naval Weapons Lab., 1965
 mod_airy a=6377340.189    b=6356034.446    Modified Airy
   andrae a=6377104.43     rf=300.0         Andrae 1876 (Den., Iclnd.)
  aust_SA a=6378160.0      rf=298.25        Australian Natl & S. Amer. 1969
    GRS67 a=6378160.0      rf=298.2471674270 GRS 67(IUGG 1967)
   bessel a=6377397.155    rf=299.1528128   Bessel 1841
 bess_nam a=6377483.865    rf=299.1528128   Bessel 1841 (Namibia)
   clrk66 a=6378206.4      b=6356583.8      Clarke 1866
   clrk80 a=6378249.145    rf=293.4663      Clarke 1880 mod.
      CPM a=6375738.7      rf=334.29        Comm. des Poids et Mesures 1799
   delmbr a=6376428.       rf=311.5         Delambre 1810 (Belgium)
  engelis a=6378136.05     rf=298.2566      Engelis 1985
  evrst30 a=6377276.345    rf=300.8017      Everest 1830
  evrst48 a=6377304.063    rf=300.8017      Everest 1948
  evrst56 a=6377301.243    rf=300.8017      Everest 1956
  evrst69 a=6377295.664    rf=300.8017      Everest 1969
  evrstSS a=6377298.556    rf=300.8017      Everest (Sabah & Sarawak)
  fschr60 a=6378166.       rf=298.3         Fischer (Mercury Datum) 1960
 fschr60m a=6378155.       rf=298.3         Modified Fischer 1960
  fschr68 a=6378150.       rf=298.3         Fischer 1968
  helmert a=6378200.       rf=298.3         Helmert 1906
    hough a=6378270.0      rf=297.          Hough
     intl a=6378388.0      rf=297.          International 1909 (Hayford)
    krass a=6378245.0      rf=298.3         Krassovsky, 1942
    kaula a=6378163.       rf=298.24        Kaula 1961
    lerch a=6378139.       rf=298.257       Lerch 1979
    mprts a=6397300.       rf=191.          Maupertius 1738
 new_intl a=6378157.5      b=6356772.2      New International 1967
  plessis a=6376523.       b=6355863.       Plessis 1817 (France)
   SEasia a=6378155.0      b=6356773.3205   Southeast Asia
  walbeck a=6376896.0      b=6355834.8467   Walbeck
    WGS60 a=6378165.0      rf=298.3         WGS 60
    WGS66 a=6378145.0      rf=298.25        WGS 66
    WGS72 a=6378135.0      rf=298.26        WGS 72
    WGS84 a=6378137.0      rf=298.257223563 WGS 84
   sphere a=6370997.0      b=6370997.0      Normal Sphere (r=6370997)

Ariki писал(а): А вот сфер всяких в разных программных продуктах используется куча

KotAlex писал(а): то получается, что шар земной вообще резиновый ...
М - да - а - а, без поллитры не разберешь ...

KotAlex писал(а):Камрады, прошу помощи.
На картах мелких масштабов ( карты стран, регионов и мира ) решаются задачи типа прямой и обратной геодезических, и ведутся разного рода геометрические построения. С целью упрощения расчетов все построения ведутся по методам сферической тригонометрии.
Интересует вопрос - насколько велики будут ошибки и искажения вследствие применения методов сферической тригонометрии для эллиптической поверхности ? имеются ли какие-либо численные выражения подобных ошибок ?
Для максимальной определенности максимально конкретизирую задачу : строится линия кратчайшего расстояния ( ортодрома ) между двумя точками земной поверхности, разнесенными друг от друга на расстояние, скажем, в четверть земной дуги. Насколько велико тут Читы будет отклонение этой линии, построенной на сфере, от этой же линии, построенной на эллипсоиде ?
Или другая задача - для заданной точки найти координаты точки, отстоящей от первой на заданных азимуте и дистанции ( дистанция, как и в первом случае - четверть земной окружности ). Насколько велики расхождения вторых точек в случае расчетов по указанным методам ?
Я знаю, что на эллипсоиде эта линия строится методом численного интегрирования дифференциальных
уравнений геодезической линии, и на основе этого можно получить численные оценки расхождений.
Но меня интересует - не имеются ли уже готовые ответы на этот вопрос ?

stout писал(а): Дан эллипсоид. Найти радиус эквивалентной сферы.
Ответ зависит от того, что понимать под эквивалентной сферой. Это сфера той же площади или того же объёма?
Или сфера с радиусом, равным среднему радиусу кривизны эллипсоида в точке?
Разные условия дают разный ответ. Насколько помню, на практике различают 5 разных радиусов. Для каждой задачи – свой радиус.

Alexa111 писал(а):
Тоже интересует данный вопрос. Как можно с вами связаться? Оставьте пожалуйста какой-нибудь контакт! Не могу послать личное сообщение..

nickleb писал(а):В R есть библиотека geosphere сферической тригонометрии:
https://cran.r-project.org/web/packages ... index.html
https://cran.r-project.org/web/packages ... sphere.pdf