Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  957 1561 (N25E19 L=270)  2110    486        3 ( 850/1343-1029/1715) -50.7 E30-E11 N12-N34
      2  531 1433 (N07E41 L=248)   255     42        1 ( 513/1412- 545/1460)  90.0 E42-E39 N05-N09
      3 1804 1185 (N27W50 L=339)  3660    569        3 (1702/1045-1929/1347) -41.2 W38-W69 N19-N34
      4  329 1032 (S20E46 L=244)   577    135        0 ( 307/ 983- 371/1089) -63.4 E47-E41 S22-S16
      5 1507  969 (N05W29 L=319)   789    198        1 (1442/ 903-1545/1067)  84.8 W26-W30 N00-N11
      6 1903  917 (N14W67 L=356)  1657    128        3 (1867/ 867-1952/ 957)  50.2 W61-W79 N11-N17
      7 1854  709 (N02W69 L=358)  1130    126        0 (1798/ 683-1906/ 744)   5.2 W60-W82 N00-N03
      8  957  634 (S28W07 L=296)   558    157        0 ( 891/ 613-1041/ 654)  32.0 W03-W12 S30-S25
      9  397  559 (S49E38 L=251)  2031    275        2 ( 300/ 486- 516/ 636)  -6.8 E50-E24 S53-S45
     10 1611  534 (S15W53 L=343)  1258    284        1 (1485/ 501-1719/ 566)  30.5 W45-W63 S20-S10
     11 1175  253 (S46W43 L=332)   329     43        0 (1157/ 229-1191/ 276)  -9.5 W39-W45 S47-S43
     12 1339  204 (S43W64 L=353)  2053    150        2 (1277/ 165-1405/ 267)  -1.8 W50-W82 S44-S39
  

    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