JWST image with scale bar

This tutorial shows how to download one of the new famous JWST release images of the Carina nebula, plot it with world coordinates (RA and DEC), and add a scalebar. Let's get started!

Packages

  • Downloads: we'll use this built-in Julia standard library to download the raw image data
  • AstroImages: we'll use this package to load and display the image
  • Plots: we'll use this package to display coordinates along the image and add the scalebar

You can install the necessary packages by running Julia, and typing ] to enter Pkg-mode. Then: add AstroImages Plots. Alternatively, you can run using Pkg; Pkg.add(["AstroImages", "Plots"]).

If you will be using these tools as part of a bigger project, it's strongly recommended to create a Julia Project to record package versions. If you're just experimenting, you can create a temporary project by running ] activate --temp.

If you're using Pluto notebooks, installing and recording package versions in a project are handled for you automatically.

Downloading the data

All JWST data is available from the MAST archive.

A mirror of the JWST initial release data is also hosted on AWS. We'll use this to directly download a specific famous image by URL.

Large Download

The size of this image is over 3GB! Ensure you have enough space on your hard drive before running the following code.

using Downloads

fname = Downloads.download(
   "https://stpubdata-jwst.stsci.edu/ero/jw02731/L3/t/"*
   "jw02731-o001_t017_nircam_clear-f187n_i2d.fits"
)

Loading the image

We can use AstroImages.jl to load the image.

using AstroImages

carina_full = load(fname); # semi-colon prevents displaying the image
[ Info: Image was loaded from HDU 2

Let's check the resolution of the image:

size(carina)
(14436, 8568)

Wow, that's over 120 megapixels!

AstroImages will fall back to downscaling it automatically before displaying it using ImageTransformations.restrict, but let's make our lives simpler and just use every 10th pixel:

carina = carina_full[begin:10:end, begin:10:end]

Carina nebula displayed as an image

Since we didn't put a semi-colon, the image was displayed. This works automatically if you use the VS Code editor, Jupyter notebooks, or Pluto Notebooks. If you're using a terminal, you can use ElectronDisplay.jl or ImageView.jl to open an interactive window.

imview is the function that was called automatically to display the image, but we can call it ourselves if we want to customize the display. We can make the image pop a little more by adjusting the display limits, or clims. These can be either a tuple of min and max values (i.e. (0.0, 100.0)) or a function to calculate them for us based on some criteria. We'll adjust the limits to remove the top and bottom 1% of all pixels:

imview(carina; clims=Percent(98))

Carina nebula displayed as an image with higher contrast

In Julia, keyword arguments come after positional arguments and they can be separated by a semi-colon.

Adding World Coordinates

The image header already contains WCS keys to specify the location of the image in the sky. Let's take a peek at one of these header keys:

println(carina["CTYPE1", Comment])
println(carina["CTYPE1"])
first axis coordinate type
RA---TAN

This shows that the coordinate type for the first axis is "RA–-TAN".

Now, we'll combine AstroImages with the Plots.jl library to display these coordinates.

using Plots

implot(carina; clims=Percent(98))

Carina nebula displayed in a plot

The image is tilted compared to the right-ascension and declination celestial coordinate system, so the grid is rotated.

The colourbar on the right displays the units of the image: Mega-Janskys per steradian: units of energy per collecting area, per spectral bandwidth, per solid-angle on the sky.

Let's turn off the grid for now:

implot(carina; grid=false, clims=Percent(98))

Carina nebula displayed in a plot without grid lines

Full documentation for the implot function is available in the AstroImages.jl docs.

Adding a Scalebar

Another way we can communicate the scale of this image is by adding a scalebar. We'll need to compute the length of a 2D vector (i.e. using the Pythagorean theorem) but we use the LinearAlgebra standard library function norm here for convenience.

using LinearAlgebra

# Start the scalebar here in pixel coordinates
start_coord_pix = [1000, 100]

# Convert to world coordinates
start_coord_world = pix_to_world(carina_full, start_coord_pix)

# Extend by 10 arcminutes of declination
stop_coord_world  = start_coord_world .+ [0, 1/60]

# Convert back to pixel coordinates
stop_coord_pix = world_to_pix(carina_full, stop_coord_world)

# Measure how long it is here in pixel coordinates
# note: in theory this depends on where in the image we make this calculation
# because the coordinate system is warped
arcmin_px = norm(stop_coord_pix .- start_coord_pix)

# Plot our image again
implot(carina; grid=false, clims=Percent(98))

# Overplot the scalebar
plot!(
    # X
    [11000, 11000+arcmin_px],
    # Y
    [1000,  1000];
    # Formatting:
    color=:white,
    linewidth=5,
    label=""
)

# add annotation
annx = 11000+arcmin_px/2
anny = 1000+100
annotate!(
    (annx, anny, text(" 1' ", 10, :white, :center, :bottom))
)

Carina nebula displayed with a scale bar

Looking good! If we only want to communiate the scale and not the location, we can hide the tick marks and colorbar for a cleaner look:


# Hide decorations
implot(
    carina[10:end-10, 10:end-10]; # crop slightly
    grid=false,
    ticks=false,
    colorbar=false,
    clims=Percent(98),
    xlabel="",
    ylabel="",
    framestyle=:none,
    background_outside=:transparent,
)

# Overplot the scalebar
plot!(
    # X
    [11000, 11000+arcmin_px],
    # Y
    [1000,  1000];
    # Formatting:
    color=:white,
    linewidth=5,
    label="",
)

# add annotation
annx = 11000+arcmin_px/2
anny = 1000+100
annotate!(
    (annx, anny, text(" 1' ", 10, :white, :center, :bottom)) 
)

Carina nebula displayed with a scale bar and no other plot decorations

To save the image, just run savefig("output.png"). pdf, svg, and other file formats are also supported.

To learn more about displaying images and world coordinates, see the AstroImages.jl documentation.

For more on plotting in general, see the Plots.jl documentation.