Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1367 1762 (N57W04 L=216)   783    101        2 (1336/1723-1393/1799)  73.3 W00-W06 N52-N62
      2 1129 1649 (N38E11 L=201)   803    203        1 (1081/1555-1179/1719) -37.6 E17-E04 N32-N42
      3 1664 1487 (N44W38 L=249)   537     52        0 (1647/1446-1686/1515) -78.7 W35-W40 N40-N46
      4 1466  784 (S06W32 L=244)   276     57        1 (1445/ 755-1483/ 817) -56.3 W31-W33 S07-S03
      5  307  730 (S38E48 L=164)   677    112        1 ( 271/ 696- 347/ 781)  35.0 E53-E42 S40-S33
      6 1361  703 (S13W28 L=240)   311     67        1 (1341/ 672-1380/ 738) -78.7 W26-W30 S15-S10
      7  365  565 (S48E44 L=168)   870     64        1 ( 335/ 528- 399/ 607)  35.0 E50-E39 S51-S44
      8  408  502 (S53E41 L=171)   464     70        1 ( 374/ 476- 447/ 523)  20.0 E45-E34 S55-S49
      9 1442  385 (S29W50 L=262)  2585    405        3 (1245/ 328-1598/ 450)  36.0 W38-W60 S38-S21
     10  499  409 (S58E30 L=182)   421     42        2 ( 473/ 386- 522/ 429)   0.0 E34-E26 S60-S55
  

    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