Вычисление области Polygon с помощью Широты и Долготы, сгенерированной от Декартова пространства и мирового файла

Учитывая серию пар координаты GPS, я должен вычислить область полигона (n-полувагон). Это является относительно маленьким (не больше, чем 50 000 кв.футов). Геокоды создаются путем применения аффинного преобразования с данными из мирового файла.

Я попытался использовать два подхода шага путем выполнения преобразования геокодов к декартовым координатам:

double xPos = (lon-lonAnchor)*( Math.toRadians( 6378137 ) )*Math.cos( latAnchor );
double yPos = (lat-latAnchor)*( Math.toRadians( 6378137 ) );

затем я использую вычисление векторного произведения для определения области.

Проблема - то, что результаты немного выключены в точности (приблизительно 1%). Есть ли что-нибудь, что я могу изучить для улучшения этого?

Спасибо.