1. A method of detecting changes between a current image of an area and a historical image of said area, wherein said current image is taken at a time later than said historical image, said method comprising: retrieving said historical image from a database based on said current image and said historical image having corresponding geographical locations; relating pixel locations of said current image and said historical image; performing fine-scale co-registration of said current image and said historical image; optimizing inter-scene phase coherence of said current image and said historical image; performing local co-registration of said current image and said historical image based on said optimizing; and performing a canonical correlation analysis based on said local co-registration; wherein performing fine-scale co-registration further comprises applying a scale invariant feature transform to said current image and said historical image; and wherein performing local co-registration further comprises: performing non-overlapping patch correlation; obtaining correlation peaks for each said patch; estimating local relative patch translations in along-track and across-track dimensions; parameterizing along-track translations into surge and heading error vectors; parameterizing across-track translations into heave and sway error vectors; re-mapping said current image onto a grid corresponding to said historical image based on interpolation functions formed from said surge, heading, sway and heave vectors; calculating a coefficient of a complex correlation between said current image and said historical image, a phase of said complex correlation forming an interferogram; decomposing said interferogram into surge and sway functions; defining an argument of a phase function based on said surge and sway functions; and multiplying said phase function by said current image.