Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  803 1652 (N52E05 L=248)   431     56        1 ( 779/1617- 826/1678)  45.0 E07-E02 N48-N54
      2 1793 1349 (N10W64 L=317)  1122    180        2 (1733/1296-1851/1424) -48.4 W60-W70 N06-N15
      3  250 1080 (N19E56 L=197)  1273    213        3 ( 185/1044- 337/1126)   6.7 E65-E47 N17-N21
      4 1282  924 (S04W14 L=267)   375    126        0 (1232/ 903-1350/ 940)   0.0 W10-W18 S04-S02
      5  324  839 (N04E51 L=202)   369     51        0 ( 300/ 812- 354/ 864)  18.4 E53-E47 N02-N05
      6 1387  539 (S31W16 L=269)  2854    599        3 (1163/ 399-1567/ 736)  25.8 E01-W31 S37-S20
      7  315  614 (S10E59 L=194)   454     44        1 ( 291/ 594- 342/ 635)  -9.5 E62-E56 S11-S08
      8  586  535 (S18E38 L=215)   435     79        0 ( 567/ 493- 610/ 575)  51.3 E40-E35 S20-S15
      9  338  461 (S21E68 L=185)   722     69        2 ( 310/ 435- 368/ 484)  10.0 E79-E62 S23-S19
     10  823  380 (S32E27 L=226)   337     87        0 ( 795/ 356- 859/ 416) -59.0 E28-E24 S34-S28
     11  580  348 (S31E49 L=204)   967    152        2 ( 553/ 285- 627/ 414)  32.0 E57-E40 S36-S26
  

    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