Let
$$I = \int_{-\pi/4}^{\pi/4} x \, \text{cosec} \, x \, dx = 2 \int_0^{\pi/4} x \, \text{cosec} \, x \, dx.$$
After integrating by parts we have
\begin{align}
I &= -\frac{\pi}{2} \ln (1 + \sqrt{2}) + 2 \int_0^{\pi/4} \ln (\text{cosec} \, x + \cot x) \, dx\\
&= -\frac{\pi}{2} \ln (1 + \sqrt{2}) + 2 \int_0^{\pi/4} \ln (1 + \cos x) \, dx - 2 \int_0^{\pi/4} \ln (\sin x) \, dx.
\end{align}
The second of these integrals is perhaps (?) reasonably well known (for an evaluation, see here). The result is:
$$\int_0^{\pi/4} \ln (\sin x) \, dx = -\frac{1}{2} \mathbf{G} - \frac{\pi}{4} \ln 2,$$
where $\mathbf{G}$ is Catalan's constant. Thus
$$I = -\frac{\pi}{2} \ln (1 + \sqrt{2}) + \mathbf{G} + \frac{\pi}{2} \ln 2 + 2I_1.$$
For the first of the integrals we will make use of the following formula, a proof of which can be found here
$$\ln (1 + \cos x) = 2 \sum_{n = 1}^\infty (-1)^{n + 1} \frac{\cos (nx)}{n} - \ln 2.$$
So
\begin{align}
I_1 &= \int_0^{\pi/4} \ln (1 + \cos x) \, dx\\
&= 2 \sum_{n = 1}^\infty \frac{(-1)^{n + 1}}{n} \int_0^{\pi/4} \cos (nx) \, dx - \ln 2 \int_0^{\pi/4} dx\\
&= -\frac{\pi}{4} \ln 2 + 2 \sum_{n = 1}^\infty \frac{(-1)^{n + 1}}{n^2} \sin \left (\frac{n \pi}{4} \right )\\
&= -\frac{\pi}{4} \ln 2 + 2 \sum_{\stackrel{{\Large{n = 1}}}{n \in \text{odd}}}^\infty \frac{(-1)^{n + 1}}{n^2} \sin \left (\frac{n \pi}{4} \right ) + + 2 \sum_{\stackrel{{\Large{n = 1}}}{n \in \text{even}}}^\infty \frac{(-1)^{n + 1}}{n^2} \sin \left (\frac{n \pi}{4} \right )\\
&= -\frac{\pi}{4} \ln 2 + 2 \sum_{n = 1}^\infty \frac{1}{(2n - 1)^2} \sin \left [\frac{\pi}{4} (2n - 1) \right ] - \frac{1}{2} \sum_{n = 1}^\infty \frac{1}{n^2} \sin \left (\frac{n \pi}{2} \right )\\
&= -\frac{\pi}{4} \ln 2 + 2 S_1 - \frac{1}{2} S_2.
\end{align}
For the second of these sums,
\begin{align}
S_2 &= \sum_{\stackrel{{\Large{n = 1}}}{n \in \text{odd}}}^\infty \frac{1}{n^2} \sin \left (\frac{\pi n}{2} \right ) + \sum_{\stackrel{{\Large{n = 1}}}{n \in \text{even}}} \frac{1}{n^2} \sin \left (\frac{\pi n}{2} \right )\\
&= \sum_{n = 0}^\infty \frac{(-1)^n}{(2n + 1)^2}\\
&= \mathbf{G},
\end{align}
after a shift of the index $n \mapsto 2n + 1$ in the odd sum has been made while the even sum is identically equal to zero.
For the first of the sums, as it converges absolutely we can split it up as folows:
$$S_1 = \sum_{\stackrel{{\Large{n = 1}}}{n \in 1,5,9,\ldots}}^\infty \frac{1}{(2n - 1)^2} \sin \left [\frac{\pi}{4} (2k - 1) \right ] + \sum_{\stackrel{{\Large{n = 1}}}{n \in 2,6,10,\ldots}}^\infty \frac{1}{(2n - 1)^2} \sin \left [\frac{\pi}{4} (2k - 1) \right ] + \sum_{\stackrel{{\Large{n = 1}}}{n \in 3,7,11,\ldots}}^\infty \frac{1}{(2n - 1)^2} \sin \left [\frac{\pi}{4} (2k - 1) \right ] + \sum_{\stackrel{{\Large{n = 1}}}{n \in 4,8,12,\ldots}}^\infty \frac{1}{(2n - 1)^2} \sin \left [\frac{\pi}{4} (2k - 1) \right ].$$
Shifting the indices as follows: $n \mapsto 4n - 3, n \mapsto 4n - 2, n \mapsto 4n - 1, n \mapsto 4n$ leads to
\begin{align}
S_1 &= \frac{1}{\sqrt{2}} \left [\sum_{n = 1}^\infty \frac{1}{(8n - 7)^2} + \sum_{n = 1}^\infty \frac{1}{(8n - 5)^2} - \sum_{n = 1}^\infty \frac{1}{(8n - 3)^2} - \sum_{n = 1}^\infty \frac{1}{(8n - 1)^2} \right ]\\
&= \frac{1}{\sqrt{2}} \left [\sum_{n = 0}^\infty \frac{1}{(8n + 1)^2} + \sum_{n = 0}^\infty \frac{1}{(8n + 3)^2} - \sum_{n = 0}^\infty \frac{1}{(8n + 5)^2} - \sum_{n = 0}^\infty \frac{1}{(8n + 7)^2} \right ]\\
&= \frac{1}{64 \sqrt{2}} \left [\sum_{n = 0}^\infty \frac{1}{(n + 1/8)^2} + \sum_{n = 0}^\infty \frac{1}{(n + 3/8)^2} - \sum_{n = 0}^\infty \frac{1}{(n + 5/8)^2} - \sum_{n = 0}^\infty \frac{1}{(n + 7/8)^2} \right ]\\
&= \frac{1}{64 \sqrt{2}} \left [\psi^{(1)} \left (\frac{1}{8} \right ) + \psi^{(1)} \left (\frac{3}{8} \right ) - \psi^{(1)} \left (\frac{5}{8} \right ) - \psi^{(1)} \left (\frac{7}{8} \right ) \right ],
\end{align}
where we have made use of the series representation for the polygamma function of order one (as known as the trigamma function). So the value for $I_1$ is:
$$I_1 = -\frac{\pi}{4} \ln 2 - \frac{1}{2} \mathbf{G} + \frac{1}{32 \sqrt{2}} \left [\psi^{(1)} \left (\frac{1}{8} \right ) + \psi^{(1)} \left (\frac{3}{8} \right ) - \psi^{(1)} \left (\frac{5}{8} \right ) - \psi^{(1)} \left (\frac{7}{8} \right ) \right ],$$
leading to a final result of
$$\int_{-\pi/4}^{\pi/4} x \, \text{cosec} \, x \, dx = -\frac{\pi}{2} \ln (1 + \sqrt{2}) + \frac{1}{16 \sqrt{2}} \left [\psi^{(1)} \left (\frac{1}{8} \right ) + \psi^{(1)} \left (\frac{3}{8} \right ) - \psi^{(1)} \left (\frac{5}{8} \right ) - \psi^{(1)} \left (\frac{7}{8} \right ) \right ].$$
Converting the trigamma functions to Clausen functions of order two
Inspired by the answer given by Zacky in terms of the Clausen function of order two, here I will show how to convert my answer in terms of the 4 trigamma functions into 2 Clausen functions of order 2.
The relation between the Clausen function of order two and the trigamma function is given by (a proof of this can be found here)
$$\text{Cl}_2 \left (\frac{q \pi}{p} \right ) = \frac{1}{(2p)^{2m} (2m - 1)!} \sum_{n = 1}^p \sin \left (\frac{qn\pi}{p} \right ) \left [\psi^{(1)} \left (\frac{n}{2p} \right ) + (-1)^q \psi^{(1)} \left (\frac{n + p}{2p} \right ) \right ].$$
Setting $m = 1, q = 1, p = 4$ gives
\begin{align}
\text{Cl}_2 \left (\frac{\pi}{4} \right ) &= \frac{1}{64} \left [\frac{1}{\sqrt{2}} \left \{\psi^{(1)} \left (\frac{1}{8} \right ) - \psi^{(1)} \left (\frac{5}{8} \right ) \right \} + \psi^{(1)} \left (\frac{1}{4} \right ) - \psi^{(1)} \left (\frac{3}{4} \right ) \right.\\
& \qquad \left. + \frac{1}{\sqrt{2}} \left \{\psi^{(1)} \left (\frac{3}{8} \right ) - \psi^{(1)} \left (\frac{7}{8} \right ) \right \} \right ], \qquad (*)
\end{align}
and setting $m = 1, q = 3, p = 4$ gives
\begin{align}
\text{Cl}_2 \left (\frac{3\pi}{4} \right ) &= \frac{1}{64} \left [\frac{1}{\sqrt{2}} \left \{\psi^{(1)} \left (\frac{1}{8} \right ) - \psi^{(1)} \left (\frac{5}{8} \right ) \right \} - \psi^{(1)} \left (\frac{1}{4} \right ) + \psi^{(1)} \left (\frac{3}{4} \right ) \right.\\
& \qquad \left. + \frac{1}{\sqrt{2}} \left \{\psi^{(1)} \left (\frac{3}{8} \right ) - \psi^{(1)} \left (\frac{7}{8} \right ) \right \} \right ]. \qquad (**)
\end{align}
On adding ($*$) to ($**$) we see that
$$\psi^{(1)} \left (\frac{1}{8} \right ) + \psi^{(1)} \left (\frac{3}{8} \right ) - \psi^{(1)} \left (\frac{5}{8} \right ) - \psi^{(1)} \left (\frac{7}{8} \right ) = 32 \sqrt{2} \left [\text{Cl}_2 \left (\frac{\pi}{4} \right ) + \text{Cl}_2 \left (\frac{3\pi}{4} \right ) \right ],$$
giving
$$\int_{-\pi/4}^{\pi/4} \frac{x}{\sin x} \, dx = -\frac{\pi}{2} \ln (1 + \sqrt{2}) + 2 \, \text{Cl}_2 \left (\frac{\pi}{4} \right ) + 2 \, \text{Cl}_2 \left (\frac{3\pi}{4} \right ).$$
\begin{align}\sum_{n=0}^\infty \frac{\binom{2n}{n}^2}{4^{2n+1}(2n+1)}=\frac{\text{G}}{\pi}\tag1\end{align}
(see p81, Deriving Forsyth-Glaisher type series for $\frac{1}{\pi}$ and Catalan's constant by an elementary method. )
From the same source,
\begin{align}\sum_{n=0}^\infty \frac{\binom{2n}{n}^2}{16^n(2n+3)}=\frac{\text{G}}{\pi}+\frac{1}{2\pi}\tag2\end{align}
ADDENDUM:
Proof for (1),
It is well known that for $n\geq 0$ integer,
\begin{align}\int_0^{\frac{\pi}{2}}\cos^{2n} x\,dx=\frac{\pi}{2}\cdot\frac{\binom{2n}{n}}{4^n}\end{align}
(Wallis formula)
Therefore for $n\geq 0$ integer,
\begin{align}\frac{\binom{2n}{n}^2\pi^2}{4^{2n+1}(2n+1)}=\int_0^1 \left(\int_0^\infty \int_0^\infty t^{2n}\cos^{2n}x \cos^{2n}y \,dx\,dy \right)\,dt\end{align}
therefore,
\begin{align}\pi^2\sum_{n=0}^{\infty}\frac{\binom{2n}{n}^2}{4^{2n+1}(2n+1)}&=\sum_{n=0}^{\infty}\left(\int_0^1 \left(\int_0^{\frac{\pi}{2}} \int_0^{\frac{\pi}{2}} t^{2n}\cos^{2n}x \cos^{2n}y \,dx\,dy \right)\,dt\right)\\
&=\int_0^1 \left(\int_0^{\frac{\pi}{2}} \int_0^{\frac{\pi}{2}} \left(\sum_{n=0}^{\infty}t^{2n}\cos^{2n}x \cos^{2n}y\right) \,dx\,dy \right)\,dt\\
&=\int_0^1 \left(\int_0^{\frac{\pi}{2}} \int_0^{\frac{\pi}{2}} \frac{1}{1-t^2\cos^2 x\cos^2 y}\,dx\,dy \right)\,dt\\
\end{align}
Perform the change of variable $u=\tan x$,$v=\tan y$,
\begin{align}\pi^2\sum_{n=0}^{\infty}\frac{\binom{2n}{n}^2}{4^{2n+1}(2n+1)}&=
\int_0^1 \left(\int_0^{\infty} \int_0^{\infty}\frac{1}{(1+u^2)(1+v^2)-t^2}\,du\,dv \right)\,dt\\
&=\int_0^1 \left(\int_0^\infty \frac{1}{\sqrt{1+v^2}}\left[\frac{\arctan\left(\frac{u\sqrt{1+v^2}}{\sqrt{1+v^2-t^2}}\right)}{\sqrt{1+v^2-t^2}}\right]_{u=0}^{u=\infty}\,dv\right)\,dt\\
&=\frac{\pi}{2}\int_0^1 \left(\int_0^\infty \frac{1}{\sqrt{1+v^2}\sqrt{1+v^2-t^2}}\,dv\right)\,dt\\
&=\frac{\pi}{2}\int_0^\infty \frac{1}{\sqrt{1+v^2}}\left[\arctan\left(\frac{t}{\sqrt{1+v^2-t^2}}\right)\right]_{t=0}^{t=1}\,dv\\
&=\frac{\pi}{2}\int_0^\infty \frac{\arctan\left(\frac{1}{v}\right)}{\sqrt{1+v^2}}\,dv\\
\end{align}
Perform the change of variable $y=\dfrac{1}{x}$,
\begin{align}\pi^2\sum_{n=0}^{\infty}\frac{\binom{2n}{n}^2}{4^{2n+1}(2n+1)}&=\frac{\pi}{2}\int_0^\infty \frac{\arctan x}{x\sqrt{1+x^2}}\,dx\\
\end{align}
Perform the change of variable $y=\arctan x$,
\begin{align}\pi^2\sum_{n=0}^{\infty}\frac{\binom{2n}{n}^2}{4^{2n+1}(2n+1)}&=\frac{\pi}{2}\int_0^{\frac{\pi}{2}}\frac{x}{\sin x}
\,dx\\
&=\frac{\pi}{2}\Big[x\ln\left(\tan\left(\frac{x}{2}\right)\right)\Big]_0^{\frac{\pi}{2}}-\frac{\pi}{2}\int_0^{\frac{\pi}{2}}\ln\left(\tan\left(\frac{x}{2}\right)\right)\,dx\\
&=-\frac{\pi}{2}\int_0^{\frac{\pi}{2}}\ln\left(\tan\left(\frac{x}{2}\right)\right)\,dx\\
\end{align}
Perform the change of variable $y=\frac{x}{2}$,
\begin{align}\pi^2\sum_{n=0}^{\infty}\frac{\binom{2n}{n}^2}{4^{2n+1}(2n+1)}&=
-\pi\int_0^{\frac{\pi}{4}}\ln(\tan x)\,dx\\
&=\pi\times \text{G}\\
\end{align}
Therefore,
\begin{align}\boxed{\sum_{n=0}^\infty \frac{\binom{2n}{n}^2}{4^{2n+1}(2n+1)}=\frac{\text{G}}{\pi}}\end{align}
Best Answer
We have that, for a sufficient function $f(x)$, $$f(x)=\prod_{f(\omega)=0}(x-\omega)$$ so that $$\frac{1}{f(x)}=\prod_{f(\omega)=0}\frac{1}{x-\omega}.$$ We then assume that we may write $$\frac{1}{f(x)}=\sum_{f(\omega)=0}\frac{b(\omega)}{x-\omega}.$$ We can show that $$b(\omega)=\prod_{f(r)=0\ ,\ r\ne\omega}\frac{1}{\omega-r}\ .$$ But on the other hand, $$\ln f(x)=\sum_{f(\omega)=0}\ln(x-\omega)$$ so that $$f'(x)=\sum_{f(\omega)=0}\frac{f(x)}{x-\omega}=\sum_{f(\omega)=0}\ \prod_{f(r)=0\ ,\ r\ne\omega}(x-r)\ .$$ Thus, for any $q$ with $f(q)=0$, we plug in $x=q$ to get $$f'(q)=\prod_{f(r)=0\ ,\ r\ne q}(q-r)=\frac{1}{b(q)}\ .$$ Which gives $$\frac{1}{f(x)}=\sum_{f(\omega)=0}\frac{1}{(x-\omega)f'(\omega)}\ .$$ Choosing $f(x)=\sin x$ implies that $$\frac{\omega}{\pi}\in\Bbb Z$$ and $$f'(\omega)=\cos\omega=(-1)^k,$$ so we end up with $$\frac{1}{\sin x}=\sum_{k\in\Bbb Z}\frac{(-1)^k}{x+\pi k}$$ for $x/\pi\not\in\Bbb Z$. This expression is equivalent to $$\frac{x}{\sin x}=\sum_{k\in \Bbb N}\frac{\pi^2k^2}{\pi^2k^2-x^2}\prod_{k\ne j\in\Bbb N}\frac{j^2}{j^2-k^2}\ .$$