Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  547 1519 (N48E28 L= 96)   372     23        0 ( 526/1498- 573/1539)  11.3 E30-E25 N46-N49
      2  651 1483 (N43E18 L=107)   488    102        0 ( 614/1456- 688/1521) -40.6 E21-E14 N40-N46
      3 1697 1388 (N12W56 L=180)   765    156        2 (1644/1338-1745/1433) -10.3 W50-W61 N09-N14
      4 1085 1175 (N15W07 L=131)   246    105        0 (1042/1148-1121/1215)  18.4 W03-W09 N13-N16
      5 1846 1136 (S06W61 L=186)   874     94        2 (1808/1106-1884/1172) -36.9 W57-W66 S09-S03
      6  215 1051 (N21E61 L= 64)   666    130        2 ( 194/1009- 270/1102) -36.9 E63-E54 N18-N24
      7  558  774 (N00E34 L= 91)   693    137        1 ( 526/ 714- 597/ 831) -90.0 E35-E32 S03-N03
      8 1552  646 (S28W28 L=152)   393     75        2 (1532/ 605-1577/ 675)  45.0 W25-W29 S30-S25
      9  773  543 (S18E26 L= 99)   254     44        0 ( 760/ 518- 789/ 573)  76.0 E27-E24 S19-S15
  

    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