Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  599 1754 (N58E45 L=277)   477     49        1 ( 568/1734- 626/1773)  45.0 E49-E40 N56-N59
      2  732 1516 (N38E19 L=303)   533    141        0 ( 664/1481- 802/1547)  15.3 E24-E13 N35-N39
      3  522 1475 (N36E37 L=285)   351     76        0 ( 482/1446- 559/1509) -24.0 E41-E32 N33-N38
      4 1826 1412 (N21W71 L= 34)  1980    142        6 (1798/1345-1862/1470) -90.0 W66-W79 N16-N24
      5 1734 1353 (N19W56 L= 18)   583     79        1 (1705/1313-1766/1388) -68.2 W53-W58 N16-N21
      6  168 1242 (N20E70 L=252)  2797    209        2 ( 135/1142- 203/1324)  70.0 E77-E66 N14-N25
      7  716 1114 (N12E19 L=303)   612    261        0 ( 636/1023- 803/1158)  31.0 E24-E13 N06-N14
      8  235 1068 (N10E57 L=265)   652    121        0 ( 207/1022- 283/1103)  51.3 E60-E52 N07-N12
      9  478  878 (S02E36 L=286)  2080    491        4 ( 408/ 710- 651/1067) -55.0 E40-E25 S12-N09
     10 1358  922 (S04W20 L=342)   393    152        0 (1310/ 873-1404/ 988)  45.0 W16-W23 S06-N00
     11  216  640 (S17E70 L=252)  4314    302        7 ( 155/ 527- 288/ 749) -70.3 E78-E63 S25-S10
     12  699  575 (S22E25 L=297)   774    280        0 ( 594/ 540- 822/ 642) -26.6 E31-E16 S24-S17
     13  931  421 (S34E11 L=311)   335     61        0 ( 906/ 390- 955/ 450)  45.0 E13-E09 S36-S31
     14  915  358 (S39E14 L=308)   438     50        1 ( 882/ 333- 950/ 374)   0.0 E16-E10 S41-S37
     15  538  278 (S45E57 L=265)   963     98        1 ( 499/ 256- 571/ 312)   8.5 E70-E50 S48-S41
     16  850  272 (S47E22 L=300)   977    141        2 ( 791/ 239- 915/ 310)  23.2 E29-E15 S49-S43
  

    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