Issue with AIA 171 Level 1.5 data

[Richard Morton]

Hi,

I am running the AIApy calibration routines to upgrade from Level 1 to Level 1.5 data.

However, wheneve I run the code the size of the array for 171 data ends up as (4094, 4094), while all other channels (correctly) end as (4096, 4096).

Has anbody experienced this?

Richard

I am running sunpy v 5.0.0 and aiapy 0.7.3

minimum working example

------------------------------------------------------

wvsrch=a.Wavelength(94*u.angstrom, 335*u.angstrom)
result = Fido.search(a.Time('2011-11-03T12:15:07', '2011-11-03T12:15:19'), a.Instrument("aia"), wvsrch)

files = Fido.fetch(result)

sdo_maps = Map(files)

aprep=[]
for m in sdo_maps:
    m_temp = update_pointing(m)
    reg = register(m_temp)
    aprep.append(reg)

[Nabil Freij]

Hi Richard,

This is a known issue with the function. Due to the nature of rotating and scaling the image, we cannot guarantee  it will return a 4096 by 4096 data array as sometimes there isn't any  data for that outer ring of 1 pixel.

Do you need them to be 4096 or you were just wondering why it was cropped by 2 pixels?

Kind regards,

Nabil

[Richard Morton]

Hi Nabil,

Thanks for the response. I actually found this information in the API just before I saw you response :D

I am doing DEM calculations and the routine I am using requires the 6 channels to be in a single array... so yep it is needed. But I can pad the edges.

Richard

[Nabil Freij]

Hi Richard,

That is a very common issue people run into when they use aiapy and then want to pass it into a DEM code.

We are still wondering how best to handle this for us and users going forward.

Sorry for the trouble!

Nabil