基于JS实现《国家基本比例尺地形图分幅和编号》标准

已于 2024-06-06 11:15:52 修改
阅读量212 收藏 2
点赞数 1
分类专栏： Cesium 文章标签： javascript 前端 服务器
于 2024-06-05 17:10:57 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/anglestar2012/article/details/139472803
版权

1、标准

GB T 13989-2012国家基本比例尺地形图分幅和编号

地址：【高清版】GB T 13989-2012国家基本比例尺地形图分幅和编号 - 道客巴巴

2、1:100万比例尺

2.1 说明

2.2 计算公式

2.3 计算代码

2.3.1 元素数据定义

由于中国只到N层，所以只定义到O.

/**
   * 100万比例尺下的行号对应的代码
   */
  static mateData100 = [
    {
      rowCode: "A",
      rowNumber: 1,
    },
    {
      rowCode: "B",
      rowNumber: 2,
    },
    {
      rowCode: "C",
      rowNumber: 3,
    },
    {
      rowCode: "D",
      rowNumber: 4,
    },
    {
      rowCode: "E",
      rowNumber: 5,
    },
    {
      rowCode: "F",
      rowNumber: 6,
    },
    {
      rowCode: "G",
      rowNumber: 7,
    },
    {
      rowCode: "H",
      rowNumber: 8,
    },
    {
      rowCode: "I",
      rowNumber: 9,
    },
    {
      rowCode: "J",
      rowNumber: 10,
    },
    {
      rowCode: "K",
      rowNumber: 11,
    },
    {
      rowCode: "L",
      rowNumber: 12,
    },
    {
      rowCode: "M",
      rowNumber: 13,
    },
    {
      rowCode: "N",
      rowNumber: 14,
    },
    {
      rowCode: "O",
      rowNumber: 15,
    },
  ];

2.3.2 代码实现

/**
   * 根据经纬度计算100万比例尺的行列号
   * @param lon
   * @param lat
   * @return {{col: number, row: number, code: string}}
   */
  static calcMapFrameNumber100 = (lon, lat) => {
    let a = Math.trunc(lat / 4.0) + 1;
    let b = Math.trunc(lon / 6.0) + 31;
    let mateInfo = MapFrameUtils.mateData100.find(
      (item) => item.rowNumber === a
    );
    return {
      row: a,
      col: b,
      code: mateInfo.rowCode,
    };
  };

3、非100万比例尺

3.1 说明

3.1.1 2000到50万比例尺格式

3.1.2 500到1000比例尺格式

3.2 计算公式

3.3 代码

3.3.1 元数据

/**
   * 100万比例尺下的行号对应的代码
   */
  static mateData100 = [
    {
      rowCode: "A",
      rowNumber: 1,
    },
    {
      rowCode: "B",
      rowNumber: 2,
    },
    {
      rowCode: "C",
      rowNumber: 3,
    },
    {
      rowCode: "D",
      rowNumber: 4,
    },
    {
      rowCode: "E",
      rowNumber: 5,
    },
    {
      rowCode: "F",
      rowNumber: 6,
    },
    {
      rowCode: "G",
      rowNumber: 7,
    },
    {
      rowCode: "H",
      rowNumber: 8,
    },
    {
      rowCode: "I",
      rowNumber: 9,
    },
    {
      rowCode: "J",
      rowNumber: 10,
    },
    {
      rowCode: "K",
      rowNumber: 11,
    },
    {
      rowCode: "L",
      rowNumber: 12,
    },
    {
      rowCode: "M",
      rowNumber: 13,
    },
    {
      rowCode: "N",
      rowNumber: 14,
    },
    {
      rowCode: "O",
      rowNumber: 15,
    },
  ];

3.3.2 实现代码

/**
   * 根据比例尺，经纬度计算地图的行列号
   * @param scale
   * @param lon
   * @param lat
   * @return {{col: number, code: *, row: number, digit: (number|*)}|null}
   */
  static calcMapFrameNumberNot100 = (scale, lon, lat) => {
    let mateData = MapFrameUtils.mateData.find((item) => item.scale === scale);
    if (mateData) {
      let cStep = ((lat % 4) / mateData.latD).toPrecision(4);
      let c = Math.round(4 / mateData.latD - Math.trunc(Number(cStep)));
      let dStep = ((lon % 6) / mateData.lonD).toPrecision(4);
      let d = Math.trunc(Number(dStep)) + 1;
      return {
        row: c,
        col: d,
        code: mateData.code,
        digit: mateData.digit,
      };
    } else {
      return null;
    }
  };

