android 使用proj4j库，转换部分解释（Java 版纠正高程没有计算）

最新推荐文章于 2024-08-18 03:11:10 发布
最新推荐文章于 2024-08-18 03:11:10 发布
阅读量1.1w 收藏 5
点赞数 1
分类专栏： Java Android 常用开发技术 文章标签： proj4
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_16064871/article/details/83117248
版权

之前写过一篇文章介绍proj4，不过后面测试发现那个库高程没有参与计算，所以自己调试源代码，封装了一个高程参与计算的。搞着也是不容易啊。下面介绍库还是proj4j-0.1.1.jar

之前的文章：android 使用proj4j库（Java版本）

1、封装坐标转换的代码

package com.mapzoom.demo.until;

import org.osgeo.proj4j.CoordinateReferenceSystem;
import org.osgeo.proj4j.CoordinateTransform;
import org.osgeo.proj4j.Proj4jException;
import org.osgeo.proj4j.ProjCoordinate;
import org.osgeo.proj4j.datum.GeocentricConverter;

/**
 * 修复proj4 高程没有转换的bug
 */
public class BLHCoordinateTransform implements CoordinateTransform {
    private CoordinateReferenceSystem srcCRS;
    private CoordinateReferenceSystem tgtCRS;
    private ProjCoordinate geoCoord = new ProjCoordinate(0.0D, 0.0D, 0.0D);
    private boolean doInverseProjection = true;
    private boolean doForwardProjection = true;
    private boolean doDatumTransform = false;
    private boolean transformViaGeocentric = false;
    private GeocentricConverter srcGeoConv;
    private GeocentricConverter tgtGeoConv;

    public BLHCoordinateTransform(CoordinateReferenceSystem srcCRS, CoordinateReferenceSystem tgtCRS) {
        this.srcCRS = srcCRS;
        this.tgtCRS = tgtCRS;
        this.doInverseProjection = srcCRS != null && srcCRS != CoordinateReferenceSystem.CS_GEO;
        this.doForwardProjection = tgtCRS != null && tgtCRS != CoordinateReferenceSystem.CS_GEO;
        this.doDatumTransform = this.doInverseProjection && this.doForwardProjection && srcCRS.getDatum() != tgtCRS.getDatum();
        if (this.doDatumTransform) {
            boolean isEllipsoidEqual = srcCRS.getDatum().getEllipsoid().isEqual(tgtCRS.getDatum().getEllipsoid());
            if (!isEllipsoidEqual) {
                this.transformViaGeocentric = true;
            }

            if (srcCRS.getDatum().hasTransformToWGS84() || tgtCRS.getDatum().hasTransformToWGS84()) {
                this.transformViaGeocentric = true;
            }

            if (this.transformViaGeocentric) {
                this.srcGeoConv = new GeocentricConverter(srcCRS.getDatum().getEllipsoid());
                this.tgtGeoConv = new GeocentricConverter(tgtCRS.getDatum().getEllipsoid());
            }
        }

    }

    public CoordinateReferenceSystem getSourceCRS() {
        return this.srcCRS;
    }

    public CoordinateReferenceSystem getTargetCRS() {
        return this.tgtCRS;
    }

    public ProjCoordinate transform(ProjCoordinate src, ProjCoordinate tgt) throws Proj4jException {
        if (this.doInverseProjection) {
            this.srcCRS.getProjection().inverseProjectRadians(src, this.geoCoord);
            //高程赋值
            this.geoCoord.z = src.z;
        } else {
            this.geoCoord.setValue(src);
        }

        if (this.doDatumTransform) {
            this.datumTransform(this.geoCoord);
        }

        if (this.doForwardProjection) {
            this.tgtCRS.getProjection().projectRadians(this.geoCoord, tgt);
            //高程赋值
            tgt.z = this.geoCoord.z;
        } else {
            tgt.setValue(this.geoCoord);
        }

        return tgt;
    }

    //涉及到七参数 或者 三参数转换
    private void datumTransform(ProjCoordinate pt) {
//        if (this.srcCRS.getDatum().isEqual(this.tgtCRS.getDatum())) {
        if (this.transformViaGeocentric) {
            this.srcGeoConv.convertGeodeticToGeocentric(pt);
            if (this.srcCRS.getDatum().hasTransformToWGS84()) {
                this.srcCRS.getDatum().transformFromGeocentricToWgs84(pt);
            }

            if (this.tgtCRS.getDatum().hasTransformToWGS84()) {
                this.tgtCRS.getDatum().transformToGeocentricFromWgs84(pt);
            }

            this.tgtGeoConv.convertGeocentricToGeodetic(pt);
        }

//        }
    }
}

