Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1710 1611 (N18W69 L=104)  2541    224        1 (1611/1545-1792/1682) -50.2 W62-W79 N12-N24
      2 1488 1618 (N25W48 L= 83)  3332    557        2 (1258/1507-1723/1678) -28.7 W31-W62 N12-N31
      3 1381 1502 (N21W34 L= 70)   259     90        0 (1338/1489-1432/1512) -26.6 W31-W37 N19-N22
      4 1740 1129 (S08W47 L= 82)  1403    311        3 (1682/1044-1786/1238)  85.2 W41-W51 S13-S01
      5 1652 1156 (S04W41 L= 76)   266     91        1 (1613/1134-1691/1176)   0.0 W37-W44 S05-S03
      6  473 1103 (N20E32 L=  3)   812    277        2 ( 376/1020- 586/1169) -42.7 E38-E25 N13-N25
      7  900  863 (S03E11 L= 25)  1948    512        4 ( 726/ 684-1043/1000) -57.7 E19-E06 S15-N05
      8 1356  930 (S09W16 L= 52)   529     95        2 (1321/ 891-1383/ 982)  45.0 W13-W18 S11-S06
      9 1214  893 (S08W08 L= 43)   407    134        2 (1182/ 830-1234/ 958)  60.3 W04-W09 S11-S04
     10 1265  925 (S08W11 L= 46)   238     36        1 (1245/ 900-1289/ 952)  63.4 W10-W12 S09-S06
     11  372  903 (N10E43 L=352)   616    146        0 ( 314/ 872- 437/ 944) -29.1 E47-E38 N07-N13
     12 1461  671 (S27W18 L= 53)  1184    241        4 (1356/ 633-1571/ 721) -11.3 W11-W26 S29-S23
     13  584  674 (S07E34 L=  1)   344     25        0 ( 558/ 649- 608/ 696) -56.3 E35-E32 S08-S05
     14 1146  402 (S37E09 L= 26)  3949    674        8 ( 996/ 185-1261/ 691)  45.9 E32-W02 S50-S18
     15  744  592 (S15E27 L=  8)   522    116        1 ( 714/ 522- 770/ 637)  76.0 E29-E24 S19-S11
     16  557  598 (S10E39 L=356)   632    101        4 ( 507/ 560- 600/ 631) -51.3 E41-E36 S13-S08
     17  536  424 (S20E49 L=346)  2876    463        4 ( 388/ 314- 647/ 569)  14.3 E77-E34 S25-S13
     18 1063  338 (S39E17 L= 18)   251     83        0 (1054/ 294-1076/ 379)  56.3 E19-E15 S42-S36
  

    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