基于Java的度分秒坐标转纯经纬度坐标的漂亮国基地信息管理的方法
作者:夜郎king
前言
众所周知,漂亮国在全球范围内部署了大量的基地,用以维持其打击能力。这里无意于讨论这些部署的事情。只是作为一个GIS开发者,使用地理信息的视角来看看其到底分布在全球的哪些地方,在之前的博客中,我们曾经介绍了使用Java去解析KMZ和KML数据,这份数据就是其它博主分享的漂亮国海外基地的分布情况。这是一份数据源,同时我们在某度百科上也能找到一些公开的资料,上面公布了一些国内或者海外基地的坐标信息。
要想实现在WebGIS中直观的显示这些位置,首先要解决数据的管理问题。之前我们在使用GDAL的KMZ数据解析中,通过在数据中查看具体的数据可以看到。其提供的数据中,应该是有重复的。我们采用基地名称和位置wkt字符串作为分组查询,可以发现有很多的数据的条数是大于2的。因此可以判定出是有数据的重复,因此要求我们对入库的数据进行去重。同时由于在互联网上公开的数据中,有一些位置信息是以度分秒的形式给出的。对于常见的坐标,我们都是采用纯数字的经纬度来展示的。那么怎么将度分秒的位置信息转换成常见的经纬度信息呢?
这就是本文的由来,本文将以java语言为例,详细介绍如何管理漂亮国的基地信息。为下一步全球的空间可视化打下坚实的基础,首先介绍如何对数据进行去重处理,然后介绍在java当中如何进行度分秒位置的转换,最后结合实现原型进行详细的说明。通过本文,您可以掌握如何管理这些矢量数据,同时知道如何进行空间数据的去重,还可以掌握如何进行度分秒数据的转换,欢迎品鉴。
一、空间表设计
这里着重介绍空间表的设计与实现,让大家了解空间表的设计。这里直接分享表的物理结构,通过SQL语句直接分享给大家,如果大家感兴趣也可以进行实践。
1、物理表结构
闲言少叙,这里将设计的表字段信息以sql语句的形式提供出来,供大家参考,表中的字段,大家可以根据实际情况进行增减。
物理表结构的sql语句如下(如果您想在本地记性测试,请一定要安装postgis扩展,否则空间数据类型将无法创建):
CREATE TABLE "public"."biz_usa_military_base" ( "id" int8 NOT NULL, "en_name" varchar(255) COLLATE "pg_catalog"."default" NOT NULL, "en_desc" varchar(1024) COLLATE "pg_catalog"."default", "cn_name" varchar(255) COLLATE "pg_catalog"."default", "remark" varchar(255) COLLATE "pg_catalog"."default", "geom" "public"."geometry", "create_by" varchar(64) COLLATE "pg_catalog"."default", "create_time" timestamp(6), "update_by" varchar(64) COLLATE "pg_catalog"."default", "update_time" timestamp(6), "type" int4, "en_country" varchar(255) COLLATE "pg_catalog"."default", "cn_country" varchar(255) COLLATE "pg_catalog"."default", "en_city" varchar(255) COLLATE "pg_catalog"."default", "cn_city" varchar(255) COLLATE "pg_catalog"."default", CONSTRAINT "pk_biz_usa_military_topics" PRIMARY KEY ("id") ); CREATE INDEX "idx_biz_usa_military_base_geom" ON "public"."biz_usa_military_base" USING gist ( "geom" "public"."gist_geometry_ops_2d" ); COMMENT ON COLUMN "public"."biz_usa_military_base"."id" IS '主键'; COMMENT ON COLUMN "public"."biz_usa_military_base"."en_name" IS '英文名称'; COMMENT ON COLUMN "public"."biz_usa_military_base"."en_desc" IS '英文描述'; COMMENT ON COLUMN "public"."biz_usa_military_base"."cn_name" IS '英文名称'; COMMENT ON COLUMN "public"."biz_usa_military_base"."remark" IS '备注'; COMMENT ON COLUMN "public"."biz_usa_military_base"."geom" IS '空间信息'; COMMENT ON COLUMN "public"."biz_usa_military_base"."create_by" IS '创建人'; COMMENT ON COLUMN "public"."biz_usa_military_base"."create_time" IS '创建时间'; COMMENT ON COLUMN "public"."biz_usa_military_base"."update_by" IS '更新人'; COMMENT ON COLUMN "public"."biz_usa_military_base"."update_time" IS '更新时间'; COMMENT ON COLUMN "public"."biz_usa_military_base"."type" IS '基地类型,1海外 0 本土'; COMMENT ON COLUMN "public"."biz_usa_military_base"."en_country" IS '部署国家英文名'; COMMENT ON COLUMN "public"."biz_usa_military_base"."cn_country" IS '部署国家中文名'; COMMENT ON COLUMN "public"."biz_usa_military_base"."en_city" IS '部署城市英文名'; COMMENT ON COLUMN "public"."biz_usa_military_base"."cn_city" IS '部署城市中文名'; COMMENT ON TABLE "public"."biz_usa_military_base" IS '美军军事基地信息表';
在设计好上述的空间数据表之后,我们就将表在PostGIS数据库中进行新建,确保没有报错。下面再来进行后台去重识别和信息转换的工作。
二、后台数据管理
后台数据的管理比较简单,这里不涉及复杂的逻辑,首先解析KMZ的数据,这部分的内容不再赘述,仅提供如何进行数据去重的处理。然后如何进行经纬度不同数据格式的转换,因此将数据同步成同一个标准。
1、数据去重
在之前的数据处理过程当中,我们将所有的基地信息都读取出来,然后统一放到一个List集合当中,最后调用MP的批量插入接口实现数据的批量插入。但是这样的数据是很容易重复的。查询数据是否重复的sql语句如下,我们采用的机制就是根据英文名称来进行检索,我们在筛选数据时,着重检查基地的英文名和位置信息是否一致,如果这两个完全一致,那么我们认为就是同一个基地,是需要将数据进行去除的。
select id,en_name,cn_name,st_asgeojson(geom) from biz_usa_military_base order by st_asgeojson(geom);
上面的数据由于我已经处理过了,因此已经将重复的数据都已经全部去掉了。默认情况下,是有蛮多重复的数据。在实际项目中,我们怎么解决这种重复数据呢?一般可以有两种思路。第一种就是循环插入,这样在插入前做一次重复查询,如果数据已经存在了就不再插入。这样做的效率是比较低的,因为没插一次还要去数据库中查询数据是否存在,消耗的时间是比较大。既然这种方式不太好,那么第二种方式就是在批量插入,但是要求是必须要在批量的集合中进行去重,这样从源头上就解决了重复插入的问题。
2、去重的具体实现
为了实现去重,我们这里选用LinkedHashMap来进行处理,选这个数据结构是因为LinkedHashMap可以实现有序性,数据会按照读取的顺序添加到集合当中。同时基于HashMap的特性可以实现重复数据的处理,map的key我们选用字符串来处理,由基地的名称+空间的WKT字符串,这样子虽然长点,但是可以进行重复的校验。在map的结构当中,如果key是重复的话,会以最后的一次put为准,这样就实现我们的去重。代码实现就更加简单了:
Map<String,UsaMilitaryBase> linkedHashMap = new LinkedHashMap<String, UsaMilitaryBase>(); for(int i=0; i<featureCount; i++){ Feature feature = layer.GetFeature(i); if(null == feature) continue; String enName = feature.GetFieldAsString("Name"); UsaMilitaryBase base = new UsaMilitaryBase(); Geometry geom = feature.GetGeometryRef(); //step 1、生成原始wkt String wkt = geom.ExportToWkt(); wkt = "SRID=" + srid +";" + wkt;//拼接srid,实现动态写入 linkedHashMap.put(enName + wkt, base); }
然后通过这种方式得到的map就是实现了去重之后的数据集合。当然,如果想实现批量插入,就需要循环map,将map的value取出来,然后添加到List中,再调用MP的批量插入的接口。
if(linkedHashMap.size() >0) { // 使用for-each循环遍历HashMap for (Map.Entry<String, UsaMilitaryBase> entry : linkedHashMap.entrySet()) { dataList.add(entry.getValue()); } usaMilitaryBaseService.saveBatch(dataList, 300);//批量插入,每300条提交一次事务 }
3、度分秒数据格式转换
在空间数据展示中,度分秒也是一种非常常见的格式,因此很多地方也是以度分秒的形式来进行展示的。其中,我们在某度上看到的漂亮国基地信息都是以度分秒的形式展示的,如下所示:
其数据的表达方式都是:
华盛顿州:
序号 基地名称、经纬度
1、 陆军刘易斯堡 47° 4'57.21"N 122°35'2.78"W
2、 惠德贝岛海军航空站 48°20'42.49"N 122°39'51.71"W
3、 埃弗雷特海军基地 47°59'20.32"N 122°13'12.91"W
4、 普吉海峡海军造船厂 47°33'41.32"N 122°37'56.17"W
像这些数据怎么保存到数据库中,怎么进行wkt类型的转换呢?这里分享java的实现方式。
4、具体的转换方法
关于如何将度分秒的数据进行转换,网上有很多例子。我也是参考了别人的代码,但是很多代码其实没有将清楚,比如南纬和西经的数据处理方法,也就是会有负数的情况出现。如果按照正常的计算逻辑,则会出现位置错误。下面分享在网友的基础之上改造的可以带符号的经纬度转换方法。在后台的控制器中增加以下代码:
@RequiresPermissions("mt:usabase:add") @Log(title = "美军军事基地", businessType = BusinessType.INSERT) @PostMapping("/add") @ResponseBody public AjaxResult addSave(UsaMilitaryBase militaraBase,String latLonStr){ if(StringUtils.isNotEmpty(latLonStr)) { //String lat1 = "34°49'36.28\"N 95°57'33.18\"W"; String lat = latLonStr.split(" ")[0]; String lon = latLonStr.split(" ")[1]; String latSuffix = lat.substring(lat.length() -1 );//取出纬度后缀N,S String lonSuffix = lon.substring(lon.length() -1 );//取出纬度后缀W,E String newLat = String.valueOf(LatLngUtil.dfm2LatLng(lat.substring(0,lat.length() -1 ))); String newLon = String.valueOf(LatLngUtil.dfm2LatLng(lon.substring(0,lon.length() -1 ))); newLat = latSuffix.equalsIgnoreCase("S") ? "-" + newLat : newLat; newLon = lonSuffix.equalsIgnoreCase("W") ? "-" + newLon : newLon; //构造wkt,默认参考为4326 String geom = "SRID=" + 4326 + ";POINT (" + newLon + " " + newLat + ")"; militaraBase.setGeom(geom); } return toAjax(mbaseService.insertEntity(militaraBase)); }
需要注意的是,在数据中,要解析出数据中的南纬和西经,也就是S和W,如果位置后面是这两个字符,需要将结果转换成负数即可。处理具体转换的代码如下,我整理成了一个工具类:
/** * 度分秒转经纬度 * * @param dms 116°25'7.85" * @return 116.418847 */ public static double dfm2LatLng(String dms) { if (dms == null) return 0; try { dms = dms.replace(" ", ""); String[] str2 = dms.split("°"); if (str2.length < 2) return 0; int d = Integer.parseInt(str2[0]); String[] str3 = str2[1].split("\'"); if (str3.length < 2) return 0; int f = Integer.parseInt(str3[0]); String str4 = str3[1].substring(0, str3[1].length() - 1); double m = Double.parseDouble(str4); double fen = f + (m / 60); double du = (fen / 60) + Math.abs(d); if (d < 0) du = -du; return Double.parseDouble(String.format("%.7f", du)); } catch (Exception e) { e.printStackTrace(); } return 0; }
5、新增界面的实现
这里简单分享新增界面的实现,界面实现采用Thymeleaf,并不是流行的vue,喜欢vue的可以自己自行改造哈。将度分秒信息一次性传给后台,后台使用上面展示的add方法接收参数,然后调用转换逻辑。
<div class="form-group"> <label class="col-sm-3 control-label">位置(度分秒):</label> <div class="col-sm-8"> <input name="latLonStr" class="form-control" type="text" required> </div> </div>
三、数据管理界面
本节将重点讲解数据管理界面以及功能介绍。
列表管理界面实现数据的查询展示,展示基地的英文名称、中文名称,数据的创建时间,支持按照基地的中英文名称进行查询。
新增基地信息页面,请注意在位置这里,输入度分秒格式的数据,同时输入基地的中英文名称、所在国家的中英文名称,所在城市的中英文名称。输入完毕后,点击确定即可。
信息编辑页面,与新建的窗口界面不同的是,新增界面可以输入位置坐标信息,而修改暂时不提供这个功能(技术上没有问题,只是从业务上禁止修改)。
总结
以上就是本文的主要内容,本文将以java语言为例,详细介绍如何管理漂亮国的基地信息。为下一步全球的空间可视化打下坚实的基础,首先介绍如何对数据进行去重处理,然后介绍在java当中如何进行度分秒位置的转换,最后结合实现原型进行详细的说明。通过本文,您可以掌握如何管理这些矢量数据,同时知道如何进行空间数据的去重,还可以掌握如何进行度分秒数据的转换,欢迎感兴趣的朋友互相交流。行文仓促,难免有不足之处,欢迎朋友们在评论区批评指正,不慎感谢。
到此这篇关于基于Java的度分秒坐标转纯经纬度坐标的漂亮国基地信息管理的文章就介绍到这了,更多相关Java度分秒坐标转纯经纬度内容请搜索脚本之家以前的文章或继续浏览下面的相关文章希望大家以后多多支持脚本之家!