如何使用htslib从VCF文件中提取每个样本的基因型信息？

#include <htslib/vcf.h>

int main() {
    // 打开VCF文件
    htsFile *vcf_file = bcf_open("path/to/vcf_file.vcf", "r");

    // 读取VCF文件头信息
    bcf_hdr_t *vcf_header = bcf_hdr_read(vcf_file);

    // 遍历每个样本
    int num_samples = bcf_hdr_nsamples(vcf_header);
    for (int i = 0; i < num_samples; i++) {
        // 读取下一条VCF记录
        bcf1_t *vcf_record = bcf_init();
        if (bcf_read(vcf_file, vcf_header, vcf_record) < 0) {
            break;
        }

        // 解析基因型信息
        int32_t *genotypes = NULL;
        int num_genotypes = 0;
        if (bcf_get_genotypes(vcf_header, vcf_record, &genotypes, &num_genotypes) > 0) {
            // 处理基因型信息
            for (int j = 0; j < num_genotypes; j++) {
                int32_t genotype = genotypes[j];
                // 进行你的处理操作
            }
        }

        // 释放内存
        free(genotypes);
        bcf_destroy(vcf_record);
    }

    // 关闭VCF文件
    bcf_hdr_destroy(vcf_header);
    bcf_close(vcf_file);

    return 0;
}