Is there a technique to exactly calculate the Hausdorff dimension of the border of this fractal? 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?
 A: $\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$.
A: You may be able to calculate a graph-directed IFS for the boundary of the tile, though the rotated tiles may make this tricky - I'm not sure whether this is classed as a self-affine tile, or whether it needs a multi-tile system.
Given a graph-directed IFS, the Hausdorff dimension of the fractal is $s$ such that the spectral radius of a special matrix $M(s)$ (constructed from the structure of the GDIFS) is $1$.
I have some numerical code (implemented in Javascript, see the function hausdorffDimension() in the source code of https://mathr.co.uk/blog/2007-10-03_graphgrow.svg) that computes this.
The freeware (closed source) software available at https://ifstile.com can

  
*
  
*extract boundary of self-affine tiles as directed graph IFS
  
*compute dimension of the boundary of self-affine tiles (numerically and analytically)
  

which seems to be what you want - I haven't tried it myself yet.  It probably uses the methods from the last reference (by the author of IFSTile).
References:


*

*R. D. Mauldin & S. C. Williams, Hausdorff dimension in graph directed constructions, Trans. Amer. Math. Soc. 309 (1988), 811–829.
https://www.ams.org/journals/tran/1988-309-02/S0002-9947-1988-0961615-4/S0002-9947-1988-0961615-4.pdf

*Robert S. Strichartz and Yang Wang, Geometry of Self–Affine Tiles I, Indiana University Mathematics Journal, Vol. 48, No. 1 (Spring, 1999), pp. 1-23 
https://www.jstor.org/stable/24900135

*Richard Kenyon, Jie Li, Robert S. Strichartz and Yang Wang, Geometry of Self–Affine Tiles II,  Indiana University Mathematics Journal, Vol. 48, No. 1 (Spring, 1999), pp. 25-42
https://www.jstor.org/stable/24900136

*Dmitry Mekhontsev An algebraic framework for finding and analyzing self-affine tiles and fractals, Dr. rer. nat. thesis (2018)
https://epub.ub.uni-greifswald.de/frontdoor/deliver/index/docId/2479/file/MekhontsevElectronic.pdf
A: You can definitely use IFStile package to define the tile, extract the boundary and calculate the boundary dimension.
Here the definition of the tile for IFStile:
########################################################
@G
$dim=2
s=$companion([1,0])
h0=[-3,0]*s
h1=[-2,0]*1
h2=[-1,0]*s
h3=[-1,2]*s
h4=[-1,1]*s^2
h5=[-1,-1]*s^2
h6=[0,1]*s^3
h7=[1,1]*s^2
h8=[1,0]*s
A=3^-1*(h0|h1|h2|h3|h4|h5|h6|h7|h8)*A
########################################################

The boundary is a graph-directed IFS

Hausdorff dimension: 
log(x)/log(3)=1.565033324743844....

where 
x^3-3x^2-11x-19=0

