hfe.py 4.6 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133
  1. # greaseweazle/image/hfe.py
  2. #
  3. # Written & released by Keir Fraser <keir.xen@gmail.com>
  4. #
  5. # This is free and unencumbered software released into the public domain.
  6. # See the file COPYING for more details, or visit <http://unlicense.org>.
  7. import struct
  8. from greaseweazle import error
  9. from greaseweazle.track import MasterTrack, RawTrack
  10. from bitarray import bitarray
  11. from .image import Image
  12. class HFE(Image):
  13. def __init__(self):
  14. self.bitrate = 250 # XXX real bitrate?
  15. # Each track is (bitlen, rawbytes).
  16. # rawbytes is a bytes() object in little-endian bit order.
  17. self.to_track = dict()
  18. @classmethod
  19. def from_file(cls, name):
  20. with open(name, "rb") as f:
  21. dat = f.read()
  22. (sig, f_rev, n_cyl, n_side, t_enc, bitrate,
  23. _, _, _, tlut_base) = struct.unpack("<8s4B2H2BH", dat[:20])
  24. error.check(sig != b"HXCHFEV3", "HFEv3 is not supported")
  25. error.check(sig == b"HXCPICFE" and f_rev <= 1, "Not a valid HFE file")
  26. error.check(0 < n_cyl, "HFE: Invalid #cyls")
  27. error.check(0 < n_side < 3, "HFE: Invalid #sides")
  28. error.check(bitrate != 0, "HFE: Invalid bitrate")
  29. hfe = cls()
  30. hfe.bitrate = bitrate
  31. tlut = dat[tlut_base*512:tlut_base*512+n_cyl*4]
  32. for cyl in range(n_cyl):
  33. for side in range(n_side):
  34. offset, length = struct.unpack("<2H", tlut[cyl*4:(cyl+1)*4])
  35. todo = length // 2
  36. tdat = bytes()
  37. while todo:
  38. d_off = offset*512 + side*256
  39. d_nr = 256 if todo > 256 else todo
  40. tdat += dat[d_off:d_off+d_nr]
  41. todo -= d_nr
  42. offset += 1
  43. hfe.to_track[cyl,side] = (len(tdat)*8, tdat)
  44. return hfe
  45. def get_track(self, cyl, side):
  46. if (cyl,side) not in self.to_track:
  47. return None
  48. bitlen, rawbytes = self.to_track[cyl,side]
  49. tdat = bitarray(endian='little')
  50. tdat.frombytes(rawbytes)
  51. track = MasterTrack(
  52. bits = tdat[:bitlen],
  53. time_per_rev = bitlen / (2000*self.bitrate))
  54. return track
  55. def emit_track(self, cyl, side, track):
  56. raw = RawTrack(clock = 5e-4 / self.bitrate, data = track)
  57. bits, _ = raw.get_revolution(0)
  58. bits.bytereverse()
  59. self.to_track[cyl,side] = (len(bits), bits.tobytes())
  60. def get_image(self):
  61. n_side = 1
  62. n_cyl = max(self.to_track.keys(), default=(0), key=lambda x:x[0])[0]
  63. n_cyl += 1
  64. # We dynamically build the Track-LUT and -Data arrays.
  65. tlut = bytearray()
  66. tdat = bytearray()
  67. # Stuff real data into the image.
  68. for i in range(n_cyl):
  69. s0 = self.to_track[i,0] if (i,0) in self.to_track else None
  70. s1 = self.to_track[i,1] if (i,1) in self.to_track else None
  71. if s0 is None and s1 is None:
  72. # Dummy data for empty cylinders. Assumes 300RPM.
  73. nr_bytes = 100 * self.bitrate
  74. tlut += struct.pack("<2H", len(tdat)//512 + 2, nr_bytes)
  75. tdat += bytes([0x88] * (nr_bytes+0x1ff & ~0x1ff))
  76. else:
  77. # At least one side of this cylinder is populated.
  78. if s1 is not None:
  79. n_side = 2
  80. bc = [s0 if s0 is not None else (0,bytes()),
  81. s1 if s1 is not None else (0,bytes())]
  82. nr_bytes = max(len(t[1]) for t in bc)
  83. nr_blocks = (nr_bytes + 0xff) // 0x100
  84. tlut += struct.pack("<2H", len(tdat)//512 + 2, 2 * nr_bytes)
  85. for b in range(nr_blocks):
  86. for t in bc:
  87. slice = t[1][b*256:(b+1)*256]
  88. tdat += slice + bytes([0x88] * (256 - len(slice)))
  89. # Construct the image header.
  90. header = struct.pack("<8s4B2H2BH",
  91. b"HXCPICFE",
  92. 0,
  93. n_cyl,
  94. n_side,
  95. 0xff, # unknown encoding
  96. self.bitrate,
  97. 0, # rpm (unused)
  98. 0xff, # unknown interface
  99. 1, # rsvd
  100. 1) # track list offset
  101. # Pad the header and TLUT to 512-byte blocks.
  102. header += bytes([0xff] * (0x200 - len(header)))
  103. tlut += bytes([0xff] * (0x200 - len(tlut)))
  104. return header + tlut + tdat
  105. # Local variables:
  106. # python-indent: 4
  107. # End: