Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1050 1759 (N48W21 L=109)   328     49        0 (1023/1737-1071/1786)  23.2 W17-W24 N46-N49
      2  747 1452 (N32E10 L= 77)   335    114        0 ( 697/1415- 796/1474)  14.0 E14-E06 N30-N32
      3  551 1405 (N32E25 L= 62)   341     87        0 ( 508/1372- 589/1435)  18.4 E28-E22 N31-N33
      4 1713 1056 (S09W45 L=132)   407     80        1 (1696/1008-1726/1098)  68.2 W43-W46 S11-S06
      5  228  988 (N12E56 L= 32)   427     60        1 ( 203/ 961- 251/1024) -76.0 E58-E53 N10-N14
      6 1297  899 (S11W14 L=101)   856    281        1 (1223/ 776-1357/1012)  42.5 W07-W19 S16-S04
      7  457  910 (N04E37 L= 50)   922    227        0 ( 359/ 872- 531/ 973)   4.1 E44-E30 N02-N06
      8  703  829 (S05E22 L= 65)   265     50        2 ( 681/ 810- 735/ 857) -45.0 E23-E20 S05-S02
      9  695  669 (S14E27 L= 61)   441     82        2 ( 668/ 641- 735/ 708) -90.0 E28-E24 S15-S11
  

    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