A system for automated metrology, measurement, and model correction of a three dimensional model for 3D printability, said system comprising: at least one processor configured to: receive or generate slices of the 3D model, the slices representing two dimensional solids of the 3D model to be printed in corresponding print layers; calculate medial axis transforms of the slices, the medial axis transforms representing the slices in terms of corresponding medial axes; and determine a local feature size at any point along a boundary of the slices as the shortest distance from the point to a corresponding medial axis, wherein the local feature size is defined in terms of the medial axis transform and used to automatically identify regions that require a differential addition of material to ensure manufacturability, wherein the at least one processor is further configured to: simulate printing of the slices to identify corresponding printability maps, each printability map partitioning a slice into at least two disjoint sets of regions, the at least two disjoint sets including a first set of regions formed of features sized smaller than a threshold feature size and a second set of regions formed of features sized larger than the threshold feature, and to: identify a largest local feature size of the second set for each of the slices; determine a padding value for each slice of the 3D model or for the entire 3D model from the largest local feature sizes of the slices; and thicken each of the slices by the corresponding padding value. at least one processor configured to: receive or generate slices of the 3D model, the slices representing two dimensional solids of the 3D model to be printed in corresponding print layers; calculate medial axis transforms of the slices, the medial axis transforms representing the slices in terms of corresponding medial axes; and determine a local feature size at any point along a boundary of the slices as the shortest distance from the point to a corresponding medial axis, wherein the local feature size is defined in terms of the medial axis transform and used to automatically identify regions that require a differential addition of material to ensure manufacturability, wherein the at least one processor is further configured to: simulate printing of the slices to identify corresponding printability maps, each printability map partitioning a slice into at least two disjoint sets of regions, the at least two disjoint sets including a first set of regions formed of features sized smaller than a threshold feature size and a second set of regions formed of features sized larger than the threshold feature, and to: receive or generate slices of the 3D model, the slices representing two dimensional solids of the 3D model to be printed in corresponding print layers; calculate medial axis transforms of the slices, the medial axis transforms representing the slices in terms of corresponding medial axes; and determine a local feature size at any point along a boundary of the slices as the shortest distance from the point to a corresponding medial axis, wherein the local feature size is defined in terms of the medial axis transform and used to automatically identify regions that require a differential addition of material to ensure manufacturability, wherein the at least one processor is further configured to: simulate printing of the slices to identify corresponding printability maps, each printability map partitioning a slice into at least two disjoint sets of regions, the at least two disjoint sets including a first set of regions formed of features sized smaller than a threshold feature size and a second set of regions formed of features sized larger than the threshold feature, and to: identify a largest local feature size of the second set for each of the slices; determine a padding value for each slice of the 3D model or for the entire 3D model from the largest local feature sizes of the slices; and thicken each of the slices by the corresponding padding value.