Расчет географических координат в Proj.4

sergeym2008 · 2012-06-28T12:10:50.000Z

Здравствуйте, не знаю по адресу ли я, но подскажите пожалуйста как при помощи библиотеки Proj.4 рассчитать географическую координату по точке (широта, долгота) и расстоянию в км.
Например к 10 гр. с.ш. и 0 гр. в.д прибавить 6600 км по широте и 6000 км по долготе и получить географическую координату.

chnav · 2012-06-28T13:07:05.000Z

По-русски это “прямая геодезическая задача (на эллипсоиде)”, по-английски “direct geodetic problem”, ищите функции со схожими названиями.

sergeym2008 · 2012-08-06T17:14:52.000Z

результат

вызывать GetShiftedCoordDegree нужно обязательно после GeodSetup

// longitude долгота в градусах
// latitude широта в градусах
// x смещение по долготе в метрах
// y смещение по широте в метрах
// result_longitude долгота результат
// result_latitude широта результат
void GetShiftedCoordDegree(double longitude, double latitude, double x, double y, double *result_longitude, double *result_latitude, int *error);

код
//------------------------------------------------------------------------------
#include “projects.h”
#include “geocent.h”

#define DEG2RAD(x) x * PI / 180.
double deg2rad(double degree)
{
while (degree > 180)
degree -= 180;
while (degree < -180)
degree += 180;

return DEG2RAD(degree);
}

#define RAD2DEG(x) x * 180. / PI;
double rad2deg(double radian)
{
while (radian > PI)
radian -= PI;
while (radian < -PI)
radian += PI;

return RAD2DEG(radian);
}

// настраиваем geod, установка проекции, эллипса и т.п.
void GeodSetup(int argc, char *argv)
{
geod_set(argc, argv); / setup projection */
}

void GetShiftedCoordDegree(double longitude, double latitude, double x, double y, double *result_longitude, double *result_latitude, int *error)
{
// process - функция в geod.c
process(deg2rad(latitude), deg2rad(longitude), deg2rad(90), x, result_latitude, result_longitude);
process(*result_latitude, *result_longitude, deg2rad(0), y, result_latitude, result_longitude);
*result_longitude = rad2deg(*result_longitude);
*result_latitude = rad2deg(*result_latitude);
}

//------------------------------------------------------------------------------