search

teuben / dc2019

Improving Image Fidelity on Astronomical Data: Radio Interferometer and Single-Dish Data Combination (workshop material)
9 stars 17 forks source link

DC_script.py fails at Step 4: Faridani short spacings combination (SSC) #20

Closed emoravec closed 3 years ago

emoravec commented 3 years ago

When running the DC_script.py, the process fails at Step 4: Faridani short spacings combination (SSC) exiting with the following error in CASA ending with TypeError: string indicies must be integers. Step4_error_150421

The terminal output is: ### Step 4 Faridani short spacings combination (SSC) ### Regridding lowres image... 2021-04-15 16:27:31 WARN imregrid::::casa You've specified that you want to regrid the spectral axis without specifying the output shape. Normally the length chosen would be that of the corresponding template axis, however, the template spectral axis is degenerate and one cannot regrid an axis such that its output length is one. So, removing axis 3 from the axes list and just copying the input spectral information to the output image 2021-04-15 16:27:31 WARN ImageRegrid::regrid You are regridding an image whose beam is not well sampled by the pixel size. Total flux can be lost when regridding such images, especially when the new pixel size is larger than the old pixel size. It is recommended to check the total flux of your input and output image, and if necessary rebin the input to have smaller pixels.

teuben commented 3 years ago

In my new one with the DC_locals.py - it also fails. Then 5 works, and 6 and 7 fail. Possibly related.

But in your next "git pull" you will need to run the configure script. I do this: ./configure --with-s4p-work=tmp and I have the GMC dataset in ../data from the scripts4paper directory, following the guidelines in ../data/README_DC2019_data

teuben commented 3 years ago

I've fixed up a dubious line in datacomb.py that caused this error. step 4 now works. 6 and 7 still fail. @LydiaMF

LydiaMF commented 3 years ago

Indeed, the line 1924 that it complains about needs to have the > ['value'] < removed.