Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1444 1712 (N52W07 L=274)   379     65        1 (1419/1683-1462/1734) -21.8 W04-W09 N50-N55
      2  904 1676 (N29E29 L=238)   357     65        0 ( 871/1652- 936/1698)   0.0 E30-E25 N27-N30
      3  631 1529 (N13E38 L=228)   617    143        1 ( 565/1485- 691/1580)  58.0 E42-E36 N08-N17
      4  981 1521 (N21E17 L=250)   554    111        0 ( 948/1468-1017/1564) -23.2 E20-E13 N19-N23
      5  528 1411 (N03E41 L=226)   286     86        0 ( 501/1383- 549/1439) -21.8 E43-E38 N02-N04
      6 1250 1352 (N19W04 L=270)   341     95        1 (1213/1307-1276/1402) -56.3 W00-W05 N15-N21
      7  452 1351 (S02E44 L=223)   300     98        1 ( 421/1328- 497/1366)   9.5 E46-E40 S03-N00
      8  877 1314 (N06E16 L=251)   389    121        0 ( 828/1260- 909/1346) -23.2 E19-E12 N03-N07
      9  255 1219 (S13E57 L=210)  1144    314        3 ( 202/1117- 339/1301)  71.6 E61-E50 S19-S07
     10  739 1238 (S01E22 L=245)   465    167        0 ( 688/1202- 818/1298)   0.0 E26-E16 S02-N01
     11 1197 1236 (N11W04 L=270)   246     67        1 (1161/1218-1229/1252)  26.6 W01-W05 N09-N12
     12  347 1083 (S19E45 L=222)  2404    542        9 ( 249/ 938- 405/1247) -54.8 E57-E38 S28-S11
     13 1008  859 (S15W04 L=270)   225     83        0 ( 967/ 837-1057/ 889)  51.3 W01-W05 S17-S12
     14  657  865 (S24E17 L=249)   207     48        0 ( 629/ 851- 688/ 882)  14.0 E19-E15 S25-S23
     15  900  747 (S24W01 L=268)   324    105        0 ( 839/ 732- 950/ 760)  26.6 E03-W03 S25-S22
     16 1407  493 (S23W41 L=308)   474    174        1 (1331/ 464-1481/ 519)  20.0 W35-W46 S24-S20
     17  735  440 (S49W00 L=267)   455     86        0 ( 706/ 418- 777/ 478) -15.9 E03-W04 S50-S47
     18  851  425 (S46W10 L=277)   402     81        0 ( 833/ 393- 873/ 458) -80.5 W07-W12 S48-S42
  

    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