Imaging a sample that includes phototransformable optical labels (“PTOLs”) with an optical system having a diffraction-limited resolution volume (DLRV), includes providing activation radiation to the PTOLs to activate a statistical subset of the PTOLs. A density of the PTOLs of the activated subset is less than an inverse of the DLRV. Excitation radiation is provided to the activated subset to excite activated PTOLs. Radiation emitted from the activated and excited PTOLs located at different focal planes of the optical system within the sample is detected with the optical system. The preceding steps are repeated one or more times, each time activating a different statistical subset of the plurality of PTOLs. Three-dimensional locations within the sample are determined, with a sub-diffraction-limited accuracy, of the activated and excited PTOLs based on the radiation emitted from the activated and excited PTOLs that is detected from the different focal planes of the optical system.