Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1668 1516 (N21W60 L=137)   565     52        1 (1638/1494-1698/1543)   6.3 W55-W64 N20-N22
      2 1309 1210 (N12W21 L= 98)   290     88        0 (1272/1192-1356/1229) -11.3 W18-W23 N11-N13
      3  218 1158 (N26E60 L= 17)   624    111        3 ( 183/1118- 260/1197) -29.1 E66-E54 N23-N28
      4 1801  938 (S15W52 L=129)   820     94        2 (1779/ 898-1832/ 982)  26.6 W49-W55 S17-S12
      5 1918  926 (S20W69 L=146)   731     71        0 (1902/ 891-1937/ 963)  26.6 W66-W74 S21-S18
      6 1462  917 (S08W25 L=102)   280     61        0 (1422/ 891-1488/ 938) -26.6 W22-W26 S09-S05
      7  446  577 (S12E48 L= 29)  2009    320        3 ( 382/ 481- 554/ 711)  56.3 E56-E40 S17-S02
      8  914  568 (S19E15 L= 62)   500    156        0 ( 894/ 489- 955/ 637)  84.3 E16-E13 S24-S14
      9 1608  583 (S33W33 L=110)   916    103        3 (1552/ 559-1656/ 614)  -6.3 W28-W38 S36-S31
     10  320  518 (S15E65 L= 12)  3346    225        7 ( 272/ 437- 361/ 603)  47.5 E76-E55 S20-S08
     11 1192  513 (S28W01 L= 78)   342     79        1 (1164/ 478-1213/ 554)  68.2 E01-W02 S30-S25
     12  638  487 (S20E36 L= 41)   542    129        2 ( 617/ 420- 663/ 542)  48.8 E41-E32 S24-S16
  

    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