Hi,
After much deliberation, I found a passable solution:
distances = np.abs(np.arange(0,resolution,1)+0.5-(resolution/2.0)) x_gradient = np.tile(distances,(resolution,1)) y_gradient = np.copy(x_gradient) y_gradient = np.swapaxes(y_gradient,0,1) distances_to_center = np.hypot(x_gradient,y_gradient) angles = np.radians( (distances_to_center/((resolution/2.0)-0.5)) * 90.0 ) lambert_weights = np.cos(angles)
This seems like it could be improved.
Ian
PS, I meant "Sinusoidal" in the title. I carried that over into my description, which should have been "cosine", instead of "sine".