2、测试代码

package com.mapzoom.demo.activity;

import android.os.Bundle;
import android.support.v7.app.AppCompatActivity;
import android.view.View;
import android.widget.TextView;

import com.mapzoom.demo.R;
import com.mapzoom.demo.until.BLHCoordinateTransform;

import org.osgeo.proj4j.CRSFactory;
import org.osgeo.proj4j.CoordinateReferenceSystem;
import org.osgeo.proj4j.CoordinateTransform;
import org.osgeo.proj4j.CoordinateTransformFactory;
import org.osgeo.proj4j.ProjCoordinate;
import org.osgeo.proj4j.io.Proj4FileReader;

import java.io.IOException;
import java.util.Locale;

public class ProjTestActivity extends AppCompatActivity implements View.OnClickListener {
    TextView showTextView, showTextView1, showTextView2, showTextView3,showTextView4,showTextView5;
    CRSFactory crsFactory = new CRSFactory();
    CoordinateTransformFactory ctf = new CoordinateTransformFactory();
    ProjCoordinate projCoordinate = new ProjCoordinate();
    BLHCoordinateTransform blhCoordinateTransform = null;
    CoordinateReferenceSystem wgs84;
    CoordinateReferenceSystem targetSystem;
    private String showResult = "";
    private CoordinateTransform transform;
    private double dNorth = 0;
    private double dEast = 0;
    private double dHigh = 0;

    @Override
    protected void onCreate(Bundle savedInstanceState) {
        super.onCreate(savedInstanceState);
        setContentView(R.layout.activity_proj_test);

        initWGS84();
        initUI();
    }

    private void initUI() {
        showTextView = (TextView) findViewById(R.id.textView);
        showTextView1 = (TextView) findViewById(R.id.textView1);
        showTextView2 = (TextView) findViewById(R.id.textView3);
        showTextView3 = (TextView) findViewById(R.id.textView4);
        showTextView4 = (TextView) findViewById(R.id.textView5);
        showTextView5 = (TextView) findViewById(R.id.textView6);
        findViewById(R.id.button).setOnClickListener(this);
        findViewById(R.id.button1).setOnClickListener(this);
        findViewById(R.id.button2).setOnClickListener(this);
        findViewById(R.id.button3).setOnClickListener(this);
        findViewById(R.id.button4).setOnClickListener(this);
        String result = String.format(Locale.ENGLISH, "%.3f", 113.0) + "\n" +
                String.format(Locale.ENGLISH, "%.3f", 23.0) + "\n" +
                String.format(Locale.ENGLISH, "%.3f", 10.0) + "\n";
        showTextView5.setText(result);
    }

    private void initWGS84() {
        //wgs84 ����4326
        Proj4FileReader proj4FileReader = new Proj4FileReader();
        String[] paramStr = new String[0];
        //Դ����ϵͳ
        try {
            paramStr = proj4FileReader.readParametersFromFile("epsg", "4326");
        } catch (IOException e) {
            e.printStackTrace();
        }
//        String wgs84_param = "+title=long/lat:WGS84 +proj=longlat +ellps=WGS84 +datum=WGS84 +units=degress";

        wgs84 = crsFactory.createFromParameters("WGS84", paramStr);
    }

    private void initMathTransform() {
        //beijng54  ����2435
        Proj4FileReader proj4FileReader = new Proj4FileReader();
        String[] paramStr = new String[0];
        try {
            paramStr = proj4FileReader.readParametersFromFile("epsg", "2435");
        } catch (IOException e) {
            e.printStackTrace();
        }
        targetSystem = crsFactory.createFromParameters("2435", paramStr);
        transform = ctf.createTransform(wgs84, targetSystem);
        blhCoordinateTransform = new BLHCoordinateTransform(wgs84, targetSystem);
    }

