In a programming method, a number of points measured on a product are roughly fitted into/matched with one or more geometrical elements, and a fitting value of the geometrical elements are output to a first data list. The method records the storing positions of the geometrical elements having fitting values that are within a predetermined tolerance value in a second data list, and obtains a best geometrical element type of the product by determining geometrical element types that correspond to each of the storing positions recorded in the second data list using a matching condition. After the best geometrical element type from the first data list is obtained, the method generates a measurement program and displays the measurement program on a display device.