Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1457 1647 (N29W51 L=250)   398     75        0 (1421/1621-1496/1672) -68.2 W49-W52 N26-N31
      2  676 1577 (N50E09 L=191)   247     42        1 ( 652/1559- 701/1602) -11.3 E11-E06 N48-N51
      3 1562 1533 (N19W51 L=250)   393     87        2 (1519/1520-1608/1550) -23.2 W47-W55 N17-N20
      4 1530 1388 (N13W41 L=240)   394     64        1 (1504/1358-1558/1412)  21.8 W38-W43 N11-N14
      5 1585 1343 (N09W43 L=242)   358    102        1 (1553/1301-1616/1383)  31.0 W39-W45 N07-N11
      6  376 1127 (N28E41 L=159)   272     43        0 ( 346/1111- 402/1146)   0.0 E43-E38 N26-N28
      7 1031 1079 (N10W02 L=201)   346    115        1 ( 986/1031-1075/1109) -63.4 W00-W03 N06-N12
      8  757 1078 (N16E14 L=185)   284    138        0 ( 701/1044- 801/1104) -45.0 E17-E12 N13-N18
      9 1629  906 (S16W33 L=233)   291     77        1 (1603/ 877-1646/ 953)  63.4 W31-W35 S17-S13
     10 1236  708 (S16W04 L=203)  1010    316        4 (1145/ 614-1341/ 797)  13.0 E02-W11 S21-S12
     11 1872  710 (S38W60 L=260)   622     63        2 (1843/ 689-1905/ 733)  -4.4 W54-W69 S40-S36
     12  331  603 (S02E58 L=141)   884    148        0 ( 253/ 563- 382/ 634) -24.4 E65-E53 S05-N00
     13  379  454 (S13E62 L=137)  1739    200        7 ( 345/ 368- 409/ 545)  47.5 E70-E56 S18-S06
     14 1601  447 (S47W27 L=226)  1232    165        3 (1535/ 414-1673/ 495)   0.0 W19-W34 S49-S44
     15  745  331 (S28E40 L=160)  2511    393        2 ( 708/ 155- 790/ 494)  38.9 E59-E28 S42-S17
     16  475  423 (S16E55 L=145)   999    142        5 ( 439/ 377- 525/ 474) -74.1 E58-E51 S19-S12
     17  548  254 (S30E65 L=134)   645     51        0 ( 524/ 235- 574/ 281) -90.0 E69-E61 S31-S27
  

    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