C ++ Конвертировать Lat Long в BNG с Proj.4
Я создаю плагин для VBS2. В программе я могу вывести сетку в LAT LONG, UTM или MGRS. Мне нужно иметь возможность конвертировать в BNG. Мне удалось создать приложение TKinter в Python, которое работает с использованием Proj.4, но теперь нужно создать его на C ++ в виде DLL.
Использование LatLong im (51.20650N 1.81906W) — известная точка, конверсия BNG составляет приблизительно SU 127 452.
#include <proj_api.h>double x = 51.20650; //atof(GetStrArgument(input, 0));
double y = 1.81906; //atof(GetStrArgument(input, 1));
char *pj_latlongc = "+proj=longlat +datum=WGS84";
char *pj_UTMc = "+proj=utm +zone=29 +ellps=WGS84";
char *osc = "+proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000 +ellps=airy +datum=OSGB36 +units=m +no_defs";
projPJ pj_OS = pj_init_plus(osc);
projPJ pj_UTM = pj_init_plus(pj_UTMc);
projPJ pj_latlong = pj_init_plus(pj_latlongc);
x *= DEG_TO_RAD;
y *= DEG_TO_RAD;
int p = pj_transform(pj_latlong, pj_OS, 1, 1, &x, &y, NULL);
К сожалению, результат совершенно неверный и не имеет смысла. Может ли кто-нибудь помочь в этом?
Решение
Переключение x для y и переключение знака x, кажется, дают лучшие результаты. Запуск этих чисел на консоли, proj выход, кажется, лучше согласуется с указанная вами точка (51.20650N 1.81906W)
echo -1.81906 51.20650 | proj -V +proj=tmerc
Longitude: 1d49'8.616"W [ -1.81906 ]
Latitude: 51d12'23.4"N [ 51.2065 ]
Так же, cs2cs возвращает, казалось бы, разумный вывод:
echo -1.81906 51.20650 | cs2cs -v +proj=longlat +datum=WGS84 +to +proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000 +ellps=airy +datum=OSGB36 +units=m +no_defs
# ---- From Coordinate System ----
#Lat/long (Geodetic alias)
#
# +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
# ---- To Coordinate System ----
#Transverse Mercator
# Cyl, Sph&Ell
# +proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000
# +ellps=airy +datum=OSGB36 +units=m +no_defs
# +towgs84=446.448,-125.157,542.060,0.1502,0.2470,0.8421,-20.4894
#--- following specified but NOT used
# +ellps=airy
412736.92 145270.05 -47.95
Это идентично выводу программы C ++
double x = -1.81906;
double y = 51.20650;
412736.924409 145270.054358
Для сравнения я проверил цифры через BNG конвертер Британской геологической службы, который, похоже, тоже согласен
Easting: 412737
Northing: 145270
Для написания этих чисел в виде Национальная Сетка Боеприпасов ссылка, просто вычтите целые кратные 100 км, чтобы получить соответствующую букву, как показано в этом Исходный код для JavaScript. Конечный результат SU 127 452, по желанию:
#include <proj_api.h>
#include <iostream>
int main() {
double x = -1.81906;
double y = 51.20650;
const char* pj_latlongc = "+proj=longlat +datum=WGS84";
const char* osc = "+proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000 +ellps=airy +datum=OSGB36 +units=m +no_defs";
projPJ pj_latlong = pj_init_plus(pj_latlongc);
projPJ pj_OS = pj_init_plus(osc);
x *= DEG_TO_RAD;
y *= DEG_TO_RAD;
int p = pj_transform(pj_latlong, pj_OS, 1, 1, &x, &y, NULL);
std::cout.setf(std::ios::fixed, std::ios::floatfield);
std::cout.setf(std::ios::showpoint);
std::cout << x << " " << y << std::endl;;
// prints 412736.924409 145270.054358// now convert to UK Grid Ref
int e1 = floor(x/100000);
int n1 = floor(y/100000);
int e2 = (int)x % 100000 / 100;
int n2 = (int)y % 100000 / 100;
char l1 = (19 - n1) - (19-n1) % 5 + ((e1 + 10)/5);
char l2 = (19 - n1) * 5 % 25 + e1 % 5;
if (l1 > 7) l1++;
if (l2 > 7) l2++;
l1 = 'A' + l1;
l2 = 'A' + l2;
std::cout << l1 << l2 << ' ' << e2 << ' ' << n2 << std::endl;
// prints SU 127 452
}
Другие решения