# Changing axis in polar seasonal plots with Cairomakie

**URL:** <https://discourse.julialang.org/t/changing-axis-in-polar-seasonal-plots-with-cairomakie/106108>\
**Category:** Visualization\
**Tags:** plotting, makie, cairomakie\
**Created:** [November 12, 2023, 8:31am UTC](https://discourse.julialang.org/t/changing-axis-in-polar-seasonal-plots-with-cairomakie/106108 "2023-11-12T08:31:17Z")\
**Posts on this page:** 7\
**Page:** 1

<div class="post-metadata">

**Author:** ![indymnv](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/indymnv/32/212594_2.png) [@indymnv](https://discourse.julialang.org/u/indymnv)\
**Post date:** [November 12, 2023, 8:31am UTC](https://discourse.julialang.org/t/changing-axis-in-polar-seasonal-plots-with-cairomakie/106108/1 "2023-11-12T08:31:17Z")

</div>

I am trying to replicate the seasonal plot with the polar axis from [this example](https://otexts.com/fpp2/seasonal-plots.html), but I couldn’t figure out how to replace the degree axis with months. This is the code:

Creating the dataframe

```julia
link = "https://vincentarelbundock.github.io/Rdatasets/csv/fpp2/a10.csv"
f = download(link)
a10 = DataFrame(CSV.File(f))
a10.year = parse.( Int64, first.(string.(a10.time),4))
a10.month .= 1
for (i, r) in enumerate(eachrow(a10))
	r.month += (i +5) % 12
end

```

should give this:

```julia
julia> first(a10, 5)
5×5 DataFrame
 Row │ rownames time value year month
     │ Int64 Float64 Float64 Int64 Int64
─────┼──────────────────────────────────────────
   1 │ 1 1991.5 3.52659 1991 7
   2 │ 2 1991.58 3.18089 1991 8
   3 │ 3 1991.67 3.25222 1991 9
   4 │ 4 1991.75 3.611 1991 10
   5 │ 5 1991.83 3.56587 1991 11

```

And now, trying to create the Plot.

```julia
psize = length(unique(a10.year))	
f = Figure()
ax =PolarAxis(f[1,1], title = "seasonal plot: " ,theta_0 = -pi/2, direction = -1, )
cbarPal = :thermal
cmap = cgrad(ColorSchemes.colorschemes[:viridis], psize, categorical = true)
x = length(unique(a10.month))
t = a10.month
tnorm = (t.-minimum(t))./(maximum(t).-minimum(t))*2*pi #normalize the axis related to X
a10.month_norm = tnorm #add to the dataframe
list_lin = [] #empty list to use in legend
for (i, value) in enumerate(unique(a10.year))
	a10_filtered = filter([:year] => (z) -> z == value , a10)
	y = a10_filtered.value
	x = a10_filtered.month_norm
	lin = lines!(ax ,x, y , label = string(value), color = cmap[i])
	push!(list_lin, lin)
end
Legend(f[1, 2], list_lin, string.(minimum(a10.year):maximum(a10.year)))

f

```

giving me this:

 ![image](https://global.discourse-cdn.com/julialang/original/3X/6/e/6ed9134111450b9401521da92fdc312c774eb608.png)

How can I replace the degrees on the circumference axes with months?

---

<div class="post-metadata">

**Author:** ![jules](https://avatars.discourse-cdn.com/v4/letter/j/41988e/32.png) [@jules](https://discourse.julialang.org/u/jules)\
**Post date:** [November 12, 2023, 11:39am UTC](https://discourse.julialang.org/t/changing-axis-in-polar-seasonal-plots-with-cairomakie/106108/2 "2023-11-12T11:39:56Z")

</div>

I’d assume you can pass a tuple `(radians, ticklabels)` to `thetaticks` [PolarAxis · Makie](https://docs.makie.org/stable/reference/blocks/polaraxis/#thetaticks)

---

<div class="post-metadata">

**Author:** ![indymnv](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/indymnv/32/212594_2.png) [@indymnv](https://discourse.julialang.org/u/indymnv)\
**Post date:** [November 12, 2023, 12:44pm UTC](https://discourse.julialang.org/t/changing-axis-in-polar-seasonal-plots-with-cairomakie/106108/3 "2023-11-12T12:44:11Z")

</div>

I have my doubts if I am doing it right because if I included in the definition of the axes considering the default value, I got just an error

```julia
julia> ax =PolarAxis(f[1,1], title = "seasonal plot: " ,theta_0 = -pi/2, direction = -1, thetathicks = (180 / pi, "°"))
ERROR: MethodError: no method matching initialize_block!(::PolarAxis; thetathicks::Tuple{Float64, String})

```

and with ` AngularTicks(180 / pi, "°")` get an error because AngularTicks is not defined

---

<div class="post-metadata">

**Author:** ![jules](https://avatars.discourse-cdn.com/v4/letter/j/41988e/32.png) [@jules](https://discourse.julialang.org/u/jules)\
**Post date:** [November 12, 2023, 12:47pm UTC](https://discourse.julialang.org/t/changing-axis-in-polar-seasonal-plots-with-cairomakie/106108/4 "2023-11-12T12:47:33Z")

</div>

`thetaticks` not `thetathicks` 🙂

And I’d guess you should pass `(range(0, 2pi, length = 13)[1:end-1], list_of_months)`

---

<div class="post-metadata">

**Author:** ![indymnv](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/indymnv/32/212594_2.png) [@indymnv](https://discourse.julialang.org/u/indymnv)\
**Post date:** [November 12, 2023, 12:58pm UTC](https://discourse.julialang.org/t/changing-axis-in-polar-seasonal-plots-with-cairomakie/106108/5 "2023-11-12T12:58:42Z")

</div>

That works fine! To sum up:

```julia
psize = length(unique(a10.year))	
f = Figure()
list_of_months = ["jan", "feb", "mar", "apr", "may", "jun", "jul", "aug", "sep", "oct", "nov", "dec"]
ax =PolarAxis(f[1,1], title = "seasonal plot Polar Axis" ,theta_0 = -pi/2, direction = -1, thetaticks =(range(0, 2pi, length = 13)[1:end-1], list_of_months))
#ax.ygridvisible = false
cbarPal = :thermal
cmap = cgrad(ColorSchemes.colorschemes[:viridis], psize, categorical = true)
x = length(unique(a10.month))
t = a10.month
tnorm = (t.-minimum(t))./(maximum(t).-minimum(t))*2*pi
a10.month_norm = tnorm
list_lin = []
for (i, value) in enumerate(unique(a10.year))
	a10_filtered = filter([:year] => (z) -> z == value , a10)
	y = a10_filtered.value
	x = a10_filtered.month_norm
	lin = lines!(ax ,x, y , label = string(value), color = cmap[i])
	push!(list_lin, lin)
end
Legend(f[1, 2], list_lin, string.(minimum(a10.year):maximum(a10.year)),)

f

```

 ![image](https://global.discourse-cdn.com/julialang/original/3X/f/8/f8fb40cedca9eeadf4890d7dfd5c0bc10fd27e2b.png)

Thank you so much 😄

---

<div class="post-metadata">

**Author:** ![jules](https://avatars.discourse-cdn.com/v4/letter/j/41988e/32.png) [@jules](https://discourse.julialang.org/u/jules)\
**Post date:** [November 12, 2023, 1:01pm UTC](https://discourse.julialang.org/t/changing-axis-in-polar-seasonal-plots-with-cairomakie/106108/6 "2023-11-12T13:01:20Z")

</div>

Nice, although something’s still wrong in your code because the data points don’t match the theta locations of the months. Probably an off-by-one error somewhere.

Ah yeah your final point is in January again, while it should be in December

---

<div class="post-metadata">

**Author:** ![indymnv](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/indymnv/32/212594_2.png) [@indymnv](https://discourse.julialang.org/u/indymnv)\
**Post date:** [November 12, 2023, 2:01pm UTC](https://discourse.julialang.org/t/changing-axis-in-polar-seasonal-plots-with-cairomakie/106108/7 "2023-11-12T14:01:43Z")

</div>

Yeah, it’s another problem that I have to face. In the example, the line is connected, so the transition becomes smooth. However, in my iteration, I am not connecting the points between December and January, which is making that part a bit more confusing for interpretation.