4、综合代码

/**
 * 采用标准：国家基本比例尺地形图分幅和编号《GB/T 13989-2012》
 * 1、根据经纬度计算图幅号
 * 2、根据图幅号图幅范围计算
 * 3、定义图幅信息
 */
export default class MapFrameUtils {
  /**
   * 根据经纬度计算100万比例尺的行列号
   * @param lon
   * @param lat
   * @return {{col: number, row: number, code: string}}
   */
  static calcMapFrameNumber100 = (lon, lat) => {
    let a = Math.trunc(lat / 4.0) + 1;
    let b = Math.trunc(lon / 6.0) + 31;
    let mateInfo = MapFrameUtils.mateData100.find(
      (item) => item.rowNumber === a
    );
    return {
      row: a,
      col: b,
      code: mateInfo.rowCode,
    };
  };
  /**
   * 根据比例尺，经纬度计算地图的行列号
   * @param scale
   * @param lon
   * @param lat
   * @return {{col: number, code: *, row: number, digit: (number|*)}|null}
   */
  static calcMapFrameNumberNot100 = (scale, lon, lat) => {
    let mateData = MapFrameUtils.mateData.find((item) => item.scale === scale);
    if (mateData) {
      let cStep = ((lat % 4) / mateData.latD).toPrecision(4);
      let c = Math.round(4 / mateData.latD - Math.trunc(Number(cStep)));
      let dStep = ((lon % 6) / mateData.lonD).toPrecision(4);
      let d = Math.trunc(Number(dStep)) + 1;
      return {
        row: c,
        col: d,
        code: mateData.code,
        digit: mateData.digit,
      };
    } else {
      return null;
    }
  };
  /**
   * 根据计算出来的行列号，比例尺，计算该格网在地图商的范围
   * @param frame100
   * @param frameNot100
   * @param scale
   * @return {{minY: number, minX: number, maxY: number, maxX: number}}
   */
  static calcMapFrameExtent = (frame100, frameNot100, scale) => {
    let mateData = MapFrameUtils.mateData.find((item) => item.scale === scale);
    let lon = (frame100.col - 31) * 6 + (frameNot100.col - 1) * mateData.lonD;
    let lat =
      (frame100.row - 1) * 4 +
      (4 / mateData.latD - frameNot100.row) * mateData.latD;

    return {
      minX: lon,
      minY: lat,
      maxX: lon + mateData.lonD,
      maxY: lat + mateData.latD,
    };
  };
  /**
   * 根据比例尺和经纬度，获取所在图幅范围的编号、行列号以及范围
   * @param scale
   * @param lon
   * @param lat
   * @return {null|{extent: {minY: number, minX: number, maxY: number, maxX: number}, col: number, code: string, row: number}}
   */
  static calcMapFrameNumber = (scale, lon, lat) => {
    let mapFrameNumber100 = MapFrameUtils.calcMapFrameNumber100(lon, lat);
    let mapFrameNumberNot100 = MapFrameUtils.calcMapFrameNumberNot100(
      scale,
      lon,
      lat
    );
    if (mapFrameNumberNot100) {
      let code =
        mapFrameNumber100.code +
        "" +
        mapFrameNumber100.col +
        "" +
        mapFrameNumberNot100.code +
        "" +
        mapFrameNumberNot100.row
          .toString()
          .padStart(mapFrameNumberNot100.digit, "0") +
        "" +
        mapFrameNumberNot100.col
          .toString()
          .padStart(mapFrameNumberNot100.digit, "0");
      let extent = MapFrameUtils.calcMapFrameExtent(
        mapFrameNumber100,
        mapFrameNumberNot100,
        scale
      );
      return {
        code: code,
        extent: extent,
        row: mapFrameNumberNot100.row,
        col: mapFrameNumberNot100.col,
      };
    } else {
      return null;
    }
  };
  /**
   * 100万比例尺下的行号对应的代码
   */
  static mateData100 = [
    {
      rowCode: "A",
      rowNumber: 1,
    },
    {
      rowCode: "B",
      rowNumber: 2,
    },
    {
      rowCode: "C",
      rowNumber: 3,
    },
    {
      rowCode: "D",
      rowNumber: 4,
    },
    {
      rowCode: "E",
      rowNumber: 5,
    },
    {
      rowCode: "F",
      rowNumber: 6,
    },
    {
      rowCode: "G",
      rowNumber: 7,
    },
    {
      rowCode: "H",
      rowNumber: 8,
    },
    {
      rowCode: "I",
      rowNumber: 9,
    },
    {
      rowCode: "J",
      rowNumber: 10,
    },
    {
      rowCode: "K",
      rowNumber: 11,
    },
    {
      rowCode: "L",
      rowNumber: 12,
    },
    {
      rowCode: "M",
      rowNumber: 13,
    },
    {
      rowCode: "N",
      rowNumber: 14,
    },
    {
      rowCode: "O",
      rowNumber: 15,
    },
  ];
  /**
   * 非100万情况下，不同比例尺的原始数据
   */
  static mateData = [
    {
      scale: 500000,
      label: "1:50万",
      lonD: 3,
      latD: 2,
      code: "B",
      digit: 3,
    },
    {
      scale: 250000,
      label: "1:25万",
      lonD: 1.5,
      latD: 1,
      code: "C",
      digit: 3,
    },
    {
      scale: 100000,
      label: "1:10万",
      lonD: 0.5,
      latD: 0.33333333,
      code: "D",
      digit: 3,
    },
    {
      scale: 50000,
      label: "1:5万",
      lonD: 0.25,
      latD: 0.16666667,
      code: "E",
      digit: 3,
    },
    {
      scale: 25000,
      label: "1:2.5万",
      lonD: 0.125,
      latD: 0.08333333,
      code: "F",
      digit: 3,
    },
    {
      scale: 10000,
      label: "1:1万",
      lonD: 0.0625,
      latD: 0.04166667,
      code: "G",
      digit: 3,
    },
    {
      scale: 5000,
      label: "1:5000",
      lonD: 0.03125,
      latD: 0.02083333,
      code: "H",
      digit: 3,
    },
    {
      scale: 2000,
      label: "1:2000",
      lonD: 0.01041667,
      latD: 0.006944444,
      code: "I",
      digit: 3,
    },
    {
      scale: 1000,
      label: "1:1000",
      lonD: 0.005208333,
      latD: 0.003472222,
      code: "J",
      digit: 4,
    },
    {
      scale: 500,
      label: "1:500",
      lonD: 0.002603889,
      latD: 0.001736111,
      code: "K",
      digit: 4,
    },
  ];
}

