Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  474 1640 (N22E52 L=270)   352     32        0 ( 449/1627- 495/1659) -11.3 E54-E49 N21-N23
      2  200 1148 (S11E60 L=261)  1264    250        2 ( 140/1083- 266/1236)  45.0 E66-E55 S15-S05
      3 1779 1060 (N12W49 L= 10)   789    239        0 (1711/ 950-1805/1133) -38.7 W41-W52 N05-N15
      4 1850 1070 (N14W56 L= 17)   453    113        1 (1822/1026-1871/1111) -26.6 W52-W58 N12-N16
      5  346  979 (S20E44 L=278)  2149    687        1 ( 166/ 836- 584/1076)  23.8 E60-E26 S30-S13
      6  191  980 (S21E59 L=263)   367     69        1 ( 175/ 949- 215/1007)  76.0 E60-E56 S23-S19
      7 1897  845 (N03W68 L= 30)   585     87        1 (1878/ 817-1914/ 890) -76.0 W65-W71 N01-N05
      8 1525  711 (S14W38 L=360)   244     29        0 (1500/ 692-1544/ 731)  56.3 W36-W39 S15-S12
      9  373  636 (S41E42 L=280)   797    213        3 ( 322/ 574- 421/ 705)  48.0 E48-E37 S46-S36
     10 1767  645 (S12W61 L= 22)   553    166        0 (1718/ 618-1800/ 661)   5.7 W55-W65 S12-S10
     11 1510  637 (S19W39 L=  1)   323    100        0 (1455/ 624-1550/ 658)  26.6 W35-W41 S20-S16
     12 1566  398 (S31W55 L= 17)  2938    489        4 (1436/ 361-1729/ 438)  25.6 W45-W68 S35-S23
  

    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