Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  418 1643 (N23E61 L= 90)  1065    233        2 ( 396/1587- 438/1697) -37.9 E67-E57 N19-N26
      2  330 1283 (S01E52 L= 99)   300     69        1 ( 315/1251- 352/1314) -36.9 E54-E49 S02-N01
      3 1586 1199 (N20W32 L=183)   658    128        1 (1546/1137-1621/1258) -50.2 W28-W35 N16-N23
      4 1832 1120 (N22W55 L=206)   423     58        1 (1807/1100-1855/1143)  21.8 W52-W57 N20-N23
      5  289  931 (S22E48 L=102)   374    102        2 ( 256/ 899- 317/ 964)  38.7 E51-E45 S25-S20
      6  625  900 (S18E22 L=129)   290    110        2 ( 589/ 873- 675/ 932)   0.0 E24-E18 S19-S15
      7  368  871 (S25E41 L=110)  1089    334        1 ( 295/ 792- 485/ 919)  36.3 E47-E31 S31-S20
      8 1143  828 (S11W11 L=162)   280     91        0 (1093/ 812-1191/ 846)  31.0 W08-W13 S12-S08
      9  916  711 (S23W00 L=151)   470    145        0 ( 853/ 679- 966/ 772)  63.4 E03-W02 S25-S17
     10  963  479 (S37W10 L=161)   466    133        0 ( 893/ 465-1034/ 499)  21.8 W04-W14 S38-S33
     11 1427  443 (S27W43 L=193)   991    203        0 (1336/ 412-1490/ 499)  45.0 W36-W46 S31-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