Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1345 1669 (N32W44 L=197)   760    268        1 (1234/1639-1449/1701) -13.2 W35-W52 N29-N35
      2 1645 1634 (N20W64 L=218)  2437    453        6 (1535/1594-1732/1670) -24.2 W56-W76 N15-N25
      3 1233 1529 (N28W29 L=183)   574    180        0 (1193/1452-1282/1604)  79.7 W27-W31 N22-N33
      4 1663 1485 (N12W56 L=210)   401     72        0 (1634/1449-1688/1516) -78.7 W55-W57 N09-N14
      5 1625 1432 (N11W50 L=204)   426    132        0 (1611/1369-1645/1492)  60.3 W48-W53 N06-N13
      6  706 1287 (N30E12 L=142)   536    139        1 ( 638/1258- 772/1320) -29.1 E16-E07 N26-N33
      7  804 1202 (N22E08 L=146)   320     90        0 ( 773/1170- 834/1245) -80.5 E08-E06 N19-N25
      8 1901 1130 (S15W64 L=218)  1592    252        3 (1858/1051-1943/1191) -70.0 W59-W67 S21-S10
      9  220  861 (N16E59 L= 94)  1566    336        3 ( 180/ 772- 249/ 974)  45.0 E67-E53 N10-N22
     10 1806  865 (S26W47 L=201)  1675    247        2 (1774/ 741-1836/ 964)  72.9 W43-W51 S35-S20
     11 1244  723 (S16W04 L=158)   912    268        2 (1144/ 594-1316/ 848)  36.3 E04-W11 S21-S10
     12  144  820 (N14E71 L= 83)   724    187        2 ( 111/ 801- 181/ 842) -11.3 E76-E66 N12-N15
     13  740  682 (S05E25 L=128)   957    205        3 ( 697/ 601- 795/ 769) -90.0 E27-E23 S11-N00
     14 1514  683 (S27W20 L=174)   281    145        0 (1494/ 609-1533/ 746)  58.0 W17-W22 S31-S23
     15 1398  557 (S31W10 L=163)   519    130        3 (1347/ 525-1458/ 595)   0.0 W05-W14 S33-S30
     16 1155  406 (S34E11 L=143)  1400    257        6 (1076/ 320-1247/ 495)  22.4 E19-E02 S38-S30
     17  557  331 (S23E54 L=100)   562    170        2 ( 539/ 292- 575/ 371)  35.5 E57-E50 S25-S20
  

    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