Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  930 1377 (N28E04 L=157)   314     85        1 ( 895/1352- 973/1393) -21.8 E05-W00 N26-N29
      2 1552 1295 (N17W38 L=199)   337    125        1 (1497/1279-1609/1315) -18.4 W34-W43 N15-N18
      3  732 1034 (N08E18 L=142)   370     77        1 ( 700/1000- 769/1066) -21.8 E20-E15 N05-N09
      4 1909 1021 (S05W71 L=231)   975     72        1 (1873/1002-1947/1043) -10.6 W65-W81 S06-S03
      5 1807  835 (S15W57 L=217)   585    124        1 (1761/ 812-1846/ 856) -20.0 W51-W62 S17-S13
      6  730  705 (S13E21 L=139)   420    152        1 ( 668/ 681- 776/ 745)   7.1 E25-E17 S14-S10
      7 1419  662 (S21W23 L=184)  1282    303        6 (1305/ 559-1527/ 745)  26.6 W15-W31 S27-S16
      8 1757  632 (S28W56 L=216)   508     62        0 (1734/ 601-1781/ 661)  53.1 W53-W58 S30-S26
      9 1239  467 (S33W11 L=171)   453     54        1 (1209/ 439-1269/ 500)  36.9 W08-W13 S35-S30
     10 1080  271 (S49E03 L=157)   513     87        1 (1053/ 251-1118/ 304) -45.0 E05-W00 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