Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  775 1581 (N22E27 L=111)  2424    588        9 ( 657/1427- 960/1796) -30.3 E47-E11 N15-N36
      2  548 1478 (N12E39 L= 99)  1473    490        4 ( 493/1268- 590/1687) -69.7 E47-E34 S01-N26
      3 1058 1663 (N33E12 L=126)   287     55        1 (1028/1639-1083/1684)  45.0 E13-E09 N32-N35
      4  498 1608 (N20E49 L= 89)   455     96        0 ( 465/1573- 540/1636)  68.2 E51-E46 N17-N22
      5 1318 1474 (N26W09 L=148)   416     93        1 (1301/1436-1341/1523) -63.4 W07-W11 N23-N29
      6 1484 1433 (N27W22 L=160)   312    107        0 (1457/1380-1511/1482) -56.3 W19-W23 N24-N30
      7 1870 1053 (N14W59 L=197)  2008    283        5 (1828/ 907-1914/1177) -51.8 W53-W67 N06-N21
      8 1061 1066 (S04W01 L=139)   317     91        1 (1028/1040-1103/1102) -21.8 E01-W04 S05-S02
      9  688 1021 (S13E20 L=118)   368    140        0 ( 641/ 968- 733/1048)  45.0 E22-E17 S16-S10
     10  764  987 (S14E15 L=123)   567    111        2 ( 709/ 950- 807/1021)   0.0 E17-E11 S15-S11
     11  475  844 (S27E33 L=105)   677    163        0 ( 404/ 786- 526/ 892) -24.4 E39-E28 S29-S24
     12 1164  818 (S16W12 L=150)   244     57        0 (1126/ 803-1192/ 835)  45.0 W10-W13 S17-S14
     13  421  720 (S35E38 L=100)   395     87        0 ( 389/ 687- 466/ 754) -23.2 E40-E33 S36-S33
     14 1201  571 (S30W22 L=160)  1327    398        0 (1018/ 507-1353/ 605)   4.4 W08-W34 S32-S27
     15  623  528 (S45E20 L=118)  1282    242        3 ( 521/ 497- 742/ 574)  20.2 E29-E10 S47-S39
     16 1713  505 (S20W65 L=203)   469    118        0 (1660/ 476-1748/ 547)  60.3 W61-W67 S23-S16
     17 1290  435 (S36W34 L=172)   665    111        0 (1247/ 399-1343/ 468)  56.3 W30-W36 S38-S32
     18  669  259 (S68E16 L=122)  1023    156        0 ( 593/ 229- 742/ 287)  -2.3 E29-E04 S69-S65
  

    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