Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1014 1569 (N36W17 L=103)   511    116        0 ( 951/1550-1072/1590) -20.6 W11-W20 N33-N38
      2 1277 1547 (N27W32 L=118)   251     83        1 (1237/1530-1319/1562)  -9.5 W28-W34 N25-N27
      3 1510 1459 (N15W42 L=128)   298     59        0 (1486/1432-1532/1478)  14.0 W40-W44 N13-N16
      4 1300 1301 (N13W23 L=109)   296    142        2 (1240/1258-1355/1324) -50.2 W20-W26 N09-N15
      5 1793 1014 (S17W49 L=135)   667    144        0 (1746/ 977-1835/1060) -45.0 W45-W52 S20-S13
      6  463 1035 (N19E34 L= 52)   259     77        1 ( 433/1004- 500/1058) -45.0 E35-E31 N16-N20
      7 1348  740 (S19W11 L= 97)   370    157        2 (1314/ 697-1383/ 777)  53.1 W08-W12 S22-S17
      8 1479  560 (S35W16 L=102)  1025    182        2 (1427/ 482-1518/ 658)  55.0 W10-W19 S39-S28
      9  915  464 (S24E23 L= 64)   771    254        2 ( 878/ 375- 963/ 555)  50.7 E26-E16 S31-S19
     10  625  354 (S24E47 L= 40)   708    173        1 ( 597/ 307- 655/ 414)  36.9 E51-E43 S26-S19
     11 1500  344 (S54W15 L=101)   446     91        0 (1453/ 323-1543/ 369)  -5.7 W09-W19 S56-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