Is there any reason why :naxis isn't included in setproperty!? For example I'm interested in a function like sub(n) in astropy.wcs that truncates the axes. It looks like the easiest way to implement this is to change naxis because everything else references that in getproperty, but I can't do that with setproperty! and had to call setfield!, so I wonder if there's some consideration behind the choice that I'm not aware of, or maybe that could be easily added? Thanks so much!