通过两点经纬度确认精度

最新推荐文章于 2023-07-29 20:58:41 发布

原创最新推荐文章于 2023-07-29 20:58:41 发布 · 444 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#其他

日常专栏收录该内容

15 篇文章

订阅专栏

vs c++ 下开发的程序：

#include"math.h"
#include"stdio.h"
#include"stdlib.h"

#define PI 3.1415926
#define EARTH_RADIUS 6378.137

double radian(double d)
{
	return d*PI/180.0;
}


int main(int argc, char *argv[])
{
	double longitudeA =atof(argv[1]);
	double latitudeA  =atof(argv[2]);

	double longitudeB =atof(argv[3]);
	double latitudeB  =atof(argv[4]);

	double distance = 0;
	double radLatA = radian(latitudeA);
	double radLatB = radian(latitudeB);
	double a = radLatA - radLatB;
	double b = radian(longitudeA) - radian(longitudeB);
	__int64 tmp =0;


    if(5 != argc)
	{
		printf("parameter count error!\n");
		exit(1);
	}

	printf("longitudeA: %16f\n",longitudeA);
	printf("latitudeA : %16f\n",latitudeA);
	printf("longitudeB: %16f\n",longitudeB);
	printf("latitudeB : %16f\n",latitudeB);



	distance = 2*asin((sqrt(pow(sin(a / 2),2) + cos(radLatA) * cos(radLatB) * pow(sin(b / 2),2) ))); 
	distance = distance * EARTH_RADIUS;

	//distance = round(distance *10000)/10000;
	tmp =  ((distance *10000) +0.5);
	distance = tmp/10000.0;

	printf("distance is %fm\n",distance*1000);

	return 0;
}

2.在ubuntu下开发的程序


$cat distancebylongitudeandlatitude.c 
#include <math.h>
#include<stdio.h>
#include<stdlib.h>

#define PI 3.1415926
#define EARTH_RADIUS 6378.137

double radian(double d)
{
	return d*PI/180.0;
}


int main(int argc, char *argv[])
{
	if(5 != argc)
	{
		printf("parameter count error!\n");
        exit(1);
	}


	double longitudeA =atof(argv[1]);
	double latitudeA  =atof(argv[2]);

	double longitudeB =atof(argv[3]);
	double latitudeB  =atof(argv[4]);

	printf("longitudeA: %16f\n",longitudeA);
	printf("latitudeA : %16f\n",latitudeA);
	printf("longitudeB: %16f\n",longitudeB);
	printf("latitudeB : %16f\n",latitudeB);


	double distance = 0;
	double radLatA = radian(latitudeA);
	double radLatB = radian(latitudeB);
	double a = radLatA - radLatB;
	double b = radian(longitudeA) - radian(longitudeB);

	distance = 2*asin((sqrt(pow(sin(a / 2),2) + cos(radLatA) * cos(radLatB) * pow(sin(b / 2),2) ))); 
	distance = distance * EARTH_RADIUS;
	distance = round(distance *10000)/10000;
	
    printf("distance is %fm\n",distance*1000);

	return 0;
}