Commit a0b58e9f authored by Mohammad Akhlaghi's avatar Mohammad Akhlaghi

UDF1 to UDF9 regions defined and cropped

The UDF1 to UDF9 regions of the MUSE UDF-Mosiac image have been
defined and cropped from the main input image. However, the input
image was not aligned with the celestial coordinates (which ImageCrop
currently needs), so it was first aligned and then the regions were
cut.
parent 44b5b5b3
......@@ -65,7 +65,7 @@ include $(foreach m, preparations download input-cutouts , \
#
# These are targets that do not correspond to a file (phony in the
# Make terminology). The following rules define what they do.
.PHONY: all
.PHONY: all clean
......@@ -75,3 +75,5 @@ include $(foreach m, preparations download input-cutouts , \
# -----------
#
# To clean the outputs if necessary.
clean:
rm -rf $(BDIR)/*
# Reproduction pipeline
onlydirconf 1
#onlyversion 0.1
# Input:
hdu 0
# Output:
maxblankfrac 0.8
\ No newline at end of file
# Configuration files for programs
imgcrop = reproduce/config/gnuastro/astimgcrop.conf
imgwarp = reproduce/config/gnuastro/astimgwarp.conf
# UDF1 region
# ============
#
# Vertices of polygon to define region to use UDF1.
udf1raa = 53.160532
udf1deca = -27.774896
udf1rab = 53.147629
udf1decb = -27.764659
udf1rac = 53.159066
udf1decc = -27.753664
udf1rad = 53.171654
udf1decd = -27.763868
# The polygon option for ImageCrop used for these vertices.
udf1polygon = --polygon=$(udf1raa),$(udf1deca):$(udf1rab),$(udf1decb):$(udf1rac),$(udf1decc):$(udf1rad),$(udf1decd)
# UDF2 region
# ============
#
# Vertices of polygon to define region to use UDF2.
udf2raa = 53.174052
udf2deca = -27.785841
udf2rab = 53.161612
udf2decb = -27.775976
udf2rac = 53.173478
udf2decc = -27.764314
udf2rad = 53.185508
udf2decd = -27.774290
# The polygon option for ImageCrop used for these vertices.
udf2polygon = --polygon=$(udf2raa),$(udf2deca):$(udf2rab),$(udf2decb):$(udf2rac),$(udf2decc):$(udf2rad),$(udf2decd)
# UDF3 region
# ============
#
# Vertices of polygon to define region to use UDF3.
udf3raa = 53.187452
udf3deca = -27.797344
udf3rab = 53.174848
udf3decb = -27.786957
udf3rac = 53.186372
udf3decc = -27.775679
udf3rad = 53.199274
udf3decd = -27.786258
# The polygon option for ImageCrop used for these vertices.
udf3polygon = --polygon=$(udf3raa),$(udf3deca):$(udf3rab),$(udf3decb):$(udf3rac),$(udf3decc):$(udf3rad),$(udf3decd)
# UDF4 region
# ============
#
# Vertices of polygon to define region to use UDF4.
udf4raa = 53.148436
udf4deca = -27.786683
udf4rab = 53.136052
udf4decb = -27.776859
udf4rac = 53.146975
udf4decc = -27.766372
udf4rad = 53.159180
udf4decd = -27.776273
# The polygon option for ImageCrop used for these vertices.
udf4polygon = --polygon=$(udf4raa),$(udf4deca):$(udf4rab),$(udf4decb):$(udf4rac),$(udf4decc):$(udf4rad),$(udf4decd)
# UDF5 region
# ============
#
# Vertices of polygon to define region to use UDF5.
udf5raa = 53.160766
udf5deca = -27.797776
udf5rab = 53.149576
udf5decb = -27.788774
udf5rac = 53.161196
udf5decc = -27.777309
udf5rad = 53.172679
udf5decd = -27.786412
# The polygon option for ImageCrop used for these vertices.
udf5polygon = --polygon=$(udf5raa),$(udf5deca):$(udf5rab),$(udf5decb):$(udf5rac),$(udf5decc):$(udf5rad),$(udf5decd)
# UDF6 region
# ============
#
# Vertices of polygon to define region to use UDF6.
udf6raa = 53.175113
udf6deca = -27.808545
udf6rab = 53.162673
udf6decb = -27.798725
udf6rac = 53.173469
udf6decc = -27.787874
udf6rad = 53.186128
udf6decd = -27.798305
# The polygon option for ImageCrop used for these vertices.
udf6polygon = --polygon=$(udf6raa),$(udf6deca):$(udf6rab),$(udf6decb):$(udf6rac),$(udf6decc):$(udf6rad),$(udf6decd)
# UDF7 region
# ============
#
# Vertices of polygon to define region to use UDF7.
udf7raa = 53.135275
udf7deca = -27.799024
udf7rab = 53.122761
udf7decb = -27.788744
udf7rac = 53.134394
udf7decc = -27.777903
udf7rad = 53.146578
udf7decd = -27.787875
# The polygon option for ImageCrop used for these vertices.
udf7polygon = --polygon=$(udf7raa),$(udf7deca):$(udf7rab),$(udf7decb):$(udf7rac),$(udf7decc):$(udf7rad),$(udf7decd)
# UDF8 region
# ============
#
# Vertices of polygon to define region to use UDF8.
udf8raa = 53.148098
udf8deca = -27.810511
udf8rab = 53.136237
udf8decb = -27.800755
udf8rac = 53.148071
udf8decc = -27.789328
udf8rad = 53.160087
udf8decd = -27.798638
# The polygon option for ImageCrop used for these vertices.
udf8polygon = --polygon=$(udf8raa),$(udf8deca):$(udf8rab),$(udf8decb):$(udf8rac),$(udf8decc):$(udf8rad),$(udf8decd)
# UDF9 region
# ============
#
# Vertices of polygon to define region to use UDF9.
udf9raa = 53.162442
udf9deca = -27.821235
udf9rab = 53.149729
udf9decb = -27.810759
udf9rac = 53.161296
udf9decc = -27.799688
udf9rad = 53.174245
udf9decd = -27.810176
# The polygon option for ImageCrop used for these vertices.
udf9polygon = --polygon=$(udf9raa),$(udf9deca):$(udf9rab),$(udf9decb):$(udf9rac),$(udf9decc):$(udf9rad),$(udf9decd)
......@@ -30,19 +30,93 @@ $(indir): $(BDIR); mkdir $@
# Align the MUSE UDF region
# -------------------------
#
# The original MUSE UDF Mosaic region is not aligned with the
# celestial coordinates. To make things easier, we are first aligning
# them to ease the processing.
udfaligned = $(foreach f, $(filters), $(indir)/udf-muse-$(f).fits)
$(udfaligned): $(indir)/udf-muse-%.fits: \
$(MUSEINPUTS)/muse-udf-%.fits $(imgwarp)
astimgwarp $< --align -h1 -o$@
# From original MUSE images
# -------------------------
#
# Cutout the udf10 and udf regions
udf1-muse-cutouts = $(foreach f, $(filters), $(indir)/udf1-muse-$(f).fits)
$(udf1-muse-cutouts): $(indir)/udf1-muse-%.fits: $(indir)/udf-muse-%.fits \
$(imgcrop) reproduce/config/internal/vertices-udf1.mk \
| $(indir)
astimgcrop $(udf1polygon) $< -o$@
udf2-muse-cutouts = $(foreach f, $(filters), $(indir)/udf2-muse-$(f).fits)
$(udf2-muse-cutouts): $(indir)/udf2-muse-%.fits: $(indir)/udf-muse-%.fits \
$(imgcrop) reproduce/config/internal/vertices-udf2.mk \
| $(indir)
astimgcrop $(udf2polygon) $< -o$@
udf3-muse-cutouts = $(foreach f, $(filters), $(indir)/udf3-muse-$(f).fits)
$(udf3-muse-cutouts): $(indir)/udf3-muse-%.fits: $(indir)/udf-muse-%.fits \
$(imgcrop) reproduce/config/internal/vertices-udf3.mk \
| $(indir)
astimgcrop $(udf3polygon) $< -o$@
udf4-muse-cutouts = $(foreach f, $(filters), $(indir)/udf4-muse-$(f).fits)
$(udf4-muse-cutouts): $(indir)/udf4-muse-%.fits: $(indir)/udf-muse-%.fits \
$(imgcrop) reproduce/config/internal/vertices-udf4.mk \
| $(indir)
astimgcrop $(udf4polygon) $< -o$@
udf5-muse-cutouts = $(foreach f, $(filters), $(indir)/udf5-muse-$(f).fits)
$(udf5-muse-cutouts): $(indir)/udf5-muse-%.fits: $(indir)/udf-muse-%.fits \
$(imgcrop) reproduce/config/internal/vertices-udf5.mk \
| $(indir)
astimgcrop $(udf5polygon) $< -o$@
udf6-muse-cutouts = $(foreach f, $(filters), $(indir)/udf6-muse-$(f).fits)
$(udf6-muse-cutouts): $(indir)/udf6-muse-%.fits: $(indir)/udf-muse-%.fits \
$(imgcrop) reproduce/config/internal/vertices-udf6.mk \
| $(indir)
astimgcrop $(udf6polygon) $< -o$@
udf7-muse-cutouts = $(foreach f, $(filters), $(indir)/udf7-muse-$(f).fits)
$(udf7-muse-cutouts): $(indir)/udf7-muse-%.fits: $(indir)/udf-muse-%.fits \
$(imgcrop) reproduce/config/internal/vertices-udf7.mk \
| $(indir)
astimgcrop $(udf7polygon) $< -o$@
udf8-muse-cutouts = $(foreach f, $(filters), $(indir)/udf8-muse-$(f).fits)
$(udf8-muse-cutouts): $(indir)/udf8-muse-%.fits: $(indir)/udf-muse-%.fits \
$(imgcrop) reproduce/config/internal/vertices-udf8.mk \
| $(indir)
astimgcrop $(udf8polygon) $< -o$@
udf9-muse-cutouts = $(foreach f, $(filters), $(indir)/udf9-muse-$(f).fits)
$(udf9-muse-cutouts): $(indir)/udf9-muse-%.fits: $(indir)/udf-muse-%.fits \
$(imgcrop) reproduce/config/internal/vertices-udf9.mk \
| $(indir)
astimgcrop $(udf9polygon) $< -o$@
udf10-muse-cutouts = $(foreach f, $(filters), $(indir)/udf10-muse-$(f).fits)
$(udf10-muse-cutouts): $(indir)/udf10-muse-%.fits: \
$(MUSEINPUTS)/muse-udf10-%.fits $(imgcrop) | $(indir)
astimgcrop $(udf10polygon) --wcsmode $< -h1 -o$@
$(MUSEINPUTS)/muse-udf10-%.fits $(imgcrop) \
reproduce/config/internal/vertices-udf10.mk | $(indir)
astimgcrop $(udf10polygon) -h1 $< -o$@
udf-muse-cutouts = $(udf1-muse-cutouts) $(udf2-muse-cutouts) \
$(udf3-muse-cutouts) $(udf4-muse-cutouts) \
$(udf5-muse-cutouts) $(udf6-muse-cutouts) \
$(udf7-muse-cutouts) $(udf8-muse-cutouts) \
$(udf9-muse-cutouts) $(udf10-muse-cutouts)
file.txt: $(udf10-muse-cutouts)
file.txt: $(udf-muse-cutouts)
echo
echo "end of input-cutouts"
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment