Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  875 1746 (N53E17 L=157)   603    110        2 ( 821/1719- 920/1763)  11.3 E22-E12 N50-N54
      2  366 1538 (N33E59 L=115)  1910    221        3 ( 300/1480- 445/1623) -12.4 E89-E47 N29-N38
      3 1425 1526 (N36W30 L=204)  1063    232        1 (1332/1458-1501/1579) -30.3 W22-W35 N30-N39
      4  641 1458 (N29E29 L=145)  1329    322        2 ( 500/1425- 785/1487)  -2.5 E40-E17 N27-N31
      5 1492 1342 (N23W32 L=206)   257     52        0 (1479/1312-1505/1377) -90.0 W30-W32 N21-N25
      6 1783 1193 (N13W55 L=229)  1053    248        1 (1697/1104-1849/1271) -40.2 W47-W62 N07-N18
      7  480 1137 (N08E36 L=138)   437    185        0 ( 434/1055- 516/1204) -48.4 E40-E32 N03-N12
      8 1856  952 (S01W62 L=236)   431     59        1 (1841/ 917-1872/ 978)  45.0 W60-W64 S03-N00
      9 1491  863 (S06W30 L=204)   313     97        1 (1466/ 822-1511/ 915) -68.2 W28-W31 S08-S03
     10  439  773 (S14E39 L=135)   331    131        0 ( 427/ 716- 467/ 832) -81.9 E40-E37 S17-S10
     11 1783  760 (S13W57 L=231)  1266    148        2 (1743/ 708-1829/ 819) -60.3 W52-W63 S16-S09
     12  523  699 (S18E33 L=141)   376     84        2 ( 496/ 650- 550/ 740)  71.6 E35-E31 S21-S15
     13  417  667 (S21E43 L=131)   375    150        0 ( 361/ 643- 477/ 710)  15.3 E48-E36 S22-S17
     14 1110  556 (S26W07 L=181)  1085    247        3 (1041/ 511-1191/ 650) -42.0 W01-W12 S29-S20
     15  898  563 (S26E08 L=166)   418    129        0 ( 838/ 539- 949/ 598)  15.9 E11-E04 S28-S23
     16 1384  215 (S56W47 L=221)   584     59        0 (1362/ 190-1398/ 246) -54.5 W41-W53 S59-S52
  

    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