Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1540 1754 (N32W68 L=230)   648     63        1 (1508/1732-1563/1778)   0.0 W63-W75 N30-N33
      2 1257 1722 (N39W39 L=201)  1002    212        0 (1162/1674-1347/1759)  -3.0 W29-W48 N37-N40
      3  481 1335 (N34E30 L=132)   736    179        2 ( 396/1315- 570/1351) -11.3 E37-E22 N31-N35
      4 1090 1267 (N15W09 L=171)   199     79        0 (1047/1254-1128/1276) -36.9 W07-W11 N13-N16
      5  679 1227 (N22E16 L=145)   303     77        0 ( 639/1207- 722/1244)  -9.5 E19-E13 N21-N23
      6 1798 1217 (S04W56 L=218)   350     77        0 (1762/1183-1821/1244) -45.0 W53-W58 S06-S02
      7 1866 1076 (S14W60 L=222)  1090    187        4 (1814/ 995-1923/1149)  16.4 W53-W71 S18-S11
      8  270 1005 (N18E51 L=110)  1545    342        1 ( 135/ 947- 387/1106)   3.8 E69-E39 N14-N21
      9 1745 1035 (S13W47 L=208)   960    203        3 (1676/ 971-1796/1098) -45.0 W41-W52 S17-S09
     10 1923  913 (S26W68 L=230)   663     59        1 (1903/ 884-1951/ 933)   5.2 W64-W75 S27-S23
     11 1554  795 (S22W28 L=189)   349     81        1 (1534/ 758-1571/ 844)  68.2 W25-W29 S24-S19
     12 1383  729 (S22W15 L=176)   662    193        0 (1319/ 647-1430/ 809)  35.0 W09-W19 S25-S18
     13 1177  662 (S21W01 L=162)   770    258        1 (1146/ 549-1239/ 775)  53.1 E03-W06 S27-S15
     14  874  661 (S14E17 L=144)   311    113        0 ( 838/ 617- 921/ 711)  29.7 E20-E13 S16-S12
     15 1408  593 (S31W14 L=176)   437    103        0 (1356/ 566-1450/ 634)   0.0 W10-W17 S32-S29
     16  818  443 (S26E28 L=133)   494    156        1 ( 783/ 373- 843/ 515)  36.0 E34-E23 S30-S22
     17 1630  448 (S50W35 L=197)   848    128        0 (1562/ 427-1687/ 476)  -6.7 W26-W43 S52-S48
  

    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