高德地图GCJ-02火星坐标系与GPS的WGS-84坐标系转换公式（Java）（可读取excel文件）

最新推荐文章于 2024-08-22 14:23:25 发布

*Major*

最新推荐文章于 2024-08-22 14:23:25 发布

阅读量2.9k

点赞数 3

本文链接：https://blog.csdn.net/qq_41375318/article/details/104399481

版权

高德地图GCJ-02火星坐标系与GPS的WGS-84坐标系转换公式（Java）

注意excel格式需要为: .xls 不能为xlsx等

额外需要读取.xls文件的jar包–jxl.jar
在这里插入图片描述
jar包在文件包中
文件包地址提取码: nanh

jdk 提取码: itwb 和 Eclipse 提取码: 13rn 的文件包

教程：
Java JDK下载与安装教程
Eclipse解压直接可用

一新建java project项目为GCJ2WGSUtils

点击File–>new–>java project(若没有java project则点击project,再选择java project)
在这里插入图片描述
输入项目名–》点击Finsh

结果：

二创建java文件

右击src–>new–>package–>取名test

在这里插入图片描述

接着右击test,选择new–>class–>取名ExcelUtil（剩余两个文件GCJ2WGS、GCJ2WGSUtils 同这个操作）
在这里插入图片描述
最后粘贴复制内容

三如果想读取excel文件，则引入jxl.jar包

右击选择项目GCJ2WGSUtils 选择Build Path -->configure Build path
在这里插入图片描述

接着选择Libraries–>Add External JSRs–>选中文件包中的jxl.jar文件
在这里插入图片描述

最后目录形式：

在这里插入图片描述
四运行

右击ExcelUtil.java–>run as -->java application

五源码

ExcelUtil.java

package test;

import jxl.Cell;
import jxl.Sheet;
import jxl.Workbook;
import jxl.read.biff.BiffException;


import java.io.FileInputStream;
import java.io.IOException;
import java.io.InputStream;

public class ExcelUtil {

  
   public static void main(String[] args) {
        try{
            jxl.Workbook wb =null;
            InputStream is = new FileInputStream("D://test1.xls");
            wb = Workbook.getWorkbook(is);

            int sheetSize = wb.getNumberOfSheets();
            Sheet sheet = wb.getSheet(0);
            int row_total = sheet.getRows();
            
            for (int j = 0; j < row_total; j++) {
                
                	
                    Cell[] cells = sheet.getRow(j);
                    
            		GCJ2WGSUtils transform = new GCJ2WGSUtils();
            		
            		double value_lat = Double.valueOf(cells[0].getContents());
            		double value_lon = Double.valueOf(cells[1].getContents());
            		
            		double lat=transform.WGSLat(value_lat,value_lon);
            		double lon=transform.WGSLon(value_lat,value_lon);
            		
            		System.out.println("WGSγ׈:"+lat+"  "+"WGSޭ׈:"+lon);
            	

             
                
            }
        }catch (IOException e) {
            // TODO Auto-generated catch block
            e.printStackTrace();
        } catch (BiffException e){
            e.printStackTrace();
        }
    }

}

GCJ2WGS.java

package test;




public class GCJ2WGS {
	
	public static void main(String[] args) {
		
		GCJ2WGSUtils transform = new GCJ2WGSUtils();	
		double lat=transform.WGSLat(122,122);
		double lon=transform.WGSLon(122,122);
		
		System.out.println("WGSγ׈:"+lat);
		System.out.println("WGSޭ׈:"+lon);


	}

}

GCJ2WGSUtils.java

package test;
public class GCJ2WGSUtils {
    //输入GCJ经纬度 转WGS纬度
    public static double WGSLat (double lat,double lon) {
        double PI = 3.14159265358979324;//圆周率
        double a = 6378245.0;//克拉索夫斯基椭球参数长半轴a
        double ee = 0.00669342162296594323;//克拉索夫斯基椭球参数第一偏心率平方
        double dLat = transformLat(lon - 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);
        return (lat - dLat);
    }

    //输入GCJ经纬度 转WGS经度
    public static double WGSLon (double lat,double lon) {
        double PI = 3.14159265358979324;//圆周率
        double a = 6378245.0;//克拉索夫斯基椭球参数长半轴a
        double ee = 0.00669342162296594323;//克拉索夫斯基椭球参数第一偏心率平方
        double dLon = transformLon(lon - 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);
        dLon = (dLon * 180.0) / (a / sqrtMagic * Math.cos(radLat) * PI);
        return (lon - dLon);
    }

    //转换经度所需
    public static double transformLon(double x, double y) {
        double PI = 3.14159265358979324;//圆周率
        double ret = 300.0 + x + 2.0 * y + 0.1 * x * x + 0.1 * x * y + 0.1 * Math.sqrt(Math.abs(x));
        ret += (20.0 * Math.sin(6.0 * x * PI) + 20.0 * Math.sin(2.0 * x * PI)) * 2.0 / 3.0;
        ret += (20.0 * Math.sin(x * PI) + 40.0 * Math.sin(x / 3.0 * PI)) * 2.0 / 3.0;
        ret += (150.0 * Math.sin(x / 12.0 * PI) + 300.0 * Math.sin(x / 30.0 * PI)) * 2.0 / 3.0;
        return ret;
    }
    //转换纬度所需
    public static double transformLat(double x, double y) {
        double PI = 3.14159265358979324;//圆周率
        double ret = -100.0 + 2.0 * x + 3.0 * y + 0.2 * y * y + 0.1 * x * y + 0.2 * Math.sqrt(Math.abs(x));
        ret += (20.0 * Math.sin(6.0 * x * PI) + 20.0 * Math.sin(2.0 * x * PI)) * 2.0 / 3.0;
        ret += (20.0 * Math.sin(y * PI) + 40.0 * Math.sin(y / 3.0 * PI)) * 2.0 / 3.0;
        ret += (160.0 * Math.sin(y / 12.0 * PI) + 320 * Math.sin(y * PI / 30.0)) * 2.0 / 3.0;
        return ret;
    }
    
    
}