An automatic and efficient edge blending method was presented to improve the imagery quality of multi-projectors display system with anomalistic and multiple overlap regions. Using the pixel-level geometric calibration,the edge blending problem was formulated as a non-linear optimization problem with the object function being the global consistency of brightness distribution. An iterative computing method accelerated by GPU was presented to solve such optimal problem. The theoretic analysis was carried out to demonstrate the importance of enlarging the region of edge blending. Real applications showed that the proposed method can make the multi-projectors jointly display a visual seamless imagery with anomalistic and multiple overlap regions and maximize the edge blending regions by using all pixels of projectors.