首页
学习
活动
专区
工具
TVP
发布
精选内容/技术社群/优惠产品,尽在小程序
立即前往

如何在使用极线(epsg:3995)投影和ggspatial时在ggplot中的66度平行北线上绘制一条线

在使用极线(epsg:3995)投影和ggspatial时,在ggplot中的66度平行北线上绘制一条线,可以按照以下步骤进行操作:

  1. 首先,确保已经安装了必要的R包,包括ggplot2、ggspatial和sf。可以使用以下命令安装这些包:
代码语言:txt
复制
install.packages("ggplot2")
install.packages("ggspatial")
install.packages("sf")
  1. 导入所需的库:
代码语言:txt
复制
library(ggplot2)
library(ggspatial)
library(sf)
  1. 创建一个包含线的数据框,其中包含线的起点和终点的经纬度坐标。假设起点坐标为(lon1, lat1),终点坐标为(lon2, lat2):
代码语言:txt
复制
data <- data.frame(lon = c(lon1, lon2), lat = c(lat1, lat2))
  1. 将数据框转换为sf对象,并设置坐标参考系统为极线(epsg:3995)投影:
代码语言:txt
复制
sf_data <- st_as_sf(data, coords = c("lon", "lat"), crs = 3995)
  1. 创建一个ggplot对象,并使用geom_sf函数将sf_data添加到图表中。同时,使用coord_sf函数设置坐标参考系统为极线(epsg:3995)投影,并使用xlim和ylim函数设置绘图范围:
代码语言:txt
复制
ggplot() +
  geom_sf(data = sf_data) +
  coord_sf(crs = 3995) +
  xlim(xmin, xmax) +
  ylim(ymin, ymax)

其中,xmin、xmax、ymin和ymax是绘图范围的边界值。

这样,就可以在使用极线(epsg:3995)投影和ggspatial时,在ggplot中的66度平行北线上绘制一条线。请注意,这里的代码示例中没有提及具体的腾讯云产品,因为与问题无关。如需了解腾讯云相关产品和产品介绍,请访问腾讯云官方网站。

页面内容是否对你有帮助?
有帮助
没帮助

相关·内容

R-ggspatial+ggplot2 实现带指北针比例尺空间地图绘制

继上次使用tmap包制作含有指北针(compass)比例尺(scale bar)以及南海小地图添加,详细内容分别见如下两篇文章:R-tmap 绘制带指北针比例尺空间地图 R-tmap+grid...虽然tmap包有着类似于ggplot2绘图语法,但对习惯使用ggplot2绘图小伙伴怎不怎么友好。那么今天本期推文就使用 ggspatial 包实现指北针比例尺添加。...主要涉及知识点如下: sf包shp文件读取及转换投影 ggplot2 + ggspatial 实现指北针比例尺添加 sf 读取shp文件并转换投影 我想经常使用R处理空间数据小伙伴们对sf包一定特别熟悉...当然指北针比例尺我们都是可以选择样式或者定制,详细内容可以查看网ggspatial官网 投影转换-高斯-克吕格投影 绘制中国地图,我们一般使用高斯-克吕格投影,具体原因如下:高斯-克吕格投影没有角度变形...总结 本期推文我们使用ggspatial包、ggplot2包以及sf包geom_sf()完成另一种指北针比例尺地图绘制,希望给大家一些绘制灵感。

2.4K20

Python地信专题 | 基于geopandas空间数据分析-坐标参考系篇

2 坐标参考系基础 2.1 CRS 一个二维平面,我们可以使用如图1所示坐标系统,通过坐标唯一确定点位置: 图1 现实世界地球作为一个球体,当我们想要用同样方式利用坐标来唯一确定地球球面上某个位置...而当我们想要在纸面或电脑屏幕上绘制平面地图,就又需要有一套将地球球面展平方法。...,但纬度经度位置没有使用统一测量单位。...;而上述两个示例中都带有towgs84=0,0,0,这是一个转换因子,需要进行数据转换使用。.../ref/epsg/查看搜索所有已知EPSG与CRS对应关系(图11): 图11 或在QGIS查看: 图12 譬如对于重庆,因为地跨东经105°11~110°11,中轴线距离108E更近,常用如下投影