    public void transform(double x, double y, double z) {
        if (transform != null) {
            projCoordinate.setValue(x, y, z);
            transform.transform(projCoordinate, projCoordinate);
            dNorth = projCoordinate.y;
            dEast = projCoordinate.x;
            dHigh = projCoordinate.z;
        }
    }

    public void blhTransform(double x, double y, double z) {
        if (blhCoordinateTransform != null) {
            projCoordinate.setValue(x, y, z);
            blhCoordinateTransform.transform(projCoordinate, projCoordinate);
            dNorth = projCoordinate.y;
            dEast = projCoordinate.x;
            dHigh = projCoordinate.z;
        }
    }

    @Override
    public void onClick(View v) {
        CRSFactory crsFactory = new CRSFactory();
        String result = "";
        String proj4 = "";
        switch (v.getId()) {
            case R.id.button:
                initWGS84();
                initMathTransform();

                transform(113, 23, 10);
                result = String.format(Locale.ENGLISH, "%.3f", dNorth) + "\n" +
                        String.format(Locale.ENGLISH, "%.3f", dEast) + "\n" +
                        String.format(Locale.ENGLISH, "%.3f", dHigh) + "\n";
                showTextView.setText(result);
                break;
            case R.id.button1:
                initWGS84();
                //tmerc 指的是高斯投影
                proj4 = "+proj=tmerc +lat_0=0 +lon_0=114 +k=1 +x_0=500000 +y_0=0 +ellps=krass +towgs84=0,0,0,0,0,0,0 +units=m +no_defs";

                targetSystem = crsFactory.createFromParameters("2435", proj4);
                blhCoordinateTransform = new BLHCoordinateTransform(wgs84, targetSystem);
                transform(113, 23, 10);

                result = String.format(Locale.ENGLISH, "%.3f", dNorth) + "\n" +
                        String.format(Locale.ENGLISH, "%.3f", dEast) + "\n" +
                        String.format(Locale.ENGLISH, "%.3f", dHigh) + "\n";
                showTextView1.setText(result);
                break;
            case R.id.button2:

                initWGS84();
                initMathTransform();

                blhTransform(113, 23, 10);
                result = String.format(Locale.ENGLISH, "%.3f", dNorth) + "\n" +
                        String.format(Locale.ENGLISH, "%.3f", dEast) + "\n" +
                        String.format(Locale.ENGLISH, "%.3f", dHigh) + "\n";
                showTextView2.setText(result);
                break;
            case R.id.button3:
                initWGS84();

                proj4 = "+proj=tmerc +lat_0=0 +lon_0=114 +k=1 +x_0=500000 +y_0=0 +ellps=krass +towgs84=0,0,0,0,0,0,0 +units=m +no_defs";
                targetSystem = crsFactory.createFromParameters("2435", proj4);
                blhCoordinateTransform = new BLHCoordinateTransform(wgs84, targetSystem);
                blhTransform(113, 23, 10);

                result = String.format(Locale.ENGLISH, "%.3f", dNorth) + "\n" +
                        String.format(Locale.ENGLISH, "%.3f", dEast) + "\n" +
                        String.format(Locale.ENGLISH, "%.3f", dHigh) + "\n";
                showTextView3.setText(result);
                break;

            case R.id.button4:
                //+proj=longlat 指的是经纬度坐标
                proj4 = "+proj=longlat +ellps=krass +no_defs";
                wgs84 = crsFactory.createFromParameters("wgs84", proj4);

                proj4 = "+proj=tmerc +lat_0=0 +lon_0=114 +k=1 +x_0=500000 +y_0=0 +ellps=krass +towgs84=0,0,0,0,0,0,0 +units=m +no_defs";
                targetSystem = crsFactory.createFromParameters("2435", proj4);
                blhCoordinateTransform = new BLHCoordinateTransform(wgs84, targetSystem);
                blhTransform(113, 23, 10);

                result = String.format(Locale.ENGLISH, "%.3f", dNorth) + "\n" +
                        String.format(Locale.ENGLISH, "%.3f", dEast) + "\n" +
                        String.format(Locale.ENGLISH, "%.3f", dHigh) + "\n";
                showTextView4.setText(result);
                break;
            default:

                break;
        }
    }


}

3、测试效果

最后面的是原始经纬度，分别转换了。用的是高斯投影。

4、部分prj语句解释

       兰勃脱投影

