Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1241 1785 (N51W44 L=173)  1197    169        2 (1172/1741-1324/1815)   0.0 W33-W55 N48-N52
      2 1477 1653 (N35W52 L=181)   624     85        2 (1442/1629-1522/1678) -45.0 W48-W55 N33-N37
      3 1720 1460 (N18W62 L=191)   416     62        0 (1689/1438-1752/1477) -33.7 W58-W64 N15-N19
      4 1379 1445 (N25W32 L=161)   230     52        1 (1347/1432-1405/1459)   0.0 W29-W34 N23-N25
      5 1778 1362 (N11W62 L=192)  1408    222        4 (1696/1307-1831/1434)   5.4 W51-W72 N08-N14
      6 1549 1257 (N10W38 L=167)   765    197        1 (1460/1221-1619/1296) -15.3 W32-W44 N08-N14
      7  223 1161 (N25E60 L= 69)   685    111        4 ( 201/1105- 245/1220)  74.1 E62-E57 N21-N28
      8  484 1092 (N18E34 L= 95)   255    112        0 ( 467/1038- 513/1142) -74.1 E35-E32 N14-N21
      9 1894 1049 (S10W66 L=195)  2252    178        5 (1847/ 986-1944/1108)  50.2 W59-W77 S14-S06
     10 1703  907 (S14W43 L=172)  1303    405        5 (1601/ 755-1766/1065)  45.0 W33-W50 S21-S04
     11 1640  935 (S10W38 L=167)   986    169        6 (1584/ 885-1691/ 985)  12.5 W33-W42 S13-S07
     12 1000  719 (S11E07 L=123)   538    138        1 ( 976/ 659-1035/ 785)  54.5 E08-E03 S14-S07
     13  602  663 (S09E33 L= 96)  1339    335        6 ( 466/ 585- 732/ 748) -40.6 E40-E25 S15-S02
     14  917  704 (S11E12 L=117)   235     54        2 ( 899/ 671- 941/ 731)  71.6 E12-E10 S12-S09
     15  784  487 (S23E26 L=104)   418     82        2 ( 767/ 445- 806/ 527)  71.6 E27-E23 S26-S20
     16  934  400 (S32E18 L=111)   920    184        2 ( 860/ 367-1026/ 436)   4.4 E24-E11 S33-S29
     17  518  387 (S27E52 L= 77)  1127    103        2 ( 482/ 346- 560/ 430)  81.9 E55-E47 S30-S23
     18 1232  262 (S51E00 L=129)   421     49        1 (1208/ 239-1254/ 284)  26.6 E03-W01 S52-S48
  

    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