Abstract: Higher-order smooth surface extraction method can overcome aliasing artifacts of marching cubes algorithm．However,it introduces extra computation burden especially with optimal embedding function calculation．To resolve the problem,a parallel implementation based on graphics processing unit was presented．The original higher-order smooth surface extraction algorithm was divided into five parts including margin region,narrow band region,embedding function margin values,optimal embedding function and triangular mesh extraction．Then they were paralleled with task assignments method．The most complexity function embedding calculation is approximated by projected Jacobin method．The experimental results show that the speedup achieves more than 9after parallelization on GeForce GT 240M GPU.

Key words: Marching cubes,Surface extraction,Graphics processing unit,Higher-order smooth