UE 基于Georeference插件 将 百度坐标系转换为 UE内部坐标

Georeference插件支持将真实世界的地理数据与UE中的场景进行对应，例如地形、建筑物、道路等，可以在虚拟环境中实现精确的地理定位和导航。

GeoRefrence插件中的 "GeographicToProjected" 函数可以将地理坐标（经纬度，一般是WGS84坐标系）转换为投影坐标（平面坐标），"ProjectedToEngine"函数用于将已经进行过投影的地理数据转换为引擎的坐标格式。

但是GeoRefrence插件不支持百度坐标系的直接转换，所以要将百度坐标系（BD-09）先转火星坐标系（GCJ-02），火星坐标系（GCJ-02）再转WGS84，具体如下：

• 百度坐标系（BD-09）转为  火星坐标系（GCJ-02）

void BD09toGCJ02(double bd_lon, double bd_lat, double& GCJ_lon, double& GCJ_lat)
{
	double x = bd_lon - 0.0065;
	double y = bd_lat - 0.006;
	double x_PI = 3.14159265358979323846 * 3000.0 / 180.0;
	double z = sqrt(x * x + y * y) - 0.00002 * sin(y * x_PI);
	double theta = atan2(y, x) - 0.000003 * cos(x * x_PI);
	GCJ_lon = z * cos(theta);
	GCJ_lat = z * sin(theta);
}

• 火星坐标系（GCJ-02）转为 WGS84

bool GCJ02ToWGS84(FString lon_Str_GCJ, FString lat_Str_GCJ, FString& LonStr_WGS, FString& LatStr_WGS)
{

	TArray<double> WGS = {1,2};

	double lon_GCJ = FCString::Atod(*lon_Str_GCJ);
	double lat_GCJ = FCString::Atod(*lat_Str_GCJ);

	if (!isInChinaBbox(lon_GCJ, lat_GCJ)) {
		WGS[0] = lon_GCJ;
		WGS[1] = lat_GCJ;
		return false;
	}
	WGS[0] = lon_GCJ; WGS[1] = lat_GCJ;
	TArray<FString> temp = WGS84ToGCJ02(WGS[0], WGS[1]);
	double dx = FCString::Atod(*temp[0]) - lon_GCJ;
	double dy = FCString::Atod(*temp[1]) - lat_GCJ;
	while (fabs(dx) > 1e-6 || fabs(dy) > 1e-6) {
		WGS[0] -= dx;
		WGS[1] -= dy;
		temp = WGS84ToGCJ02(WGS[0], WGS[1]);
		dx = FCString::Atod(*temp[0]) - lon_GCJ;
		dy = FCString::Atod(*temp[1]) - lat_GCJ;
	}

	std::stringstream lonss;// = std::to_string(WGS[0]);
	lonss << std::setprecision(15) << WGS[0];
	std::string lonstr = lonss.str();

	std::stringstream latss;// = std::to_string(WGS[0]);
	latss << std::setprecision(15) << WGS[1];
	std::string lanstr = latss.str();

	LonStr_WGS = lonstr.c_str();
	LatStr_WGS = lanstr.c_str();
	return true;
}


int isInChinaBbox(double lon, double lat) {
	return lon >= 72.004 && lon <= 137.8347 && lat >= 0.8293 && lat <= 55.8271;
}


TArray<FString> WGS84ToGCJ02(double lon_WGS, double lat_WGS) {
	TArray<FString> GCJ ;
	if (!UPSBlueprintFunctionLibrary::isInChinaBbox(lon_WGS, lat_WGS)) {
		GCJ.Insert(FString::SanitizeFloat(lon_WGS), 0);
		GCJ.Insert(FString::SanitizeFloat(lat_WGS), 1);
		return GCJ;
	}
	double* d = delta(lon_WGS, lat_WGS);
	GCJ.Insert(FString::SanitizeFloat(lon_WGS + d[0]), 0);
	GCJ.Insert(FString::SanitizeFloat(lat_WGS + d[1]), 1);
	return GCJ;
}

double* delta(double lon, double lat) {
	double* d = (double*)malloc(2 * sizeof(double));
	double dLon = transformLon(lon - 105.0, lat - 35.0);
	double dLat = transformLat(lon - 105.0, lat - 35.0);

	double radLat = (lat / 180.0) * PI;
	double magic = sin(radLat);

	magic = 1.0 - ee * magic * magic;

	double sqrtMagic = sqrt(magic);
	d[0] = (dLon * 180.0) / (a / sqrtMagic * cos(radLat) * PI);
	d[1] = (dLat * 180.0) / ((a * (1.0 - ee)) / (magic * sqrtMagic) * PI);
	return d;
}

设置原点位置的参数：可以通过超图桌面的“坐标点转换”工具   

例如本项目：118.52803879715°，31.2563464533562°（84坐标系下的 经纬度）的转换  ：

ProjectedCRS 填的参数是：UTM坐标系下场景区域对应的EPSG编号

UTM（通用横轴墨卡托投影）是一种将地球表面划分为60个纵向区域的投影方法，每个区域覆盖6度经度范围，每个纵向区域都用一个数字标识，从1到60。我的项目是UTM Zone 51N表示北纬0度到6度之间的区域，对应的EPSG就是32651。

中国UTM分区可以参考：中国UTM分区_utm分度带103e-CSDN博客

坐标简单介绍：

1. WGS84（World Geodetic System 1984）：WGS84是一种全球通用的地理坐标系统，常用于全球定位系统（GPS）和卫星导航系统。在中国范围内，WGS84坐标可以用于大多数定位和地图应用。

2. GCJ-02（国测局坐标）：GCJ-02是中国国家测绘局（NGCC）使用的一种加密坐标系统，也称为火星坐标系。GCJ-02坐标通常用于中国地图服务、导航和定位应用。这是因为中国政府对地理信息的保密政策，对原始WGS84坐标进行了加密处理。

3. BD-09（百度坐标）：百度坐标是百度地图使用的一种自有坐标系统，主要用于中国地区。BD-09坐标是基于GCJ-02坐标进行了特殊的偏移计算。百度坐标主要用于百度地图及相关应用中。

参考：

1. 百度坐标(BD09)、GCJ02、与WGS84之间的转换 C/C++ UI实现_gcj java转c++-CSDN博客
2. 【精选】UE4 GeoReferencing Plugin-CSDN博客