标准分幅下的图幅号转换成经纬度坐标【原理+源代码】

本文主要是介绍标准分幅下的图幅号转换成经纬度坐标【原理+源代码】,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

最近要批量的把标准分幅下的图幅号转换成经纬度坐标,所以这两天写了个程序来搞定这件事情。
先举个例子说明一下这个程序的作用。
例如:计算出图幅号I50G021040的经纬度范围,即最大经度、最小经度、最大纬度、最小纬度。
运用我编写的这个程序,可以直接算出来,这个图幅号的经纬度范围,最大经度为115.3125°,最小经度为115.25°,最大纬度为31.167°,最小纬度为31.125°。
说一下转换的原理吧。

原理

我国的基本比例尺有7种:1:1万,1:2.5万,1:5万,1:10万,1:25万,1:50万,1:100万。
根据国家标准规定,我国基本比例尺地形图均以1:100万地形图位基础,按规定的经差和纬差划分图幅。其中,1:100万地形图的分幅采用国际1:100万地图分幅标准。每幅1:100万地形图的范围是经差6°、纬差4°;纬度60°-76°为经差12°、纬差4°;纬度76°-88°之间经差24°、纬差4°。(中国在1:100万分幅中都是按经差6°、纬差4°分幅的)
100万地形图的分幅和编号(北半球)
1幅1:100万地形图可以划分成4(2*2)幅1:50万地形图、16(4*4)幅1:25万地形图、144(12*12)幅1:10万地形图、576(24*24)幅1:5万地形图、2304(48*48)幅1:2.5万地形图、9216(96*96)幅1:1万地形图、36864(192*192)幅1:0.5万地形图。
这里写图片描述
由于历史原因,我国地形图的编号在20世纪90年代以前很不统一。20世纪90年代以后,1:1 000 000-1:5000地形图的编号均以1:1 000 000地形图编号为基础,采用行列编号的方法。将1:1 000 000地形图按所含各比例尺地形图的经差和纬差划分成若干行和列,横行从上到下、纵列从左到右按顺序分别用3位阿拉伯数字(数字码)表示,不足3位者前面补零,取行号在前、列号在后的排列形式标记;各种比例尺地形图分别采用不同的字符作为其比例尺代码。
比例尺代码


接下来,把图幅号分解开来讲解。以上面提到的I50G021040图幅号为例。

  1. I,表示的是1:100万行号;50,表示的是1:100万列号。用于表示此范围在1:100万分幅地形图下的经纬度范围,即最大纬度是32°N,最小纬度是28°N,最大经度是120°,最小经度是114°。(查看前面的分幅图可知道:I,表示28°-32°N范围;50,表示114°-120°E)
  2. G,表示的是比例尺代码。用于记录地形图的实际比例尺。根据前面的比例尺代码表格可知,G表示1:10000比例尺,即1幅1:1 000 000地形图划分成96*96幅1:10000地形图。据此,我们可以算出1幅1:10000地形图的纬差为:6/96°,经差为:4/96°。
  3. 021,表示的是图幅行号数字码;040,表示的是图幅列号数字码。用于确定该地形图确切的经纬度范围。由于此图幅号为96*96,所以图幅行号和列号数字码都不能超过096。
    据此我们可以算出此图幅号的最大纬度、最小纬度、最大经度、最小经度。
    最大纬度=1:1 000 000地形图的最大纬度-(图幅行号数字码-1)*1:10000地形图的纬差=32°N - (21-1)*6/96°=31.1667°N
    最小纬度=最大纬度-纬度差=31.1667°N-6/96°=31.125°N
    最小经度=1:1 000 000地形图的最小经度+(图幅列数数字码-1)*1:10000地形图的经差=114°E+(40-1)*4/96°=115.25°E
    最大经度=最小经度+1:10000地形图的经差=115.25°E+4/96°=115.3125°E

程序代码

核心算法部分代码:

        #region

这篇关于标准分幅下的图幅号转换成经纬度坐标【原理+源代码】的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1089785

相关文章

Java五子棋之坐标校正

上篇针对了Java项目中的解构思维,在这篇内容中我们不妨从整体项目中拆解拿出一个非常重要的五子棋逻辑实现:坐标校正,我们如何使漫无目的鼠标点击变得有序化和可控化呢? 目录 一、从鼠标监听到获取坐标 1.MouseListener和MouseAdapter 2.mousePressed方法 二、坐标校正的具体实现方法 1.关于fillOval方法 2.坐标获取 3.坐标转换 4.坐

计算绕原点旋转某角度后的点的坐标

