Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  642 1759 (N33E48 L=286)  1213    212        2 ( 580/1683- 682/1808) -32.7 E54-E39 N27-N36
      2  568 1678 (N25E48 L=286)  1103    158        3 ( 523/1626- 623/1725)  63.4 E51-E43 N21-N28
      3  592 1454 (N11E36 L=298)  1101    261        2 ( 532/1368- 670/1569)  78.7 E38-E33 N04-N19
      4  756 1451 (N13E25 L=309)   357    115        1 ( 734/1402- 785/1494) -80.5 E26-E23 N10-N16
      5  955 1329 (N09E10 L=324)   397    143        1 ( 912/1271-1004/1397)  83.7 E11-E07 N05-N14
      6 1850 1223 (N24W56 L= 29)   593     69        0 (1831/1185-1872/1262) -53.1 W53-W58 N22-N26
      7  683 1218 (S02E24 L=310)   197     77        0 ( 655/1187- 710/1246) -31.0 E26-E21 S03-N00
      8 1567 1199 (N15W30 L=  3)   417    145        0 (1521/1130-1590/1240) -56.3 W26-W31 N10-N16
      9 1397 1189 (N10W19 L=352)   199     69        0 (1385/1154-1409/1227) -63.4 W17-W19 N08-N12
     10  207 1036 (S19E59 L=274)   947    145        1 ( 183/ 970- 237/1103) -63.4 E62-E55 S22-S14
     11 1790 1002 (N09W51 L= 25)   495    104        1 (1775/ 947-1810/1053)  90.0 W49-W52 N05-N12
     12 1656 1021 (N07W39 L= 13)   244     71        1 (1619/1001-1688/1043)   0.0 W36-W41 N05-N07
     13  686  920 (S19E19 L=315)   400    148        1 ( 658/ 853- 713/ 987)  77.5 E19-E17 S23-S14
     14 1029  895 (S14W03 L=336)   202     97        0 ( 980/ 874-1068/ 911)  36.9 W00-W04 S15-S12
     15  514  708 (S35E29 L=304)   710    265        1 ( 407/ 679- 621/ 772)  26.6 E39-E21 S38-S29
     16 1798  669 (S09W64 L= 37)  1273    155        3 (1750/ 611-1836/ 736) -74.1 W57-W70 S11-S04
     17 1719  602 (S14W59 L= 33)   511    119        2 (1698/ 544-1736/ 661) -50.2 W55-W62 S17-S11
     18 1599  489 (S24W55 L= 29)  2023    344        5 (1538/ 323-1661/ 629) -49.9 W50-W66 S34-S15
  

    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