Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1114 1781 (N43E19 L=172)   886    196        1 (1067/1716-1163/1847) -29.7 E26-E12 N38-N48
      2  519 1688 (N21E54 L=137)   414     24        0 ( 497/1673- 545/1705)  33.7 E56-E51 N19-N21
      3  555 1595 (N15E46 L=145)   603    178        0 ( 536/1534- 575/1661) -41.2 E50-E41 N11-N19
      4 1252 1102 (N03W11 L=202)   244     51        1 (1222/1080-1292/1120)  14.0 W08-W12 N01-N04
      5 1777 1030 (N15W47 L=238)   366     83        0 (1764/ 984-1789/1076) -68.2 W45-W48 N12-N17
      6  273 1051 (S22E50 L=141)   476    155        1 ( 224/1007- 322/1072)  32.0 E54-E46 S25-S20
      7 1410 1028 (N03W22 L=213)   276     59        0 (1386/ 992-1430/1062)  90.0 W20-W22 N00-N05
      8 1438  851 (S06W28 L=219)   810    132        0 (1413/ 768-1464/ 923) -63.4 W25-W31 S09-S01
      9 1175  783 (S16W15 L=206)  1079    224        1 (1050/ 741-1276/ 839)  42.3 W08-W19 S21-S11
     10 1854  813 (N06W62 L=254)   470    166        0 (1818/ 798-1894/ 830)  16.7 W58-W68 N04-N07
     11  875  523 (S40W06 L=197)  1711    334        3 ( 712/ 490-1056/ 567)  15.6 E07-W18 S44-S35
     12 1480  462 (S25W47 L=238)  1706    398        0 (1314/ 433-1626/ 505)  24.4 W36-W58 S29-S19
  

    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