Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1078 1769 (N48W09 L=282)  1051    143        4 (1031/1723-1141/1828) -41.2 W04-W13 N43-N54
      2 1543 1656 (N37W46 L=319)  1612    344        3 (1360/1619-1686/1683)  -5.0 W29-W63 N34-N39
      3 1301 1527 (N28W21 L=294)  3398    752        2 ( 960/1374-1590/1609) -20.1 E02-W39 N17-N34
      4 1001 1396 (N20E00 L=273)   289     83        0 ( 969/1358-1036/1430) -59.0 E02-W01 N17-N22
      5 1794 1326 (N14W58 L=331)  1540    204        3 (1731/1252-1859/1407) -48.0 W52-W65 N09-N19
      6 1336 1360 (N17W21 L=294)   230     90        0 (1326/1319-1352/1400)  78.7 W20-W22 N14-N19
      7  603 1352 (N19E26 L=247)   347     70        2 ( 569/1336- 646/1375) -21.8 E28-E23 N18-N20
      8 1718 1294 (N12W49 L=322)   631    120        1 (1659/1254-1757/1329) -29.7 W44-W52 N09-N14
      9 1577 1242 (N09W36 L=310)   268     62        2 (1563/1206-1595/1272)  63.4 W35-W38 N06-N11
     10  755 1167 (N07E16 L=257)   460    119        0 ( 706/1124- 801/1215) -45.0 E18-E13 N04-N10
     11 1765  890 (S12W51 L=324)  2298    329        6 (1691/ 746-1839/1052)  68.7 W45-W58 S20-S02
     12 1143  871 (S12W07 L=280)   216     60        1 (1122/ 840-1166/ 897)  56.3 W05-W08 S13-S10
     13 1418  769 (S19W25 L=298)   710    245        0 (1386/ 662-1481/ 869)  67.4 W22-W28 S25-S13
     14 1878  786 (S18W67 L=341)   660     92        2 (1853/ 755-1903/ 830) -32.0 W63-W72 S20-S15
     15  301  690 (S19E55 L=218)   476    102        1 ( 277/ 638- 324/ 730)  71.6 E57-E52 S22-S16
     16 1676  542 (S34W52 L=325)  1775    265        2 (1559/ 481-1762/ 627)  30.5 W41-W58 S38-S28
     17 1103  589 (S30W04 L=277)   229     67        0 (1077/ 559-1136/ 622)  53.1 W01-W06 S31-S27
     18 1310  461 (S39W20 L=293)   423     90        1 (1282/ 435-1355/ 500) -39.8 W17-W24 S41-S36
     19  968  421 (S41E07 L=266)   402     55        1 ( 926/ 404-1002/ 441)   0.0 E10-E04 S42-S39
  

    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