在python3中绘制散点图,其中x轴是以公里为单位的纬度/经度,y轴是深度

如何解决在python3中绘制散点图,其中x轴是以公里为单位的纬度/经度,y轴是深度

我试图找到绘制一些数据的最佳方法。基本上我有一个包含纬度、经度、深度、sample_ID、Group_ID 列的数据文件。我想生成一个二维散点图,其中 y 是深度,x 是从北到南的公里距离(或相对于在指定方向采样的第一个站计算横断面距离),类似于 ODV 样式地图下面一个:

enter image description here

更新

我想在我最初的问题中添加更多信息。经过更多的搜索和测试,我在 R 中找到了一个可能的解决方案,使用 geosphere 包和 distGEO 函数将我的坐标转换为以公里为单位的距离,然后可以映射。 (https://www.rdocumentation.org/packages/geosphere/versions/1.5-10/topics/distGeo)

如果有人知道 python 方法来做到这一点,那就太好了!

更新

ODV 不允许我进行我需要的定制。我想生成一个这样的图,我可以在其中指定元数据变量来为点着色。更具体地说,我的数据文件中的 group_ID 列在下面我的文件示例中看到。

Latitude    Longitude   Depth_m Sample_ID   Group_ID
49.7225 -42.4467    10  S1  1
49.7225 -42.4467    50  S2  1
49.7225 -42.4467    75  S3  1
49.7225 -42.4467    101 S4  1
49.7225 -42.4467    152 S5  1
49.7225 -42.4467    199 S6  1
46.312  -39.658 10  S7  2
46.312  -39.658 49  S8  2
46.312  -39.658 73  S9  2
46.312  -39.658 100 S10 2
46.312  -39.658 153 S11 2
46.312  -39.658 198 S12 2

不过,它给我带来了很多麻烦。我已经使用半正弦计算计算了坐标之间的距离,但是一旦我到达那里,我不确定如何使用这些距离合并到散点图中。这是我目前所拥有的:

import matplotlib.pyplot as plt
import pandas as pd
import numpy as np
#import haversine as hs
from math import radians
from sklearn.neighbors import DistanceMetric
df=pd.read_csv("locations.csv",sep="\t",index_col="Sample_ID")
#plt.scatter(df['Latitude'],df['Depth_m'])
#plt.show()
df['Latitude'] = np.radians(df['Latitude'])
df['Longitude'] = np.radians(df['Longitude'])
dist = DistanceMetric.get_metric('haversine')
x = dist.pairwise(np.unique(df[['Latitude','Longitude']].to_numpy(),axis=0))*6373
print(x)

这段代码为我提供了一个坐标距离矩阵,但老实说,我无法弄清楚如何将其拉入散点图,该散点图将 x 轴设置为从北到南。特别是因为必须考虑具有相同坐标的多个深度。非常感谢任何帮助绘图!

解决方法

对于距离计算,您可以使用 geopy 包,特别是 geopy.distance.geodesic(),通过假设特定椭球来计算沿弧的距离(例如WGS84)。

要生成类似于您所描述的图,您可以使用 ma​​tplotlib 库的散点图功能,特别是 ma​​tplotlib.pyplot.scatter()

下面的代码示例将引导您完成距离计算(从某个参考纬度/经度到另一个纬度/经度的距离……这不一定是 N-S 分量,但计算起来很容易)。以及如何使用 Group_ID 字段生成散点图以使用两种方法为点着色。

time.Time

和输出数字... enter image description here

,

我不确定你想用距离做什么,但从概念上讲,你需要像我所做的那样将你的 x 输出作为一个新列输入到你的数据框中。就组的不同颜色而言,我会使用 { {1}},因为它们有一个 seaborn 参数。请查看下面第一个散点图的输出,并尝试使用第二个散点图执行的操作:

hue

enter image description here

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。

相关推荐


依赖报错 idea导入项目后依赖报错,解决方案:https://blog.csdn.net/weixin_42420249/article/details/81191861 依赖版本报错:更换其他版本 无法下载依赖可参考:https://blog.csdn.net/weixin_42628809/a
错误1:代码生成器依赖和mybatis依赖冲突 启动项目时报错如下 2021-12-03 13:33:33.927 ERROR 7228 [ main] o.s.b.d.LoggingFailureAnalysisReporter : *************************** APPL
错误1:gradle项目控制台输出为乱码 # 解决方案:https://blog.csdn.net/weixin_43501566/article/details/112482302 # 在gradle-wrapper.properties 添加以下内容 org.gradle.jvmargs=-Df
错误还原:在查询的过程中,传入的workType为0时,该条件不起作用 <select id="xxx"> SELECT di.id, di.name, di.work_type, di.updated... <where> <if test=&qu
报错如下,gcc版本太低 ^ server.c:5346:31: 错误:‘struct redisServer’没有名为‘server_cpulist’的成员 redisSetCpuAffinity(server.server_cpulist); ^ server.c: 在函数‘hasActiveC
解决方案1 1、改项目中.idea/workspace.xml配置文件,增加dynamic.classpath参数 2、搜索PropertiesComponent,添加如下 <property name="dynamic.classpath" value="tru
删除根组件app.vue中的默认代码后报错:Module Error (from ./node_modules/eslint-loader/index.js): 解决方案:关闭ESlint代码检测,在项目根目录创建vue.config.js,在文件中添加 module.exports = { lin
查看spark默认的python版本 [root@master day27]# pyspark /home/software/spark-2.3.4-bin-hadoop2.7/conf/spark-env.sh: line 2: /usr/local/hadoop/bin/hadoop: No s
使用本地python环境可以成功执行 import pandas as pd import matplotlib.pyplot as plt # 设置字体 plt.rcParams['font.sans-serif'] = ['SimHei'] # 能正确显示负号 p
错误1:Request method ‘DELETE‘ not supported 错误还原:controller层有一个接口,访问该接口时报错:Request method ‘DELETE‘ not supported 错误原因:没有接收到前端传入的参数,修改为如下 参考 错误2:cannot r
错误1:启动docker镜像时报错:Error response from daemon: driver failed programming external connectivity on endpoint quirky_allen 解决方法:重启docker -> systemctl r
错误1:private field ‘xxx‘ is never assigned 按Altʾnter快捷键,选择第2项 参考:https://blog.csdn.net/shi_hong_fei_hei/article/details/88814070 错误2:启动时报错,不能找到主启动类 #
报错如下,通过源不能下载,最后警告pip需升级版本 Requirement already satisfied: pip in c:\users\ychen\appdata\local\programs\python\python310\lib\site-packages (22.0.4) Coll
错误1:maven打包报错 错误还原:使用maven打包项目时报错如下 [ERROR] Failed to execute goal org.apache.maven.plugins:maven-resources-plugin:3.2.0:resources (default-resources)
错误1:服务调用时报错 服务消费者模块assess通过openFeign调用服务提供者模块hires 如下为服务提供者模块hires的控制层接口 @RestController @RequestMapping("/hires") public class FeignControl
错误1:运行项目后报如下错误 解决方案 报错2:Failed to execute goal org.apache.maven.plugins:maven-compiler-plugin:3.8.1:compile (default-compile) on project sb 解决方案:在pom.
参考 错误原因 过滤器或拦截器在生效时,redisTemplate还没有注入 解决方案:在注入容器时就生效 @Component //项目运行时就注入Spring容器 public class RedisBean { @Resource private RedisTemplate<String
使用vite构建项目报错 C:\Users\ychen\work>npm init @vitejs/app @vitejs/create-app is deprecated, use npm init vite instead C:\Users\ychen\AppData\Local\npm-