Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  684 1768 (N31E47 L=214)  1608    169        4 ( 615/1710- 748/1823)  78.7 E51-E44 N25-N36
      2 1508 1793 (N64W20 L=281)   614     50        0 (1491/1767-1524/1819)  47.5 W16-W27 N58-N70
      3  806 1688 (N27E33 L=228)  1125    236        2 ( 747/1645- 884/1757)  26.6 E38-E28 N24-N32
      4  721 1640 (N22E36 L=225)   262     49        0 ( 693/1620- 748/1657)   0.0 E38-E33 N21-N23
      5  416 1572 (N12E56 L=205)   579    145        0 ( 381/1548- 470/1608)   0.0 E60-E50 N10-N14
      6 1544 1446 (N32W23 L=284)  1667    417        0 (1417/1303-1692/1603) -24.8 W10-W36 N26-N39
      7  497 1539 (N11E47 L=214)   619    132        2 ( 459/1503- 526/1595) -74.1 E50-E44 N08-N15
      8 1646 1467 (N37W32 L=293)   371     85        0 (1607/1434-1681/1501) -23.2 W27-W35 N35-N38
      9 1806 1327 (N33W49 L=310)   552     59        0 (1779/1296-1831/1352)  -9.5 W45-W52 N31-N34
     10 1216 1156 (N05W07 L=269)   363     85        0 (1171/1128-1255/1186)  45.0 W05-W09 N02-N07
     11  524 1128 (S12E33 L=229)   594    131        1 ( 474/1078- 574/1168)  45.0 E35-E29 S15-S08
     12  160 1007 (S25E66 L=195)  2311    173        6 ( 123/ 931- 195/1093)  84.3 E72-E60 S29-S19
     13  370  941 (S26E42 L=220)   423    120        1 ( 334/ 910- 424/ 980)  40.6 E44-E37 S28-S22
     14 1194  876 (S11W13 L=274)   430    157        0 (1102/ 854-1248/ 897)  20.6 W08-W16 S14-S09
     15  553  779 (S32E25 L=237)  2359    513        4 ( 365/ 653- 776/ 862)  29.9 E42-E09 S44-S24
     16 1601  507 (S20W55 L=316)   777    175        0 (1532/ 461-1681/ 561)  56.3 W51-W59 S23-S14
     17 1309  522 (S28W33 L=294)   532    180        0 (1226/ 504-1388/ 546)  14.0 W26-W38 S29-S25
     18 1486  356 (S31W57 L=319)  2370    283        4 (1349/ 328-1609/ 391)  26.6 W45-W68 S36-S25
  

    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