The most obvious improvement to the techniques used here is to weight the data to account for the missing stations and for the missing data within each station. The resulting spectra show an imbalance in the imaging directions, and each set of spectra show noise that is probably caused by the muting and tapering.
The response of the earth to gravitational waves of various frequencies could be better defined than that of the crude assumptions used here by modelling using the known structure of the earth. The effects of the earth's rotation should also be included in the response.
The earthquake removal and muting of noise-dominated stations could be improved. Smaller earthquakes might be ignored. Large earthquakes might have better time constraints on the mutes. The small spikes remaining after the earthquake removal should be removed, perhaps with a median filter. The effects of the taper and gain should be examined with more care.
To provide some expected responses, a numerical connection between a given binary system and the earth's response should be made. It appears the the most difficult part of this problem is defining the details of the earth's response.
The output spectra might be better displayed as an amplitude chart similar to the dartboard form of Figure . These results might also be deconvolved to produce better resolution of the source position, making identification of weak sources easier over several years.