Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1036 1687 (N50W12 L=140)   532     85        1 (1005/1645-1059/1720)  50.2 W09-W15 N46-N52
      2 1263 1667 (N44W31 L=158)   345     66        0 (1246/1636-1285/1690)  38.7 W28-W33 N41-N45
      3  959 1606 (N44W03 L=131)   488     94        2 ( 929/1573- 989/1641) -59.0 W01-W06 N41-N47
      4  577 1456 (N38E29 L= 99)   287     58        1 ( 561/1433- 604/1483) -36.9 E30-E26 N36-N39
      5 1248 1046 (N04W14 L=142)   483    161        1 (1173/1020-1328/1091) -24.0 W09-W18 N02-N08
      6 1908  948 (S13W68 L=196)   666    128        1 (1882/ 918-1928/ 985)  45.0 W64-W72 S14-S10
      7 1116  509 (S28E00 L=127)   664    201        1 (1035/ 445-1184/ 586)  38.7 E06-W04 S32-S23
      8  354  486 (S23E59 L= 69)  1067    322        0 ( 314/ 440- 406/ 524)  21.8 E66-E51 S27-S20
  

    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