Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1790 1437 (N45W48 L=139)  1116    125        0 (1731/1403-1847/1459)  11.9 W39-W59 N41-N47
      2  547 1194 (S09E33 L= 58)   261     58        1 ( 515/1162- 569/1220)  76.0 E34-E30 S11-S07
      3 1744 1074 (N18W44 L=135)   608    151        1 (1702/1006-1796/1128) -15.3 W39-W50 N16-N21
      4  589  998 (S19E25 L= 66)   453    153        0 ( 532/ 952- 642/1034)  58.0 E28-E22 S23-S15
      5  429  995 (S24E37 L= 54)   432    155        0 ( 359/ 961- 490/1014)  10.3 E42-E31 S25-S21
      6  299  758 (S41E48 L= 43)  1125    182        1 ( 257/ 698- 354/ 837)  40.6 E55-E41 S46-S34
      7 1307  528 (S25W34 L=125)  1230    201        2 (1238/ 468-1403/ 571)  55.0 W27-W39 S30-S20
      8 1149  476 (S33W27 L=118)   272    102        1 (1125/ 431-1183/ 523) -63.4 W25-W29 S35-S29
      9 1013  337 (S47W28 L=119)  1314    211        1 ( 962/ 283-1087/ 436) -21.3 W19-W37 S50-S40
     10 1184  382 (S37W35 L=126)   522     82        1 (1147/ 354-1223/ 413)  15.9 W30-W38 S38-S35
  

    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