注:本文是我为某历史地图项目组撰写的扫描地图配准方案,现分享出来,希望能为大家提供帮助。阅读本文需要一定的地图学基础。

在学习和工作中,由于无法获取矢量数据,经常需要以图片格式地图为工作底图,并且需要对底图进行配准。这些地图包括早期纸质地图的扫描图、各种软件绘制的图片地图等,为方便表述本文统一简称为扫描地图

不同的扫描地图由于制图区域的不同,往往会根据特定的区域,选择最合适的地图投影,使得变形相对较小且分布较为均匀。在地图配准过程中,根据地图投影来选择合适的配准方法,可以提高配准的精度和效率。

1. 扫描地图中常用地图投影介绍

识别地图投影最主要的方法是根据经纬网的形态,如果没有经纬网,则只能通过和已知投影的地图进行对比来进行估计。

扫描地图中常用的地图投影有无投影坐标系(地理坐标系)、圆柱投影、圆锥投影、方位投影等四种。少数地图会使用伪圆柱投影、伪圆锥投影等,如国外世界地图常用的摩尔威德投影、罗宾森投影,国内亚洲地图常用的彭纳投影,国内世界地图常用的等差分纬线多圆锥投影,由于数量较少本文暂不涉及。

1. 无投影坐标系,即地理坐标系,并且几乎均WGS-84地理坐标系,该坐标系以实际的经纬度数值为坐标。WGS-84地理坐标系在历史地图十分常见,尤其是布哈林、史图馆、Worldhistorymaps等网络渠道发布的图片。WGS-84坐标系下,经线和纬线均为直线且互相垂直,等间隔的经纬线所形成的网格均为正方形。与其他投影相比,该投影下世界地图和中国地图均略扁一些。

图1 无投影坐标系

2.圆柱投影是一类投影的总称,简单来说,可以视为将地球表面投影到圆柱面上再展开为平面。扫描地图中一般只会用到墨卡托投影,又称正轴等角圆柱投影。目前,主流的各大电子地图均采用墨卡托投影(略有修改)作为投影坐标系。墨卡托投影下,经线和纬线均为直线且互相垂直,等间隔的经线之间距离相等,等间隔的纬线之间的距离随纬度增大而增大。与其他投影相比,该投影下世界地图和中国地图均略高一些。

图2 圆柱投影(墨卡托投影)

3.圆锥投影是一类投影的总称,与圆柱投影类似,可以视为将地球表面投影到圆锥面上再展开为平面。扫描地图中,以等面积圆锥投影为主。目前,圆锥投影多用于纸质出版的地图,尤其是分国图、区域图,全球地图极少使用该投影。圆锥投影下,经线为直线且从极点呈放射状,纬线为圆且互为同心圆,经纬线互相垂直(指在相交处正交)。与其他投影相比,该投影下的中国地图与实际最为接近。

图3 圆锥投影

4.方位投影也是一类投影的总称,可以视为将地球表面直接投影到平面上。方位投影与圆锥投影极为相似,在扫描地图中,也广泛地应用于区域图中。方位投影和圆锥投影最大也是最明显的区别在于,方位投影的经线多为弧线(极地除外)。方位投影下的中国地图也与实际十分接近。方位投影多用于极地地图、东西半球图,近年来出版的竖版中国地图也使用了方位投影。

图4 方位投影(斜轴等面积方位投影,中心点为110°E、30°N)

简单总结一下就是:

  • 无投影坐标系:经纬线互相垂直,经纬网呈正方形;
  • 墨卡托投影:经纬线互相垂直,经纬网呈矩形,随着纬度升高纬线间距越大;
  • 圆锥投影:经纬线互相垂直,经线为同一点放出的射线,纬线为圆弧且互为同心圆;
  • 方位投影:经纬网形态与圆锥投影相近,但经线多为弧线。

2. 不同地图投影下扫描地图的配准

WGS-84地理坐标系由于直接将经纬度作为地图坐标,虽然地图有一定的变形,但是依旧得到了广泛使用,适宜作为底图存储的坐标系。本文也将介绍将各类扫描地图配准到WGS-84地理坐标系下的方法,示例软件为ArcGIS 10.2。

2.1 WGS-84地理坐标系

WGS-84地理坐标系,理论上WGS-84地理坐标系下的扫描地图只需三个配准点即可完成配准,一般会在图片四个角上选取配准点。根据实际情况以及验证配准精度,还会适当地增加配准点的数量。

由于配准方法较为简单,具体案例不再赘述。

2.2 墨卡托投影

