Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1258 1709 (N45W37 L=352)  1218    199        1 (1202/1646-1344/1747)   6.3 W28-W46 N41-N46
      2  621 1703 (N60E26 L=289)   826     96        1 ( 579/1661- 654/1741)  29.7 E32-E21 N55-N63
      3 1466 1574 (N31W45 L=359)   464     89        2 (1427/1554-1512/1591) -38.7 W41-W48 N28-N32
      4 1728 1513 (N20W66 L= 21)   947    144        0 (1667/1475-1782/1567) -41.2 W62-W72 N17-N24
      5  999 1443 (N32W06 L=321)  1964    519        1 ( 780/1367-1230/1527) -33.7 E08-W19 N23-N42
      6  251 1402 (N40E64 L=251)   861     92        2 ( 214/1366- 288/1437) -12.5 E74-E56 N37-N42
      7 1594 1332 (N14W43 L=358)   433    133        0 (1543/1278-1623/1369) -63.4 W40-W45 N10-N16
      8  177 1093 (N21E65 L=250)   886    110        2 ( 161/1039- 196/1147)  66.8 E67-E61 N17-N24
      9  270 1072 (N19E53 L=261)   491    148        0 ( 255/1015- 296/1131)  81.9 E54-E50 N15-N22
     10  601 1015 (N12E26 L=288)  1152    310        1 ( 482/ 909- 707/1094) -42.9 E34-E20 N04-N18
     11 1575  901 (S11W32 L=347)   535    134        0 (1558/ 830-1596/ 964)  60.3 W30-W34 S14-S07
     12 1696  938 (S11W42 L=357)   292     89        0 (1657/ 924-1741/ 955) -15.9 W39-W46 S12-S09
     13 1748  865 (S17W46 L=  1)   310     68        0 (1727/ 850-1782/ 894) -31.0 W44-W49 S18-S14
     14  916  683 (S12E12 L=303)   768    325        1 ( 852/ 571- 975/ 795) -90.0 E14-E09 S19-S04
     15 1028  695 (S13E05 L=310)   294    112        0 (1012/ 639-1047/ 744)  81.9 E05-E04 S16-S09
     16 1225  470 (S32W03 L=318)  1435    306        3 (1112/ 404-1338/ 589)  27.9 E05-W12 S34-S25
     17  741  224 (S43E45 L=270)  1162    171        2 ( 674/ 189- 817/ 262)   3.4 E55-E34 S46-S40
  

    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