问题: A点(x, y)按顺时针旋转 theta 角度后点的坐标为A1点(x1,y1)  ,求x1 y1坐标用(x,y)和 theta 来表示 方法一: 设 OA 向量和x轴的角度为 alpha , 那么顺时针转过 theta后 ,OA1 向量和x轴的角度为 (alpha - theta) 。 使用圆的参数方程来表示点坐标。A的坐标可以表示为: \[\left\{ {\begin{ar

数据库原理与安全复习笔记(未完待续)

1 概念 产生与发展:人工管理阶段 → \to → 文件系统阶段 → \to → 数据库系统阶段。 数据库系统特点:数据的管理者(DBMS);数据结构化;数据共享性高,冗余度低,易于扩充;数据独立性高。DBMS 对数据的控制功能:数据的安全性保护;数据的完整性检查;并发控制;数据库恢复。 数据库技术研究领域:数据库管理系统软件的研发;数据库设计;数据库理论。数据模型要素 数据结构:描述数据库

计算机组成原理——RECORD

第一章 概论 1.固件  将部分操作系统固化——即把软件永恒存于只读存储器中。 2.多级层次结构的计算机系统 3.冯*诺依曼计算机的特点 4.现代计算机的组成:CPU、I/O设备、主存储器(MM) 5.细化的计算机组成框图 6.指令操作的三个阶段:取指、分析、执行 第二章 计算机的发展 1.第一台由电子管组成的电子数字积分和计算机(ENIAC) 第三章 系统总线

GaussDB关键技术原理:高性能(二)

GaussDB关键技术原理:高性能(一)从数据库性能优化系统概述对GaussDB的高性能技术进行了解读,本篇将从查询处理综述方面继续分享GaussDB的高性能技术的精彩内容。 2 查询处理综述 内容概要:本章节介绍查询端到端处理的执行流程,首先让读者对查询在数据库内部如何执行有一个初步的认识,充分理解查询处理各阶段主要瓶颈点以及对应的解决方案,本章以GaussDB为例讲解查询执行的几个主要阶段

陀螺仪LSM6DSV16X与AI集成(8)----MotionFX库解析空间坐标

陀螺仪LSM6DSV16X与AI集成.8--MotionFX库解析空间坐标 概述视频教学样品申请源码下载开启CRC串口设置开启X-CUBE-MEMS1设置加速度和角速度量程速率选择设置FIFO速率设置FIFO时间戳批处理速率配置过滤链初始化定义MotionFX文件卡尔曼滤波算法主程序执行流程lsm6dsv16x_motion_fx_determin欧拉角简介演示 概述 本文将探讨

【计算机组成原理】部分题目汇总

计算机组成原理 部分题目汇总 一. 简答题 RISC和CICS 简要说明,比较异同 RISC(精简指令集)注重简单快速的指令执行,使用少量通用寄存器,固定长度指令,优化硬件性能,依赖软件(如编译器)来提升效率。 CISC(复杂指令集)包含多样复杂的指令,能一条指令完成多步操作,采用变长指令,减少指令数但可能增加执行时间,倾向于硬件直接支持复杂功能减轻软件负担。 两者均追求高性能,但RISC

MySQL数据库锁的实现原理

MySQL数据库的锁实现原理主要涉及到如何确保在多用户并发访问数据库时,保证数据的完整性和一致性。以下是MySQL数据库锁实现原理的详细解释: 锁的基本概念和目的 锁的概念:在数据库中,锁是用于管理对公共资源的并发控制的机制。当多个用户或事务试图同时访问或修改同一数据时,数据库系统通过加锁来确保数据的一致性和完整性。 锁的目的:解决多用户环境下保证数据库完整性和一致性的问题。在并发的情况下,会

OSG数学基础:坐标系统

坐标系是一个精确定位对象位置的框架,所有的图形变换都是基于一定的坐标系进行的。三维坐标系总体上可以分为两大类:左手坐标系和右手坐标系。常用的坐标系:世界坐标系、物体坐标系和摄像机坐标系。 世界坐标系 世界坐标系是一个特殊的坐标系,它建立了描述其他坐标系所需要的参考框架。从另一方面说,能够用世界坐标系来描述其他坐标系的位置,而不能用更大的、外部的坐标系来描述世界坐标系。世界坐标系也被广泛地

线性回归(Linear Regression)原理详解及Python代码示例

一、线性回归原理详解         线性回归是一种基本的统计方法,用于预测因变量(目标变量)与一个或多个自变量(特征变量)之间的线性关系。线性回归模型通过拟合一条直线(在多变量情况下是一条超平面)来最小化预测值与真实值之间的误差。 1. 线性回归模型         对于单变量线性回归,模型的表达式为:         其中: y是目标变量。x是特征变量。β0是截距项(偏置)。β1