Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1369 1764 (N36W55 L=148)  1515    158        3 (1299/1722-1445/1793) -14.0 W47-W65 N34-N38
      2  683 1716 (N60E01 L= 92)   440     85        1 ( 640/1693- 713/1738)   6.3 E06-W03 N57-N61
      3 1775  863 (S26W45 L=138)  4372    758       12 (1690/ 628-1943/1156)  31.2 W35-W74 S39-S15
      4  238 1103 (N29E54 L= 38)  1668    212        5 ( 142/1036- 325/1140) -18.9 E68-E45 N23-N32
      5  344  817 (N10E49 L= 44)   670    159        1 ( 319/ 757- 377/ 884)  66.8 E50-E44 N06-N13
      6  587  809 (N04E31 L= 62)   254    107        0 ( 543/ 793- 633/ 843) -51.3 E32-E28 N02-N07
      7  222  770 (N10E63 L= 30)  1609    153        3 ( 156/ 702- 261/ 829)  66.8 E72-E58 N05-N12
      8 1468  460 (S42W13 L=105)   606    155        3 (1405/ 419-1536/ 498)   0.0 W07-W19 S43-S40
      9  312  472 (S10E70 L= 22)   671     62        0 ( 281/ 456- 332/ 496)   4.8 E77-E65 S10-S08
     10  688  206 (S36E59 L= 34)  1224    110        0 ( 668/ 145- 711/ 255)  19.1 E76-E50 S41-S32
  

    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