Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1196 1775 (N57W15 L=107)   353     64        1 (1167/1752-1222/1794)  14.0 W11-W19 N54-N59
      2  392 1510 (N31E55 L= 37)   521     84        2 ( 350/1493- 435/1527)   0.0 E60-E49 N29-N32
      3  718 1378 (N24E22 L= 69)   613    202        1 ( 631/1341- 800/1417)  -8.7 E29-E16 N21-N26
      4  865 1251 (N16E11 L= 81)   248    104        0 ( 807/1236- 915/1262)  -8.1 E14-E07 N15-N17
      5 1294 1199 (N14W17 L=108)   242    107        0 (1242/1175-1336/1229) -18.4 W13-W19 N12-N16
      6 1495  729 (S15W32 L=124)   504    224        1 (1452/ 639-1546/ 802)  65.6 W29-W35 S20-S09
      7  729  711 (S18E18 L= 73)   271     97        0 ( 690/ 696- 779/ 724)   0.0 E20-E14 S18-S16
      8  914  460 (S35E06 L= 86)   991    259        3 ( 859/ 354- 981/ 556) -66.4 E09-W00 S43-S27
      9 1479  504 (S30W36 L=128)   361     68        0 (1449/ 479-1517/ 530)  45.0 W33-W38 S31-S27
  

    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