"+proj=lcc  +a=6378137  +b=6356752.3142  +x_0=500000  +y_0=0  +lon_0=114  +lat_1=35  +lat_2=35  +k=1 +lat_0=0"

       该字符串定义的就是兰勃脱投影（+proj=lcc），椭球长轴为6378137m（+a=6378137），短轴为6356752.3142m（+b=6356752.3142），向东偏移500000m（+x_0=500000），向北偏移0m（+y_0=0），中央子午线是东经114°（+lon_0=114），第一标准纬线是北纬35°（+lat_1=35），第二标准纬线是北纬35°（+lat_2=35），比例为1（+k=1），起始纬度为0°（+lat_0=0，赤道）。

       需要说明的是，第一标准纬线和第二标准纬线一致是为切投影，或者不输入第二标准纬线，两者不一样时，为割投影。

       墨卡托投影

"+proj=merc  +a=6378137  +b=6356752.3142  +x_0=500000  +y_0=0  +lon_0=114  +k=1  +lat_ts=60"

       该字符串定义的是墨卡托投影（+proj=merc），椭球长轴为6378137m（+a=6378137），短轴为6356752.3142m（+b=6356752.3142），向东偏移500000m（+x_0=500000），向北偏移0m（+y_0=0），中央子午线是东经114°（+lon_0=114），比例为1（+k=1），真纬度为北纬60°（+lat_ts=60，也就是旧库里说的割纬线半径）。

       高斯—克吕格投影

"+proj=tmerc  +a=6378137  +b=6356752.3142  +x_0=500000  +y_0=0  +lon_0=114  +k=1  +lat_0=0"

       该字符串定义的高斯投影（+proj=tmerc），椭球长轴为6378137m（+a=6378137），短轴为6356752.3142m（+b=6356752.3142），向东偏移500000m（+x_0=500000），向北偏移0m（+y_0=0），中央子午线是东经114°（+lon_0=114），比例为1（+k=1，当+k=0.996时可变为通用横轴墨卡托投影），起始纬度为0°（+lat_0=0，赤道）。

 

其中（+proj=longlat）代表是经纬度坐标

其中（+proj=geocent）代表是空间坐标坐标

一条语句代表一个椭球下面的坐标。有什么什么椭球下的经纬度坐标（longlat）、空间坐标|（geocent）、投影坐标（平面坐标）。所以投影是有很多类型的，就是不同的平面坐标了。转换出来的。

