Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  217 1203 (S15E60 L=177)   901    141        1 ( 177/1155- 260/1254)  66.8 E63-E56 S18-S11
      2 1518 1187 (N17W24 L=261)   226     60        0 (1497/1146-1533/1220) -45.0 W22-W25 N14-N18
      3 1399 1078 (N07W19 L=256)   529    163        0 (1335/1019-1470/1144)  -9.5 W13-W25 N05-N09
      4 1759  946 (N12W48 L=285)  2232    353        3 (1717/ 778-1782/1108) -66.0 W41-W53 N02-N20
      5  377 1053 (S21E41 L=195)   572    176        3 ( 308/1025- 436/1078)  21.8 E47-E37 S22-S18
      6  753 1018 (S14E15 L=221)   381    100        0 ( 704/1001- 809/1035)  26.6 E17-E11 S15-S12
      7  235  782 (S41E52 L=185)   976    320        3 ( 167/ 745- 319/ 822)  23.6 E59-E43 S44-S36
      8  689  747 (S31E12 L=224)   201     68        0 ( 654/ 724- 723/ 778)  -9.5 E15-E09 S31-S29
      9 1169  695 (S20W18 L=255)   433    130        0 (1105/ 672-1230/ 725)   6.3 W13-W22 S20-S18
     10  638  531 (S47E11 L=226)   625    118        1 ( 586/ 502- 683/ 573)  40.6 E15-E07 S49-S42
  

    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