5、测试结果

5.1 测试代码

MapFrameUtils.mateData.forEach((item) => {
  let mapFrame = MapFrameUtils.calcMapFrameNumber(item.scale, 114.5625, 39.375);
  console.log(item.scale, mapFrame);
});

5.2 测试结果

与标准中结果一致

生活真难
关注
专栏目录
国家基本比例尺地形图分幅和编号实现
xyt1172228177的博客
12-07 1918
国家基本比例尺地形图分幅和编号，通过经纬度坐标查询对应图幅号并生成图幅四至。
GIS开发入门，OpenLayers、Leaflet、Maplibre-gl和Cesiumjs地图引擎介绍以及几种地图服务vms、vmts、TMS和XYZ介绍
热门推荐
GIS地图技术分享(GIS入门、Openlayers教程、Leaflet教程)，做最好的GIS地图开发教程
04-10 1万+
OpenLayer、Leaflet、Maplibre-gl和Cesiumjs（3d）OpenLayer 对各种地图服务支持较好，比如各种vms,vmts等，动态效果较差。Leaflet 对vue支持较好，与前端交互方便，对地图服务支持不好Maplibre-gl:继承自Mapbox-2d的开源分支，Mapbox商业化后Maplibre-gl项目从Mapbox-2d的开源分支开始继续维护更新，号称地图服务中动态效果最好的库，插件较多，个人感觉像是地图引擎里的游戏引擎，各种动态效果极好。
新图幅号计算和图幅框计算
LongJ_Sir的博客
09-14 1642
本文的图幅号计算直接引用了博客园地址：使用Java实现的新版地图图幅编号计算公式（新版图幅规则(GB/T 13989—2012)） - JuCheap - 博客园 的代码，微小的调整了一下，计算图幅框四至是参照 地图学第02章：地图的分幅和编号来实现，地图学第02章 地图的分幅和编号 （新） - 道客巴巴 /** * <p> * 地图坐标对象 * </p> * * @Author jucheap * @Email jucheap@jucheap.com */ p.
arcgis for js4二维地图比例尺
qq_37320631的博客
04-15 602
arcgis for js 4
根据经纬度求算所在的地形图分幅编号
weixin_44605704的博客
07-31 1905
【代码】根据经纬度求算所在的地形图分幅编号。
uniapp 地图分幅网格生成 小程序基于map组件
weixin_44309299的博客
04-09 423
// 获取小数部分 const fractional = function(x) { x = Math.abs(x); return x - Math.floor(x); } const formatInt = function(x, len) { let result = '' + x; len = len - result.length; while (len > 0) { result = '0' + result; len--; } return result; }
OGC标准WMTS服务概念与地图商的瓦片编号流派-web地图切片加载
周陆军的博客，分享it技术。
07-26 1897
还不知道地图栅格化切片等相关GIS原理的，推荐阅读《webGIS底图栅格化与实时数据合成处理原理,地图API设计,xyz加载》 OGC概念 OGC全称——开放地理空间信息联盟(Open Geospatial Consortium), 它的主要目的就是制定与空间信息、基于位置服务相关的标准。而这些所谓的标准其实就是一些接口或编码的技术文档，不同的厂商、各种GIS产品都可以对照这些文档来定义开放服务的接口、空间数据存储的编码、空间操作的方法。 OGC已经是一个比较“官方”的标准化机构了，它不但包括了..
我国常用的地理和投影坐标系的WKID.docx
11-10
- **应用**：适用于中国全境的地形图和大比例尺地图制作。其中，数字表示不同的投影带编号，如“Xian_1980_GK_Zone_13”代表第13号投影带。 2. **Xian_1980_3_Degree_GK_Zone_25 至 Xian_1980_3_Degree_GK_Zone_...
地图瓦片资源下载器使用
05-09
每个瓦片覆盖一定的地理范围，具有固定的比例尺和分辨率。当用户在浏览器中浏览地图时，地图服务会根据用户的缩放级别和位置动态加载相应的瓦片，以形成连续且平滑的视觉效果。常见的瓦片格式有JPEG、PNG等，它们...
通过CAD软件绘制土木工程的道路平面图
# 1. CAD软件简介 ## 1.1 CAD软件的定义和作用 在土木工程中，CAD（计算机辅助设计）软件是一类...CAD软件可以帮助工程师创建精确的工程图纸，实现设计方案的可视化和优化。 ## 1.3 适用于绘制道路平面图的常见CAD
火星坐标系转北京54坐标方法
dernan Thinker
05-29 6335
实现火星坐标系转北京54 1. App获得高德地图坐标系下的经纬度，请求参数给服务器； 2. 服务器解析Json数据，将高德地图坐标的点经纬度经过算法处理成WGS84坐标系下的点坐标； 3. WGS84坐标系下的点坐标经ArcGIS REST API调用投影处理成北京54坐标系； 4. 处理偏移问题，进行identify查询，获取土壤养分数据，返回App； 5. App获得数据，显示到用户界面。
国家基本比例尺地形图分幅和编号
阿童林的专栏
02-19 1790
原文首发地址：http://www.atolin.net/285/ 国家基本比例尺地形图分幅和编号 1、定义详见标准文档:《国家基本比例尺地形图分幅和编号》 2、设计及实现: 国家基本比例尺地形图分幅和编号功能设计实现 目录及文件概述: MapSubdivision——》国家基本比例尺地形图分幅和编号功能设计实现； PyMesh——》Python接口封装； Te
【Vue】a-radio-group单选框数字与字符串回显问题
最新发布
叶为正的博客
10-08 130
根据后端接口返回数字或字符串，来配置Ant Design of Vue中a-radio-group单选框的回显问题，注意value中的取值是不一样的。
JavaScript ArrayBuffer的读写与类型转换
白瑕 的博客
10-01 495
所有视图的基本单位是ArrayBuffer, ArrayBuffer中只存储二进制格式的数据.ArrayBuffer不可直接读写, 必须通过视图(比如DataView)暴露的API.
vue-live2d看板娘集成方案设计使用教程
每个人都是独一无二的，把握好自己的节奏，跟着自己的心走。
10-05 1279
vue-live2d看板娘集成方案设计使用教程
D3.js(v3) + React 教程：构建带坐标轴和比例尺的散点图
"这篇文章主要介绍了如何使用D3.js(v3)和React来创建一个带有坐标轴和比例尺的散点图。散点图是一种用于显示数据分布和关系的图表，通常包含水平和垂直两个坐标轴。本文提供的示例数据是一组二维坐标，用于在图表上...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

生活真难

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值