Nextflow中如何实现变量类型转换？

Nextflow 基于 Groovy 语言，因此可以直接使用 Groovy 的原生类型转换方法处理变量类型转换问题，以下是针对你场景的具体解决方案：

基础类型转换 #

从 Bash 脚本输出（或文件读取）的数值通常是字符串格式，你可以通过以下方法转换为数值类型：

• 转整数：使用 toInteger() 方法（适用于整数结果）
• 转浮点数：使用 toFloat() 或 toDouble() 方法（适用于带小数的覆盖度结果）

注意：转换前建议先用 trim() 去除字符串中的换行符、空格等冗余字符，避免转换失败。

示例代码：

// 读取存储覆盖度结果的文件内容
def coverage_str = file('coverage.txt').text.trim()
// 转换为双精度浮点数
def coverage = coverage_str.toDouble()

逻辑判断实现 #

转换为数值类型后，即可直接进行逻辑判断：

if (coverage < 10.0) {
    println "警告：样本覆盖度仅为 ${coverage}x，低于阈值"
} else if (coverage >= 30.0) {
    println "样本覆盖度良好（${coverage}x）"
} else {
    println "样本覆盖度合格（${coverage}x）"
}

控制小数位数 #

如果需要格式化输出小数位数，有两种常用方式：

1. 字符串格式化：使用 String.format() 精确控制小数位数
2. 数值取整：通过 round() 方法结合数值运算实现

示例代码：

// 方法1：保留两位小数的字符串格式化
def formatted_coverage = String.format("%.2f", coverage)
println "格式化后覆盖度：${formatted_coverage}x"

// 方法2：通过数值运算保留两位小数
def rounded_coverage = (coverage * 100).round() / 100.0
println "取整后覆盖度：${rounded_coverage}x"

完整场景示例 #

结合你用 samtools + awk 计算覆盖度的场景，以下是一个完整的 Nextflow 流程示例：

process calculateCoverage {
    input:
    path bam_file

    output:
    path 'coverage.txt' into coverage_ch

    script:
    """
    samtools depth -a ${bam_file} | awk '{sum+=\$3; count++} END {print sum/count}' > coverage.txt
    """
}

workflow {
    // 输入BAM文件
    bam_files = Channel.fromPath('data/*.bam')
    
    calculateCoverage(bam_files)
        .map { cov_file ->
            // 读取并转换覆盖度数值
            def cov_str = cov_file.text.trim()
            def cov = cov_str.toDouble()
            
            // 格式化小数位数
            def formatted_cov = String.format("%.2f", cov)
            
            // 状态判断
            def status = cov < 10.0 ? "低覆盖度" : "正常"
            
            // 返回结果
            [cov_file.baseName, formatted_cov, status]
        }
        .view { sample, coverage, status ->
            println "样本 ${sample}：覆盖度 ${coverage}x，状态：${status}"
        }
}

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