Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1581 1527 (N27W50 L= 84)   596     79        0 (1543/1481-1612/1568) -81.9 W48-W51 N23-N30
      2  965 1414 (N31W02 L= 36)   349    136        0 ( 901/1399-1031/1440)  -6.3 E02-W07 N30-N32
      3 1807 1386 (N13W66 L=100)   702    174        0 (1772/1364-1837/1417)   5.7 W61-W72 N12-N14
      4  556 1314 (N30E28 L=  6)  1742    379        4 ( 404/1267- 712/1371)   2.3 E41-E16 N27-N32
      5 1138 1180 (N14W09 L= 43)   658    149        4 (1074/1146-1211/1220) -32.0 W05-W13 N10-N16
      6  528  709 (S07E37 L=357)  1965    388        2 ( 471/ 575- 570/ 871)  71.6 E44-E33 S15-N03
      7  773  761 (S06E19 L= 15)   650    147        2 ( 698/ 736- 851/ 786) -18.4 E23-E14 S08-S04
      8 1625  706 (S24W36 L= 70)   779    114        2 (1573/ 659-1674/ 740)  23.2 W31-W40 S26-S21
      9  716  475 (S24E30 L=  4)   461    109        0 ( 683/ 426- 747/ 513)  35.5 E33-E26 S27-S21
     10  559  303 (S36E51 L=343)  1214    180        4 ( 510/ 252- 599/ 372)  70.0 E55-E49 S41-S30
  

    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