Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1245 1829 (N57E16 L=353)  1472    179        2 (1171/1787-1328/1857)  13.2 E24-E06 N53-N60
      2  617 1701 (N25E50 L=318)  2956    391        6 ( 572/1553- 688/1835) -24.1 E74-E36 N16-N34
      3 1665 1625 (N57W40 L= 49)   890     89        1 (1642/1576-1683/1666)  38.7 W35-W52 N52-N64
      4 1706 1229 (N27W39 L= 48)   298     72        1 (1677/1197-1731/1260)  -9.5 W36-W42 N25-N28
      5 1861  909 (N13W62 L= 71)   748    122        4 (1831/ 854-1884/ 964)  69.4 W59-W65 N09-N17
      6  502  880 (S26E29 L=340)   691    185        2 ( 452/ 800- 535/ 953) -45.0 E34-E25 S30-S22
      7  671  734 (S30E14 L=355)   262    111        0 ( 620/ 718- 723/ 759)  39.8 E17-E10 S32-S27
      8  801  653 (S32E03 L=  6)   263     88        1 ( 762/ 635- 837/ 678)  51.3 E05-E01 S33-S28
      9 1355  592 (S19W33 L= 42)   351    115        0 (1293/ 578-1409/ 607)  33.7 W29-W36 S21-S16
     10  624  333 (S62E09 L=360)   630     75        2 ( 595/ 302- 669/ 363)  26.6 E13-E03 S66-S58
     11 1406  284 (S35W59 L= 68)  1365    141        3 (1349/ 261-1465/ 317)  18.4 W52-W65 S38-S31
  

    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