Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  743 1826 (N62E42 L= 11)   867     99        2 ( 697/1806- 798/1850)  -9.1 E56-E31 N59-N64
      2  553 1615 (N42E44 L=  9)   385     66        0 ( 542/1575- 569/1648) -56.3 E46-E41 N38-N44
      3  414 1471 (N30E50 L=  3)   484     88        1 ( 383/1439- 446/1512) -32.0 E55-E46 N28-N33
      4 1701 1398 (N26W54 L=106)   425     65        1 (1674/1373-1737/1413)   0.0 W50-W58 N24-N27
      5  964 1358 (N24E04 L= 48)   474    196        0 ( 886/1299-1035/1398) -33.7 E09-W00 N20-N27
      6 1743  855 (S08W51 L=104)   432     30        1 (1718/ 831-1767/ 876)  45.0 W49-W53 S09-S06
      7 1633  845 (S08W41 L= 94)   426    115        0 (1586/ 827-1675/ 870)   8.1 W37-W45 S09-S06
      8 1814  832 (S10W60 L=112)   484     64        1 (1795/ 805-1834/ 859)  53.1 W57-W62 S11-S07
      9 1206  733 (S15W12 L= 65)   947    289        2 (1120/ 681-1284/ 806) -15.3 W06-W17 S18-S09
     10  596  688 (S18E29 L= 24)   979    183        1 ( 532/ 627- 677/ 728)  26.6 E34-E22 S22-S15
     11 1212  568 (S26W13 L= 66)   936    374        0 (1064/ 535-1377/ 604)   5.2 W03-W25 S28-S23
     12 1566  544 (S28W42 L= 94)   433     73        0 (1529/ 515-1601/ 574)  31.0 W39-W44 S30-S25
  

    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