Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1644 1474 (N13W54 L=340)   592    111        1 (1601/1453-1676/1500) -15.9 W50-W57 N11-N14
      2 1674 1139 (S05W43 L=329)   442    166        1 (1614/1071-1727/1177) -53.1 W38-W45 S10-S02
      3 1205 1114 (N07W12 L=299)   218     78        0 (1170/1082-1226/1156)  36.9 W09-W14 N05-N09
      4  947 1074 (N12E03 L=283)   488    162        1 ( 869/1053-1026/1097) -18.4 E07-W02 N10-N13
      5  751 1053 (N15E15 L=271)   410     86        2 ( 718/1014- 783/1083)  36.9 E17-E12 N13-N16
      6  301  846 (N13E52 L=234)  3020    606       10 ( 158/ 772- 483/ 930) -26.6 E66-E39 N05-N19
      7 1478  833 (S16W21 L=308)   338     55        1 (1448/ 806-1512/ 858) -36.9 W19-W23 S17-S14
      8 1344  634 (S24W08 L=295)   607    185        2 (1293/ 558-1374/ 735)  45.0 W03-W12 S27-S18
      9  558  545 (S10E42 L=244)   396     98        0 ( 539/ 504- 578/ 587)  63.4 E43-E40 S12-S06
  

    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