Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1185 1815 (N49E15 L= 53)   564     65        0 (1158/1782-1218/1838)  74.1 E16-E12 N45-N52
      2 1481 1725 (N53W15 L= 83)   472     51        0 (1462/1694-1505/1749)  60.3 W13-W18 N49-N56
      3  508 1709 (N23E58 L=  9)   654     72        0 ( 485/1670- 533/1748) -28.6 E64-E52 N20-N26
      4 1123 1625 (N31E10 L= 58)   672    146        0 (1075/1548-1162/1677) -33.7 E14-E05 N26-N34
      5 1153 1456 (N21E03 L= 64)   207     47        0 (1127/1442-1186/1469)  33.7 E04-E01 N19-N21
      6  131 1096 (S19E72 L=356)   649     58        0 (  98/1081- 162/1112)   0.0 E78-E66 S20-S18
      7 1219  458 (S34W30 L= 98)   256     65        0 (1191/ 445-1267/ 473)  31.0 W28-W33 S35-S32
      8 1068  207 (S56W42 L=109)   436     57        0 (1047/ 178-1087/ 239) -29.1 W38-W47 S58-S53
      9 1117  145 (S59W59 L=126)  1125    111        0 (1097/  87-1137/ 197)  -7.1 W47-W87 S61-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