Форумы GIS-Lab.info

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

static const double DEG_TO_RAD = 0.017453292519943295769236907684886;
static const double EARTH_RADIUS_IN_METERS = 6372797.560856;
double latitudeArc  = (latitude - right.latitude) * DEG_TO_RAD;
double longitudeArc = (longitude - right.longitude) * DEG_TO_RAD;
double latitudeH = sin(latitudeArc * 0.5);
latitudeH *= latitudeH;
double lontitudeH = sin(longitudeArc * 0.5);
lontitudeH *= lontitudeH;
double tmp = cos(latitude * DEG_TO_RAD) * cos(right.latitude * DEG_TO_RAD);
double distanceArc =  2.0 * asin(sqrt(latitudeH + tmp * lontitudeH));
return EARTH_RADIUS_IN_METERS * distanceArc;

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

bool inside = false;
for (int i = 0, j = numPoints - 1; i < numPoints; j = i++)
{
	if (((points[i].y > testPoint.y) != (points[j].y > testPoint.y)) &&
		(testPoint.x < (points[j].x - points[i].x) * (testPoint.y - points[i].y) / (points[j].y - points[i].y) + points[i].x))
	{
		inside = !inside;
	}
}
return inside;