如何把空间数据从CGCS2000转换到WGS84和BD09 |
您所在的位置：网站首页 › wgs坐标转换为cgcs2000 › 如何把空间数据从CGCS2000转换到WGS84和BD09 |
如何把空间数据从CGCS2000转换到WGS84和BD09
|
目录
1 背景说明1.1 空间数据1.2 CGCS20001.3 WGS841.4 BD09
2 读取数据库中的geometry数据3 代码实现3.1 CGCS2000转WGS843.2 BD09坐标系中点的表示3.3 WGS84转BD093.4测试结果
4 依赖和环境5 提示6 结语
1 背景说明
最近的项目中涉及到空间数据的转换，主要是把CGCS200坐标转换为WGS84坐标和BD09坐标。我在网上找了许久，只找到了WGS84坐标转BD09坐标的工具类，但没有找到CGCS200转换到WGS84坐标的工具类。我实在没有办法，就只有自己写。但我没有GIS行业的教育背景，不了解坐标转换的专业知识，不过我断定GIS行业一定有现成的相关方法。通过查阅相关资料，最后在ArcGIS提供的 ArcGIS Runtime API for Java中找到了相关API。我借此实现了个工具类，测试效果也还不错。现在写出这篇文章，共享我的代码，希望能帮到有同样问题的开发者。 为了方便GIS行业外的朋友更流畅地阅读本文，先阐释下相关名词概念。 1.1 空间数据从数据库的角度来看，是指在表中以geometry类型存在的数据。 从java语言的角度来看，是指com.esri.arcgisruntime.geometry包下各种类的对象,包括Point、Multipoint、Polyline、Polygon等等. 1.2 CGCS2000CGCS2000：China Geodetic Coordinate System 2000，2000国家大地坐标系，是我国当前最新的国家大地坐标系。该坐标系在ArcGIS中有众多子坐标系。如下图所示。本文选用的是CGCS2000_3_Degree_GK_CM_105E。下方的坐标系详情中有WKID，您需要记录下这个值，后续会用到。 WGS84：World Geodetic System 1984，是为GPS全球定位系统使用而建立的坐标系统，在全球有广泛的使用。 1.4 BD09BD09：百度地图坐标系，它是在标准Web墨卡托的基础上进行GCJ-02加偏之后，再加上百度自身的加偏算法。 2 读取数据库中的geometry数据我采用的方式是在数据库用SQL语句查询出geometry字段的String表达形式，再做后续处理。代码片段如下： @Query( value = "select SHAPE.STAsText() FROM LAND WHERE ID=?", nativeQuery = true) String getShapeString (String id);下面这个字符串是本文通过此方法在数据库中查询出来的，此字符串表示一个地块的坐标，逗号分隔出每个坐标点，多个坐标点连成一个面。(此数据是测试数据，由笔者随意编造，没有描述任何真实意义) POLYGON ((370907.330099999904697 3418398.496099999174576, 370899.395100000314615 3418410.75789999861853, 370909.574500000104318 3418415.690899999812245)) 3 代码实现 3.1 CGCS2000转WGS84上面的这个字符串表示的是CGCS2000坐标系下的数据。下面的代码实现了将多边形地块坐标转换为WGS84，再从WGS84转换为BD09。 package cn.wja.utils; import cn.wja.model.BD09Point; import cn.wja.model.PolygonResult; import com.esri.arcgisruntime.geometry.*; import java.util.ArrayList; import java.util.List; /** * 用于坐标转换的工具类 * * @author wang.jingan * @since 2019-07-20 10:39:16 */ public class CoordinateSystemTransform { /** * 通过String字符串得到三个坐标系的空间数据对象 * * @param shape 转化为String的Geometry数据 * @return PolygonResult对象，内含三个坐标系的空间数据对象 */ public static PolygonResult shapeStringToPolygons(String shape) { //获取两个坐标系 SpatialReference cgcs2000 = SpatialReference.create(4544); SpatialReference wgs84 = SpatialReference.create(4326); //创建一个CGCS2000的点收集器，为创建Polygon对象做准备 PointCollection land = new PointCollection(cgcs2000); //准备返回的百度坐标点集合 List pointsBD09 = new ArrayList(); //获取需要的字符串 String shapeSubstring = shape.substring(10, shape.length() - 2); //得到每个点的原始数据 String[] points = shapeSubstring.split(","); //遍历每个点 for (String pointString : points) { //去前后的空格 pointString = pointString.trim(); //把经纬度分开成两个数据 String[] xy = pointString.split(" "); Double x = Double.valueOf(xy[0]); Double y = Double.valueOf(xy[1]); //把点添加到点收集器 land.add(x, y); //得到CGCS2000点数据 Point pointCgcs2000 = new Point(x, y, cgcs2000); //把CGCS2000点数据转化为WGS84点数据 Point pointwgs84 = (Point) GeometryEngine.project(pointCgcs2000, wgs84); //获取点数据中的经纬度,X为经度,Y为纬度 double pointwgs84X = pointwgs84.getX(); double pointwgs84Y = pointwgs84.getY(); //把wgs84点坐标转化为bdo9点坐标 double[] bd09 = GPSUtil.gps84_To_bd09(pointwgs84Y, pointwgs84X); BD09Point coordinatesDTO = new BD09Point(); coordinatesDTO.setLat(bd09[0]); coordinatesDTO.setLon(bd09[1]); pointsBD09.add(coordinatesDTO); } //从点收集器中创建Polygon，得到一个cgcs2000的Polygon Polygon polygonCgcs2000 = new Polygon(land); //转换为wgs84的Polygon Polygon polygonWgs84 = (Polygon) GeometryEngine.project(polygonCgcs2000, wgs84); //装入返回数据 PolygonResult result = new PolygonResult(); result.setCgcs2000(polygonCgcs2000); result.setWgs84(polygonWgs84); result.setBd09(pointsBD09); return result; } } 3.2 BD09坐标系中点的表示在ArcGIS中wkid用于表示相关坐标系的代码， 并用SpatialReference.create（int wkid） 进行创建对应的坐标系。但BD09坐标系 在ArcGIS中没有相关的wkid，所以BD09坐标是通过点的集合表示，而WGS84的坐标是通过com.esri.arcgisruntime.geometry 包下的Polygon类的对象表示。 BD09坐标系的点是笔者自己定义的一个类。代码如下： package cn.wja.model; public class BD09Point { private Double lat; private Double lon; public Double getLat() { return lat; } public void setLat(Double lat) { this.lat = lat; } public Double getLon() { return lon; } public void setLon(Double lon) { this.lon = lon; } @Override public String toString() { return "BD09Point{" + "lat=" + lat + ", lon=" + lon + '}'; } } 3.3 WGS84转BD09此工具类网上有很多，同时涉及到很多测量学和数学知识，我也不太擅长那些方面，所以就没有自己写，直接用了 高德，百度，Google地图定位偏移以及坐标系转换 这篇文章的代码，您也可以自行在网上寻找类似工具。 3.4测试结果如下图所示，测试结果未发现异常。 这篇文章首次发表是2019年7月，当时我用使用的ArcGIS Runtime SDK for Java版本为100.4.0。 com.esri.arcgisruntime arcgis-java 100.4.0当时是只需要在pom.xml中加入上面的依赖即可，默认配置的maven仓库或者阿里云仓库都可以自动下载jar包，也不需要安装任何ArcGIS环境，就可以使用运行。 但文章发布后大约半年，就陆续有开发者留言反馈说，不能正确加载jar包，报错提示“程序包com.esri.arcgisruntime.geometry不存在”。开始我以为是他们的Maven配置出了问题，我就手动下载100.8.0的jar文件，放在了我的CSDN资源中，供需要的朋友下载。之所以下载100.8.0版本，是因为原来的100.4.0的jar包我也没有刻意保存，也找不到了。但解决完jar包加载问题后，又有不少开发者反馈出现了如下报错： Caused by: java.lang.RuntimeException: Could not find runtime in any of: - A directory specified by calling ArcGISRuntimeEnvironment.setInstallDirectory() - The current directory F:\Code\Intimate\other\gis-demo - A location specified by the environment variable ARCGISRUNTIMESDKJAVA_100_10_0 - Within the “.arcgis” directory in the user’s home path C:\Users\admin.arcgis 我在自己的电脑上在测试，也重现了这个问题。开始我以为是100.8.0作为更高的版本，有其新的要求。我就花了不少功夫到arcgis的官网中找到并下载到了100.4.0的jar包，但是运行后也会报一样的错。我看这个jar包上传时间是2020年10月30日，估计是新下载的100.4.0，相对于2019年的文件，也做了些更新吧。 其实上述问题的实质就是ArcGIS Runtime API for Java的运行环境没有成功构建。 那么解决办法就只有两个了： 第一个方法是找到老版本的100.4.0.jar。我尝试了很久，都没有找到相关资源，所以就暂时放弃了，如果有读者找到，可以联系我。第二个方法是使用新版本的jar包，并在此基础上构建ArcGIS Runtime API for Java的运行环境。我进行了一些探索后，成功实现了。 解决办法请参见我新写的另一篇博客：《如何构建ArcGIS Runtime API for Java的运行环境》（左侧博客标题有超链接，请点击跳转）。 您也可以直接下载示例代码：ArcGIS Runtime API for Java使用示例：CGCS2000转WGS84和BD09 （左侧资源标题含有超链接，请点击跳转 ）。 5 提示虽然我是通过公开途径下载的 ArcGIS Runtime API for Java 的相关jar包和资源文件，并未曾接收到任何关于收费的提示，但我也未查询到 ArcGIS Runtime API for Java 的供应商Esri 有表明这款软件是可以免费使用的开源软件。因此如需在实际项目开发中使用ArcGIS Runtime API for Java，请先向Esri公司咨询授权的相关事宜，避免不必要的纠纷和损失。 我在撰写本文时就职的公司是购买了ArcGIS 的软件使用权，在使用 ArcGIS Runtime API for Java时，也得到了Esri的技术支持。 如果您希望使用开源软件完成类似功能，建议查询一下geotools。 6 结语本文的参考文献已经在正文中注明了出处。如果转载本文，也请注明出处和作者。这篇文章是我在CSDN上写的第一篇博客，时间仓促，技术有限，也不太了解CSDN上写博客的技巧，难免有些纰漏，欢迎读者指出其中的不足，也希望大家在博客下的留言区积极讨论。 |
【本文地址】
今日新闻 |
推荐新闻 |