A computer-implemented method for 3D reconstruction, the method comprising: - providing 2D images which all represent a same real object which is fixed, and, for each 2D image, camera parameters which define a perspective projection, the camera parameters including camera extrinsic parameters; - providing, for each 2D image, a smooth map, the smooth map having pixel values each representing a measurement of contour presence; and - determining a 3D modeled object that represents the real object, the determining iteratively optimizing an energy that rewards, for each smooth map, projections of silhouette vertices of the 3D modeled object having pixel values, on the smooth map, representing a high measurement of contour presence; wherein the 3D modeled object is a control mesh of a subdivision surface, the control mesh being determined from an initial dense mesh approximating the real object. wherein the method further comprises, at each iteration of the optimizing: - providing a current first mesh, which is the control mesh; - providing a current second mesh, the current second mesh being obtained by subdividing the current first mesh; - for each 2D image: ▪ projecting the current second mesh on the 2D image; ▪ computing a 2D silhouette of the 3D modeled object in the 2D image; and ▪ identifying points of the current second mesh corresponding to the 2D silhouette, as the silhouette vertices; and - modifying the current first mesh to reduce the energy. - providing 2D images which all represent a same real object which is fixed, and, for each 2D image, camera parameters which define a perspective projection, the camera parameters including camera extrinsic parameters; - providing, for each 2D image, a smooth map, the smooth map having pixel values each representing a measurement of contour presence; and - determining a 3D modeled object that represents the real object, the determining iteratively optimizing an energy that rewards, for each smooth map, projections of silhouette vertices of the 3D modeled object having pixel values, on the smooth map, representing a high measurement of contour presence; wherein the 3D modeled object is a control mesh of a subdivision surface, the control mesh being determined from an initial dense mesh approximating the real object. wherein the method further comprises, at each iteration of the optimizing: - providing a current first mesh, which is the control mesh; - providing a current second mesh, the current second mesh being obtained by subdividing the current first mesh; - for each 2D image: ▪ projecting the current second mesh on the 2D image; ▪ computing a 2D silhouette of the 3D modeled object in the 2D image; and ▪ identifying points of the current second mesh corresponding to the 2D silhouette, as the silhouette vertices; and - modifying the current first mesh to reduce the energy. - providing a current first mesh, which is the control mesh; - providing a current second mesh, the current second mesh being obtained by subdividing the current first mesh; - for each 2D image: ▪ projecting the current second mesh on the 2D image; ▪ computing a 2D silhouette of the 3D modeled object in the 2D image; and ▪ identifying points of the current second mesh corresponding to the 2D silhouette, as the silhouette vertices; and ▪ projecting the current second mesh on the 2D image; ▪ computing a 2D silhouette of the 3D modeled object in the 2D image; and ▪ identifying points of the current second mesh corresponding to the 2D silhouette, as the silhouette vertices; and - modifying the current first mesh to reduce the energy.