Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  320 1466 (N06E59 L=134)  1994    366        9 ( 249/1372- 386/1572) -72.9 E66-E53 N00-N13
      2  857 1378 (N10E17 L=175)   406    136        0 ( 835/1309- 880/1434) -60.3 E19-E14 N06-N13
      3  111 1155 (S14E72 L=120)  2832    253        6 (  84/1087- 141/1220) -76.0 E77-E68 S17-S09
      4  217 1163 (S13E59 L=134)   502    117        1 ( 194/1114- 237/1211)  90.0 E61-E56 S16-S10
      5  320  941 (S25E45 L=147)  2340    483        4 ( 213/ 772- 432/1102)  56.9 E55-E37 S36-S13
      6  800 1053 (S10E13 L=179)   233     72        2 ( 771/1017- 818/1094)  78.7 E14-E12 S12-S07
      7 1262  901 (S09W16 L=209)   359     70        2 (1224/ 878-1304/ 927)   9.5 W13-W19 S10-S07
      8 1258  745 (S17W20 L=213)  1074    281        5 (1118/ 656-1353/ 865)  59.0 W13-W24 S25-S08
      9  900  748 (S26E01 L=191)   517     97        1 ( 843/ 707- 945/ 788)  66.8 E04-W00 S29-S22
     10 1569  481 (S24W52 L=244)  3065    749        2 (1265/ 368-1846/ 583)  36.9 W36-W76 S38-S08
     11  541  548 (S47E25 L=167)   453     65        1 ( 503/ 521- 576/ 570) -15.9 E28-E21 S48-S44
     12 1300  519 (S29W31 L=224)   525    119        1 (1230/ 502-1356/ 543)  20.6 W26-W34 S31-S26
  

    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