Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  524 1624 (N58E22 L= 98)   410     68        0 ( 495/1606- 563/1644) -12.5 E26-E17 N56-N60
      2 1133 1601 (N35W27 L=147)  1115    231        2 (1018/1569-1240/1627) -34.7 W19-W32 N30-N39
      3  288 1317 (N42E49 L= 71)   555     76        1 ( 250/1284- 320/1348) -29.1 E55-E45 N39-N44
      4 1715 1153 (S07W46 L=166)   309     49        0 (1690/1123-1731/1175)  71.6 W44-W47 S08-S05
      5 1686  962 (S17W39 L=159)  1307    375        1 (1569/ 843-1794/1093)  14.0 W27-W51 S20-S12
      6 1357 1004 (S04W18 L=138)   439    119        0 (1325/ 949-1392/1060) -82.9 W16-W18 S08-N00
      7  228  766 (N10E62 L= 58)  1438    215        2 ( 167/ 685- 292/ 835) -78.7 E72-E58 N03-N13
      8 1425  531 (S35W11 L=131)  1170    223        2 (1380/ 445-1456/ 627)  69.0 W05-W14 S40-S27
      9  551  324 (S24E57 L= 63)   435     48        0 ( 521/ 309- 582/ 342)  -9.5 E60-E54 S24-S22
  

    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