Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1403 1752 (N41W58 L=292)   616    109        2 (1374/1712-1446/1778)   8.7 W51-W64 N39-N42
      2  502 1399 (N40E31 L=204)   691    113        1 ( 454/1378- 562/1427) -24.4 E36-E25 N37-N42
      3  852  977 (N08E11 L=224)   252     52        0 ( 835/ 953- 879/1004) -90.0 E12-E09 N06-N09
      4 1827  862 (S22W55 L=289)  1733    163        4 (1787/ 783-1864/ 938)  48.8 W49-W59 S27-S17
      5 1702  614 (S35W42 L=276)   501    118        0 (1668/ 567-1731/ 666)  39.8 W38-W44 S38-S32
      6  765  618 (S12E24 L=210)   501    154        0 ( 728/ 560- 813/ 656)  39.8 E27-E20 S15-S10
      7  635  610 (S11E33 L=202)   284     30        0 ( 610/ 587- 659/ 633)  14.0 E34-E30 S11-S09
      8  773  490 (S21E28 L=207)   321    119        0 ( 713/ 477- 824/ 513)  -6.3 E32-E23 S21-S19
      9 1499  403 (S46W22 L=257)   416     74        0 (1471/ 370-1526/ 434)  38.7 W19-W24 S48-S43
  

    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