Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  425 1546 (N30E54 L=144)  1464    150        2 ( 363/1499- 501/1590) -10.8 E65-E44 N27-N32
      2  282  882 (S13E52 L=145)   560     82        0 ( 254/ 843- 307/ 924)  50.2 E55-E49 S15-S09
      3  220  830 (S17E61 L=137)   441     59        0 ( 204/ 803- 235/ 865) -63.4 E62-E58 S18-S14
      4  487  750 (S19E35 L=163)   415     92        0 ( 453/ 720- 516/ 795)  45.0 E37-E32 S21-S16
      5 1226  739 (S14W15 L=213)   632    242        0 (1116/ 703-1318/ 771)  -4.1 W07-W21 S15-S12
      6 1708  644 (S18W54 L=252)  1173    199        3 (1634/ 599-1754/ 677)  26.6 W47-W60 S20-S15
      7  779  362 (S45E15 L=183)   495     69        1 ( 746/ 337- 810/ 386)  39.8 E18-E11 S48-S42
      8  546  341 (S50E42 L=156)   447     68        0 ( 507/ 324- 577/ 363) -15.9 E45-E38 S52-S48
  

    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