Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  806 1805 (N46E30 L=260)  1140    102        2 ( 755/1777- 866/1832) -13.0 E37-E24 N42-N48
      2 1169 1699 (N41W03 L=293)   415    106        0 (1121/1678-1222/1717)  15.9 W00-W07 N38-N42
      3 1625 1617 (N42W44 L=335)   693    124        1 (1577/1591-1680/1640)   0.0 W38-W50 N39-N43
      4  482 1475 (N18E42 L=249)  1549    295        3 ( 397/1383- 592/1580) -26.6 E53-E31 N12-N25
      5 1401 1229 (N11W21 L=311)   451    208        0 (1348/1170-1471/1312) -41.2 W17-W26 N07-N15
      6 1749 1269 (N18W48 L=338)   605    155        0 (1686/1235-1816/1305)  -8.7 W42-W55 N17-N19
      7 1686 1155 (N10W42 L=332)   364     70        0 (1657/1130-1720/1188) -18.4 W39-W45 N08-N11
      8  684  716 (S27E20 L=271)  1379    425        0 ( 580/ 541- 803/ 873)  53.5 E29-E12 S39-S16
      9 1665  790 (S12W45 L=336)   249     50        0 (1652/ 767-1681/ 825) -45.0 W43-W46 S12-S09
     10 1792  590 (S21W68 L=358)  1209     97        3 (1754/ 539-1834/ 621) -21.8 W61-W75 S23-S19
     11  955  582 (S33W00 L=290)   203     38        0 ( 930/ 568- 982/ 597)  45.0 E01-W01 S34-S32
  

    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