Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1289 1764 (N55W32 L=168)  1731    334        1 (1153/1723-1426/1802)   6.6 W14-W57 N52-N58
      2  419 1452 (N31E48 L= 87)  2059    367        2 ( 327/1328- 545/1583) -24.0 E69-E33 N23-N39
      3 1151 1555 (N38W11 L=146)   482     71        0 (1116/1538-1202/1576) -18.4 W08-W14 N37-N40
      4 1349 1505 (N34W26 L=161)   345     67        0 (1313/1485-1387/1522)  15.9 W22-W29 N32-N35
      5  897 1254 (N18E08 L=128)   902    284        1 ( 824/1172- 998/1305)  24.4 E12-E01 N13-N21
      6  161 1222 (N15E73 L= 63)  1733    145        4 ( 127/1169- 206/1267)  31.0 E81-E66 N11-N17
      7 1002 1225 (N16E01 L=134)   227     39        1 ( 974/1206-1027/1242)   0.0 E02-W00 N15-N17
      8  573 1025 (N04E29 L=106)  1072    239        0 ( 499/ 939- 667/1097)  33.7 E34-E22 S01-N08
      9  503  861 (S06E35 L=101)   416     65        1 ( 469/ 836- 535/ 885)  26.6 E37-E31 S07-S04
     10  453  642 (S21E42 L= 94)  1723    281        2 ( 397/ 508- 514/ 739)  52.6 E51-E35 S30-S13
     11 1208  488 (S32W12 L=148)  4399    720       10 ( 948/ 331-1437/ 692)  36.3 E07-W27 S44-S17
     12  858  518 (S29E13 L=123)  1922    550        2 ( 685/ 382-1068/ 688)  35.0 E28-W02 S39-S17
  

    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