1.8K21

(数据科学学习手札75)基于geopandas空间数据分析——坐标参考系篇

而当我们想要在纸面或电脑屏幕上绘制平面地图,就又需要有一套将地球球面展平方法,上述这些用于不同情况下定义对象位置信息坐标系统,就称为坐标参考系统(Coordinate Reference System...2.1.2 投影坐标系 地理坐标系虽然解决了我们地球球面上定位问题,但纬度经度位置没有使用统一测量单位,因为经度不变情况下,纬度每变化1单位因为是对固定弧长映射,所以真实距离是固定不变...11:声明对应11区(因为这里是横轴墨卡托所以拥有独立分区,但并不是所有CRS都有分区,且Proj4区号加S才为南半球分区11S,否则默认为北半球分区) datum=WGS84:声明基准面为WGS84...  与投影坐标系相比,没有单位units信息,因为地理坐标系通常单位为十进制度数;而上述两个示例中都带有towgs84=0,0,0,这是一个转换因子,需要进行数据转换使用。.../ref/epsg/查看搜索所有已知EPSG与CRS对应关系(图10): 图10   或在QGIS查看: 图11   譬如对于重庆,因为地跨东经105°11~110°

1.6K30

GIS坐标系测绘原理:大地水准面基准面参考椭球体EPSGSRIWKT

历史上出现了很多不同参考椭球体,很多还仍然使用。...现实生活,从一大堆统计点中拟合出一条有规律线,就需要用最小二乘法,其实就是正交投影思路。对应数学描述为:当W平面Ax = y无解,转换为Px= y形式,使其有解。...同一条经线上,纬度越低变形越大,最大值位于赤道上。 等角投影,无角度变形,面积比为长度比平方。 长度比等变形线平行于中央轴子午线。 优点:长度和面积变形是最小(比起其它投影)。...正轴投影,纬线为同心圆,其间隔由投影中心向外逐渐缩小,经线为同心圆半径。横轴投影,中央经线赤道为相互垂直直线,其他经线纬线分别为对称于中央经线赤道曲线。...优点:没有角度变形,由每一点向各方向长度比相等,它经纬线都是平行直线,且相交成直角。 缺点:长度和面积变形明显,纬线间隔从基准纬线处向两逐渐增大。

3.8K11

Web墨卡托公开小秘密

由于地球是一个赤道略宽两略扁不规则梨形球体,故其表面是一个不可展平曲面,所以运用任何数学方法进行这种转换都会产生误差变形,为按照不同需求缩小误差,就产生了各种投影方法,eg:墨卡托投影,高斯克吕格投影...墨卡托投影今天对于地图展示仍起着极为重要作用,目前世界各国绘制海洋地图仍广泛使用墨卡托投影.我们平时看到谷歌地图,百度地图,包括我们高德地图,都是使用墨卡托投影....首先墨卡托投影“等角”特性,保证了对象形状不变,以及方向相互位置正确,墨卡托投影“圆柱”特性,也保证了纬线经线平行相互垂直,而且经线间隔相同,但是纬线间隔从赤道向两级逐渐增大,“等角”不可避免带来面积巨大变形...投影过程,将表示地球参考椭球体近似的作为正球体处理(正球体半径 R = 椭球体半长轴 a)。...Auxiliary Sphere 就是告知你,这个坐标投影过程,将椭球体近似为正球体做投影变换,虽然基准面是WGS 1984 椭球面。

2.6K10

CARVE:2001-2013 年阿拉斯加火灾排放数据库(AKFED)

","9001"]]] 数据文件坐标系范围: : 2495174.3196 南: 866174.3196 西: -842644.8858 东: 843055.1142 应用推导 该产品时间跨度...包含所有不确定性来源 1000 次模拟,总碳消耗全区域标准偏差为 0.50 千克碳/平方米。 不确定性主要来源是黑云杉碳消耗模型,其次是土地覆盖分类。...黑云杉林主要分布寒冷、排水不畅、向北或低地地方,而白云杉落叶树种(主要是杨树桦树)则主要分布温暖、排水良好、向南且没有永久冻土地方(Viereck,1973 年;Bonan,1989 年)。...草地灌木丛生态系统出现在早期演替林分、排水不良地点、陡峭山坡以及林线线以上地方。阿拉斯加内陆植被马赛克因火灾火灾后演替而不断重塑。...所有分析均在阿拉斯加阿尔伯斯等面积投影范围内进行:中央子午线 = 西经 154 度;标准平行线 1 = 北纬 55 度;标准平行线 2 = 北纬 65 度;原点纬度 = 北纬 50 度;1983 年北美基准

5800

如何处理地图投影转换

之前在学习ggplot2geom_polygon()图层制作地图图形,从来没有苦恼过投影问题,因为coord_map()中直接给出投影转换参数,如果要制作基于国家地图,直接赋值为polyconic...需要花大把时间导入这些数据模型,并从模型抽取出geom_polygon所支持点、线、多边形数据,才能按照ggplot2所规范可视化语法进行制图。...R语言中支持GIS数据模型包一共有两个:sp包sf包,旧版ggplot2,geom_polygon高度依赖从sp导入数据对象(虽然也可以从sf获取)。...由于投影投影坐标系已经被投影算法转换,所以使用geom_text等图层函数,务必要使用与几何对象投影一致经纬度点,这里使用sf点中心计算函数最为快捷。...使用sf模型,导入素材通常要检查模型是否包含默认投影,如果有则可以直接进行转换,没有则最好先转化为WGS84(4236),然后再往其他投影坐标系进行转换。

1.7K30

R-ggplot2 标准中国地图制作

由于基础图表绘制系列推文还在加紧准备,所以这期推文还是对一些感兴趣图表进行绘制,这期涉及空间图表绘制,主要涉及知识点如下: ggplot2 + sf 中国标准地图绘制 cowplot 包实现地图子图插入...接下来我们进行探索性绘图,主要使用geom_sf() 进行绘制,代码如下: library(ggspatial) ggplot() + geom_sf(data = china,fill="NA"...比例尺指北针添加,则是依靠ggspatial包。 散点数据读取 由于只获得数据经纬度信息,其他数据均为虚构,数据预览如下: ?...注意红色框内容 接下来将数据地图中显示,具体代码如下: ggplot() + geom_sf(data = china,fill="NA",size=.5,color="black") +...这里需注意是xlimylim设置,具体值都是 +proj=laea +lat_0=40 +lon_0=104 投影坐标系下值,而常规坐标转不同坐标下操作也很简单,可使用st_transform

2.6K10

实用帖之R-ggplot2 标准中国地图制作

这期涉及空间图表绘制,主要涉及知识点如下: ggplot2 + sf 中国标准地图绘制 cowplot 包实现地图子图插入 ggplot2 + sf 中国标准地图绘制 本期绘图我们是按照如下图表进行绘制...接下来我们进行探索性绘图,主要使用geom_sf() 进行绘制,代码如下: library(ggspatial) ggplot() + geom_sf(data = china,fill="NA"...比例尺指北针添加,则是依靠ggspatial包。 散点数据读取 由于只获得数据经纬度信息,其他数据均为虚构,数据预览如下: ?...注意红色框内容 接下来将数据地图中显示,具体代码如下: ggplot() + geom_sf(data = china,fill="NA",size=.5,color="black") +...这里需注意是xlimylim设置,具体值都是 +proj=laea +lat_0=40 +lon_0=104 投影坐标系下值,而常规坐标转不同坐标下操作也很简单,可使用st_transform

1.4K10

Mastercam9.1

例如限定某一层,则绘制该层图素才能被选择,完成诸如分析,删除等操作。设置OFF,则系统可以认得出任何一个图层图素 WCS    世界坐标系         设置系统视角管理。...常用在图形文件转换肘,当有些构图面视角与Mastercam软件不兼容,可将其图素转正。...)或生成通过投影点沿着曲面法向及给定长度一矢量线         Prep/Dist 法向/距离        生成与一直线、圆弧或曲线法线上相距给定距离点         Grid 网格点 生成一系列网状点...生成二组相交曲面间线         Project 投影线        生成曲线曲面上投影线。...Baseline 基准标注 选一条线性尺寸线作为基准,以后生成尺寸线,均以该基准线一端点引出尺寸线                 Chained 串联标注 选一条线性尺寸线,以后生成尺寸线,均以该基准线一端点引出尺寸线有缘学习交流关注桃报

2.5K20

R-ggplot2 标准中国地图制作

由于基础图表绘制系列推文还在加紧准备,所以这期推文还是对一些感兴趣图表进行绘制,这期涉及空间图表绘制,主要涉及知识点如下: ggplot2 + sf 中国标准地图绘制 cowplot 包实现地图子图插入...接下来我们进行探索性绘图,主要使用geom_sf() 进行绘制,代码如下: library(ggspatial) ggplot() + geom_sf(data = china,fill="NA"...比例尺指北针添加,则是依靠ggspatial包。 散点数据读取 由于只获得数据经纬度信息,其他数据均为虚构,数据预览如下: ?...注意红色框内容 接下来将数据地图中显示,具体代码如下: ggplot() + geom_sf(data = china,fill="NA",size=.5,color="black") +...这里需注意是xlimylim设置,具体值都是 +proj=laea +lat_0=40 +lon_0=104 投影坐标系下值,而常规坐标转不同坐标下操作也很简单,可使用st_transform

1.3K20

R语言 | ggplot2绘制标准中国地图

由于基础图表绘制系列推文还在加紧准备,所以这期推文还是对一些感兴趣图表进行绘制,这期涉及空间图表绘制,主要涉及知识点如下: ggplot2 + sf 中国标准地图绘制 cowplot 包实现地图子图插入...ggplot2 + sf 中国标准地图绘制 本期绘图我们是按照如下图表进行绘制(该图来自朋友公号): 其实我很早就想试着用R进行此类图表绘制了,但由于强大Arcgis绘图进度等原因,导致迟迟没有进行...(部分): 接下来我们进行探索性绘图,主要使用geom_sf() 进行绘制,代码如下: library(ggspatial) ggplot() + geom_sf(data = china,fill...(scatter_df,coords = c("lon", "lat"),crs = 4326) 结果如下: 注意红色框内容 接下来将数据地图中显示,具体代码如下: ggplot() +...xlimylim设置,具体值都是 +proj=laea +lat_0=40 +lon_0=104 投影坐标系下值,而常规坐标转不同坐标下操作也很简单,可使用st_transform进行操作。

4.8K31

R可视乎|空间地理数据可视化(1)

最常用投影方式是墨卡托投影(Universal Transverse Mercator,UTM),这种投影方式将地球划分为60个经度为6度区域,每个区域都使用横向墨卡托投影绘制出一个南北方向范围...半球(南或)以及区东经北纬坐标(以米为单位)给出。...3.3 R 设置坐标参考系统 地球形状可以用一个扁椭球形模型来近似,它在赤道上隆起,扁平,目前世界上有很多不同参考椭球体来使用,最常用是全球定位系统(GPS)所使用世界大地测量系统...除此之外,还有欧洲石油调查组(EPSG)所制定地图,由于坐标系不同,各地地图也会不同,例如中国:以地球几何球心为中心EPSG 代码为 4479;以地球椭球焦点为中心EPSG 代码为 4480...WGS84 EPSG 代码为 4326。 R 语言中,CRS 是用 proj4 字符串指定,这些字符串指定了投影、椭球体基准点属性。

3.4K30

(模仿精进数据可视化01) 全国38城居住自由指数可视化

因此与其matplotlib中极坐标系基础上想方法隐藏部分参考线,不如逆向思维,从构造参考线角度出发,自己组织构造参考线,会更加自由灵活。...2.2.1 构建坐标系统   因为极坐标系参考线非常类似俯视南北极点所看到经纬线,因此我们可以利用地图学坐标参考系里正射投影(Orthographic),可以理解为纯粹半球: ?...图4   我们只需要设定中心点参数南极点或北极点,再配合简单经纬度相关知识就可以伪造出任意经纬线,再利用geopandas投影变换向设定好正射投影进行转换,再作为平面坐标进行绘图即可。   ...转换到设置好正射投影crs上,再作为不同图层进行叠加绘制: ?...图12   那么接下来我们要做事就so easy了,只需要分别得到两者去除重叠面后,剩余部分,以对应填充色彩叠加绘制图11图像上就可以啦~,利用geopandasdifference即可轻松实现

79110

模仿精进数据可视化01:国内38城居住自由指数

,发布了「2020 新一线城市居住报告」: 图1 而在这个报告中有几张数据可视化作品还是比较可圈可点,作为(模仿精进数据可视化)系列文章开篇之作,我将基于我观察原始数据可视化作品进而构思出方式...因此与其matplotlib中极坐标系基础上想方法隐藏部分参考线,不如逆向思维,从构造参考线角度出发,自己组织构造参考线,会更加自由灵活。...2.2.1 构建坐标系统 因为极坐标系参考线非常类似俯视南北极点所看到经纬线,因此我们可以利用地图学坐标参考系里「正射投影」(Orthographic),可以理解为纯粹半球: 图4 我们只需要设定中心点参数南极点或北极点...转换到设置好「正射投影」crs上,再作为不同图层进行叠加绘制: 图5 嘿嘿,是不是底层参考线已经有内味了~ 2.2.2 绘制指标折线 坐标系以及参考线逻辑定了下来之后,接下来我们需要将原作品中所展现...图12 那么接下来我们要做事就so easy了,只需要分别得到两者去除重叠面后,剩余部分,以对应填充色彩叠加绘制图11图像上就可以啦~,利用geopandasdifference即可轻松实现

84130

R-gstat+ggplot2-克里金(Kriging)插值计算及可视化绘制

,也可以参考之前文档(R-gstat-ggplot2 IDW计算及空间插值可视化绘制)。...我们根据数据分布选择 model="Exp",使用如下代码进行拟合线绘制: model.variog<-vgm(psill=125, model="Exp", nugget=45, range=.6)...接下来我们就这两种情况进行可视化绘制。 插值结果可视化绘制 我们有了规整好df类型数据,这就可以方便使用ggplot2进行可视化绘制。...注意:这里裁剪方法存在一些问题:面对较大面积,裁剪处不够明确(出现过多范围,这个问题在本系列结束给出完整文档中会给出解决方法)。...针对model="Gau"结果 由于步骤类似,我们这里直接给出代码绘图结果即可: library(sf) library(tidyverse) library(ggspatial) library(

3.7K32

墨卡托投影坐标系(Mercator Projection)原理及实现C代码

墨卡托投影今天对于航海事业起着极为重要作用,目前世界各国绘制海洋地图仍广泛使用墨卡托投影,国际水路局(IHB)规定:“除特殊情况外,各国都要用墨卡托投影绘制海图”。...国际水路局发行《大洋水深总图》是把全世界分成24幅编辑南北纬72度之间就是使用墨卡托投影绘成。...墨卡托投影有一个特别的特性:所有罗盘等角线,或称斜航线(就是与所经过所有经线形成相同角度航线,也称恒向航线)墨卡托投影下都是直线。这使得航海领域这个投影非常重要。   ...墨卡托投影对透视圆筒投影改造点:要使圆筒投影称为等角性质,必须使由赤道向两经线逐渐伸长倍数与经线上各点相应纬度扩大倍数相同。 ? ? 透视圆筒投影 ?...,文献基础自己理解上实现,仅略知皮毛,疏漏在所难免,承蒙诸君不吝告知,将不胜感激!

4K50

音视频知识图谱 2022.06

全景图像及视频,等距柱状投影实现思路是用相同数量采样点保存每条纬线上数据,从而得到对应二维平面上矩形视频。...特点:全景视频,连续投影主要为 ERP EAP,很明显,他们两者主要区别于平面坐标 y 与球面坐标 φ 映射关系。...ERP 纬线被投影为等间隔水平线,而 EAP 不是,因为 EAP 保证了投影等面积。但两者都存在拉伸问题,引入了较大冗余。这种投影方式赤道处无失真,即赤道为标准纬线,失真向两迅速递增。...像任何圆柱投影一样,远离赤道拉伸越严重。两将一点拉伸为一条线,产生了无限失真。 AEP,Adjusted Equal-area Projection。...事实上,EAP 是 AEP 投影参数 β=1 一种特定情形。

57230
领券