使用SeuratDisk加载h5seurat时遇'dims'需包含所有(i,j)对错误求助

问题场景 #

使用SeuratDisk将h5ad文件转为h5seurat格式后，执行加载操作触发如下错误：

Validating h5Seurat file
Initializing RNA with data
Error in sparseMatrix(i = x[["indices"]][] + 1, p = x[["indptr"]][], x = x[["data"]][], :
'dims' must contain all (i,j) pairs

完整复现代码：

library(SeuratDisk)
Convert("xxx.h5ad","h5seurat",overwrite = TRUE)
seuratobj <- LoadH5Seurat("xxxx.h5seurat",assays ='RNA')

解决建议 #

1. 验证原h5ad文件的稀疏矩阵完整性 #

该错误大概率源于原h5ad中的稀疏矩阵存在索引不匹配（如indptr/indices与data维度不对应）。可在Python中用scanpy验证：

import scanpy as sc
adata = sc.read_h5ad("xxx.h5ad")
# 检查矩阵基本信息
print(adata.X.shape)
print(adata.X.nnz)
# 小数据集可尝试转稠密矩阵验证完整性
dense_X = adata.X.toarray()

若此步骤报错，说明原h5ad文件本身存在问题，需重新生成或修复原文件。

2. 转换时明确指定assay参数 #

转换过程中手动指定assay，避免自动识别偏差：

Convert("xxx.h5ad", dest = "h5seurat", overwrite = TRUE, assay = "RNA")

3. 更新Seurat及SeuratDisk至最新版本 #

旧版本工具可能存在转换逻辑bug，更新后重试：

install.packages("Seurat")
remotes::install_github("mojaveazure/seurat-disk")

4. 手动修正h5seurat的矩阵维度信息 #

若上述方法无效，可直接修改h5seurat文件中的维度配置：

library(hdf5r)
# 以读写模式打开h5seurat文件
h5 <- H5File$new("xxxx.h5seurat", mode = "r+")
# 查看当前错误的维度设置
h5[["assays/RNA/data/dims"]][]
# 替换为原h5ad对应的基因数、细胞数（需自行获取正确数值）
h5[["assays/RNA/data/dims"]][] <- c(基因数, 细胞数)
# 保存并关闭文件
h5$close()

修改完成后重新执行LoadH5Seurat加载。

5. 临时方案：转为稠密矩阵后再转换 #

若需应急查看数据，可先将h5ad转为稠密矩阵再进行格式转换：

import scanpy as sc
adata = sc.read_h5ad("xxx.h5ad")
adata.X = adata.X.toarray()
adata.write_h5ad("xxx_dense.h5ad")

再用SeuratDisk转换这个稠密版的h5ad文件。

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