Форумы GIS-Lab.info

Код: Выделить всё

echo "10dN 0dE 84.85453 6000" | geod +ellps=WGS84 +units=km 

Код: Выделить всё

echo "10dN 54d47'38.321\"E 0 6600" | geod +ellps=WGS84 +units=km

Код: Выделить всё

// 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);	
}