Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  676 1613 (N18E39 L=167)  1021    311        1 ( 637/1469- 703/1753) -38.7 E51-E31 N10-N27
      2  317 1553 (N07E67 L=140)  6001    310        7 ( 229/1428- 394/1681) -76.8 E76-E59 S01-N16
      3  933 1622 (N26E24 L=183)   290    100        0 ( 914/1572- 957/1664) -66.8 E25-E21 N22-N29
      4  274 1134 (S18E53 L=154)  1868    292        3 ( 187/ 999- 350/1234)  50.9 E60-E47 S27-S11
      5  395  922 (S28E38 L=169)  2661    571        3 ( 148/ 775- 627/1015)  26.1 E71-E22 S41-S16
      6 1269  887 (S07W17 L=224)   442    246        0 (1210/ 820-1351/ 933)  45.0 W14-W22 S11-S03
      7  977  705 (S25W07 L=213)   245     79        0 ( 937/ 687-1015/ 720)  45.0 W04-W08 S26-S22
      8  397  595 (S49E36 L=170)   386     50        1 ( 377/ 575- 423/ 624) -26.6 E38-E33 S50-S47
      9 1663  519 (S13W59 L=266)  1039    140        2 (1600/ 490-1734/ 546)  35.0 W54-W65 S16-S09
     10 1701  455 (S15W69 L=275)  3019    185        2 (1629/ 384-1748/ 519) -84.3 W62-W79 S20-S10
     11  782  404 (S50W06 L=212)   306     75        0 ( 762/ 368- 805/ 440) -29.7 W02-W09 S51-S47
     12  973  331 (S48W24 L=231)  1480    255        4 ( 860/ 270-1080/ 378)  10.8 W13-W35 S51-S45
  

    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