Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  800 1855 (N39E52 L=151)   510     91        0 ( 772/1817- 822/1895) -18.4 E61-E46 N35-N42
      2 1294 1695 (N43E05 L=198)   566    185        0 (1233/1617-1343/1759) -28.3 E12-W01 N38-N46
      3  435 1402 (S00E48 L=155)   779    130        1 ( 405/1335- 469/1443) -45.0 E50-E43 S03-N02
      4 1356 1286 (N18W11 L=214)   701    169        0 (1294/1211-1409/1338) -16.7 W06-W16 N14-N19
      5  937  839 (S19W00 L=203)  1062    474        0 ( 721/ 749-1126/ 881)  43.5 E10-W09 S30-S11
      6  507  711 (S39E26 L=177)   912    305        1 ( 471/ 583- 577/ 813)  86.4 E28-E21 S46-S30
      7 1039  677 (S25W11 L=214)   517    174        0 ( 970/ 641-1123/ 706)  48.8 W08-W15 S28-S20
      8 1133  556 (S29W22 L=225)   429    116        1 (1076/ 542-1194/ 574)  29.7 W18-W25 S30-S26
      9  442  429 (S63E35 L=168)  1299    162        2 ( 376/ 396- 516/ 483)  -5.2 E46-E22 S65-S59
     10 1386  431 (S27W45 L=248)   671    198        1 (1314/ 396-1467/ 473)   3.6 W37-W53 S29-S25
     11 1617  399 (S20W67 L=270)   929    159        2 (1555/ 358-1680/ 448)  66.0 W62-W71 S24-S15
  

    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