Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1315 1740 (N40W50 L=241)   509     75        1 (1275/1725-1352/1757) -45.0 W47-W53 N38-N42
      2 1498 1624 (N27W53 L=244)  1107    260        0 (1377/1551-1595/1692) -57.3 W48-W57 N20-N34
      3  988 1664 (N47W20 L=211)   304     84        0 ( 941/1653-1030/1675) -23.2 W16-W23 N44-N48
      4 1493 1466 (N19W43 L=234)  1545    451        1 (1311/1336-1653/1543) -54.5 W33-W49 N06-N27
      5 1860 1086 (S14W60 L=251)   419     30        1 (1841/1067-1883/1109)  45.0 W57-W63 S14-S12
      6  151  962 (N21E70 L=121)   912    138        1 ( 117/ 911- 200/1015) -26.6 E76-E63 N17-N23
      7 1904  917 (S26W64 L=255)  1194     91        1 (1872/ 871-1937/ 953)  11.3 W58-W72 S27-S23
      8 1304  876 (S08W12 L=203)   416    173        0 (1283/ 807-1329/ 951)  69.4 W10-W14 S12-S03
      9 1635  528 (S42W31 L=222)  1561    211        2 (1545/ 462-1701/ 607)  14.9 W21-W38 S45-S37
     10  460  419 (S17E58 L=133)   948    171        2 ( 428/ 347- 509/ 495)  65.6 E60-E53 S22-S11
     11  730  258 (S34E48 L=143)   460     53        1 ( 704/ 235- 756/ 279)  78.7 E49-E44 S36-S31
  

    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