Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1619 1247 (N15W42 L=336)   431    125        0 (1562/1206-1659/1278) -39.8 W37-W45 N12-N17
      2 1319 1143 (N10W19 L=313)   331     58        1 (1297/1115-1344/1180)  76.0 W17-W20 N08-N12
      3 1414 1101 (N07W25 L=319)   276    112        1 (1378/1072-1471/1124)  23.2 W22-W29 N05-N08
      4  399 1042 (N06E42 L=252)   323     83        0 ( 370/1001- 429/1074) -38.7 E44-E39 N03-N07
      5 1869  854 (S11W66 L=  0)  1333    156        4 (1829/ 777-1901/ 919) -42.3 W59-W74 S16-S06
      6  451  849 (S06E39 L=255)   923    210        0 ( 361/ 789- 523/ 887) -22.6 E45-E33 S09-S03
      7  256  653 (S19E63 L=231)   469     82        0 ( 245/ 619- 274/ 701)  59.0 E65-E61 S21-S16
      8  509  562 (S25E39 L=255)   665    232        0 ( 461/ 500- 542/ 648)  51.3 E45-E35 S29-S18
      9  757  516 (S28E21 L=274)   472    115        0 ( 716/ 474- 791/ 561)  40.6 E24-E17 S31-S24
     10 1630  502 (S33W48 L=343)   766    118        0 (1592/ 472-1684/ 541) -14.0 W43-W55 S35-S30
  

    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