• cesium加载高层级离线影像地图瓦片(天地图、19级Arcgis)


    实际加载效果如图:

    1、下载离线地图瓦片方式(多种任选其一,个人倾向于Qgis工具):

    方式1、采用第三方下载工具如:91卫图、水经注、全能电子地图下载器、bigemap等等。(这些有的下载层级不够,有的需要付费,结合个人情况选择)

    方式2爬虫技术实现,网上贴的相关代码如:Java为例

    1. import java.io.BufferedOutputStream;
    2. import java.io.File;
    3. import java.io.FileOutputStream;
    4. import java.io.InputStream;
    5. import java.io.OutputStream;
    6. import java.util.concurrent.ExecutorService;
    7. import java.util.concurrent.Executors;
    8. import org.apache.http.HttpEntity;
    9. import org.apache.http.HttpResponse;
    10. import org.apache.http.HttpStatus;
    11. import org.apache.http.client.methods.HttpGet;
    12. import org.apache.http.impl.client.CloseableHttpClient;
    13. import org.apache.http.impl.client.HttpClients;
    14. public class map{
    15. //矢量(行政) - 等经纬度
    16. public static String vec_c = "http://{server}.tianditu.gov.cn/vec_c/wmts?SERVICE=WMTS&REQUEST=GetTile&VERSION=1.0.0&LAYER=vec&STYLE=default&TILEMATRIXSET=c&FORMAT=tiles&TILEMATRIX={z}&TILEROW={y}&TILECOL={x}&tk={tk}";
    17. //矢量(行政) - 墨卡托
    18. public static String vec_w = "http://{server}.tianditu.gov.cn/vec_w/wmts?SERVICE=WMTS&REQUEST=GetTile&VERSION=1.0.0&LAYER=vec&STYLE=default&TILEMATRIXSET=w&FORMAT=tiles&TILEMATRIX={z}&TILEROW={y}&TILECOL={x}&tk={tk}";
    19. //矢量注记(行政) - 等经纬度
    20. public static String cva_c = "http://{server}.tianditu.gov.cn/cva_c/wmts?SERVICE=WMTS&REQUEST=GetTile&VERSION=1.0.0&LAYER=cva&STYLE=default&TILEMATRIXSET=c&FORMAT=tiles&TILEMATRIX={z}&TILEROW={y}&TILECOL={x}&tk={tk}";
    21. //矢量注记(行政) - 墨卡托
    22. public static String cva_w = "http://{server}.tianditu.gov.cn/cva_w/wmts?SERVICE=WMTS&REQUEST=GetTile&VERSION=1.0.0&LAYER=cva&STYLE=default&TILEMATRIXSET=w&FORMAT=tiles&TILEMATRIX={z}&TILEROW={y}&TILECOL={x}&tk={tk}";
    23. //影像 - 等经纬度
    24. public static String img_c = "http://{server}.tianditu.gov.cn/img_c/wmts?SERVICE=WMTS&REQUEST=GetTile&VERSION=1.0.0&LAYER=img&STYLE=default&TILEMATRIXSET=c&FORMAT=tiles&TILEMATRIX={z}&TILEROW={y}&TILECOL={x}&tk={tk}";
    25. //影像 - 墨卡托
    26. public static String img_w = "http://{server}.tianditu.gov.cn/img_w/wmts?SERVICE=WMTS&REQUEST=GetTile&VERSION=1.0.0&LAYER=img&STYLE=default&TILEMATRIXSET=w&FORMAT=tiles&TILEMATRIX={z}&TILEROW={y}&TILECOL={x}&tk={tk}";
    27. //影像注记 - 等经纬度
    28. public static String cia_c = "http://{server}.tianditu.gov.cn/cia_c/wmts?SERVICE=WMTS&REQUEST=GetTile&VERSION=1.0.0&LAYER=cia&STYLE=default&TILEMATRIXSET=c&FORMAT=tiles&TILEMATRIX={z}&TILEROW={y}&TILECOL={x}&tk={tk}";
    29. //影像注记 - 墨卡托
    30. public static String cia_w = "http://{server}.tianditu.gov.cn/cia_w/wmts?SERVICE=WMTS&REQUEST=GetTile&VERSION=1.0.0&LAYER=cia&STYLE=default&TILEMATRIXSET=w&FORMAT=tiles&TILEMATRIX={z}&TILEROW={y}&TILECOL={x}&tk={tk}";
    31. //地形 - 等经纬度
    32. public static String ter_c = "http://{server}.tianditu.gov.cn/ter_c/wmts?SERVICE=WMTS&REQUEST=GetTile&VERSION=1.0.0&LAYER=ter&STYLE=default&TILEMATRIXSET=c&FORMAT=tiles&TILEMATRIX={z}&TILEROW={y}&TILECOL={x}&tk={tk}";
    33. //地形 - 墨卡托
    34. public static String ter_w = "http://{server}.tianditu.gov.cn/ter_w/wmts?SERVICE=WMTS&REQUEST=GetTile&VERSION=1.0.0&LAYER=ter&STYLE=default&TILEMATRIXSET=w&FORMAT=tiles&TILEMATRIX={z}&TILEROW={y}&TILECOL={x}&tk={tk}";
    35. //地形注记 - 等经纬度
    36. public static String cta_c = "http://{server}.tianditu.gov.cn/cta_c/wmts?SERVICE=WMTS&REQUEST=GetTile&VERSION=1.0.0&LAYER=cta&STYLE=default&TILEMATRIXSET=c&FORMAT=tiles&TILEMATRIX={z}&TILEROW={y}&TILECOL={x}&tk={tk}";
    37. //地形注记 - 墨卡托
    38. public static String cta_w = "http://{server}.tianditu.gov.cn/cta_w/wmts?SERVICE=WMTS&REQUEST=GetTile&VERSION=1.0.0&LAYER=cta&STYLE=default&TILEMATRIXSET=w&FORMAT=tiles&TILEMATRIX={z}&TILEROW={y}&TILECOL={x}&tk={tk}";
    39. public static String[] servers = {"t0", "t1", "t2", "t3", "t4", "t5", "t6", "t7"};
    40. public static void main(String[] args){
    41. String basePath = "E:/天地图";
    42. //这里放你的天地图开发者秘钥,注意天地图API访问次数限制
    43. String tk = "75f0434f240669f4a2df635927XXXXXXX";
    44. String[] urlArr = {img_w};//要下载的图层
    45. int minZoom = 17;
    46. int maxZoom = 18;
    47. double startLat = 34.29;//开始纬度(从北到南)
    48. double endLat = 33.67;//结束纬度(从北到南)
    49. double startLon = 107.99;//开始经度(从西到东)
    50. double endLon = 109.14;//结束经度(从西到东)
    51. ExecutorService exe = Executors.newFixedThreadPool(6);
    52. //等经纬度第一层是1x2,纬度数量是2^0,经度数量是2^1
    53. //墨卡托投影第一层是2x2,纬度数量是2^1,经度数量是2^1
    54. for(int i=0; i<urlArr.length; i++){
    55. String url = urlArr[i].replace("{tk}", tk);
    56. System.out.println(url);
    57. String layerName = url.split("tianditu.gov.cn/")[1].split("/wmts?")[0];
    58. System.out.println(layerName);
    59. if(layerName.endsWith("c")){
    60. //等经纬度
    61. for(int z=minZoom; z<=maxZoom; z++){
    62. if(z > 17){//天地图10级以上只有中国区域
    63. startLat = startLat > 34.29 ? 34.29 : startLat;//开始纬度(从北到南)
    64. endLat = endLat < 33.67 ? 33.67 : endLat;//结束纬度(从北到南)
    65. startLon = startLon < 107.99 ? 107.99 : startLon;//开始经度(从西到东)
    66. endLon = endLon > 109.14 ? 109.14 : endLon;//结束经度(从西到东)
    67. }
    68. double deg = 360.0 / Math.pow(2, z) / 256;//一个像素点代表多少度
    69. int startX = (int)((startLon + 180) / deg / 256);//减数取整
    70. int endX = (int)((endLon + 180) / deg / 256);//加数取整
    71. int startY = (int)((90 - startLat) / deg / 256);
    72. int endY = (int)((90 - endLat) / deg / 256);
    73. for(int y=startY; y<=endY; y++){
    74. for(int x=startX; x<=endX; x++){
    75. final String newUrl = url.replace("{server}", servers[(int) (Math.random()*servers.length)]).replace("{z}", z+"").replace("{x}", x+"").replace("{y}", y+"");
    76. final String filePath = basePath + "/" + layerName + "/" + z + "/" + x + "/" + y + ".png";
    77. exe.execute(new Runnable() {
    78. @Override
    79. public void run() {
    80. File file = new File(filePath);
    81. if(!file.exists()){
    82. if(!file.getParentFile().exists()){
    83. file.getParentFile().mkdirs();
    84. }
    85. boolean loop = true;
    86. int count = 0;
    87. while(loop && count<5){//下载出错进行重试,最多5
    88. count++;
    89. try {
    90. InputStream in = getFileInputStream(newUrl);
    91. OutputStream out = new BufferedOutputStream(new FileOutputStream(file));
    92. byte[] b = new byte[8192];
    93. int len = 0;
    94. while((len = in.read(b)) > -1){
    95. out.write(b, 0, len);
    96. out.flush();
    97. }
    98. out.close();
    99. in.close();
    100. loop = false;
    101. } catch (Exception e) {
    102. loop = true;
    103. }
    104. }
    105. if(loop){
    106. System.out.println("下载失败:"+newUrl);
    107. }
    108. }
    109. }
    110. });
    111. }
    112. }
    113. }
    114. }else{
    115. //墨卡托
    116. /**
    117. * earthRadius = 6378137;//地球赤道半径6378137
    118. * aa = earthRadius * (Math.PI - 0);//赤道周长的一半,20037508.3427892
    119. * 85.05112877980659 = (Math.atan(Math.exp(aa / earthRadius))-Math.PI/4)*2 * 180 / Math.PI;//墨卡托最大有效纬度
    120. */
    121. if(startLat > 85.051128){
    122. System.out.println("墨卡托投影起始纬度最大为:85.051128,下载区域将被限定。");
    123. startLat = 85.051128;
    124. }
    125. if(endLat <- 85.051128){
    126. System.out.println("墨卡托投影终止纬度最小为:-85.051128,下载区域将被限定。");
    127. endLat = -85.051128;
    128. }
    129. for(int z=minZoom; z<=maxZoom; z++){
    130. if(z > 17){//天地图10级以上只有中国区域
    131. startLat = startLat > 34.29 ? 34.29 : startLat;//开始纬度(从北到南)
    132. endLat = endLat < 33.67 ? 33.67 : endLat;//结束纬度(从北到南)
    133. startLon = startLon < 107.99 ? 107.99 : startLon;//开始经度(从西到东)
    134. endLon = endLon > 109.14 ? 109.14 : endLon;//结束经度(从西到东)
    135. }
    136. double deg = 360.0 / Math.pow(2, z) / 256;//一个像素点代表多少度
    137. int startX = (int)((startLon + 180) / deg / 256);
    138. int endX = (int)((endLon + 180) / deg / 256);
    139. /**
    140. * 这里是基于等经纬度坐标转墨卡托坐标的公式得到的算法。()
    141. * 先算出当前地图级别下纬度方向用该有多少个像素点(一个瓦片是256个点,一共有Math.pow(2, z)个瓦片)
    142. * 再算出指定纬度坐标与赤道的距离(北半球为正南半球为负)(这里抵消了带入赤道周长的计算,因为它是一个固定值,把赤道周长定为1就抵消了):Math.log(Math.tan((90 + startLat) * Math.PI / 360)) / (Math.PI / 180)
    143. * 就可以得到纬度坐标距赤道有多少个点:(Math.log(Math.tan((90 + startLat) * Math.PI / 360)) / (Math.PI / 180)) / (360/Math.pow(2, z)/256) + 0.5)
    144. * 然后就可以计算出其瓦片坐标。
    145. */
    146. int startY = (((int)Math.pow(2, z) * 256 / 2) - (int)((Math.log(Math.tan((90 + startLat) * Math.PI / 360)) / (Math.PI / 180)) / (360/Math.pow(2, z)/256) + 0.5)) / 256;
    147. int endY = (((int)Math.pow(2, z) * 256 / 2) - (int)((Math.log(Math.tan((90 + endLat) * Math.PI / 360)) / (Math.PI / 180)) / (360/Math.pow(2, z)/256) + 0.5)) / 256;
    148. for(int y=startY; y<=endY; y++){//加入判断是等经纬度还是墨卡托
    149. for(int x=startX; x<=endX; x++){
    150. final String newUrl = url.replace("{server}", servers[(int) (Math.random()*servers.length)]).replace("{z}", z+"").replace("{y}", y+"").replace("{x}", x+"");
    151. //System.out.println(newUrl);
    152. final String filePath = basePath + "/" + layerName + "/" + z + "/" + y + "/" + x + ".png";
    153. exe.execute(new Runnable() {
    154. @Override
    155. public void run() {
    156. File file = new File(filePath);
    157. if(!file.exists()){
    158. if(!file.getParentFile().exists()){
    159. file.getParentFile().mkdirs();
    160. }
    161. boolean loop = true;
    162. int count = 0;
    163. while(loop && count<5){//下载出错进行重试,最多5
    164. count++;
    165. try {
    166. InputStream in = getFileInputStream(newUrl);
    167. OutputStream out = new BufferedOutputStream(new FileOutputStream(file));
    168. byte[] b = new byte[8192];
    169. int len = 0;
    170. while((len = in.read(b)) > -1){
    171. out.write(b, 0, len);
    172. out.flush();
    173. }
    174. out.close();
    175. in.close();
    176. loop = false;
    177. } catch (Exception e) {
    178. loop = true;
    179. }
    180. }
    181. if(loop){
    182. System.out.println("下载失败:"+newUrl);
    183. }
    184. }
    185. }
    186. });
    187. }
    188. }
    189. }
    190. }
    191. }
    192. exe.shutdown();
    193. while (true) {
    194. try {
    195. Thread.sleep(1000L);//主线程休眠1秒,等待线程池运行结束,同时避免一直死循环造成CPU浪费
    196. } catch (InterruptedException e) {
    197. }
    198. if (exe.isTerminated()) {//线程池所有线程都结束运行
    199. break;
    200. }
    201. }
    202. }
    203. //获取文件下载流
    204. public static InputStream getFileInputStream(String url) throws Exception{
    205. InputStream is = null;
    206. CloseableHttpClient httpclient = HttpClients.createDefault();
    207. HttpGet request = new HttpGet(url);
    208. request.setHeader("User-Agent", "Mozilla/4.0 (compatible; MSIE 5.0; Windows NT; DigExt)");
    209. HttpResponse response = httpclient.execute(request);
    210. response.setHeader("Content-Type", "application/octet-stream");
    211. int statusCode = response.getStatusLine().getStatusCode();
    212. if (statusCode == HttpStatus.SC_OK) {
    213. HttpEntity entity = response.getEntity();
    214. is = entity.getContent();
    215. }
    216. return is;
    217. }
    218. }

    这个亲测确实能够下载,但是存在下载的瓦片出现nginx发布后,引用实际cesium贴图位置错乱。应该是生成的XYZ目录存在一些问题。没有细究。

    方式3、通过QGIS下载(极力推荐

    QGIS下载地址:Download QGIS

    操作界面切换到中文界面如下:

    自定义下载瓦片区域,如下图第4步,会出现下拉选择列表,自行选择即可

    层级可以批处理分层下载结果如下图。

    注:根目录是层级Z,下一级目录是X,下下一级目录.png是Y

    2、Nginx发布离线瓦片数据

    nginx配置nginx.conf

    1. server {
    2. listen 9006;
    3. server_name 127.0.0.1;
    4. # 这里需要写相对路径 绝对路径会报错404
    5. root ../tianmap;
    6. autoindex on;
    7. location / {
    8. # 支持跨域
    9. add_header Access-Control-Allow-Origin *;
    10. add_header 'Access-Control-Allow-Methods' 'GET, POST, OPTIONS';
    11. add_header 'Access-Control-Allow-Headers' 'DNT,User-Agent,X-Requested-With,If-Modified-Since,Cache-Control,Content-Type,Range';
    12. add_header 'Access-Control-Expose-Headers' 'Content-Length,Content-Range';
    13. }
    14. # 不存缓存 每次重新验证
    15. add_header Cache-Control "no-cache,must-revalidate";
    16. }

    发布访问结果如下:

    3、cesium工程加载离线影像瓦片

    1. var MapModel = new Cesium.UrlTemplateImageryProvider({
    2. url:"http://127.0.0.1:9006/{z}/{x}/{y}.png",
    3. tilingScheme: new Cesium.WebMercatorTilingScheme(),
    4. fileExtension: 'png',
    5. minimumLevel: 1,
    6. maximumLevel: 18
    7. });
    8. viewer.imageryLayers.addImageryProvider(MapModel);

    到此结束

  • 相关阅读:
    Kubernetes面试题分享
    修改this.$notify通知的样式
    优化 C++ 字符串拼接:高效方法与代码示例
    CTFHUB ICS(3)
    【从零开始学习 SystemVerilog】5.3、SystemVerilog 通信—— Mailbox(邮箱)
    JOSEF信号继电器 JX-18A/2 电压 220VAC辅助电源 板后接线
    数字孪生技术助力城市轨道交通智慧升级
    Libuv源码解析 - uv_loop整个初始化模块
    使用 Learner Lab - 使用 API Gateway 触发 AWS Lambda
    操作系统的分类
  • 原文地址:https://blog.csdn.net/OuNuo5280/article/details/137857732