Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1110 1890 (N57E35 L=  2)   759     60        1 (1091/1863-1132/1917) -15.9 E43-E29 N53-N60
      2 1475 1718 (N56W12 L= 49)   675    104        1 (1447/1668-1516/1744)  67.4 W08-W18 N50-N62
      3 1679 1580 (N53W39 L= 76)   728     73        1 (1657/1541-1700/1616)  74.1 W35-W43 N49-N56
      4 1026 1437 (N19E11 L= 26)   284     97        0 ( 977/1414-1072/1456)   9.5 E14-E08 N17-N19
      5  351 1407 (N01E55 L=342)   347     44        1 ( 325/1388- 375/1423)  11.3 E58-E53 N00-N02
      6 1614 1275 (N27W31 L= 68)   439    161        0 (1568/1230-1635/1335)  82.9 W27-W32 N22-N30
      7 1620 1107 (N17W34 L= 71)   645    191        1 (1594/1022-1667/1170)  90.0 W32-W37 N11-N21
      8 1733 1013 (N15W45 L= 82)  1133    276        1 (1653/ 951-1789/1076)  39.3 W39-W50 N10-N19
      9 1889  973 (N18W65 L=102)  2133    156        3 (1848/ 924-1932/1025)  45.0 W58-W74 N14-N22
     10  233  945 (S29E56 L=341)   868    168        1 ( 195/ 852- 257/1010) -66.0 E60-E52 S33-S24
     11 1825  905 (N12W58 L= 95)  1030    120        2 (1787/ 863-1853/ 973)  76.0 W54-W61 N08-N16
     12  452  818 (S31E32 L=  5)  1136    303        3 ( 344/ 759- 578/ 916)  39.8 E42-E23 S37-S22
     13  345  862 (S31E43 L=354)   569    161        0 ( 280/ 822- 420/ 893)  26.6 E49-E35 S35-S27
     14  180  837 (S36E65 L=332)   932     77        5 ( 156/ 799- 210/ 881) -78.7 E69-E59 S39-S33
     15 1184  417 (S35W31 L= 68)   563     59        1 (1151/ 397-1225/ 445)  36.9 W28-W34 S36-S32
     16  974  422 (S41W18 L= 55)   345    103        1 ( 920/ 406-1026/ 443)  20.6 W13-W21 S43-S39
     17 1467  280 (S33W66 L=103)   488     31        2 (1444/ 264-1488/ 299)   0.0 W61-W69 S34-S31
  

    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