Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1232 1784 (N56E05 L= 35)  1779    236        3 (1134/1741-1330/1816)  35.2 E14-W05 N48-N62
      2 1148 1406 (N24E00 L= 40)   371    112        1 (1093/1379-1196/1436)  -8.1 E04-W03 N22-N24
      3 1792 1077 (N17W53 L= 93)  4303    535       10 (1680/ 937-1914/1232) -55.5 W41-W71 N10-N27
      4 1855 1101 (N20W61 L=101)   842     71        1 (1821/1071-1885/1126)  -6.3 W56-W65 N18-N22
      5  315 1066 (S13E49 L=351)   381    110        0 ( 286/1016- 355/1102)  50.2 E51-E45 S15-S09
      6  909  775 (S18E02 L= 38)   453    172        0 ( 837/ 750- 955/ 820)  18.4 E06-W00 S20-S15
      7  177  753 (S35E70 L=330)  1048     85        2 ( 146/ 713- 203/ 785) -16.7 E77-E65 S37-S32
      8 1351  556 (S23W32 L= 72)   553    160        0 (1279/ 529-1410/ 595)  49.4 W28-W34 S25-S18
      9 1113  484 (S33W19 L= 59)   244     53        0 (1094/ 457-1132/ 511) -26.6 W17-W21 S34-S30
     10 1511  384 (S30W55 L= 95)  2358    376        2 (1352/ 322-1650/ 458)  41.6 W44-W65 S37-S21
  

    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