Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  784 1648 (N52E07 L=249)   338     82        0 ( 762/1616- 815/1678)  33.7 E09-E03 N49-N54
      2 1751 1365 (N12W60 L=316)  1211    243        2 (1678/1287-1830/1462) -60.3 W56-W66 N05-N19
      3  212 1075 (N19E61 L=195)   979    163        0 ( 164/1047- 282/1108)   3.8 E68-E52 N17-N20
      4 1235  905 (S04W11 L=267)   224     86        1 (1188/ 886-1269/ 916) -26.6 W08-W12 S05-S02
      5 1595  822 (S16W34 L=290)   248    108        0 (1585/ 767-1605/ 876)  81.9 W33-W34 S19-S12
      6 1337  516 (S31W12 L=269)  2650    562        8 (1121/ 383-1520/ 700)  23.6 E05-W27 S37-S22
      7  289  614 (S10E62 L=194)   477     24        0 ( 272/ 589- 307/ 635)  36.9 E65-E59 S11-S08
      8  971  583 (S20E11 L=245)   224     81        1 ( 939/ 547-1002/ 616) -68.2 E12-E09 S22-S17
      9  777  409 (S29E29 L=227)   516    155        0 ( 763/ 344- 805/ 484)  65.6 E31-E25 S34-S23
     10  558  360 (S30E50 L=206)   745    142        0 ( 536/ 299- 602/ 418)  29.4 E58-E42 S34-S25
  

    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