mmsx
关注
专栏目录
proj4j库（Java版本）
04-26
proj4j库在在业界很出名的。有开源的c代码，这个如果要编译到android 的so完整调用就好了。不过有很多问题。我在Java上找到纯Java写的proj4j库。使用这个库跟epsg就可以很好对应起来。就是每个投影都是有世界通用的编号。我的博客：https://blog.csdn.net/qq_16064871
android 使用proj4j库（Java版本）
mmsx
03-23 6063
proj4j库在在业界很出名的。有开源的c代码，这个如果要编译到android 的so完整调用就好了。不过有很多问题。我在Java上找到纯Java写的proj4j库。使用这个库跟epsg就可以很好对应起来。就是每个投影都是有世界通用的编号。1、库把epsg定义都放在jar包下面，就是打包进行了。可以查看，查找你想要的epsg代号，以及里面内容epsg：2345 大概是广州用的范围。你可以做找到你想...
JAVA引用proj4j
最新发布
weixin_30673943的博客
08-18 56
我整理的一些关于【Java】的项目学习资料（附讲解～～）和大家一起分享、学习一下：https://d.51cto.com/bLN8S1如何在Java项目中引用Proj4j库 在GIS（地理信息系统）工作中，proj4j是一个常用的库，能够帮助我们处理坐标转换。本文将指导你如何在Java项目中引用proj4j库。以下是实现...
proj4 android,Proj4
weixin_29446845的博客
05-28 254
Proj.4是开源GIS最著名的地图投影库，许多GIS开源软件的投影都直接使用Proj.4的库文件。功能主要有经纬度坐标与地理坐标的转换，坐标系的转换，包括基准变换等，下面以命令行方式和编程方式来说明经纬度坐标与地理坐标转换功能的使用。1.球面坐标与投影坐标proj [-control] [+control] [files]-：表示输出控制。+：表示设定投影参数。CTRL+Z退出。正常情况：...
Proj4j.jar
05-15
Proj4是开源GIS最著名的地图投影库，GRASS GIS, MapServer, PostGIS, Thuban, OGDI, Mapnik, TopoCad, GDAL/OGR等软件的投影都直接或间接的使用Proj4。 Proj4的功能主要有经纬度坐标与地理坐标的转换，坐标系的转换，包括基准变换等，下面以命令行方式和编程方式来说明经纬度坐标与地理坐标转换功能的使用。 Java环境下使用Proj4所需要的jar包
编译Android版本的Proj.4库
zzzzllll2006的博客
11-16 1268
Proj4（http://proj4.org/）是一个强大的坐标转换库，官方支持CMake编译，操作非常简单，但却没有编译成Android库的指导。   在Linux系统下编译Android版本的proj库操作如下（别问我Windows下怎么办，我也不知道，不过可以试试msys2）：   1.创建独立工具链 创建一个独立的NDK工具链，具体参考http://blog.csdn.net/s...
GIS开发之- 空间坐标转换需要用到的 proj4j库（Java版本）
03-07
在Java环境下，由于可用的资源相对较少，proj4j库是一个非常宝贵的工具，它为Java开发者提供了进行这种转换的能力，而无需涉及C语言编译的复杂过程。 proj4j是Proj.4库的一个Java实现，Proj.4是一个广泛使用的开源...
proj4j:Proj.4库的Java端口，用于坐标重投影
05-15
Proj4J Proj4J是一个Java库，用于在不同地理空间坐标参考系统之间转换坐标。 它被设计为与proj.4参数兼容，并从proj.4源代码中导出了其某些实现。 在本地构建，测试和安装 mvn clean install
proj4j-0.1.1.jar GIS开发空间坐标转换、投影可以使用
11-01
maven依赖 <groupId>io.jeo <artifactId>proj4j <version>0.1.1 </dependency>
proj4j的 jar 包 roj4j-0.1.0 proj4j-0.1.0-javadoc sources等
03-26
proj4j是一个专门为Java开发的开源库，用于执行地理坐标系统的转换。这个库实现了PROJ.4库的功能，使得在Java环境中处理地理坐标和地图投影变得简单。标题中提到的"roj4j-0.1.0"、"proj4j-0.1.0-javadoc"和"sources...
Android开发之Proj4j坐标转换
gjgfjgy的博客
04-07 367
【代码】Android开发之Proj4j坐标转换。
proj4j_proj4_
09-30
上传文件说明proj4j源码上传主要是dui投影坐标解析
编译后的Proj4开源库
06-19
它的功能主要有经纬度坐标与地理坐标的转换，坐标系的转换，包括基准变换等。
proj4j:OSGeo Proj4j 项目的镜像，有一些更改以上传到 Maven Central
06-23
项目4j OSGeo Proj4j 项目的镜像，有一些更改以上传到 Maven Central。
proj4j-0.1.rar
07-31
proj4 android java库及文档
proj.4 java_Proj4J
weixin_39516865的博客
02-28 766
Welcome to Proj4JProj4J is a Java library to transform point coordinates from one geographic coordinate system to another, including datum transformations.The core of this library is a port of the ​PR...
安卓开发--proj4j坐标转换快速上手
韦_恩的博客
01-08 1181
本节专门讲解proj4j的使用，关于坐标转换原理，等我后面有时间了再讲解下。 现在最新版本的1.3.0版本使用起来很方面，不需要你手动导包，只要配置下gradle就行了。
proj4 Android交叉编译
chuwuxi0805的博客
09-21 358
和我那篇GDAL Android交叉编译博客基本一致: 首先要配置工具链位置，然后是设置sysroot的位置 关键是configure的配置： 1、工具链配置 2、configure ./configure --prefix=/home/root/proj-4.9...
Proj4J库 0.2.0版本maven地址
06-02
Proj4J库的Maven中央仓库地址为： ```xml <dependency> <groupId>org.osgeo.proj4j</groupId> <artifactId>proj4j</artifactId> <version>0.2.0</version> </dependency> ``` 你可以将上述代码添加到Maven项目的POM文件中，然后使用Maven或其他构建工具下载和安装Proj4J库的0.2.0版本。注意，这需要你的项目配置了Maven中央仓库或其他可以访问该库的仓库。
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值