基于Java的度分秒坐标转纯经纬度坐标的漂亮国基地信息管理

目录

前言

一、空间表设计

1、物理表结构

二、后台数据管理

1、数据去重

2、去重的具体实现

3、度分秒数据格式转换

4、具体的转换方法

5、新增界面的实现

三、数据管理界面

总结

---

前言

        众所周知，漂亮国在全球范围内部署了大量的基地，用以维持其打击能力。这里无意于讨论这些部署的事情。只是作为一个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 linkedHashMap = new LinkedHashMap();
for(int i=0; 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 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方法接收参数，然后调用转换逻辑。

"form-group">    
    "col-sm-3 control-label">位置（度分秒）：
    
"col-sm-8">
            "latLonStr" class="form-control" type="text" required>