Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  847 1817 (N62W23 L=356)   344     76        0 ( 826/1779- 863/1842)  29.1 W18-W27 N58-N64
      2 1340 1774 (N39W54 L= 26)   662     84        0 (1291/1753-1382/1789) -20.6 W49-W58 N36-N40
      3  647 1713 (N61E06 L=326)   357     57        2 ( 618/1694- 676/1733) -15.9 E09-E02 N59-N63
      4 1313 1629 (N31W39 L= 12)   495    135        0 (1253/1598-1371/1664)   0.0 W33-W45 N30-N32
      5 1372 1546 (N25W38 L= 10)   632    189        0 (1284/1527-1464/1564) -22.6 W32-W44 N22-N27
      6  864 1534 (N40W06 L=338)   469    104        2 ( 824/1504- 913/1554)  -8.1 W02-W09 N37-N41
      7 1404 1433 (N18W34 L=  7)   226     94        0 (1383/1387-1426/1476)  71.6 W33-W35 N14-N20
      8  423 1383 (N43E34 L=299)   632    110        1 ( 373/1348- 479/1415) -28.6 E39-E28 N39-N45
      9  988 1329 (N24W06 L=339)   209     63        0 ( 969/1294-1011/1354)  26.6 W04-W08 N22-N25
     10 1074 1321 (N21W11 L=344)   358    111        0 (1031/1268-1108/1354) -71.6 W09-W12 N17-N24
     11 1815 1261 (S04W59 L= 31)   779     94        1 (1787/1210-1848/1303)  11.3 W53-W64 S06-S02
     12  143 1037 (N26E68 L=265)   840    131        1 ( 108/ 984- 186/1096) -29.7 E75-E61 N22-N30
     13 1601  949 (S14W32 L=  4)   304     48        1 (1576/ 927-1625/ 978)  14.0 W29-W33 S15-S12
     14 1692  888 (S20W37 L= 10)   320     87        0 (1652/ 871-1734/ 911) -29.7 W34-W41 S21-S17
     15 1346  740 (S18W11 L=343)   496    131        0 (1319/ 674-1368/ 795)  60.3 W07-W12 S21-S14
     16  773  647 (S08E24 L=308)   668    209        0 ( 747/ 554- 813/ 744)  48.0 E28-E19 S13-S03
     17  770  438 (S21E33 L=300)   248     80        0 ( 757/ 397- 791/ 478)  45.0 E34-E29 S23-S18
     18  852  270 (S35E37 L=295)  2533    382        4 ( 699/ 190-1016/ 360)   6.3 E57-E21 S38-S30
  

    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