Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1119 1314 (N18W00 L=282)  3112    755        2 ( 851/1018-1344/1584)  15.8 E23-W19 N04-N30
      2 1712 1564 (N50W49 L=331)   579    120        0 (1675/1550-1740/1581)   7.1 W44-W53 N48-N51
      3 1450 1553 (N41W20 L=302)   288     71        0 (1425/1536-1481/1569)   0.0 W18-W23 N39-N42
      4 1868 1092 (N19W60 L=342)  1593    315        3 (1800/1008-1926/1166) -21.3 W52-W70 N15-N22
      5  368 1031 (S13E43 L=239)   338     91        3 ( 356/ 988- 387/1086) -80.5 E44-E41 S15-S09
      6  147 1005 (S18E66 L=216)   703    153        1 ( 126/ 943- 166/1062)  69.4 E68-E64 S22-S14
      7 1380  935 (N00W23 L=305)   393    121        1 (1336/ 893-1433/ 986) -20.6 W19-W27 S01-N02
      8  416  828 (S24E37 L=245)   593    174        0 ( 348/ 793- 483/ 860)   0.0 E42-E31 S25-S22
      9  958  789 (S16W00 L=282)  1740    501        5 ( 729/ 743-1185/ 848)  25.6 E13-W13 S21-S08
     10 1396  370 (S34W44 L=325)  5884    759       14 (1082/ 324-1648/ 450)  25.4 W21-W61 S43-S23
     11 1018  367 (S43W16 L=298)   453     94        0 ( 967/ 349-1058/ 386)  25.5 W11-W19 S45-S41
  

    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