台湾20m DEM高程瓦片数据介绍
一、背景
读过之前文章的读者都知道，中国台湾地区12.5m的高程DEM瓦片已经发布,满足地形可视化没有问题，但是如果要拿来做地形分析以及精确可视化，数据方面还有一些问题，比如台湾周边海面高度是16米，而海岸地形高度为1米，很显然是不太合理，如果需要使用，需要对数据进行裁剪修整。不过台湾地区推出了20m的高程数据，因此，我将这份数据下载后，通过投影变换、瓦片切分等过程，生成了今天要介绍的这份数据。
二、数据介绍
这次带来的是分辨率为20m的地形瓦片数据，虽然分辨率不及12.5m，但是它有以下优点：
1.水平精度和垂直精度有保证，中国台湾地区官方测绘部门生产的数据，精度上相对12.5m数据、ALOS及哥白尼DEM一众卫星遥感生产出来的高程数据更优。
2.现势性非常好，中国台湾测绘部门每年都会更新发布高程数据，本次生产地形瓦片使用的是最新2022年的数据，热腾腾刚出炉的。
3.数据细节更优，相对于12.5m的分辨率，分辨率差了点，但是数据细节更丰富。
以中山高速以及汐止-杨梅高架道路为例，下图是12.5m数据晕渲图截图：
可以看到，道路无法在12.5m数据中显现；下面看看该地区20m最新高程数据晕渲图截图：
在20m数据中能清晰看到道路的轮廓，同时桃园市龟山区地势相对平坦。由于一部分区域属于军事管制区没有对外公开，导致这份数据苗栗县缺失了一部分,通过原始数据晕渲图效果可以清晰看到缺失部分区域范围。
下面介绍这份瓦片数据的一些基本信息，这份数据分辨率是20m，和30m一样，最高只能切到14级。
可以看一下layer.json配置文件信息
三、数据效果
将数据使用Nginx进行发布，具体发布方法可以参考我之前的地形瓦片发布文章，这里给出nginx中nginx.conf的配置server方法
server
{
listen 8802;
server_name localhost;
proxy_set_header Host $host;
proxy_set_header X-Real-IP $remote_addr;
proxy_set_header X-Forwarded-For $proxy_add_x_forwarded_for;
add_header Access-Control-Allow-Origin "$http_origin";
add_header Access-Control-Allow-Credentials true;
add_header Access-Control-Allow-Methods 'GET, POST, OPTIONS';
add_header Access-Control-Allow-Headers 'DNT,X-Mx-ReqToken,Keep-Alive,User-Agent,X-Requested-With,If-Modified-Since,Cache-Control,Content-Type,Authorization';
if ($request_method = 'OPTIONS')
{
return 204;
}
location / {
autoindex on;
root D:/Data/DEM/DEM-20M-WGS84-CHINA-TAIWAN-MESH/;
#index index.html index.htm;
}
error_page 500 502 503 504 /50x.html;
location = /50x.html {
root html;
}
}
前端使用Cesium进行接入，具体接入关键代码如下
//创建viewer
const viewer = new Cesium.Viewer('cesium-app', {
selectionIndicator: false,
baseLayerPicker: false,
homeButton: false,
infoBox: false,
timeline: false,
navigationHelpButton: false,
navigationInstructionsInitiallyVisible: false,
vrButton: false,
fullscreenButton: false,
geocoder: false,
animation: false,
sceneModePicker: false,
shadows:true,
imageryProvider: new Cesium.ArcGisMapServerImageryProvider({
url: 'https://services.arcgisonline.com/ArcGIS/rest/services/World_Imagery/MapServer'
})
});
viewer.cesiumWidget.creditContainer.style.display = "none";
viewer.scene.globe.enableLighting = false;
viewer.scene.highDynamicRange = true;
viewer.scene.msaaSamples = 16;
viewer.scene.fog.enabled = true;
viewer.scene.fog.enableLighting = true;
//加载地形
viewer.terrainProvider = new Cesium.CesiumTerrainProvider({
url:"http://192.168.1.77:8802/",
tilingScheme: new Cesium.GeographicTilingScheme(),
requestVertexNormals:true
});
首先创建viewer对象，设置大气、光照、抗锯齿等参数增强地形表达效果；再通过CesiumTerrainProvider创建viewer的terrainProvider对象，这样数据就加载到场景当中了。在这里说明一下,我的后端nginx和前端cesium不是一台机器，后端机器ip是192.168.1.77,开放的端口号是8802，因此我们接入的url全称是"http:\192.168.1.77::8802"，也需要设置tilingScheme,我们的数据都是wgs84投影，所以设置的是GeographicTilingScheme对象；同时设置requestVertexNormals为true，因为地形瓦片是Mesh格式，在生成Mesh格式terrain数据时，也生成了每个点的法向量，这样开启光照，地形会有光照明暗阴影，效果更加自然。
下面看一下加载到Cesium上的截图
通过曾文水库大坝截图可以看到，大坝轮廓比较清晰，大坝里的水体与坝底由于高程落差，在三维可视化中表现很直观。
四、结语
总体而言，虽然这份数据缺失了一部分，但是瑕不掩瑜，不失是一份质量很高的高程数据。对于这份数据，大家有问题可以后台私信，或者在评论区留言，我将及时回复大家的问题。