Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  902 1844 (N55E37 L=203)   946    108        2 ( 848/1815- 957/1864)   0.0 E45-E29 N52-N56
      2 1286 1713 (N52W09 L=249)   317     39        1 (1261/1692-1312/1729)  21.8 W06-W11 N49-N53
      3  561 1662 (N33E51 L=190)   959    164        1 ( 497/1614- 619/1712) -13.4 E63-E42 N30-N35
      4 1456 1598 (N45W25 L=265)   608     73        0 (1425/1567-1488/1623)  -9.5 W21-W28 N42-N47
      5  934 1516 (N29E15 L=225)   301     91        0 ( 890/1498- 985/1533)  26.6 E17-E11 N27-N30
      6 1740 1210 (N22W49 L=289)  1002    183        1 (1670/1142-1800/1275) -26.6 W42-W55 N18-N25
      7  242  773 (S28E58 L=182)   739    138        2 ( 200/ 720- 266/ 840)  39.3 E66-E53 S32-S23
      8 1711  651 (S12W56 L=296)   380     33        1 (1694/ 628-1733/ 673) -45.0 W53-W58 S13-S10
  

    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