Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1186 1634 (N35W09 L=293)   521     82        1 (1164/1594-1209/1684)  81.9 W07-W10 N32-N39
      2 1374 1623 (N36W23 L=308)   360     39        1 (1347/1594-1403/1642)   0.0 W21-W25 N33-N37
      3 1594 1469 (N26W38 L=323)   857    183        1 (1521/1423-1680/1513) -18.4 W32-W45 N23-N28
      4 1842 1265 (N15W59 L=343)  3543    281        9 (1745/1153-1904/1356) -85.6 W48-W70 N08-N21
      5 1794  959 (S03W53 L=337)   298     49        1 (1775/ 932-1811/ 989) -71.6 W51-W54 S04-S01
      6  176  881 (S14E63 L=222)  1138    116        1 ( 141/ 834- 199/ 946)  50.2 E68-E59 S16-S10
      7 1206  615 (S28W14 L=299)  1773    463        1 (1117/ 488-1331/ 733)  39.8 W08-W21 S36-S20
      8  352  570 (S34E52 L=232)  1063    223        1 ( 274/ 509- 433/ 603)   9.5 E67-E43 S37-S31
      9 1360  542 (S32W26 L=311)   547     98        0 (1310/ 517-1410/ 579)  33.7 W22-W29 S33-S29
     10 1768  558 (S27W64 L=348)  1012     98        0 (1707/ 541-1817/ 575)   8.1 W55-W71 S28-S25
     11 1559  350 (S44W55 L=340)  2883    275        3 (1427/ 324-1689/ 384)   3.0 W39-W77 S45-S40
     12  759  336 (S51E21 L=264)  1989    379        4 ( 573/ 310- 915/ 381)   1.7 E40-E05 S52-S46
  

    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