Compensation for rendering device non-uniformities is provided for halftoned images. A spatially dependent tone reproduction curve (TRC) provides compensation values. Pixel location information is used to access TRC values. For example, the values are modification values. The modification values are added to the pixel values to generate combined values. Quantization is applied to the combined values to prepare compensated image data for rendering. For example, Rank Ordered Error Diffusion is applied to the combined values. The combined values may include diffused error from previously processed pixels. Gray values may be estimated for the respective pixels. The estimated gray values may be used to access compensation information from a TRC that is both spatially and gray value dependent. Mathematical basis decomposition is used to reduce TRC memory requirements. For example, Discrete Cosine Transformation, Singular Value Decomposition or Principal Component Analysis is used to determine a compact form for the TRC.