BJ54坐标或者XA80坐标到WGS84坐标的转换（小工具附代码）

前几天自我总结时候发现了之前给我姐写的一个winform的小程序，用来简单进行BJ54坐标到WGS84坐标的转换。

觉得挺有意思的，就贴出来分享一下。

PS：只用了简单的高斯变换，没有四、七参数的变化，所以转出来还是有较大误差。不过可以判断出大概位置

具体说明如下：

1.打开程序，主界面如下，可以进行BJ54坐标系的XY坐标到WGS84经纬度坐标的单点或者多点转换。

2.输入带号和X,Y坐标，点击单点转换，即可进行坐标的单点转换。

3.新建一个低版本（.xls格式）的excel表，将X，Y坐标输入excel表，保存。

4.在程序里先输入带号，然后点击多点转换，选择对应的excel表点击打开。

5.程序会自动计算坐标，完成会提示另存为，此时选择将计算后新生成表保存即可。

6.找到新生成的表，里面就是对应的经纬度坐标了。

具体的代码如下：

namespace WindowsFormsApplication1
{public partial class Form1 : Form{public Form1(){InitializeComponent();}private void button1_Click(object sender, EventArgs e){double a = double.Parse(textBox1.Text);double b = double.Parse(textBox2.Text);double f = double.Parse(textBox5.Text);double c=0;double d = 0;GaussProjInvCal(a, b, f,ref c, ref d);textBox3.Text = c.ToString();textBox4.Text = d.ToString();}void GaussProjInvCal(double X, double Y, double daihao,ref double longitude, ref double latitude){int ProjNo; 带宽double longitude1, latitude1, longitude0, X0, Y0, xval, yval;double e1, e2, f, a, ee, NN, T, C, M, D, R, u, fai, iPI;iPI = 0.0174532925199433; 3.1415926535898/180.0;a = 6378245.0; f = 1.0 / 298.3; //54年北京坐标系参数a=6378140.0; f=1/298.257; //80年西安坐标系参数ProjNo = (int)(X / 1000000L); //查找带号longitude0 = daihao*3;longitude0 = longitude0 * iPI; //中央经线X0 = ProjNo * 1000000L + 500000L;Y0 = 0;xval = X - X0; yval = Y - Y0; //带内大地坐标e2 = 2 * f - f * f;e1 = (1.0 - Math.Sqrt(1 - e2)) / (1.0 + Math.Sqrt(1 - e2));ee = e2 / (1 - e2);M = yval;u = M / (a * (1 - e2 / 4 - 3 * e2 * e2 / 64 - 5 * e2 * e2 * e2 / 256));fai = u + (3 * e1 / 2 - 27 * e1 * e1 * e1 / 32) * Math.Sin(2 * u) + (21 * e1 * e1 / 16 - 55 * e1 * e1 * e1 * e1 / 32) * Math.Sin(4 * u)+ (151 * e1 * e1 * e1 / 96) * Math.Sin(6 * u) + (1097 * e1 * e1 * e1 * e1 / 512) * Math.Sin(8 * u);C = ee * Math.Cos(fai) * Math.Cos(fai);T = Math.Tan(fai) * Math.Tan(fai);NN = a / Math.Sqrt(1.0 - e2 * Math.Sin(fai) * Math.Sin(fai));R = a * (1 - e2) / Math.Sqrt((1 - e2 * Math.Sin(fai) * Math.Sin(fai)) * (1 - e2 * Math.Sin(fai) * Math.Sin(fai)) * (1 - e2 * Math.Sin(fai) * Math.Sin(fai)));D = xval / NN;//计算经度(Longitude) 纬度(Latitude)longitude1 = longitude0 + (D - (1 + 2 * T + C) * D * D * D / 6 + (5 - 2 * C + 28 * T - 3 * C * C + 8 * ee + 24 * T * T) * D* D * D * D * D / 120) / Math.Cos(fai);latitude1 = fai - (NN * Math.Tan(fai) / R) * (D * D / 2 - (5 + 3 * T + 10 * C - 4 * C * C - 9 * ee) * D * D * D * D / 24+ (61 + 90 * T + 298 * C + 45 * T * T - 256 * ee - 3 * C * C) * D * D * D * D * D * D / 720);//转换为度 DDlongitude = longitude1 / iPI;latitude = latitude1 / iPI;}private void label6_Click(object sender, EventArgs e){}private void button2_Click(object sender, EventArgs e){try{string filepath = "";OpenFileDialog opf = new OpenFileDialog();if (opf.ShowDialog() == DialogResult.OK)filepath = opf.FileName;DataSet ds = new DataSet();ds = ExcelToDS(filepath);DataTable dt = ds.Tables["table1"];foreach (DataRow dr in dt.Rows){double a = double.Parse(dr[0].ToString());double b = double.Parse(dr[1].ToString());double f = double.Parse(textBox5.Text);double c = 0;double d = 0;GaussProjInvCal(a, b, f, ref c, ref d);dr[0] = c.ToString();dr[1] = d.ToString();}DataGridViewExportToExcel(ds, "坐标excel");}catch (Exception ex){MessageBox.Show(ex.Message, "软件提示");}}public DataSet ExcelToDS(string Path){string strConn = "Provider=Microsoft.Jet.OLEDB.4.0;" + "Data Source=" + Path + ";" + "Extended Properties=Excel 8.0;";OleDbConnection conn = new OleDbConnection(strConn);conn.Open();string strExcel = "";OleDbDataAdapter myCommand = null;DataSet ds = null;strExcel = "select * from [sheet1$]";myCommand = new OleDbDataAdapter(strExcel, strConn);ds = new DataSet();myCommand.Fill(ds, "table1");return ds;} public void DataGridViewExportToExcel(DataSet ds, string strTitle){SaveFileDialog saveFileDialog = new SaveFileDialog();saveFileDialog.Filter = "Excel files (*.xls)|*.xls";saveFileDialog.FilterIndex = 0;saveFileDialog.RestoreDirectory = true;saveFileDialog.CreatePrompt = false;saveFileDialog.FileName = strTitle + ".xls";if (saveFileDialog.ShowDialog() == DialogResult.Cancel) //导出时，点击【取消】按钮{return;}Stream myStream = saveFileDialog.OpenFile();StreamWriter sw = new StreamWriter(myStream, 　　　　　　System.Text.Encoding.GetEncoding(-0));string strHeaderText = "";try{//写标题for (int i = 0; i < ds.Tables[0].Columns.Count; i++){if (i > 0){strHeaderText += "\t";}strHeaderText += ds.Tables[0].Columns[i].ToString();}sw.WriteLine(strHeaderText);//写内容string strItemValue = "";for (int j = 0; j < ds.Tables[0].Rows.Count; j++){strItemValue = "";for (int k = 0; k < ds.Tables[0].Columns.Count; k++){if (k > 0){strItemValue += "\t";}strItemValue += ds.Tables[0].Rows[j][k].ToString();}sw.WriteLine(strItemValue); //把dgv的每一行的信息写为sw的每一行}}catch (Exception ex){MessageBox.Show(ex.Message, "软件提示");}finally{sw.Close();myStream.Close();}}}}

BJ54坐标或者XA80坐标到WGS84坐标的转换（小工具附代码）相关推荐

WGS84经纬度坐标转北京54高斯投影地方坐标流程
坐标系统之间的坐标转换既包括不同的参心坐标之间的转换,或者不同的地心坐标系之间的转换,也包括参心坐标系与地心坐标系之间的转换以及相同坐标系的直角坐标与大地坐标之间的坐标转换,还有大地坐标与高斯平面坐标 ...
wgs84坐标格式转换度分秒_使用ArcGIS实现WGS84经纬度坐标到北京54高斯投影坐标的转换...
[摘要] 本文针对从事测绘工作者普遍遇到的坐标转换问题,简要介绍ArcGIS实现WGS84经纬度坐标到北京54高斯投影坐标转换原理和步骤. [关键词] ArcGIS 坐标转换投影变换 1 坐标转换 ...
火星坐标、百度坐标、WGS-84坐标相互转换及墨卡托投影坐标转经纬度JavaScript版...
转自:https://www.cnblogs.com/fwc1994/p/5884115.html 火星坐标.百度坐标.WGS-84坐标相互转换及墨卡托投影坐标转经纬度JavaScript版火星坐标 ...
wgs84坐标格式转换度分秒_一起爬山吗？寻找GIS坐标系统中“隐秘的角落”
随着网剧<隐秘的角落>的爆红,笛卡尔与他的心形线的传说再次重回大众视野,r＝a(1-sinθ) 就是"心形线"方程式. 笛卡尔法国著名数学家,哲学家,创建了数学坐标系, ...
【转】火星坐标、百度坐标、WGS-84坐标相互转换及墨卡托投影坐标转经纬度JavaScript版...
原文地址:https://www.cnblogs.com/fwc1994/p/5884115.html 火星坐标火星坐标是国家测绘局为了国家安全在原始坐标的基础上进行偏移得到的坐标,基本国内的电子地 ...
使用QGIS插件转换火星坐标、百度坐标和WGS84坐标
最近因工作需要,分别从高德和百度下载了一批POI点数据.由于高德地图采用国家测绘地理信息局GCJ02坐标系(即俗称火星坐标系),百度采用自己的BD09坐标系,而国际来源地图大多采用WGS84坐标系,导 ...
gps wgs84坐标与高德gcj02坐标互转
vue中 wgs84坐标转为高德gcj02坐标,详见高德官方文档:其他坐标转高德坐标-坐标变换-教程-地图 JS API | 高德地图API //设置gps wgs84坐标转为高德gcj02conve ...
微信地图wgs84坐标，gcj02坐标，bd09坐标转换
微信小程序三种常见经纬度坐标系的转化遇见问题:在其它端点位显示正常,在小程序上点位发生偏移,微信小程序是gcj02坐标,就是xxToGcj02,具体情况具体分析我们常用的地图api坐标系有wgs8 ...
火星坐标, WGS84坐标, 百度坐标
地图坐标系目前包括: 地球坐标 (WGS84) WGS84:World Geodetic System 1984,是为GPS全球定位系统使用而建立的坐标系统. 国际标准,从 GPS 设备中取出的数据的 ...
WGS84坐标转BD-09坐标（wgs84转百度）
推荐使用方法2 !!! Web服务API | 百度地图API SDK 方法1 定义一个待转换的坐标数组. // 待转换的坐标数组 var points = [new BMap.Point(116.4 ...

BJ54坐标或者XA80坐标到WGS84坐标的转换（小工具附代码）

BJ54坐标或者XA80坐标到WGS84坐标的转换（小工具附代码）相关推荐

最新文章

热门文章