Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1689 1301 (N22W41 L=199)  1550    287        1 (1613/1179-1762/1405) -47.5 W34-W47 N15-N27
      2  914 1290 (N08E10 L=148)   255     82        0 ( 880/1257- 950/1327) -31.0 E12-E07 N06-N10
      3 1113 1136 (N02W04 L=162)   371    195        1 (1047/1096-1165/1175)  18.4 W00-W07 S01-N03
      4 1556 1072 (N05W32 L=190)   212     60        0 (1531/1047-1583/1098) -14.0 W30-W34 N04-N06
      5 1156 1008 (S05W08 L=166)   335    167        0 (1076/ 975-1229/1032)   0.0 W03-W12 S06-S04
      6  225  885 (S23E58 L=100)   735     96        4 ( 197/ 849- 261/ 932) -59.0 E61-E53 S25-S20
      7  166  869 (S24E66 L= 92)   912    120        2 ( 143/ 822- 191/ 920)  63.4 E70-E62 S26-S20
      8 1850  745 (S07W66 L=225)   536     32        2 (1826/ 722-1868/ 772)  90.0 W63-W69 S08-S05
      9 1246  683 (S23W19 L=178)   331     94        1 (1195/ 668-1288/ 702)   0.0 W15-W22 S23-S21
     10 1755  613 (S17W62 L=220)   791    136        0 (1696/ 587-1816/ 635)  22.6 W56-W68 S19-S14
  

    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