Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1592 1576 (N45W25 L=142)   420     54        0 (1564/1551-1618/1604) -31.0 W21-W27 N42-N46
      2  896 1484 (N16E20 L= 98)   469    178        1 ( 871/1403- 929/1564) -52.1 E23-E15 N11-N20
      3 1019 1448 (N17E12 L=106)   259     72        0 ( 982/1435-1057/1463)  11.3 E14-E09 N15-N18
      4 1789 1420 (N41W47 L=165)   448    116        0 (1748/1372-1825/1461) -18.4 W42-W52 N39-N43
      5 1120 1347 (N13E03 L=115)   216     83        0 (1101/1302-1144/1390) -38.7 E05-W00 N11-N15
      6 1811 1107 (N21W50 L=167)   789    175        0 (1797/1028-1825/1193) -68.2 W47-W52 N16-N26
      7 1575  887 (N01W36 L=154)   446     86        0 (1538/ 847-1608/ 920)   0.0 W33-W39 N00-N02
      8  846  597 (S36W01 L=119)  1082    288        2 ( 727/ 564- 987/ 648)  42.7 E06-W08 S41-S29
      9 1039  490 (S37W19 L=136)   377     78        0 (1023/ 450-1055/ 541) -63.4 W16-W20 S39-S33
     10 1477  445 (S25W50 L=167)   525    147        0 (1428/ 394-1521/ 510) -77.5 W47-W51 S29-S20
     11  971  151 (S64W51 L=168)   646     74        0 ( 951/ 118- 995/ 188) -24.0 W42-W63 S67-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