One aspect of the invention provides a method for drift correction to correct a 3D point collection dataset to compensate for drift over time. The method includes: (a) separating the 3D dataset into n segments, wherein n>1; (b) for each of the n segments, reconstructing a volume image as a 3D histogram in which a count for each voxel in the histogram equals a number of localization estimates falling within the voxel; (c) performing 3D cross-correlation between pairs of the n segments; (d) identifying a correlation peak in a result of the 3D cross-correlation to determine a shift distance between pairs of the n segments; (e) solving an overdetermined system of shift distances to determine independent shifts; and (f) offsetting positions from a plurality of segments in the 3D point collection dataset with the independent shifts calculated in step (e) to correct for drift.