Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1475 1736 (N32W62 L=137)  1359    127        2 (1420/1694-1522/1770) -56.3 W59-W67 N29-N35
      2  462 1463 (N48E31 L= 43)   413     71        1 ( 429/1435- 488/1488)  18.4 E35-E28 N45-N49
      3 1828 1150 (S10W57 L=132)  1060    147        3 (1786/1096-1878/1214) -66.0 W53-W62 S13-S04
      4 1761  985 (S17W46 L=121)  1759    267        3 (1736/ 855-1784/1108)  71.6 W42-W49 S24-S09
      5  610 1036 (N18E24 L= 50)   302     81        2 ( 573/1010- 641/1070) -59.0 E26-E22 N15-N20
      6  375  853 (N12E46 L= 29)   825    216        3 ( 334/ 767- 433/ 901) -90.0 E48-E40 N06-N15
      7  176  749 (N09E72 L=  3)  1824    185        1 ( 134/ 692- 231/ 830) -45.0 E78-E64 N04-N13
      8 1720  791 (S27W39 L=114)   314     75        0 (1682/ 767-1748/ 823)   9.5 W35-W42 S27-S25
      9 1559  424 (S48W22 L= 96)   578     90        0 (1526/ 392-1594/ 462)  23.2 W18-W25 S50-S46
     10  918  293 (S35E31 L= 43)   878    161        0 ( 861/ 243- 962/ 366)  26.6 E38-E24 S37-S30
     11 1419  309 (S53W06 L= 80)   789    116        1 (1359/ 275-1470/ 341)  -5.2 W00-W11 S54-S49
  

    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