墨卡托投影的经线间隔保持均匀,但是纬线间隔随着纬度升高而增大,呈现不均匀的变化,如果直接配准,即使配准十几个点也会有极大的误差。因此,首先将工作空间的投影转换为墨卡托投影,再进行配准。具体步骤如下:

1)更改工作空间投影。在内容列表中,右键点击图层,在弹出框中选择属性。在数据框属性中,打开坐标系选项卡,在其中依次选择投影坐标系-World-Mercator(World),亦可在搜索框中直接输入名称进行搜索,然后点击确定。

如果待配准图片是世界地图,且中央经线(即地图位于正中央的经线)不是0°,则需要对投影进行修改。

在数据框属性-坐标系中,右键点击Mercator(World),选择复制并修改,在弹出的投影坐标系属性框中,将Central_Meridian修改为相应的数值,建议修改名称,然后点击确定。

图5.1 在坐标系中选择墨卡托投影
图5.2 修改墨卡托投影的参数

2)按照2.1的方案进行配准。此时,扫描地图和工作空间的投影坐标系已保持一致,参考2.1的流程进行配准即可。理论上只需三个点即可完成配准,可根据实际效果适当增加配准点,变换方式选择一阶多项式(仿射)即可。对于有明显扫描变形的地图,必须增加较多的配准点,变换方式可以选择二阶多项式,甚至三次多项式。完成后,将配准好的图片导出。

3)投影变换。步骤2配准好的图片是墨卡托投影下的,需要将其转换为WGS-84坐标系下的图片。打开工具箱,依次打开数据管理工具-投影和变换-栅格-投影栅格这一工具。在输入栅格中输入配准好的图片,输出栅格数据集中输入其保存路径,输出坐标系选择WGS-84坐标系(GCS_WGS_84),然后点击确定。

配准示例如下图所示:

6.1 原始扫描地图(墨卡托投影,中央经线为100°W)
图6.2 配准后的地图(WGS-84地理坐标系,中央经线为0°)

2.3 圆锥投影和方位投影

圆锥投影和方位投影其流程较为类似,因此统一进行说明。其基本步骤如下:

1)更改工作空间投影。圆锥投影的参数包括中央经线和标准纬线,方位投影的参数中央经线和原点纬线。圆锥投影一般基于Asia_North_Albers_Equal_Area_Conic进行修改,方位投影一般基于North_Pole_Lambert_Azimuthal_Equal_Area进行修改。中央经线的确定方式如下,一般将扫描地图中居于中部(不一定是正中央)且和图片底边垂直的经线视为中央经线,如果无恰好垂直的,可根据两侧的经线进行估算;而标准纬线和原点纬线则根据中央经线确定,一般将中央经线上,居于正中的纬线视为标准纬线或原点纬线。

在修改圆锥投影的过程中,投影参数会显示有两个标准纬线,均设置为根据扫描地图获得的标准纬线数值即可,同时圆锥投影也有原点纬线这一参数,如果不为0,请设置为0。

2)后续步骤参考2.2的步骤2和3即可。此处不再赘述。

配准示例如下图:

图7.1 原始地图(中高纬度、方位投影,估测其中心点为10°E、50°N)
图7.2 配准后的地图(中高纬度、WGS-84坐标系)
图8.1 原始地图(低纬度、方位投影,估测其中心点为20°E、0°)
图8.2 配准后的地图(低纬度、WGS-84坐标系)
图9.1 原始地图(中低纬度、圆锥投影,估测其中央经线为114°E,标准纬线为30°N)
图9.2 配准后的地图(中低纬度、WGS-84坐标系)
图10.1 原始地图(中高纬度、圆锥投影,估测其中央经线为102.5°E,标准纬线为50°N)
图10.2 配准后的地图(中高纬度、WGS-84坐标系)

