AskOverflow.Dev

AskOverflow.Dev Logo AskOverflow.Dev Logo

AskOverflow.Dev Navigation

  • 主页
  • 系统&网络
  • Ubuntu
  • Unix
  • DBA
  • Computer
  • Coding
  • LangChain

Mobile menu

Close
  • 主页
  • 系统&网络
    • 最新
    • 热门
    • 标签
  • Ubuntu
    • 最新
    • 热门
    • 标签
  • Unix
    • 最新
    • 标签
  • DBA
    • 最新
    • 标签
  • Computer
    • 最新
    • 标签
  • Coding
    • 最新
    • 标签
主页 / coding / 问题 / 79247551
Accepted
Dr. Hillier Dániel
Dr. Hillier Dániel
Asked: 2024-12-03 20:13:12 +0800 CST2024-12-03 20:13:12 +0800 CST 2024-12-03 20:13:12 +0800 CST

Matplotlib 无法使用 plt.xscale('log') 一致地转换绘制元素

  • 772

绘制数据时,垂直线被 matplotlib 错误地转换。使用线性 x 轴,我的曲线和指向曲线上特定位置的垂直线完全匹配。 垂直线如预期接触曲线 在 plt.xscale('log') 之后,垂直线的终点不再在我的曲线上。 对数变换后的第二张截图

import numpy as np
import matplotlib.pyplot as plt

# Generate example data: x values from 1 to 16 (log scale)
x = np.array([1, 2, 4, 8, 16])
y = np.array([0.1, 0.2, 0.3, 0.6, 0.9])

# Perform linear interpolation to find where y crosses 0.5
crossings = np.where(np.diff(np.sign(y - 0.5)))[0]

if len(crossings) > 0:
    crossing_index = crossings[0]
    x1, x2 = x[crossing_index], x[crossing_index + 1]
    y1, y2 = y[crossing_index], y[crossing_index + 1]
    
    # Linear interpolation to find exact x for y = 0.5
    x_nd50 = x1 + (0.5 - y1) * (x2 - x1) / (y2 - y1)
    
    print(f"Interpolated x (ND50) = {x_nd50}")

# Plot the data with a logarithmic x scale
plt.plot(x, y, label="Data", color="blue")

# Plot the vertical line at the interpolated ND50 value
plt.plot([x_nd50, x_nd50], [0, 0.5], color='red', linestyle="--", alpha=0.7)
plt.scatter(x_nd50, 0.5, color='red', marker='x', alpha=0.7)

# First screenshot taken at this point!

# Set x-axis to log scale (log2)
plt.xscale('log')

# Second screenshot taken at this point!

# Show the plot
plt.xlabel('Log-scaled X')
plt.ylabel('Y')
plt.legend()
plt.show()

matplotlib
  • 1 1 个回答
  • 18 Views

1 个回答

  • Voted
  1. Best Answer
    9769953
    2024-12-03T20:30:17+08:002024-12-03T20:30:17+08:00

    您可以使用线性插值计算 x 位置。但在对数图上,线性插值不起作用。

    事实上,Matplotlib 在对数图上显示直线,但这些直线与线性图上的直线不同:如果它们相同,它们将是曲线,反之亦然。Matplotlib 只是尽力而为,没有关于点之间曲线形状的额外信息。(并且它假设您选择了对数空间,以便曲线将显示为直线,因此它默认为点之间的直线。)

    因此,您需要使用对数 x 数据计算 x 位置。您可以按如下方式解决:

    [...]
    
    xlog = True
    # Perform linear interpolation to find where y crosses 0.5
    crossings = np.where(np.diff(np.sign(y - 0.5)))[0]
    
    if len(crossings) > 0:
        crossing_index = crossings[0]
        x1, x2 = x[crossing_index], x[crossing_index + 1]
        y1, y2 = y[crossing_index], y[crossing_index + 1]
    
        if xlog:
            x1 = np.log10(x1)
            x2 = np.log10(x2)
        # Linear interpolation to find exact x for y = 0.5
        x_nd50 = x1 + (0.5 - y1) * (x2 - x1) / (y2 - y1)
        if xlog:
            x_nd50 = 10**x_nd50
        print(f"Interpolated x (ND50) = {x_nd50}")
    
    [...]
    
    # Set x-axis to log scale (log2)
    if xlog:
        plt.xscale('log')
    
    [...]
    
    • 1

