I settled on 0 based indexing as the right way to go when I was doing this calculation: https://madhadron.com/posts/2009-07-17-determining-affine-tr...

When you start calculating offsets and transformations of coordinates, starting at 1 forces you to carry around extra additive constants everywhere.