Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  831 1812 (N36E42 L=166)   735    203        1 ( 800/1752- 858/1860) -32.5 E47-E36 N32-N40
      2  779 1230 (S01E19 L=189)   207     98        0 ( 752/1198- 819/1272) -11.3 E21-E16 S02-N00
      3  519 1085 (S16E32 L=177)   846    196        1 ( 438/1012- 578/1154)  61.4 E35-E28 S21-S10
      4 1820 1020 (N18W51 L=260)   597    104        0 (1794/ 972-1841/1070) -39.8 W47-W54 N15-N20
      5  972  981 (S10E02 L=207)   245     88        0 ( 938/ 957-1013/1004)  45.0 E03-W01 S12-S08
      6  653  923 (S22E19 L=189)   351    134        0 ( 620/ 892- 716/ 977)  -7.1 E22-E14 S23-S19
      7 1739  883 (N08W48 L=257)  1251    234        0 (1690/ 793-1800/ 964) -18.4 W41-W57 N04-N10
      8 1253  859 (S09W17 L=226)   604    132        2 (1228/ 790-1270/ 933) -63.4 W15-W19 S12-S04
      9 1185  802 (S14W15 L=224)   242    100        1 (1124/ 771-1222/ 824)   0.0 W11-W18 S14-S12
     10  438  692 (S42E31 L=177)   989    307        0 ( 334/ 669- 556/ 713)  21.8 E41-E21 S45-S37
     11 1258  535 (S26W30 L=238)   642    217        0 (1163/ 504-1383/ 562)  33.7 W24-W37 S30-S21
     12 1616  479 (S17W57 L=266)   627    189        1 (1555/ 449-1667/ 503)  14.0 W51-W63 S18-S15
     13 1051  343 (S45W29 L=237)   295     68        0 (1016/ 331-1093/ 355)  45.0 W26-W31 S46-S42
  

    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