RLM produces an unbiased estimator of the mean or mean function for symmetric distribution and is calibrated for the normal distribution. I don't know how well this is approximated by the log of an exponentially distributed variable, but it won't exactly satisfy the assumptions.
There should be a more direct way of estimating the parameter for the exponential distribution in a robust way, but I never tried.(one idea would be to estimate a trimmed mean and use the estimated distribution to correct for the trimming. scipy.stats.distributions have an `expect` method that can be used to calculate the mean of a trimmed distribution, i.e. conditional on lower and upper bounds)
What's your sample size?(for very large sample sizes one approach that is sometimes used, is to fit a distribution to the central part of a histogram)