Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1400 1740 (N56W08 L=219)   607     66        1 (1369/1719-1429/1765)  39.8 W05-W11 N52-N59
      2 1157 1610 (N36E08 L=203)  1039    297        1 (1094/1475-1243/1713) -35.8 E16-W02 N28-N42
      3 1661 1496 (N44W37 L=248)   900    119        1 (1618/1451-1693/1530) -26.6 W32-W41 N40-N46
      4 1481  778 (S06W34 L=245)   260     17        0 (1465/ 756-1498/ 803) -33.7 W32-W35 S07-S04
      5  311  713 (S39E48 L=163)   577    111        1 ( 290/ 662- 353/ 739)  50.2 E50-E42 S42-S36
      6 1384  702 (S12W30 L=241)   293     90        0 (1356/ 661-1414/ 735)  90.0 W28-W31 S15-S09
      7  372  559 (S49E44 L=167)   823     75        0 ( 337/ 522- 402/ 598)  35.0 E49-E38 S51-S44
      8  410  494 (S53E41 L=170)   429     71        0 ( 382/ 473- 452/ 514)  22.6 E45-E33 S55-S50
      9 1468  387 (S28W53 L=263)  2890    461        3 (1253/ 323-1645/ 454)  35.8 W39-W65 S38-S19
     10  503  406 (S58E30 L=181)   421     30        1 ( 480/ 384- 525/ 423)   0.0 E33-E26 S60-S55
  

    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