Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  662 1840 (N34E60 L=134)   698    102        1 ( 643/1797- 677/1894) -27.9 E72-E54 N30-N39
      2  727 1867 (N38E59 L=135)   645     52        1 ( 697/1846- 760/1884)   0.0 E64-E53 N36-N39
      3  806 1623 (N22E32 L=162)   763    250        1 ( 765/1510- 862/1718) -38.7 E38-E23 N16-N28
      4  465 1593 (N12E54 L=140)   716     94        2 ( 425/1563- 507/1630)  68.2 E58-E51 N10-N15
      5 1100 1542 (N26E11 L=183)   313    104        0 (1084/1487-1116/1588) -56.3 E12-E08 N22-N29
      6  362 1481 (N04E57 L=137)  1489    180        4 ( 319/1388- 416/1543) -60.9 E62-E52 S02-N07
      7 1194 1310 (N14W02 L=196)   363    117        0 (1162/1255-1210/1357) -51.3 E01-W03 N11-N16
      8  399 1067 (S19E40 L=154)  1645    312        3 ( 280/ 952- 501/1169)  48.6 E49-E34 S28-S11
      9  875 1119 (S05E11 L=183)   190     81        0 ( 842/1083- 896/1145) -26.6 E13-E08 S06-S03
     10  505  844 (S29E28 L=166)  2913    671        4 ( 214/ 722- 806/ 934)  27.9 E60-E09 S44-S16
     11 1411  803 (S07W28 L=222)   268     98        0 (1381/ 764-1455/ 831)  68.2 W26-W29 S09-S04
     12 1144  623 (S24W19 L=213)   421    151        0 (1074/ 590-1208/ 653)  58.0 W15-W22 S28-S20
     13 1617  533 (S14W54 L=248)   531     70        0 (1593/ 495-1643/ 569) -78.7 W52-W56 S16-S11
     14  554  479 (S53E18 L=176)   269     84        0 ( 525/ 446- 585/ 516) -33.7 E21-E14 S54-S50
     15 1633  438 (S18W63 L=257)   772    105        2 (1589/ 400-1672/ 468)  71.6 W58-W65 S21-S15
     16 1017  297 (S49W30 L=224)  1848    338        4 ( 859/ 229-1138/ 389)   1.7 W12-W45 S50-S45
  

    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