根据经纬度坐标配准_扫描地图投影的识别及其配准相关推荐

  1. python点云快速配准_分享一个V-SLAM中点云配准算法改进的方法

    近年来,随着自主导航的广泛应用,视觉同时定位与地图构建(Visual simultaneous localization and mapping, V-SLAM)通过自身携带的视觉传感器对周围环境进行 ...

  2. python坐标定位_如何利用Python识别并定位图片中某一个色块的坐标?

    依赖python包|opencv.numpy.aircv 第一步:查找图片在原始图片上的坐标点 import aircv as ac def matchImg(imgsrc,imgobj,confid ...

  3. python坐标表示_已知经纬度坐标求两点间距离,用python表示

    已知经纬度坐标求两点间距离,用python表示 已知地球上任意两点(lon1, lat1),( lon2, lat2)的经纬度坐标,求两点间的距离用haversine表示: 1.首先先将经纬度坐标的角 ...

  4. 经纬度5位数和6位数差多少_各大在线地图,使用的经纬度坐标为什么都精确到小数点之后第六位...

    版权归作者所有,任何形式转载请联系作者. 作者:openGPS(来自豆瓣) 来源:https://www.douban.com/note/707875309/ 疑问由来 做过在线地图开发的同学,可能都 ...

  5. 知道经纬度坐标怎么计算两点间距离_【我的时间拣屎】亚里士多德:地球是圆的,我计算了地球的圆周...

    人们对宇宙的认识是一步步发展的,时间简史这本书最开始就是从宇宙图像的发展引入,今天这个故事就是关于我们对地球的认识. 人类对地球的认识是发展的,人类很早就认识到地球不是平的.这里的表述特意没有说地球是 ...

  6. 根据经纬度确定行政区域_基于JavaScript实现高德地图和百度地图提取行政区边界经纬度坐标...

    前言 近来由于工作需要,需要提取某些城市的经纬度坐标,稍微搜索了一下,发现百度地图和高德地图都提供了相关的函数和例子.那么剩余的工作也就比较简单了,保存坐标,然后转换为WGS坐标,这样才能和现有的GP ...

  7. python取省边界_提取行政区边界经纬度坐标(高德+百度)

    前言 近来由于工作需要,需要提取某些城市的经纬度坐标,稍微搜索了一下,发现百度地图和高德地图都提供了相关的函数和例子.那么剩余的工作也就比较简单了,保存坐标,然后转换为WGS坐标,这样才能和现有的GP ...

  8. 利用python查询电脑配置_干货|利用Python将地址转换为经纬度坐标

    本文主要讲述利用Python将文本格式的地址转换为数字格式的经纬度坐标数据,主要步骤有: 注册高德地图API账号,申请Web服务的Key 了解并测试地理编码API服务 利用Python实现地址转坐标 ...

  9. wgs84坐标格式转换度分秒_使用ArcGIS实现WGS84经纬度坐标到北京54高斯投影坐标的转换...

    [摘 要] 本文针对从事测绘工作者普遍遇到的坐标转换问题,简要介绍ArcGIS实现WGS84经纬度坐标到北京54高斯投影坐标转换原理和步骤. [关键词] ArcGIS 坐标转换 投影变换 1 坐标转换 ...

最新文章

  1. springboot tomcat配置_告诉你,Spring Boot 真是个牛逼货!
  2. html5中加亮文本,html实现高亮关键字
  3. 最全的jquery datatables api 使用详解
  4. 优先队列priority_queue的使用方式
  5. 纵观计算机网络发展历程,人工智能在计算机网络技术中的应用分析
  6. 大家都认为烧烤店是“暴利”行业
  7. tsp的理论和实践(10)tsp系列暂停一下
  8. log4j的详细介绍
  9. Redis源码分析系列十一:createClient后面内容
  10. VScode中文注释乱码问题解决
  11. 《全面软件质量管理》核心观点摘录
  12. 如何去掉图片上的文字?
  13. 打开程序员心理B面,这些黑红话题他们亲自回应丨1024特辑
  14. Android实现头像上传至数据库与保存 简易新闻(十七 下)
  15. Java | MySQL Emoji表情写入数据库时报错问题解决,Incorrect string value: ‘\xF0\x9F\x98\x84\xF0\x9F...‘ for column
  16. 树莓派+PHP+Mairadb数据库读取DHT11
  17. Java base64转化成 inputStream
  18. FAT32 学习及注意事项
  19. Android 来电秀总结
  20. 软件开发团队成员分工_分析软件开发人员的能力–选择合适的团队成员

热门文章

  1. 获取IE (控件)的所有链接(包括Frameset, iframe)zz
  2. UDLD是一个Cisco私有的二层协议
  3. 敬业福和花花卡算啥?这次不来,你亏了
  4. php [ext/gd/gd.lo] error 1,关于PHP一些编译安装出现的错误总结
  5. Perspective Mockups mac(PS透视模型动作插件)支持ps2021
  6. 22021年江苏高考成绩查询,江苏高考成绩查询系统入口
  7. 【chorme插件开发】第一节:chrome内核浏览器插件开发的说明
  8. editor.md国产markdown编辑器使用实例
  9. nginx下虚拟目录配置301域名重定向
  10. html语言table,html中的table详解