Abstract: The bathymetry data combined with intensity image are presented,and the feature points are obtained through the minimal angles of lines,then the corresponding image points are decided by the maximal correlate coefficient in searching space.Finally,the second order polynomial is used to simulate the deformation model.After imagery has been co-registered,wavelet is used to fuse imagery.It is shown that this algorithm can be used in the flat seafloor or the isotropic seabed.Some examples have been given.