相关问题

  • 二维散点图 - mcolors.Normalize 与 mcolors.LogNorm

  • 在 Streamlit 中渲染人口普查地图的正确方法

  • 使用 python 通过 matplotlib 在图表中添加自定义信息

  • 如何将 3D 直方图转换为热图

  • 为什么使用 freqz 并在 Octave 和 SciPy 之间进行比较时,绘图会翻转?

Sidebar

Stats

  • 问题 205573
  • 回答 270741
  • 最佳答案 135370
  • 用户 68524
  • 热门
  • 回答
  • Marko Smith

    Vue 3:创建时出错“预期标识符但发现‘导入’”[重复]

    • 1 个回答
  • Marko Smith

    为什么这个简单而小的 Java 代码在所有 Graal JVM 上的运行速度都快 30 倍,但在任何 Oracle JVM 上却不行?

    • 1 个回答
  • Marko Smith

    具有指定基础类型但没有枚举器的“枚举类”的用途是什么?

    • 1 个回答
  • Marko Smith

    如何修复未手动导入的模块的 MODULE_NOT_FOUND 错误?

    • 6 个回答
  • Marko Smith

    `(表达式,左值) = 右值` 在 C 或 C++ 中是有效的赋值吗?为什么有些编译器会接受/拒绝它?

    • 3 个回答
  • Marko Smith

    何时应使用 std::inplace_vector 而不是 std::vector?

    • 3 个回答
  • Marko Smith

    在 C++ 中,一个不执行任何操作的空程序需要 204KB 的堆,但在 C 中则不需要

    • 1 个回答
  • Marko Smith

    PowerBI 目前与 BigQuery 不兼容:Simba 驱动程序与 Windows 更新有关

    • 2 个回答
  • Marko Smith

    AdMob:MobileAds.initialize() - 对于某些设备,“java.lang.Integer 无法转换为 java.lang.String”

    • 1 个回答
  • Marko Smith

    我正在尝试仅使用海龟随机和数学模块来制作吃豆人游戏

    • 1 个回答
  • Martin Hope
    Aleksandr Dubinsky 为什么 InetAddress 上的 switch 模式匹配会失败,并出现“未涵盖所有可能的输入值”? 2024-12-23 06:56:21 +0800 CST
  • Martin Hope
    Phillip Borge 为什么这个简单而小的 Java 代码在所有 Graal JVM 上的运行速度都快 30 倍,但在任何 Oracle JVM 上却不行? 2024-12-12 20:46:46 +0800 CST
  • Martin Hope
    Oodini 具有指定基础类型但没有枚举器的“枚举类”的用途是什么? 2024-12-12 06:27:11 +0800 CST
  • Martin Hope
    sleeptightAnsiC `(表达式,左值) = 右值` 在 C 或 C++ 中是有效的赋值吗?为什么有些编译器会接受/拒绝它? 2024-11-09 07:18:53 +0800 CST
  • Martin Hope
    The Mad Gamer 何时应使用 std::inplace_vector 而不是 std::vector? 2024-10-29 23:01:00 +0800 CST
  • Martin Hope
    Chad Feller 在 5.2 版中,bash 条件语句中的 [[ .. ]] 中的分号现在是可选的吗? 2024-10-21 05:50:33 +0800 CST
  • Martin Hope
    Wrench 为什么双破折号 (--) 会导致此 MariaDB 子句评估为 true? 2024-05-05 13:37:20 +0800 CST
  • Martin Hope
    Waket Zheng 为什么 `dict(id=1, **{'id': 2})` 有时会引发 `KeyError: 'id'` 而不是 TypeError? 2024-05-04 14:19:19 +0800 CST
  • Martin Hope
    user924 AdMob:MobileAds.initialize() - 对于某些设备,“java.lang.Integer 无法转换为 java.lang.String” 2024-03-20 03:12:31 +0800 CST
  • Martin Hope
    MarkB 为什么 GCC 生成有条件执行 SIMD 实现的代码? 2024-02-17 06:17:14 +0800 CST

热门标签

python javascript c++ c# java typescript sql reactjs html

Explore

  • 主页
  • 问题
    • 最新
    • 热门
  • 标签
  • 帮助

Footer

AskOverflow.Dev

关于我们

  • 关于我们
  • 联系我们

Legal Stuff

  • Privacy Policy

Language

  • Pt
  • Server
  • Unix

© 2023 AskOverflow.DEV All Rights Reserve