Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1122 1317 (N18W00 L=281)  3097    765        6 ( 855/1021-1353/1595)  16.0 E23-W19 N04-N31
      2 1713 1564 (N50W49 L=331)   551    112        1 (1675/1550-1740/1581)  14.0 W44-W53 N48-N51
      3 1450 1554 (N41W20 L=302)   275     65        0 (1431/1533-1483/1570)   0.0 W18-W23 N39-N42
      4 1870 1091 (N19W61 L=342)  1829    324        4 (1799/1005-1928/1165) -21.3 W52-W70 N15-N22
      5  151 1007 (S18E65 L=216)   662    162        1 ( 130/ 942- 174/1061)  63.4 E67-E63 S22-S14
      6  374 1020 (S13E43 L=239)   296     49        1 ( 359/ 987- 393/1047)  76.0 E43-E40 S15-S11
      7 1387  933 (N00W23 L=305)   391    118        1 (1339/ 891-1439/ 984) -23.2 W19-W27 S01-N02
      8  419  826 (S24E37 L=244)   564    173        1 ( 348/ 793- 488/ 858)   0.0 E42-E31 S25-S22
      9  963  788 (S16W01 L=282)  1745    474        4 ( 739/ 742-1194/ 850)  24.8 E13-W13 S21-S08
     10 1368  367 (S35W42 L=323)  5926    957       14 ( 979/ 322-1637/ 459)  21.6 W12-W61 S44-S22
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y