10/02/23 Meeting notes - Sun subtraction follow Up (pre-peeling) #22
Replies: 10 comments 6 replies
-
No, wsclean will only subtract MODEL_DATA (and predict into MODEL_DATA). To keep things clean, if you want to manipulate extra columns, I suggest you use the https://github.com/o-smirnov/jovial/blob/reduce/jove-pol.yml#L335 So: image and deconvolve the Sun per-scan. Wsclean will leave a model for it in MODEL_DATA. Now you need to rephase that back to the original field centre (does chgcentre also update MODEL_DATA?) After rephasing, copy MODEL_DATA to SUN_MODEL_DATA. Then: Experiment 1: subtract SUN_MODEL_DATA from CORRECTED_DATA (into a new column) and re-image that column -- hopefully the aliases are reduced/gone Experiment 2: repeat selfcal on the field (without dEs), using MODEL_DATA for the field. Your QuartiCal model should look like MODEL_DATA:SUN_MODEL_DATA, and you should ask it to subtract direction 1. Hopefully, you get a better image of the field than in Experiment 1, because now selfcal is taking the Sun into account, and also subtracting it from the corrected data. Experiment 3: repeat selfcal on the field with solvable dEs -- so now you're solving for a dE towards the Sun. Subtract direction 1. Hopefully, you get an even better image of the field. @IanHeywood does that sound like a reasonable plan? |
Beta Was this translation helpful? Give feedback.
-
Sounds good to me.
I believe so, but I'm not sure if it will handle custom columns (e.g.
I think it's worth emphasising that the aliasing is an unexpected additional problem that we have encountered with this observation. Although it is definitely another aspect by which the Sun can have a detrimental effect on the science target, I think it's best if this is considered to be separate from the main issue, which is the removal of sidelobes from the Sun that are encroaching on the map centre. If the imager is setup correctly then aliasing shouldn't be a problem, whereas the main problem is where the R&D is at.
Won't you have to re-predict the field centre model into |
Beta Was this translation helpful? Give feedback.
-
Yes, form the Wsclean documentation, as soon we are not specifying -datacolumn in chgcentre to only phase-rotate the visibilities in the given column, the columns DATA, MODEL_DATA and CORRECTED_DATA will all be processed if they exist. When a measurement set contains multiple data columns (e.g., DATA, MODEL_DATA, CORRECTED_DATA), each column will be updated (as long as they have a standard name). |
Beta Was this translation helpful? Give feedback.
-
I'm pretty sure you can use wsclean to image any custom-named column. It's just the predict step that's hardwired for MODEL_DATA.
Hmmm, but when you -predict from a model image, does it not assume that the centre of the image corresponds to the current phase centre? I have a suspicion this is the case. If so, then we have to -predict the Sun using the rephased MSs, into MODEL_DATA, then phase-rotate them back to the field, then copy MODEL_DATA to SUN_MODEL_DATA, then -predict the field model into MODEL_DATA.
Completely agree. But in this case the aliases are actually a useful diagnostic for us, because their (hoped-for) disappearance will be an additional indication of the Sun being subtracted correctly.
Yep, see above.
Absolutely, I've tended to make the predict step separate in all my recipes now, precisely for that reason. Also makes it easier to resume a workflow from the predict step. |
Beta Was this translation helpful? Give feedback.
-
Indeed, but the question/answer here was about chgcentre!
Ah maybe. Is what you are pointing out here is an unavoidable issue when predicting using degridding? I guess if you're using crystalball / DFTs then you could predict the Sun into an MS that was shifted back to the phase centre. We should tread carefully here and make sure the software is doing what we want. |
Beta Was this translation helpful? Give feedback.
-
It's avoidable if the degridder knows to apply a phaseshift term. I guess we could check if Andre does that, but there's no design reason why he should (since he refuses to do it for gridding, providing chgcentre instead...) Crystallball, on the other hand, will predict correctly regardless. I guess @Victoria-Samboco should just check how many components she gets in the source list for the Sun -- if not too many, then just crystallball away and don't worry about image models. |
Beta Was this translation helpful? Give feedback.
-
Hi @Victoria-Samboco, here are some thoughts on your post above.
|
Beta Was this translation helpful? Give feedback.
-
This is looking great now, excellent work. :)
Did you mean radius here? |
Beta Was this translation helpful? Give feedback.
-
Beta Was this translation helpful? Give feedback.
-
Hoiw exactly are you subtracting it? I suggest you upload the recipe to a gist and link to it. There must be something wrong with the recipe -- if everything was down correctly, you should be getting results as good as your residual image in #22 (reply in thread). |
Beta Was this translation helpful? Give feedback.
-
According the meeting we had today, I will be pointing the progress of the Sun subtraction follow up for discussion, following the steps below:
-niter
set to 0-niter
different to zero. Deconvolution with wsclean in other to find a good model of the SUN.Q: Is there a way to specify that we want to subtract SUN_MODEL_DATA from CORRECTED_DATA if the -subtract-model parameter in wsclean use the default MODEL_DATA?
Beta Was this translation helpful? Give feedback.
All reactions