Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1518 1670 (N53W37 L=308)   409     39        1 (1503/1647-1538/1692)  51.3 W34-W40 N50-N55
      2  380 1603 (N27E66 L=205)  1143     97        2 ( 355/1563- 414/1651) -33.7 E73-E60 N24-N30
      3  735 1556 (N30E30 L=242)  1072    226        2 ( 695/1455- 799/1647) -86.2 E31-E27 N22-N37
      4 1658 1535 (N43W49 L=320)   685     63        0 (1639/1498-1678/1573)  80.5 W45-W51 N40-N46
      5  359 1450 (N17E57 L=215)   406     80        1 ( 331/1414- 386/1487)  78.7 E58-E54 N14-N19
      6 1024 1416 (N25E06 L=265)   374     76        1 ( 984/1390-1056/1436) -21.8 E08-E03 N23-N25
      7  225 1052 (S09E59 L=213)  1081    124        3 ( 189/1001- 249/1107)  90.0 E63-E56 S12-S05
      8  335  938 (S14E46 L=225)   803    251        2 ( 291/ 852- 389/1038)  58.4 E50-E42 S20-S07
      9  316 1008 (S10E49 L=223)   409     89        1 ( 285/ 987- 341/1034) -11.3 E51-E46 S11-S08
     10  248  943 (S15E55 L=216)   488     69        0 ( 225/ 915- 276/ 973) -18.4 E58-E52 S16-S13
     11 1810  741 (S06W64 L=335)   757    144        1 (1788/ 701-1851/ 795) -22.6 W59-W71 S08-S03
     12  195  751 (S29E67 L=204)  1268    101        3 ( 152/ 717- 240/ 779)  12.3 E82-E58 S31-S26
     13 1005  640 (S24W04 L=276)   271     88        0 ( 985/ 619-1033/ 686) -53.1 W02-W06 S25-S20
     14 1787  615 (S14W68 L=339)  1398    104        1 (1753/ 569-1826/ 652)  35.5 W62-W73 S17-S11
     15 1682  598 (S17W56 L=327)   378     97        1 (1640/ 574-1720/ 624)  33.7 W52-W58 S18-S14
     16 1311  483 (S30W30 L=301)  1381    376        0 (1149/ 444-1480/ 526)  24.4 W18-W40 S34-S24
     17 1184  251 (S50W33 L=305)   504     74        1 (1146/ 237-1224/ 272)   7.1 W29-W37 S51-S48
     18 1369  229 (S49W57 L=328)   600     61        1 (1342/ 207-1400/ 250)  -4.4 W49-W62 S50-S47
  

    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