Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1225 1816 (N51E11 L= 83)   457     41        0 (1196/1801-1254/1838)  45.0 E12-E08 N48-N53
      2  797 1749 (N32E37 L= 57)   507     84        1 ( 769/1706- 824/1787) -29.1 E41-E32 N29-N34
      3  431 1648 (N18E60 L= 33)  1337    122        5 ( 388/1613- 492/1696)  18.4 E67-E53 N15-N21
      4  775 1600 (N21E30 L= 63)   385     51        0 ( 750/1571- 803/1633)  90.0 E31-E29 N19-N23
      5 1147 1442 (N20E03 L= 91)   213    101        0 (1127/1412-1179/1494) -38.7 E05-W00 N18-N22
      6 1207  685 (S22W20 L=113)   263     82        0 (1169/ 673-1255/ 705)  21.8 W17-W22 S22-S20
      7 1551  371 (S29W62 L=155)   889    140        0 (1498/ 314-1607/ 417)  90.0 W59-W64 S33-S24
      8  979  226 (S59W32 L=125)  1477    238        2 ( 907/ 138-1056/ 312) -11.9 W17-W55 S62-S53
  

    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