When a lot of input data are not distributed uniformly n a chord-span direction or when the given shape is complicated, it is very difficult to obtain an inverse matrix which represents the smooth Bi-cubic B-spline surface of the initial shape. To overcome this problem, we suggest image Surface Expansion Method(ISE Method) which is suggested for vertex creation of B-spline curves and surfaces. Its basic concept, convergency and verification are shown. Also B-spline curves and Surfaces represented by ISE Method were compared with those represented by the existing method which is based on the inverse matrix method, the pseudoinverse matrix method and the chord length approximation method for vertex yielding. Ship Hull Forms which have Knuckle, Bulbous Bow, Transom and Stern frame were represented by the ISE Method.