Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  894 1723 (N41E28 L=190)   382     54        0 ( 874/1692- 914/1754) -45.0 E31-E25 N38-N43
      2 1294 1421 (N28W10 L=229)  1928    541        4 (1152/1189-1431/1557) -36.9 E02-W22 N16-N35
      3 1168 1203 (N12W05 L=224)   685    236        1 (1068/1172-1272/1243)   0.0 E01-W11 N10-N13
      4 1767 1075 (N16W51 L=270)   665     86        2 (1748/1029-1789/1120)  71.6 W48-W52 N13-N19
      5  255  937 (S20E54 L=165)  2849    416        4 ( 134/ 807- 371/1079)  31.0 E78-E43 S30-S09
      6 1783  876 (N05W56 L=274)   533    116        1 (1756/ 818-1812/ 926)  76.0 W53-W58 N00-N08
      7 1106  721 (S18W11 L=230)   318     41        1 (1081/ 698-1131/ 746) -26.6 W08-W12 S19-S16
      8 1632  663 (S11W49 L=267)  1046    135        2 (1553/ 633-1686/ 696)  29.1 W42-W53 S13-S08
      9  848  296 (S54W05 L=224)   693    101        1 ( 797/ 269- 894/ 327)   5.7 W00-W10 S56-S52
  

    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