Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  759 1354 (N12E23 L=103)   347     74        2 ( 728/1324- 800/1385) -11.3 E25-E20 N10-N13
      2  666 1132 (S03E24 L=103)   686    190        2 ( 604/1047- 737/1225) -32.5 E29-E17 S06-N01
      3 1863 1117 (N22W59 L=185)  1411    150        2 (1820/1071-1894/1169) -39.8 W54-W62 N19-N24
      4  237 1102 (S12E57 L= 70)   636    149        0 ( 221/1030- 251/1164) -63.4 E59-E55 S16-S08
      5  891 1047 (S03E08 L=118)   297     86        1 ( 880/ 999- 911/1096) -68.2 E09-E07 S05-N00
      6  607  881 (S19E23 L=103)  1140    395        0 ( 442/ 778- 739/1006)  46.5 E34-E16 S28-S08
      7 1358  894 (S03W23 L=149)   213     62        0 (1345/ 864-1372/ 933) -63.4 W21-W23 S04-N00
      8 1488  698 (S11W36 L=163)   279     60        1 (1434/ 671-1525/ 727)  51.3 W33-W38 S13-S08
      9 1697  539 (S16W60 L=187)   722     93        2 (1663/ 508-1733/ 569)  45.0 W57-W63 S17-S13
     10 1192  377 (S38W30 L=156)  1068    260        0 (1053/ 358-1309/ 396)  22.4 W20-W37 S42-S34
  

    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