Higher Order Deviation from Ellipse

[eckerl]

Hi,

I am using Photutils extensivly to perform photometric analysis and i am quite happy with your program. I recently encountered a really boxy galaxy, so the normal isophotes left a "cross" pattern in the center. By setting "high_harmonics= True" i thought this will change since i am measuring significant a4,b4 parameters. Unfortunatly visually nothing seemed to change. I subtracted the normal elliptical_isophotes from the ones by setting "high_harmonics=True" and the Residuals are minimal. For me it seems like the measurement of the higher harmonics is correct but the build_ellipse_model does not cover it. I wanted to ask if someone/you runs into the same problems by fitting boxy isophotes or not. Any help would be super nice.

Best regards