I'm trying to adjust batch effect using deseq2 limma::removeBatchEffect like below:

###### Batch Correction with limma removeBatchEffect #######

dds <- DESeqDataSetFromMatrix(countData = data,
                              colData = coldata,
                              design = ~ Samplebatch + cond)

dds <- DESeq(dds)
dds

dds$batch <- as.numeric(dds$Samplebatch)
dds$cond
dds$batch

## vst after adding batch information
vsd_abc <- vst(dds, blind = FALSE)
vsd_abc$Samplebatch
vsd_abc$batch

design0 <- model.matrix(~cond, colData(vsd_abc))
assay(vsd_abc) <- limma::removeBatchEffect(assay(vsd_abc), vsd_abc$batch, design=design0) ### Batch Correction

After batch effect correction I would like to get TPM / FPKM for performing ssGSEA and various other analyses which are normalized for the length of the gene.

Is it possible to get TPM / FPKM after batch effect removal using limma:removeBatchEffect?

Any help is appreciated. Thank you !!

Use DESeq2::fpkm(), convert to log2 and correct these values with limma. Or correct any other normalized counts on log2 scale as you like.