We develop a methodology for estimating unobservable characteristics of the individual natural history of metastatic cancer from the volume of the primary tumor and site-specific volumes of metastases measured before, or shortly after, the start of treatment. In particular, we address the question as to what information about natural history of cancer can and cannot be gained from this type of data. Estimation of the natural history of cancer is based on parameterization of a very general mathematical model of cancer progression accounting for primary tumor growth, shedding of metastases, their selection, latency and growth in a given secondary site. This parameterization assumes Gompertz (and, as a limiting case, exponential) growth of the primary tumor, exponential growth of metastases, and exponential distribution of metastasis latency times. We find identifiable parameters of this model and give a rigorous proof of their identifiability. As an illustration, we analyze a clinical case of renal cancer patient who developed 55 lung metastases whose volumes were measured through laborious reading of CT images. The model with maximum likelihood parameters provided an excellent fit to this data. We uncovered many aspects of this patient's cancer natural history and showed that, according to the model, onset of metastasis occurred long before primary tumor became clinically detectable.