Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  581 1767 (N29E57 L=164)   453     85        1 ( 569/1729- 600/1803) -36.9 E62-E52 N26-N32
      2  373 1554 (N11E60 L=162)  2033    289        5 ( 290/1465- 469/1675)  79.4 E65-E55 N04-N20
      3  504 1631 (N18E52 L=169)   432     46        1 ( 485/1596- 525/1657) -53.1 E54-E50 N15-N19
      4  677 1586 (N18E37 L=185)   630    178        1 ( 638/1538- 738/1656) -21.0 E43-E30 N15-N22
      5  560 1496 (N10E41 L=181)   838    213        1 ( 528/1399- 609/1584) -39.8 E45-E33 N05-N15
      6  448 1239 (S07E41 L=181)   248     30        1 ( 429/1216- 468/1262) -14.0 E42-E38 S08-S05
      7  360 1185 (S12E47 L=175)   542    156        0 ( 318/1148- 417/1211)  14.0 E49-E41 S14-S10
      8  968 1133 (S02E06 L=216)   218     41        1 ( 940/1112- 988/1151)   0.0 E07-E04 S03-S01
      9  536  925 (S24E28 L=193)   322     97        0 ( 502/ 899- 580/ 966) -23.2 E31-E24 S24-S21
     10  257  861 (S33E54 L=168)   360     63        0 ( 241/ 829- 281/ 884) -90.0 E55-E50 S34-S31
     11  236  779 (S38E58 L=164)   813     77        1 ( 193/ 757- 276/ 805)   4.8 E64-E52 S39-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