Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  351 1462 (N50E43 L=  9)   485    118        1 ( 304/1440- 403/1478)  -8.7 E50-E36 N47-N51
      2  656 1320 (N32E15 L= 37)  1596    400        6 ( 466/1262- 823/1359) -17.1 E30-E04 N24-N35
      3  481 1143 (N26E31 L= 21)   320     74        0 ( 450/1120- 518/1168) -45.0 E32-E28 N23-N27
      4  680 1054 (N16E19 L= 33)   223    111        0 ( 624/1036- 727/1065) -33.7 E22-E16 N13-N17
      5  270 1026 (N23E51 L=  1)   456     81        1 ( 241/ 992- 303/1055) -45.0 E54-E48 N20-N25
      6 1656  879 (S21W35 L= 87)   469     73        1 (1635/ 832-1675/ 921)  68.2 W32-W36 S23-S18
      7  728  606 (S11E29 L= 24)   313     94        1 ( 687/ 581- 778/ 627)  -9.5 E31-E25 S11-S09
      8  542  449 (S16E48 L=  4)   330     72        1 ( 517/ 409- 563/ 473)  78.7 E49-E46 S18-S13
      9 1062  331 (S37E20 L= 32)  1548    273        3 ( 959/ 255-1153/ 430)  18.4 E31-E10 S40-S32
     10  634  284 (S29E54 L=359)  1352    173        1 ( 604/ 208- 661/ 353)  26.6 E66-E46 S34-S24
     11  889  237 (S39E39 L= 14)   632     83        2 ( 855/ 219- 939/ 263) -38.7 E41-E34 S41-S36
  

    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