摘要：如需转载，请在显著位置注明个人微信公众号为什么要做地图投影简而言之，地球表面是一个三维的曲面，在曲面上进行测量是非常困难的。拓展同一个世界，不同的面孔链接在不同投影下的这个世界。。。
(如需转载，请在显著位置注明个人微信公众号stdrei)
为什么要做地图投影简而言之，地球表面是一个三维的曲面，在曲面上进行测量是非常困难的。不信你拿个地球仪量一下两点的距离或者计算个夹角试试。将三维的曲面投影到二维平面，这样我们学的平面几何才有用武之地。
什么是投影如果要了解地图投影，首先需要知道地理坐标系和投影坐标系。简单的说，地理坐标系是参考平面为椭球面的球面坐标，最常见的是以经纬度来量算的球面坐标系统。投影坐标系是定义在一个二维平面的坐标系统，其地图单位通常为米。将球面坐标转化为平面坐标的过程便称为投影。 因此，投影坐标系总是基于地理坐标系的（连坐标系都要讲出身的）。
因为地球是一个赤道略宽两极略扁的不规则的梨形球体，其表面是一个不可展平的曲面，所以运用任何数学方法进行这种投影转换都会产生误差和变形，为按照不同的需求缩小误差，就产生了各种投影方法。比如这篇很有意思的文章你所不知的有趣投影方法。也可以看文末的几幅地球不同投影的图。
EPSG Code如果经常跟地图打交道，你一定会被各种花式投影变换弄的眼花缭乱。EPSG：European Petroleum Survey Group (EPSG)成立于1986年，并在2005年重组为OGP(Internation Association of Oil & Gas Producers)，它负责维护并发布坐标参照系统的数据集参数，以及坐标转换描述，该数据集被广泛接受并使用，通过一个Web发布平台进行分发，同时提供了微软Acess数据库的存储文件，通过SQL 脚本文件，MySQL, Oracle 和PostgreSQL等数据库也可使用。目前已有的椭球体，投影坐标系等不同组合都对应着不同的ID号，这个号在EPSG中被称为EPSG code，它代表特定的椭球体、单位、地理坐标系或投影坐标系等信息。
开源的QGIS软件中就直接采用了EPSG。
The CRSs available in QGIS are based on those defined by the European Petroleum Search Group (EPSG) and the Institut Geographique National de France (IGNF) and are largely abstracted from the spatial reference tables used in GDAL. EPSG identifiers are present in the database and can be used to specify a CRS in QGIS.
一个查询EPSG code的网站 epsg.io
pyproj小试牛刀pyproj是PROJ4的python接口封装，直接看一个官网的例子吧。直接利用epsg code来定义投影参数。
from pyproj import Proj,transform # The Proj class can convert from geographic (longitude,latitude) # to native map projection (x,y) coordinates and vice versa, # or from one map projection coordinate system directly to another. p1 = Proj(init="epsg:26915") p2 = Proj(init="epsg:26715") x1, y1 = p1(-92.199881,38.56694) x2, y2 = transform(p1,p2,x1,y1) print "%9.3f %11.3f" % (x1,y1) print "%9.3f %11.3f" % (x2,y2)
输出为
569704.566 4269024.671 569722.342 4268814.027基于geopandas的矢量地图投影
import shapely, geopandas, fiona import seaborn as sns from fiona.crs import from_epsg,from_string # Data shp = "E:NationalGISdataProvince.shp" shp_df = geopandas.GeoDataFrame.from_file(shp) # #IndexError报错的话，用arcgis将shapefile文件重新导出一遍试试 shp_df.head()
# 根据当前的兰伯特投影绘制 shp_df.plot(column="GDP_1994",colormap="Set1")
# 转换到经纬度坐标 shp_df_wgs84 = shp_df.to_crs(from_epsg(4326)) shp_df_wgs84.plot(column="GDP_1994",colormap="Set1")
# 国家2000坐标系 # EPSG:4508 CGCS2000 / Gauss-Kruger CM 111E shp_df_4508 = shp_df.to_crs(from_epsg(4508)) shp_df_4508.plot(column="GDP_1994",colormap="Set1")
# 除了直接用ESPG code，也可以自己定义投影参数 ESRI_54024 = """ +proj=bonne +lon_0=0 +lat_1=60 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs shp_df_3408 = shp_df.to_crs(from_string(ESRI_54024)) shp_df_3408.plot(column="GDP_1994",colormap="Set1")代码和矢量地图数据下载
百度网盘下载地址： http://pan.baidu.com/s/1c1BExH2
密码请在关注微信公众号stdrei后，输入projection直接获取。
链接 http://www.viewsoftheworld.ne...
在不同投影下的这个世界。。。
文章版权归作者所有，未经允许请勿转载,若此文章存在违规行为，您可以联系管理员删除。
转载请注明本文地址：https://www.ucloud.cn/yun/38272.html
摘要：有许多地图的示例，包括地铁图，基于数据的地图，结合的地图等，每个示例都是单独的实现，代码也各有不同，于是我们想，是否能将这些方案统一起来，实现一种地图解决方案呢于是我们开始整合，从世界地图到国内各市县的地图，再到各种地铁图，都用技术呈现出来 Qunee for HTML5有许多地图的示例，包括地铁图，基于SVG数据的地图，结合leaflet的地图等，每个示例都是单独的实现，代码也各有...
摘要：它的处理单位是顶点，也就是每个顶点，都会调用一次顶点着色器。这里给一个的代码，把顶点着色器和片元着色器的代码，放到了一个文件中，不过引擎会解析文件，转换个着色器代码段，传递给。通常，顶点着色器只处理顶点坐标，进行空间转换。 ...
摘要：最后，为了使规划模块获得更平滑的姿态，提高定位系统的鲁棒性，采用了带有滑动窗口的姿态图，优化窗口中包含跟踪良好的帧数据，如果滑动窗口的大小超过阈值，历史记录中的一帧将根据车辆状态从滑动窗口中剔除。 点云PCL免费知识星球，点云论文速读。 文章：Coarse-to-fine Semanti...
摘要：代表第一章下的代表各章下的小细节内容。描述了李群的局部性质，是单位元附近的正交空间。否则，一个群可能存在多个李代数注意点旋转矩阵的导数，可以由旋转向量指定。误差项是将点的投影位置与观测位置做差，所以称为重投影误差。 1）代表第一章 1）下的1.2代表各章下的小细节内容。层次按1. ——1.1...
摘要：本篇文章已授权微信公众号郭霖独家发布老规矩先上图最近没有什么时间，后面项目再补上详细说明百度地图新增点聚合功能。百度地图是把整个地球是按照一个平面来展开，并且通过墨卡托投影投射到坐标轴上面。上图很明显墨卡托投影把整张世界地图投影成。 本篇文章已授权微信公众号 guolin_blog （郭霖）独家发布 老规矩先上图最近 没有什么时间，后面项目再补上详细说明 showImg(https:/...
阅读 1481·2021-11-24 09:39
阅读 3376·2021-11-24 09:39
阅读 2133·2021-11-15 11:37
阅读 1944·2021-11-11 16:55
阅读 4362·2021-10-14 09:43
阅读 3518·2021-10-08 10:05
阅读 2782·2021-09-13 10:26
阅读 2114·2021-09-08 09:35
