数轴图表特定位置标记需求及Python、R实现方案问询

你需要的自定义序列区间着色图不需要依赖现成的timeline、chromomap类封装库，直接用Python的matplotlib或者R的ggplot2基础绘图能力就可以实现，自定义程度可以完全匹配你的需求。

Python 实现方案（基于matplotlib） #

• 首先把你的数据整理为[序列名称, 区间起始位置, 区间结束位置, 对应颜色]的结构化格式，同时提前记录4个序列各自的总长度
• 用matplotlib.patches.Rectangle逐个绘制每个序列的底色，再叠加绘制指定区间的着色块即可
  示例代码如下，可直接替换其中的示例数据为你自己的实际数据：

import matplotlib.pyplot as plt
from matplotlib.patches import Rectangle

# 示例绘图数据，替换为你的实际数据
plot_data = [
    ["Length_1", 1, 2, "#2ecc71"],
    ["Length_1", 3, 7, "#ffffff"],
    ["Length_1", 8, 15, "#e74c3c"],
    ["Length_2", 1, 4, "#3498db"],
    ["Length_2", 5, 9, "#ffffff"],
    ["Length_2", 10, 12, "#2ecc71"],
    # 补充Length_3、Length_4的所有区间配置
]

# 4个序列的总长度，替换为你的实际数据
seq_meta = {
    "Length_1": 15,
    "Length_2": 12,
    "Length_3": 20,
    "Length_4": 18
}
seq_order = list(seq_meta.keys())
bar_height = 0.6

fig, ax = plt.subplots(figsize=(10, 4))
for seq_idx, seq_name in enumerate(seq_order):
    # 绘制序列完整底色
    ax.add_patch(Rectangle((1, seq_idx - bar_height/2), seq_meta[seq_name], bar_height, color="#f0f0f0"))
    # 绘制指定区间的着色块
    for item in plot_data:
        if item[0] == seq_name:
            s, e, color = item[1], item[2], item[3]
            ax.add_patch(Rectangle((s, seq_idx - bar_height/2), e - s + 1, bar_height, color=color, ec="#333"))

# 轴配置可根据需求调整
ax.set_yticks(range(len(seq_order)))
ax.set_yticklabels(seq_order)
ax.set_xlabel("序列位置")
ax.set_xlim(0, max(seq_meta.values()) + 1)
ax.set_ylim(-0.8, len(seq_order) - 0.2)
plt.tight_layout()
plt.show()

R 实现方案（基于ggplot2） #

和Python逻辑一致，用geom_rect绘制矩形色块即可，示例代码如下：

library(ggplot2)

# 示例绘图数据，替换为你的实际数据
plot_data <- data.frame(
  seq_name = c("Length_1", "Length_1", "Length_1", "Length_2", "Length_2", "Length_2"),
  start = c(1,3,8,1,5,10),
  end = c(2,7,15,4,9,12),
  color = c("#2ecc71", "#ffffff", "#e74c3c", "#3498db", "#ffffff", "#2ecc71")
)

# 4个序列的总长度，替换为你的实际数据
seq_meta <- data.frame(
  seq_name = c("Length_1", "Length_2", "Length_3", "Length_4"),
  total_len = c(15,12,20,18),
  y_pos = 1:4
)

plot_data <- merge(plot_data, seq_meta, by = "seq_name")

ggplot() +
  # 绘制序列完整底色
  geom_rect(data = seq_meta, aes(xmin = 1, xmax = total_len, ymin = y_pos - 0.3, ymax = y_pos + 0.3), fill = "#f0f0f0") +
  # 绘制指定区间的着色块
  geom_rect(data = plot_data, aes(xmin = start, xmax = end + 1, ymin = y_pos - 0.3, ymax = y_pos + 0.3, fill = I(color)), color = "#333") +
  scale_y_continuous(breaks = seq_meta$y_pos, labels = seq_meta$seq_name) +
  labs(x = "序列位置", y = "") +
  theme_minimal()

两种方案都支持自定义边框样式、刻度格式、图例添加，适配任意长度的序列和任意数量的着色区间。

内容的提问来源于stack exchange，提问作者user16600616