Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1245 1665 (N39E07 L=339)   396     90        1 (1206/1643-1270/1697)  80.5 E09-E04 N36-N42
      2  361 1626 (N13E68 L=279)   817    157        1 ( 323/1595- 411/1659)  26.6 E71-E62 N10-N15
      3  756 1633 (N22E35 L=311)   284     63        1 ( 727/1615- 788/1651)  11.3 E37-E32 N20-N22
      4  283 1283 (S10E56 L=290)   341     62        0 ( 254/1256- 306/1311)  53.1 E58-E54 S11-S07
      5  392 1140 (S16E42 L=304)   490    163        4 ( 332/1096- 415/1204)  76.0 E46-E40 S19-S11
      6 1332 1127 (N08W14 L=  0)   237     89        0 (1315/1082-1348/1179) -80.5 W13-W15 N05-N11
      7 1875 1052 (N22W56 L= 43)  1079    189        2 (1832/1001-1916/1106) -10.3 W51-W62 N20-N24
      8  370 1020 (S23E41 L=305)   670    193        2 ( 305/ 994- 450/1059)   4.4 E47-E34 S24-S21
      9 1125  511 (S32W23 L=  9)  1502    430        4 ( 955/ 433-1309/ 598)   4.1 W09-W37 S33-S28
  

    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