Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1426 1624 (N31W48 L=348)   691    169        0 (1365/1575-1496/1673) -73.3 W45-W50 N25-N35
      2 1694 1568 (N19W68 L=  8)  1295    154        2 (1628/1517-1762/1607) -42.0 W62-W74 N13-N22
      3 1016 1511 (N36W13 L=313)   515    105        2 ( 965/1483-1077/1528) -20.6 W09-W17 N33-N37
      4 1464 1436 (N19W39 L=339)   562    132        1 (1406/1411-1521/1457) -20.6 W36-W44 N17-N21
      5 1274 1171 (N09W18 L=318)   312    109        0 (1234/1125-1310/1214) -80.5 W16-W19 N06-N12
      6  510 1110 (N23E31 L=269)   885    318        2 ( 376/1073- 663/1142) -20.0 E42-E20 N18-N27
      7 1713 1081 (S08W45 L=345)   613    112        1 (1687/1028-1746/1135)  39.8 W42-W49 S10-S05
      8 1318 1031 (N01W17 L=317)   202     68        1 (1299/ 994-1343/1061)  56.3 W15-W18 S01-N02
      9  350 1010 (N20E46 L=254)   717     91        1 ( 306/ 970- 385/1056)  33.7 E50-E42 N17-N22
     10  154 1002 (N22E70 L=230)   955    111        2 ( 122/ 956- 184/1047) -33.7 E76-E64 N18-N24
     11 1815  934 (S20W52 L=352)   961    174        1 (1749/ 899-1863/ 975) -17.1 W45-W58 S22-S16
     12 1255  914 (S04W11 L=310)   294     77        2 (1217/ 885-1286/ 935) -45.0 W08-W12 S06-S02
     13  389  885 (N12E44 L=256)   354    109        0 ( 332/ 864- 439/ 898) -12.5 E49-E40 N09-N13
     14 1319  630 (S23W08 L=308)  1210    343        1 (1232/ 510-1430/ 777)  33.7 W00-W18 S29-S16
     15 1526  548 (S36W22 L=322)   533    127        1 (1507/ 483-1541/ 612) -84.3 W20-W23 S40-S30
     16  729  471 (S19E33 L=267)   346     84        0 ( 704/ 435- 755/ 516)  80.5 E33-E31 S22-S16
     17 1127  292 (S43E15 L=285)  3138    392        4 ( 979/ 168-1222/ 426)  26.6 E36-E03 S51-S34
     18  711  399 (S24E38 L=262)   301    107        1 ( 652/ 380- 746/ 421)   0.0 E43-E34 S25-S22
     19  938  156 (S51E43 L=257)   549     76        0 ( 907/ 141- 971/ 177) -36.9 E46-E38 S52-S47
  

    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