Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1356 1726 (N46W17 L=258)  1337    178        2 (1274/1685-1431/1756)  -8.1 W09-W24 N43-N49
      2  478 1555 (N22E46 L=195)   357     63        0 ( 461/1524- 496/1586)  90.0 E47-E44 N20-N24
      3 1199 1409 (N20W06 L=247)  1109    330        2 (1083/1294-1312/1507) -33.7 E01-W14 N13-N25
      4 1746 1389 (N27W49 L=290)   917    139        2 (1702/1331-1789/1451) -45.0 W45-W53 N24-N31
      5 1842 1275 (N22W59 L=299)   992    202        1 (1822/1175-1864/1358)  71.6 W55-W62 N15-N27
      6 1614 1317 (N20W36 L=277)   486     76        0 (1580/1283-1657/1347) -26.6 W33-W39 N18-N21
      7 1742 1310 (N22W47 L=288)   314     57        1 (1724/1287-1767/1340) -45.0 W45-W49 N20-N23
      8  290  980 (S14E50 L=190)   719    178        2 ( 246/ 912- 337/1061)  45.0 E55-E46 S18-S09
      9  647  869 (S19E22 L=219)   222     84        0 ( 610/ 841- 683/ 887)  31.0 E24-E19 S20-S17
     10  199  845 (S23E62 L=178)   453     66        0 ( 163/ 829- 232/ 864)  11.3 E67-E57 S23-S21
     11  473  704 (S30E36 L=205)   517     86        0 ( 441/ 662- 507/ 739)  35.5 E39-E32 S33-S28
     12 1340  585 (S28W28 L=268)  1514    387        1 (1185/ 478-1492/ 676)  46.7 W19-W35 S37-S20
  

    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