Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  639 1719 (N29E48 L=139)  1176    185        2 ( 571/1671- 695/1746)  16.7 E53-E43 N25-N30
      2  364 1582 (N14E64 L=123)  2827    325        6 ( 305/1503- 431/1658)  70.0 E69-E58 N08-N19
      3  479 1580 (N16E52 L=135)   426    111        2 ( 443/1551- 504/1605)  26.6 E55-E50 N14-N18
      4 1020  829 (S15W05 L=192)   266     81        2 ( 990/ 797-1053/ 853)  63.4 W03-W06 S16-S12
      5 1236  703 (S16W21 L=208)   640    222        1 (1136/ 629-1321/ 793)  65.2 W16-W23 S22-S09
      6 1614  730 (S05W45 L=232)   705    140        2 (1588/ 654-1637/ 789) -69.4 W41-W46 S09-S01
      7  372  626 (S44E39 L=148)   391    117        1 ( 326/ 596- 418/ 657) -12.5 E43-E34 S45-S42
      8  734  563 (S38E07 L=181)  1121    330        1 ( 583/ 536- 898/ 588)  16.9 E18-W05 S42-S33
      9 1442  521 (S21W41 L=228)   441    134        1 (1372/ 494-1496/ 558)  49.4 W37-W43 S24-S17
     10  318  532 (S53E49 L=139)   641    235        1 ( 263/ 513- 369/ 554)  14.9 E56-E41 S55-S50
     11 1460  405 (S27W49 L=236)  2836    806        4 (1171/ 336-1658/ 447)  12.8 W27-W71 S35-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