I came up with a self-similar fractal that fits into itself like a jigsaw puzzle. While the surface area is clearly 2-dimensional, it is unclear to me how to compute the Hausdorff dimension of the border. I could use the box counting method to find an approximation. But since I know how to construct the fractal, I was wondering if it were possible to exactly calculate it.
Above you can see the fractal I am talking about. I have divided it into colored areas to make its construction more clear. The red area is just a scaled down version (with scale factor 1/3) of the whole thing. The other colored areas are scaled down versions (with scale factor 1/3) as well, but are also rotated. The blue areas are rotated 90 degrees anti-clockwise. The green areas are rotated 180 degrees. The yellow area is rotated 90 degrees clockwise. The center of the fractal is inside the big gap area, which could be filled with the red area. If this gap area gets filled recursively with smaller and smaller red areas, then the center will remain outside the fractal after any finite amount of iterations. It may not look like it, but the 9 colored areas are placed in such a way that their centers form a grid pattern. To make the grid pattern more clear I've added the image below.
Unfortunately, I don't know how to give a formal mathematical definition of the fractal, so I hope the drawings I made and the description are clear enough. I want to make it clear that it's not important to me to know the Hausdorff dimension for this specific case, but rather how to generally calculate the Hausdorff dimension of the border of well defined closed areas. Is there some sort of technique for this?



$\DeclareMathOperator\unit{unit}$Let's start with the dimension $\Delta$ of the area, which we already know to be $\Delta=2$. Let each of the colored squares have an area of $1\text{ unit}^\Delta$. Then the area of the template is: $$A_0=9\cdot\unit^\Delta$$ Now let's get more fine grained measuring equipment that is $\eta$ times more precise. That is, we will be measuring in $(\frac 1\eta\cdot\unit)^\Delta$ instead of $\unit^\Delta$. As it is we will find $9$ smaller squares inside each of the $9$ squares, so the total area is: $$A_1 = 9^2\cdot\left(\frac 1\eta\cdot\unit\right)^\Delta=9^2\cdot\eta^{-\Delta}\cdot\unit^\Delta$$
With the correct dimension $\Delta=2$, we must have $A=A_0=A_1$, and therefore: $$9=9^2\cdot\eta^{-\Delta}\quad\Rightarrow\quad \eta^\Delta=9\quad\Rightarrow\quad \eta^2=9\quad\Rightarrow\quad\eta=3$$ We measure on a $3$ times smaller scale in $1$ iteration.
Let $D$ be the dimension of the circumference. The circumference with the same $\unit$ is: $$C_0 = 20\cdot\unit^D$$ When we measure again with $\frac 1\eta\cdot\unit$ we find: $$C_1 = 112\cdot\left(\frac 1\eta\cdot\unit\right)^D = 112\cdot\eta^{-D}\cdot\unit^D$$
We get with $C=C_0=C_1$: $$20=112\cdot\eta^{-D}= 112\cdot 3^{-D}\quad\Rightarrow\quad 3^D=\frac{112}{20}\quad\Rightarrow\quad D=\log_3 \frac{112}{20}\approx 1.568$$ As expected this is a bit higher than Koch's curve that has $D=1.268$ and lower than $2$.