The above four-element sets of black subsquares of a 4×4 square array
are 15 of the 60 Göpel tetrads , and 20 of the 80 Rosenhain tetrads , defined
by R. W. H. T. Hudson in his 1905 classic Kummer's Quartic Surface .
Hudson did not view these 35 tetrads as planes through the origin in a finite
affine 4-space (or, equivalently, as lines in the corresponding finite projective
3-space).
In order to view them in this way, one can view the tetrads as derived,
via the 15 two-element subsets of a six-element set, from the 16 elements
of the binary Galois affine space pictured above at top left.
This space is formed by taking symmetric-difference (Galois binary)
sums of the 15 two-element subsets, and identifying any resulting four-
element (or, summing three disjoint two-element subsets, six-element)
subsets with their complements. This process was described in my note
"The 2-subsets of a 6-set are the points of a PG(3,2)" of May 26, 1986.
The space was later described in the following —