Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1326 1758 (N49E01 L= 51)   548     85        1 (1288/1717-1355/1784) -26.6 E04-W02 N46-N52
      2  842 1695 (N29E31 L= 20)   400     53        0 ( 819/1662- 861/1727) -45.0 E33-E28 N26-N31
      3 1341 1518 (N30W07 L= 58)   990    259        1 (1266/1421-1399/1602) -30.3 W00-W12 N24-N34
      4  989 1384 (N12E11 L= 41)   461    149        2 ( 972/1307-1009/1451) -71.6 E12-E08 N07-N16
      5 1342 1186 (N10W14 L= 65)   273     72        1 (1324/1145-1367/1220)  90.0 W12-W14 N07-N12
      6  232 1035 (S22E56 L=355)  1976    208        3 ( 171/ 958- 315/1103)  35.5 E64-E47 S27-S17
      7  115 1046 (S23E74 L=337)   773     95        3 (  95/1004- 132/1094)  56.3 E79-E70 S25-S19
      8  566  806 (S30E24 L= 27)   269     94        0 ( 525/ 780- 602/ 841)  45.0 E27-E21 S32-S27
      9 1373  376 (S34W46 L= 98)   398     57        1 (1346/ 357-1405/ 399)   0.0 W42-W49 S35-S32
     10 1168  131 (S57W68 L=119)   813     72        0 (1152/  90-1185/ 167)  -7.6 W58-W88 S58-S54
  

    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