If we label the regions as follows:

then the dual graph (not counting the exterior region) is shown here. The number of proper colorings is the value of the chromatic polynomial of the graph at 4. This software gives the chromatic polynomial as

Its value at x = 4 is 576.