Hi,
I do not see the original data (the JPEG looks clipped to me). However, in general for vignetting and "hot spots" due to flat errors- a simple circular symmetry (axial) will do the trick (likely much better/easier than placing samples in the corners- although it might work OK). In addition why choose subtraction as the method? A flat error is certainly multiplicative in nature (and not "additive" like additional light). So division seems like the best way.
I plan on demonstrating some of these things, including a nice example of Local Normalization, during my AIC workshop.
-adam