charmap decoding error

[ElricPi]

Hello,

I faced another error while running this code to get the observation for the .fits file as well as the .tbl file. I have attached the screenshot - can you kindly look into it? I'll be grateful once more. I have also attached screenshot for the modules which i entered. The data was taken from the following archive: https://sha.ipac.caltech.edu/applications/Spitzer/SHA/#id=SearchByPosition&RequestClass=ServerRequest&DoSearch=true&SearchByPosition.field.radius=0.13888889000000001&UserTargetWorldPt=11.888;-25.28822;EQ_J2000&TargetPanel.field.targetName=ngc%20253&SimpleTargetPanel.field.resolvedBy=nedthensimbad&MoreOptions.field.prodtype=pbcd&shortDesc=Position&isBookmarkAble=true&isDrillDownRoot=true&isSearchResult=true

Also, the sample files for this code are running perfectly. Specifically the sample data for ngc 7023. Both the .fits and .tbl files.

Thank You.

[alexmaragko]

The charmap decoding is a warning and not the actual error. The "Format not recognized" is the problem here, and this is because you are trying to load an IRAC map observation. pyPAHdb however will only handle IRS spectra and spectral cubes.

[ElricPi]

Hello, Thank you for pointing the error but unfortunately I'm still not able to work with it and I think I'm having some trouble in understanding the kind of data I need to select. I would be grateful if you could kindly point me towards a Spitzer data which would avoid the error. Also I would really appreciate it if the data could be provided within the following bandpass: 5.2 - 14.5 microns.

Thank you once more.

Edit: Specifically I'm interested in the following objects: NGC 7027, NGC 7023, HD 97048, HD 44179, IRAS 02575+6017, NGC 253.

[PAHdb]

pyPAHdb works on infrared spectra. You provided pyPAHdb with an IRAC image. Spitzer spectra can be directly obtained from the SHA and/or, for example, CASSIS. Note that these data would still require some pre-processing before feeding it to pyPAHdb. For example, you might need to splice the different orders and certainly would want to remove any continuum and emission line components. For a somewhat comprehensive overview for the steps involved you could take a look at Boersma+ (2018).

Unfortunately, we are not in a position to provide you with the pre-processed data; this is the researcher's own responsibility. For some of the objects you list no Spitzer data are available. Though, you can perhaps look at the Sloan ISO SWS Atlas for those spectra. Note that a Spitzer spectral cube for NGC 7023 is actual included with pyPAHdb, which is used by example_fits.py. These data are from Boersma+ (2018).

Since pyPAHdb shows expected behavior, I'm closing the issue.