Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1258 1719 (N41W44 L=238)   406     63        0 (1216/1706-1295/1730) -31.0 W40-W47 N39-N42
      2  964 1639 (N46W17 L=212)   243     50        1 ( 939/1621- 995/1651) -45.0 W14-W19 N44-N47
      3 1520 1574 (N24W51 L=246)   410     91        1 (1480/1544-1555/1603) -68.2 W50-W53 N21-N26
      4 1451 1442 (N19W39 L=234)  1402    406        0 (1282/1336-1603/1520) -53.6 W30-W45 N08-N27
      5 1147 1350 (N22W16 L=211)  1994    493        4 ( 921/1264-1347/1408) -11.7 W00-W29 N18-N25
      6 1404 1322 (N14W31 L=226)   867    156        5 (1339/1254-1446/1374) -76.0 W28-W33 N09-N17
      7 1475 1222 (N06W32 L=227)   424     68        0 (1456/1179-1496/1257)  36.9 W29-W33 N04-N08
      8  909 1067 (N12E06 L=189)   360    166        0 ( 842/1021- 982/1096)   5.7 E10-W00 N11-N13
      9 1881  900 (S26W61 L=255)  1101     97        1 (1846/ 859-1911/ 935)  11.3 W55-W66 S27-S23
     10 1596  505 (S42W27 L=222)  1231    188        1 (1511/ 449-1658/ 569)   7.6 W18-W33 S45-S38
     11  436  402 (S18E62 L=133)   719    125        1 ( 398/ 361- 474/ 461)  82.9 E64-E59 S21-S13
  

    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