Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  190 1165 (S10E62 L=263)   850    184        0 ( 134/1117- 234/1240)  45.0 E68-E58 S13-S05
      2 1760 1057 (N11W47 L= 12)   806    233        1 (1702/ 948-1786/1138) -48.0 W40-W49 N04-N14
      3  494 1076 (S12E34 L=291)   217     52        2 ( 466/1057- 524/1091)   0.0 E35-E31 S13-S11
      4  425 1030 (S16E38 L=287)   852    291        0 ( 312/ 995- 536/1078)   0.0 E47-E30 S17-S14
      5  223  942 (S23E55 L=270)  1429    365        0 ( 149/ 837- 302/1059)  47.0 E62-E48 S30-S15
      6 1891  896 (N06W65 L= 30)  1908    368        3 (1862/ 826-1939/ 978) -66.0 W58-W74 N00-N09
      7  812  742 (S27E08 L=317)   243     81        1 ( 778/ 708- 840/ 782)  59.0 E09-E06 S29-S24
      8  351  642 (S41E44 L=281)   944    278        4 ( 301/ 568- 404/ 723)  42.5 E50-E38 S46-S34
      9 1482  654 (S18W37 L=  2)   391    127        0 (1415/ 637-1536/ 690)  35.5 W32-W39 S20-S15
     10 1786  623 (S12W64 L= 28)  1675    374        3 (1709/ 564-1848/ 671)  -3.2 W53-W71 S14-S08
     11 1553  408 (S30W53 L= 18)  2945    486        6 (1419/ 370-1729/ 456)  28.6 W42-W67 S35-S22
  

    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