Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  902 1874 (N56E23 L=288)   904    123        1 ( 864/1831- 961/1897)  66.0 E27-E16 N50-N59
      2  919 1728 (N41E14 L=297)   828    138        4 ( 859/1685- 977/1778) -28.6 E19-E08 N38-N45
      3 1133 1574 (N31W04 L=315)   357    151        0 (1114/1499-1151/1636) -78.7 W02-W05 N25-N35
      4  386 1588 (N29E55 L=256)   546     76        0 ( 355/1559- 416/1614)  45.0 E57-E51 N26-N30
      5  519 1554 (N26E40 L=271)   393     89        0 ( 502/1506- 548/1595)  90.0 E41-E38 N22-N29
      6  653 1561 (N27E30 L=281)   534    114        2 ( 600/1514- 690/1590) -26.6 E34-E26 N24-N29
      7 1217 1372 (N18W10 L=321)   213     65        0 (1184/1358-1251/1386)  11.3 W07-W12 N16-N18
      8  695  830 (S18E20 L=291)   326     99        0 ( 678/ 784- 717/ 886)  71.6 E20-E18 S20-S14
      9  478  760 (S23E36 L=275)   917    219        1 ( 402/ 685- 579/ 812)  23.6 E43-E27 S27-S19
     10  857  733 (S23E09 L=302)   206     67        0 ( 817/ 723- 896/ 744)   0.0 E11-E06 S23-S22
     11  543  526 (S38E35 L=276)   277     49        1 ( 516/ 509- 563/ 547)  31.0 E37-E32 S39-S36
     12  823  494 (S39E11 L=299)   617    105        1 ( 781/ 465- 875/ 535)  15.9 E14-E07 S41-S36
  

    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