Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1605 1711 (N28W73 L= 43)   931     90        0 (1575/1677-1636/1743)  51.3 W69-W76 N25-N30
      2 1527 1657 (N28W58 L= 28)   666     69        1 (1502/1621-1561/1691)  68.2 W54-W61 N25-N30
      3 1064 1574 (N38W19 L=350)   458    100        1 (1014/1533-1099/1597) -63.4 W16-W22 N34-N40
      4 1644 1546 (N19W59 L= 30)  2143    341        4 (1480/1508-1767/1583) -21.4 W44-W74 N13-N24
      5  294 1382 (N45E53 L=278)   859     93        0 ( 252/1362- 340/1402) -15.9 E59-E45 N42-N46
      6  485 1314 (N36E30 L=300)   362     67        2 ( 453/1284- 518/1339)  36.9 E33-E27 N34-N37
      7 1501 1230 (N07W33 L=  4)   222     79        0 (1489/1191-1511/1273)  78.7 W31-W34 N04-N09
      8 1381  864 (S10W17 L=347)   873    228        3 (1353/ 760-1407/ 988)  61.7 W13-W20 S16-S03
      9  720  709 (S04E25 L=306)   389    103        1 ( 696/ 656- 746/ 766)  45.0 E27-E22 S06-S01
     10 1055  376 (S33E16 L=315)  1823    349        7 ( 912/ 286-1135/ 496)  21.4 E29-E06 S37-S26
     11  556  411 (S19E48 L=283)   543     35        1 ( 528/ 386- 581/ 439)  18.4 E50-E44 S21-S17
     12  767  395 (S25E34 L=297)   454    112        1 ( 732/ 351- 797/ 439)  32.0 E37-E29 S27-S22
  

    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