目录

前言

一、空间表设计

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<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当中如何进行度分秒位置的转换，最后结合实现原型进行详细的说明。通过本文，您可以掌握如何管理这些矢量数据，同时知道如何进行空间数据的去重，还可以掌握如何进行度分秒数据的转换，欢迎感兴趣的朋友互相交流。行文仓促，难免有不足之处，欢迎朋友们在评论区批评指正，不慎感谢。