Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  903 1555 (N20E23 L=282)   333     77        0 ( 876/1517- 920/1595) -38.7 E26-E20 N17-N22
      2 1693 1272 (N28W36 L=340)  2245    460        1 (1663/1114-1766/1446) -68.2 W31-W44 N18-N39
      3 1579 1305 (N26W26 L=330)   353    100        0 (1558/1256-1604/1346)  90.0 W24-W27 N22-N29
      4 1877 1244 (N34W58 L=  3)   922    145        1 (1831/1177-1909/1296) -18.4 W51-W63 N31-N36
      5  427 1066 (S19E38 L=267)   251     71        0 ( 398/1039- 455/1097)  45.0 E40-E36 S21-S17
      6  177  993 (S29E64 L=241)  1669    203        3 (  87/ 974- 289/1016)   4.6 E86-E49 S30-S25
      7  857  829 (S22E04 L=301)   310    120        0 ( 835/ 778- 885/ 871) -59.0 E05-E01 S24-S19
      8 1722  839 (N05W48 L=353)   625     82        1 (1678/ 821-1771/ 858)  18.4 W44-W53 N02-N06
      9 1143  639 (S24W19 L=324)   567    244        0 (1033/ 614-1255/ 669)  36.9 W12-W24 S28-S19
     10 1441  481 (S23W45 L=350)  1250    300        3 (1368/ 377-1504/ 586) -76.0 W41-W50 S31-S15
     11 1565  494 (S18W53 L=358)   336     56        0 (1536/ 477-1593/ 513)  45.0 W50-W55 S19-S16
     12 1550  413 (S23W58 L=  3)  1785    191        3 (1445/ 372-1620/ 457)  42.0 W51-W64 S28-S18
     13  663  348 (S61E02 L=302)   310     54        0 ( 634/ 336- 694/ 359)  31.0 E05-W01 S62-S59
  

    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