Transform coordinate system

        # TODO: needs expert attention

def set_coordsys(self, coordsys):
        """
        Transform coordinate system

        # TODO: needs expert attention
        """
        if coordsys in self.coordsys_mapping:
            self.coordsys = self.coordsys_mapping[coordsys]
        else:
            self.coordsys = coordsys