Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  968 1870 (N46E38 L=  8)  1091    119        2 ( 929/1840-1013/1901)  53.1 E41-E33 N41-N49
      2 1072 1688 (N34E18 L= 27)  2269    447        4 (1011/1543-1117/1840) -42.1 E31-E09 N24-N44
      3  383 1653 (N15E66 L=339)   976    159        3 ( 357/1624- 415/1691)  68.2 E69-E62 N12-N17
      4  354 1491 (N04E58 L=347)   499    127        2 ( 314/1459- 405/1522)  26.6 E61-E53 N01-N06
      5  698 1400 (N06E29 L= 17)   930    225        0 ( 627/1342- 796/1480) -13.0 E35-E21 N03-N09
      6  926 1319 (N07E13 L= 32)   225     49        0 ( 909/1292- 947/1345) -71.6 E14-E11 N05-N08
      7 1059 1288 (N08E05 L= 40)   524    170        2 ( 981/1217-1115/1334) -18.4 E10-W00 N05-N10
      8  502 1313 (S04E39 L=  7)   368     96        0 ( 456/1288- 541/1333)   0.0 E42-E35 S04-S02
      9 1014 1104 (S03E03 L= 43)  1081    273        1 ( 877/ 992-1125/1184)  63.4 E07-W01 S12-N04
     10  553 1006 (S20E28 L= 18)   252     49        0 ( 521/ 987- 584/1026)   0.0 E30-E25 S20-S18
     11  525  916 (S26E28 L= 18)   378     52        1 ( 505/ 883- 550/ 942)  90.0 E29-E26 S27-S23
  

    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