4D imaging of the ionosphere allows the progress of an ionospheric storm to be studied across large geographical regions near-continuously in latitude, longitude, altitude and time. We review our data-assimilative imaging algorithm, Ionospheric Data Assimilation Three-Dimensional (IDA3D), and present our latest results. IDA3D permits high fidelity tracking of ionospheric F-region patches, revealing their development and transport, and could be used for ionospheric corrections to radar and other applications. A study of the October -November 2003 storm period illustrates the global coverage that can be achieved with ionospheric imaging, and permits study of the development of positive and negative ionospheric storms, and important features such as Storm Enhanced Density (SED) that may have application in the FAA.s Wide Area Augmentation System (WAAS). Finally the future development envisaged for IDA3D is discussed.