实例介绍
【实例简介】
【实例截图】
【核心代码】
using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
using System.Threading.Tasks;
namespace ConsoleApp1
{
public struct CTPoint {
public double lng;
public double lat;
}
public class CoordTransform
{
static double x_PI = 3.14159265358979324 * 3000.0 / 180.0;
static double PI = 3.1415926535897932384626;
static double a = 6378245.0;
static double ee = 0.00669342162296594323;
/**
* 百度坐标系 (BD-09) 与 火星坐标系 (GCJ-02)的转换
*
* @param lng 经度
* @param lat 纬度
*/
public static CTPoint bd09togcj02(double bd_lon, double bd_lat)
{
double x = bd_lon - 0.0065;
double y = bd_lat - 0.006;
double z = Math.Sqrt(x * x y * y) - 0.00002 * Math.Sin(y * x_PI);
double theta = Math.Atan2(y, x) - 0.000003 * Math.Cos(x * x_PI);
double gg_lng = z * Math.Cos(theta);
double gg_lat = z * Math.Sin(theta);
CTPoint cTPoint;
cTPoint.lat = gg_lat;
cTPoint.lng = gg_lng;
return cTPoint;
}
/**
* 火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换
*
* @param lng 经度
* @param lat 纬度
*/
public static CTPoint gcj02tobd09(double lng, double lat)
{
double z = Math.Sqrt(lng * lng lat * lat) 0.00002 * Math.Sin(lat * x_PI);
double theta = Math.Atan2(lat, lng) 0.000003 * Math.Cos(lng * x_PI);
double bd_lng = z * Math.Cos(theta) 0.0065;
double bd_lat = z * Math.Sin(theta) 0.006;
CTPoint cTPoint;
cTPoint.lng = bd_lng;
cTPoint.lat = bd_lat;
return cTPoint;
}
/**
* WGS84转GCj02
*
* @param lng 经度
* @param lat 纬度
* @return
*/
public static CTPoint wgs84togcj02(double lng, double lat)
{
double nLng = lng;
double nLat = lat;
if (!out_of_china(lng, lat))
{
double dlat = transformlat(lng - 105.0, lat - 35.0);
double dlng = transformlng(lng - 105.0, lat - 35.0);
double radlat = lat / 180.0 * PI;
double magic = Math.Sin(radlat);
magic = 1 - ee * magic * magic;
double Sqrtmagic = Math.Sqrt(magic);
dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * Sqrtmagic) * PI);
dlng = (dlng * 180.0) / (a / Sqrtmagic * Math.Cos(radlat) * PI);
nLng = lat dlat;
nLat = lng dlng;
}
CTPoint cTPoint;
cTPoint.lng = nLng;
cTPoint.lat = nLat;
return cTPoint;
}
public static double transformlat(double lng, double lat)
{
double ret = -100.0 2.0 * lng 3.0 * lat 0.2 * lat * lat 0.1 * lng * lat
0.2 * Math.Sqrt(Math.Abs(lng));
ret = (20.0 * Math.Sin(6.0 * lng * PI) 20.0 * Math.Sin(2.0 * lng * PI)) * 2.0 / 3.0;
ret = (20.0 * Math.Sin(lat * PI) 40.0 * Math.Sin(lat / 3.0 * PI)) * 2.0 / 3.0;
ret = (160.0 * Math.Sin(lat / 12.0 * PI) 320 * Math.Sin(lat * PI / 30.0)) * 2.0 / 3.0;
return ret;
}
/**
* GCJ02 转换为 WGS84
*
* @param lng 经度
* @param lat 纬度
* @return
*/
public static CTPoint gcj02towgs84(double lng, double lat)
{
double nLng = lng;
double nLat = lat;
if (!out_of_china(lng, lat))
{
double dlat = transformlat(lng - 105.0, lat - 35.0);
double dlng = transformlng(lng - 105.0, lat - 35.0);
double radlat = lat / 180.0 * PI;
double magic = Math.Sin(radlat);
magic = 1 - ee * magic * magic;
double sqrtmagic = Math.Sqrt(magic);
dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * PI);
dlng = (dlng * 180.0) / (a / sqrtmagic * Math.Cos(radlat) * PI);
double mglat = lat dlat;
double mglng = lng dlng;
nLng = lng * 2 - mglng;
nLat = lat * 2 - mglat;
}
CTPoint cTPoint;
cTPoint.lng = nLng;
cTPoint.lat = nLat;
return cTPoint;
}
public static double transformlng(double lng, double lat)
{
double ret = 300.0 lng 2.0 * lat 0.1 * lng * lng 0.1 * lng * lat 0.1 * Math.Sqrt(Math.Abs(lng));
ret = (20.0 * Math.Sin(6.0 * lng * PI) 20.0 * Math.Sin(2.0 * lng * PI)) * 2.0 / 3.0;
ret = (20.0 * Math.Sin(lng * PI) 40.0 * Math.Sin(lng / 3.0 * PI)) * 2.0 / 3.0;
ret = (150.0 * Math.Sin(lng / 12.0 * PI) 300.0 * Math.Sin(lng / 30.0 * PI)) * 2.0 / 3.0;
return ret;
}
/**
* 判断是否在国内，不在国内则不做偏移
*
* @param lng 经度
* @param lat 纬度
*/
static bool out_of_china(double lng, double lat)
{
// 纬度3.86~53.55,经度73.66~135.05
return !(lng > 73.66 && lng < 135.05 && lat > 3.86 && lat < 53.55);
}
}
}
好例子网口号：伸出你的我的手 — 分享！
小贴士
感谢您为本站写下的评论，您的评论对其它用户来说具有重要的参考价值，所以请认真填写。
- 类似“顶”、“沙发”之类没有营养的文字，对勤劳贡献的楼主来说是令人沮丧的反馈信息。
- 相信您也不想看到一排文字/表情墙，所以请不要反馈意义不大的重复字符，也请尽量不要纯表情的回复。
- 提问之前请再仔细看一遍楼主的说明，或许是您遗漏了。
- 请勿到处挖坑绊人、招贴广告。既占空间让人厌烦，又没人会搭理，于人于己都无利。
关于好例子网
本站旨在为广大IT学习爱好者提供一个非营利性互相学习交流分享平台。本站所有资源都可以被免费获取学习研究。本站资源来自网友分享，对搜索内容的合法性不具有预见性、识别性、控制性，仅供学习研究，请务必在下载后24小时内给予删除，不得用于其他任何用途，否则后果自负。基于互联网的特殊性，平台无法对用户传输的作品、信息、内容的权属或合法性、安全性、合规性、真实性、科学性、完整权、有效性等进行实质审查；无论平台是否已进行审查，用户均应自行承担因其传输的作品、信息、内容而可能或已经产生的侵权或权属纠纷等法律责任。本站所有资源不代表本站的观点或立场，基于网友分享，根据中国法律《信息网络传播权保护条例》第二十二与二十三条之规定，若资源存在侵权或相关问题请联系本站客服人员，点此联系我们。关于更多版权及免责申明参见 版权及免责申明
网友评论
我要评论