Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1710 1351 (N04W52 L=148)   303     96        0 (1688/1309-1731/1394) -90.0 W51-W53 N01-N06
      2  706 1205 (N23E14 L= 81)   520     97        0 ( 653/1184- 764/1223) -23.2 E17-E10 N20-N24
      3 1738  947 (S19W43 L=138)   365     89        0 (1706/ 910-1766/ 987)  18.4 W39-W45 S20-S16
      4  212  919 (N17E59 L= 37)  1710    270        4 ( 110/ 881- 347/ 985) -20.3 E73-E45 N12-N22
      5 1602  905 (S17W31 L=126)   407    127        0 (1562/ 852-1645/ 948)  20.6 W27-W35 S18-S15
      6  555  859 (N07E31 L= 64)   308     91        1 ( 532/ 828- 585/ 900) -90.0 E32-E29 N04-N09
      7 1631  708 (S30W30 L=126)   859    269        0 (1563/ 607-1694/ 826)  39.8 W23-W37 S34-S24
      8  414  770 (N05E44 L= 52)   334    100        0 ( 372/ 743- 465/ 798) -33.7 E46-E40 N02-N06
      9  835  720 (S07E18 L= 77)   311     98        0 ( 788/ 687- 881/ 752)   8.1 E21-E14 S08-S06
     10  405  609 (S04E50 L= 45)   740    152        1 ( 353/ 557- 457/ 675)  22.6 E55-E43 S06-S01
     11 1489  426 (S45W15 L=111)   461     99        1 (1440/ 399-1542/ 445) -37.9 W11-W20 S49-S41
  

    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