Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1250 1812 (N52E15 L=219)   451     46        0 (1228/1785-1274/1835) -81.9 E16-E12 N48-N55
      2  714 1735 (N27E45 L=189)   527    100        0 ( 687/1690- 738/1780) -54.5 E47-E42 N23-N30
      3 1422 1724 (N51W04 L=238)   693     83        0 (1375/1694-1477/1753)   5.7 W00-W10 N49-N53
      4  408 1648 (N15E64 L=170)   981    152        2 ( 384/1575- 434/1717) -45.0 E69-E58 N09-N19
      5  370 1394 (S01E52 L=182)   591    144        1 ( 305/1357- 413/1422)   5.2 E57-E46 S03-N00
      6  236 1131 (S19E57 L=177)  1120    134        2 ( 180/1092- 306/1162)  19.7 E63-E49 S21-S16
      7  826 1081 (S08E12 L=222)   279    116        0 ( 807/1030- 857/1131)  90.0 E13-E11 S11-S04
      8 1521 1085 (N11W27 L=260)   224     53        0 (1497/1055-1545/1108) -14.0 W24-W28 N10-N12
      9 1217  967 (S04W12 L=246)   422    142        0 (1189/ 911-1260/1032) -39.8 W09-W15 S06-N00
     10  598  854 (S27E21 L=213)   632    228        0 ( 566/ 752- 658/ 925)  90.0 E22-E18 S32-S21
     11 1212  731 (S17W19 L=253)   205     58        0 (1185/ 715-1240/ 756)   0.0 W16-W21 S17-S15
     12 1423  687 (S13W33 L=267)   333    121        0 (1366/ 670-1488/ 702)  29.7 W29-W36 S14-S10
     13 1528  583 (S15W44 L=278)   262     47        0 (1513/ 556-1546/ 613) -63.4 W43-W46 S16-S12
     14  538  528 (S50E20 L=214)   265     49        0 ( 519/ 497- 559/ 557) -53.1 E22-E17 S51-S47
     15  684  475 (S49E05 L=228)   752    150        3 ( 605/ 452- 760/ 503)  17.1 E12-W01 S50-S45
     16 1354  440 (S28W41 L=275)   471     98        1 (1311/ 424-1405/ 466)  38.7 W37-W44 S29-S25
     17 1542  353 (S26W63 L=297)   562     52        2 (1513/ 333-1574/ 375)  36.9 W59-W65 S27-S24
     18 1166  224 (S48W49 L=282)  1038    128        3 (1110/ 176-1223/ 268)  -3.0 W40-W60 S49-S45
  

    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