Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1468 1432 (N27W35 L=301)  1273    285        3 (1404/1352-1568/1536) -65.2 W30-W40 N21-N34
      2  267 1429 (N32E66 L=200)   762     78        0 ( 236/1400- 297/1459) -11.3 E75-E60 N29-N33
      3  719 1157 (N14E19 L=247)   343    118        2 ( 686/1109- 768/1183)  38.7 E21-E15 N11-N15
      4 1477  854 (S09W28 L=295)   632    169        3 (1452/ 783-1516/ 929)  58.0 W26-W31 S12-S04
      5  778  891 (S03E16 L=250)   207     93        0 ( 725/ 876- 806/ 914)  21.8 E19-E14 S03-S01
      6 1747  833 (S12W51 L=317)   992    196        2 (1700/ 737-1771/ 894)  73.3 W47-W53 S18-S08
      7 1332  812 (S10W18 L=285)   237     79        1 (1300/ 782-1367/ 839)  26.6 W16-W20 S11-S08
      8 1870  746 (S20W70 L=337)  1669    147        5 (1823/ 679-1901/ 794) -53.1 W61-W81 S24-S16
      9  420  588 (S22E47 L=219)   674    147        1 ( 399/ 527- 446/ 658)  78.7 E49-E44 S26-S16
     10  284  561 (S24E66 L=200)   810     98        1 ( 256/ 526- 314/ 596)  39.8 E71-E58 S26-S21
     11  616  559 (S24E31 L=235)   277     78        0 ( 587/ 527- 644/ 594)  33.7 E34-E28 S26-S21
     12  701  517 (S27E26 L=240)   334    106        2 ( 652/ 488- 736/ 556)  26.6 E30-E22 S29-S24
  

    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