An automatic Polycube construction and parameterization method is proposed in this paper. This approach first decomposes the input mesh into a set of feature regions. Each region is approximated by a simple Polycube which is further split into several patches. Since each patch corresponds to a rectangular sub-surface of the Polycube, they can be parameterized independently. A smoothing procedure between patches is performed to reduce the overall parametric distortion. This method makes Polycube parameterization more convenient in texture mapping and other applications.