Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  881 1844 (N41E39 L=177)  1077    301        1 ( 852/1775- 920/1916) -32.9 E48-E31 N36-N47
      2  782 1536 (N16E28 L=188)   502    137        0 ( 740/1477- 827/1595)  83.7 E29-E26 N11-N20
      3  375 1439 (N02E53 L=163)   424    112        0 ( 364/1381- 396/1484) -71.6 E54-E51 S01-N05
      4  730 1338 (N03E25 L=191)   467    127        0 ( 677/1293- 758/1402)  81.9 E26-E23 N00-N07
      5 1680 1304 (N27W34 L=250)   725    139        0 (1642/1222-1720/1371) -49.4 W30-W38 N23-N31
      6  864 1266 (N02E15 L=201)   904    286        2 ( 709/1231- 996/1287)  13.2 E24-E07 N00-N05
      7 1023 1182 (N01E04 L=212)   287     70        0 ( 997/1137-1041/1216) -33.7 E06-E02 N00-N02
      8  952 1003 (S10E03 L=212)  1431    358        6 ( 791/ 950-1061/1091)  47.5 E12-W03 S16-S03
      9 1914  847 (N10W69 L=285)   642    136        0 (1885/ 814-1933/ 871)  -7.1 W64-W74 N08-N10
     10  639  645 (S39E15 L=200)  1684    438        5 ( 463/ 625- 849/ 676)  20.1 E30-W00 S43-S32
     11 1220  618 (S24W23 L=239)  1093    343        1 (1055/ 590-1385/ 659)  26.6 W13-W33 S29-S18
  

    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