QGIS любой версии прекрасно умеет это всё пересчитывать, параметры же обеих СК есть, вопрос именно в том, как это организовать конкретно в калькуляторе атрибутов, где доступна одна встроенная фукнция transform.
Потому что нужно для текущего слоя динамически считать координаты центроидов сразу в двух СК, одна из которых - пользовательская. Ключевые моменты - слой один, и считать динамически (объект изменился - сразу пересчиталась координата в атрибутах). Иначе всё было бы просто.