Здравствуйте, не знаю по адресу ли я, но подскажите пожалуйста как при помощи библиотеки Proj.4 рассчитать географическую координату по точке (широта, долгота) и расстоянию в км.
Например к 10 гр. с.ш. и 0 гр. в.д прибавить 6600 км по широте и 6000 км по долготе и получить географическую координату.
По-русски это “прямая геодезическая задача (на эллипсоиде)”, по-английски “direct geodetic problem”, ищите функции со схожими названиями.
результат
вызывать 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);
}
//------------------------------------------------------------------------------