Good Mooning.... should be said in a French accent...
Trying out a new DMK51 imager with 1600x1200 resolution. Very nice looking at the Moon through a Ed80MM telescope, used a 2x Barlow for the closer imagers. Just summed 100 frames from a recorded avi and used sharpen once in ImageJ. I tried a few stars but the noise level went through the roof, which is a problem as I want to get the current comets, Linear and Lemmon, may be pons-G if I get up early enough...not likely....8) .
Saturday, January 19, 2013
Saturday, December 8, 2012
STEREO HI2B Comet C/2012 V4
Update - 20121214 - C/2012 V4 has crossed the HI1B imager in images of 20121211. Looks elongated due to is close proximity to the imager and its orbital speed. each HI1B image is a sum of several short duration images.
Currently in STEREO HI2B imager is a comet showing a dynamic streaming tail. This is due to interactions with the solar wind, standby for more interactions if a CME passes by.
This is Comet C/2012 V4, it has been said that this is a return of a lost comet D/1827 M1 (Pons-Gambart). But this has yet to be officially confirmed.
Sunday, December 2, 2012
Fun with Geminids Meteors 2012
The Geminds are coming, you may be lucky to have no clouds and be in dark enough skies to visually count them. Or even capture a few images. But if the clouds come you can still hear then and now also capture images of then using these examples of DIY meteor counters.
Get Processing programing environment from this location. Install it on to your computer. Try some of the examples, it very easy to get started.
The first example is a simple test to see that you can hear and see the audio trace,. You may need to adjust your audio settings to make this work. The programs are trying to listen to the audio line on your PC. Try it without any USB headsets if they are connected. changing the audio mixer settings may help.
Get the text code here for RadMCmin.txt. Copy the text into a new sketch and save. When you have the space weather Radio stream active, run the RADMCmin program , You will see a small window open to display a spectrogram and a signal trace of the signal level. A good meteor will look this..
This one doesn't count or save anything, your challenge would be to make changes to this version to make your own meteor counter.
Yes, here's one I made earlier, this one saves a time stamped screen image of the event and also counts each event. It has a trigger control to set sensitivity, look at the code for all the control keys.
Get the text code here for RadMCpingo.txt ....
You will need to make a small directory tree under the pingo directory, like this...
So when you save a few good images then make a timeline of the events with your favorite image software, like this... I used ImageJ...
Since you now have a UTC time line of events you could then try to correlate these events to other sources of info like this one, Fireball Network in the USA. The Space weather radio signal is being detected in New Mexico. The fireball network has a few cameras there also.
Get Processing programing environment from this location. Install it on to your computer. Try some of the examples, it very easy to get started.
The first example is a simple test to see that you can hear and see the audio trace,. You may need to adjust your audio settings to make this work. The programs are trying to listen to the audio line on your PC. Try it without any USB headsets if they are connected. changing the audio mixer settings may help.
Get the text code here for RadMCmin.txt. Copy the text into a new sketch and save. When you have the space weather Radio stream active, run the RADMCmin program , You will see a small window open to display a spectrogram and a signal trace of the signal level. A good meteor will look this..
This one doesn't count or save anything, your challenge would be to make changes to this version to make your own meteor counter.
Yes, here's one I made earlier, this one saves a time stamped screen image of the event and also counts each event. It has a trigger control to set sensitivity, look at the code for all the control keys.
Get the text code here for RadMCpingo.txt ....
You will need to make a small directory tree under the pingo directory, like this...
So when you save a few good images then make a timeline of the events with your favorite image software, like this... I used ImageJ...
Since you now have a UTC time line of events you could then try to correlate these events to other sources of info like this one, Fireball Network in the USA. The Space weather radio signal is being detected in New Mexico. The fireball network has a few cameras there also.
Labels:
Comet,
Geminids,
meteors,
phaethon,
processing.org,
radio,
Radio Meteors
Friday, November 30, 2012
Computing evolutions
Mark this point in time as a turning point in computer evolution. Look at these links...
http://www.kickstarter.com/projects/adapteva/parallella-a-supercomputer-for-everyone
http://www.adapteva.com/
These Adapteva guys have achieved thumb nail sized parallel processing chips with 16, 64, and soon 1024 CPU on board. Add one to your Mobile or add a few to your standard computer. I want one now, but the OS and generalised software is yet to catch up. Expect full desktop, laptop, and tablets to be using these in the next two years or sooner... check out the cost/performance ...
http://www.kickstarter.com/projects/adapteva/parallella-a-supercomputer-for-everyone
http://www.adapteva.com/
These Adapteva guys have achieved thumb nail sized parallel processing chips with 16, 64, and soon 1024 CPU on board. Add one to your Mobile or add a few to your standard computer. I want one now, but the OS and generalised software is yet to catch up. Expect full desktop, laptop, and tablets to be using these in the next two years or sooner... check out the cost/performance ...
Tuesday, November 13, 2012
20121114 Eclipse from Melbourne
Sunday, November 11, 2012
Comet 168P via Itelescope T5
I am not sure what the little flash is about. It is part of the raw image data, and still visible in calibrated images on the right hand side. I don't think it is part of the comet. Could be a cosmic hit on the sensor or perhaps a micro meteor burst? or a strange internal reflection?
Friday, October 19, 2012
Fun with Radio Meteors
Update -20121116- New Features in my software.
Well, The dynamic level setting, was to complicated and unnecessary. The better way to detect events that helps with event visualisation is more sensitive. Even the small ones can be counted .
This has also simplified the interface controls.
Update -20121114- New Features in my software.
The interface now has dynamic level setting so that only high level events are triggered. Found a better way to detect events, helps with event visualisation . Added a Pingo Meter to the right hand side. ...
Update -20121112- New Features in my software.
The interface now has dynamic level setting so that only high level events are triggered. Also added a colourgram counter. A few more features in mind....
Update -20121025 - Mostly small and a few large meteors
Update 20121023 - If you don't want to do the programing with my little script.
Then this is the software that the big guys are using , ARGO can be found here
http://www.sdradio.eu/weaksignals/argo/index.html
It works with the Spaceweather audio stream just as mine does. Have fun......
UPDate: 20121022 - While watching the Audio trace from Space Weather Radio, this large burst happened by. Not sure what type of burst this was, usually they are spikes.
Original post -
Remember my post "Fun With Geminids 2010". It was a bit complicated to get working.
I have been using Processing.org software lately, lots of fun. There are a few examples of how to show an audio spectra gram in real time.
With this little gem of code, LiveSpectrogram - Takes successive FFTs and renders them onto the screen as gray scale, scrolling left. From Dan Ellis dpwe@ee.columbia.edu 2010-01-15.
I wrapped a simple interface and triggered image saver around it.
I have it running now, but the big meteors just don't want to happen, so here's a few I have captured during testing.
Now, it seems to work only if you have your mic line on or a headset mic on.But, if you now go to http://spaceweatherradio.com/ and
turn on the http://topaz.streamguys.tv/~spaceweather/ audio stream.
It will trigger on the larger Meteors, if you adjust the sensitivity correctly.
Just beware it will create a lot of images if the sensitivity is to low.
And if you make a lot of noise, whistle or speak, it will trigger also.
You will need to be vwary vwary quite when hunting radio meteors....
How to make for your self - Download Processing.org software, load the following text into a new sketch and save. Make a data directory under this sketch folder. Then run it.
The images appear in the sketch folder.
Have fun...... try to catch a few ..... Orionids Oct 21-22, 2012
.......hat tip to http://astroblogger.blogspot.com.au/...
---------------------\/\/\/\/\/\/ cut here \/\/\/\/\/\/ -----------------
/**
*** Radio Meteor Counter **** agw 20121021 ****
*** modified from
*
* LiveSpectrogram
* Takes successive FFTs and renders them onto the screen as grayscale, scrolling left.
*
* Dan Ellis dpwe@ee.columbia.edu 2010-01-15
*/
import ddf.minim.analysis.*;
import ddf.minim.*;
Minim minim;
AudioInput in;
FFT fft;
int colmax = 400;
int rowmax = 256;
int[][] sgram = new int[rowmax][colmax];
int[][] sgramb = new int[rowmax][colmax+1];
int col;
int leftedge;
int count = 0 ;
int sgwt = 0 ;
int sens = 155;
int lvl = 0;
int lvlm = 0;
int trig = 99 ;
String mode = "AUTO";
//
PFont fontA = createFont("Verdana", 12);
//
void setup()
{
size(colmax+ 150, rowmax+125, P3D);
// colorMode(RGB, 255);
background(0);
noStroke();
frameRate(100);
//
fontA = loadFont("CourierNew-12.vlw");
//
textFont(fontA, 12);
textMode(SCREEN);
textAlign(LEFT);
//
//
minim = new Minim(this);
in = minim.getLineIn(Minim.STEREO,2048);
fft = new FFT(in.bufferSize(), in.sampleRate());
fft.window(FFT.HAMMING);
}
void draw()
{
// background(0);
// perform a forward FFT on the samples in the input buffer
fft.forward(in.mix);
for(int i = 0; i < rowmax /* fft.specSize() */; i++)
{
// fill in the new column of spectral values
sgram[i][col] = (int)Math.round(Math.max(0,2*20*Math.log10(1000*fft.getBand(i))));
sgramb[i][col+1] = (int)Math.round(Math.max(0,2*20*Math.log10(1000*fft.getBand(i))));
//
if ( sgramb[i][col] >= sens) {
//
count++;
trig = 10;
//
}
fill(0, sgram[i][col], 255);
rect(460, 381, 5, -((sgram[i][col]*2)));
//
//
}
// next time will be the next column
col = col + 1;
// wrap back to the first column when we get to the end
if (col == colmax) { col = 0; }
// Draw points.
// leftedge is the column in the ring-filled array that is drawn at the extreme left
// start from there, and draw to the end of the array
for (int i = 0; i < colmax-leftedge; i++) {
for (int j = 0; j < rowmax; j++) {
//
// sgwt =(sgram[j][i+leftedge])-(sgramb[j][255]);
//stroke(sgwt);
stroke(sgram[j][i+leftedge]);
//
point(i,height-j);
//rect(i,height-j,2,2);
}
}
// Draw the rest of the image as the beginning of the array (up to leftedge)
for (int i = 0; i < leftedge; i++) {
for (int j = 0; j < rowmax; j++) {
//
stroke(sgram[j][i]);
point(i+colmax-leftedge,height-j);
//rect(i+colmax-leftedge,height-j,2,2);
}
}
// Next time around, we move the left edge over by one, to have the whole thing
// scroll left
leftedge = leftedge + 1;
// Make sure it wraps around
if (leftedge == colmax) { leftedge = 0; }
// Add frequency axis labels
int x = colmax + 5; // to right of spectrogram display
stroke(255);
line(x,0,x,height); // vertical line
// Make text appear centered relative to specified x,y point
textAlign(LEFT,CENTER);
for (float freq = 0.0; freq < in.sampleRate()/2; freq += 500.0) {
int y = height - fft.freqToIndex(freq); // which bin holds this frequency?
line(x,y,x+3,y); // add tick mark
text(Math.round(freq)+" Hz", x+5, y); // add text label
}
// zoom around mouse
//copy((mouseX-25), (mouseY-25), 50, 50, 1, 1, 100, 100);
fill(0, 0, 0);
rect(1, 1, 465, 120);
fill(255, 255, 255);
//
lvlm = mouseY-125;
if(lvlm < 0){ lvlm = 0 ;}
if (lvlm > 255) { lvlm =255;}
lvl = sgramb[lvlm][395] ;
//
text(mouseX+" "+mouseY+" "+lvl, 10,10 );
userinfo() ;
//text( date()+" "+time());
// 395 338 main ffequency location.
if (trig <= 10) {
copy(335, 125, 75, 256, 475, 125, 75, 256);
fill(0, 0, 0);
rect(475, 1, 72, 120);
fill(255, 255, 255);
text(mode+"\nSENS\n"+sens+"\nMtR\n"+count+"\n"+trig, 480,50);
trig-=1 ;
//
if (trig == 1) {
saveimg();
trig = 0 ;
}
}
if (trig <= 0) {
trig = 99 ;
}
}
void keyPressed() {
//
if (key =='x') {
exit();
}
if (key =='c') {
count = 0;
background(0);
text(mode+"\nSENS\n"+sens+"\nMtR\n"+count, 480,50);
}
if (key =='-') {
sens-=1 ;
background(0);
text(mode+"\nSENS\n"+sens+"\nMtR\n"+count, 480,50);
}
if (key =='=') {
sens+=1 ;
background(0);
text(mode+"\nSENS\n"+sens+"\nMtR\n"+count, 480,50);
}
if (key =='m') {
if (mode == "AUTO") { mode = "MANU"; }
else { mode = "AUTO"; }
background(0);
text(mode+"\nSENS\n"+sens+"\nMtR\n"+count, 480,50);
}
}
void saveimg() {
saveFrame("Meteor-####.jpg");
}
void userinfo() {
//
textFont(fontA, 14);
text("RADIO METEOR COUNTER\n"+"[E(x)it][Sensitivity -(-)/+(=)]\n"+"[(C)lear counter][(M)ode AUTO/MANU ]", 50, 50 );
textFont(fontA, 12);
//
}
void stop()
{
// always close Minim audio classes when you finish with them
in.close();
minim.stop();
super.stop();
}
---------------------/\/\/\/\/\ cut here /\/\/\/\/\-----------------
Well, The dynamic level setting, was to complicated and unnecessary. The better way to detect events that helps with event visualisation is more sensitive. Even the small ones can be counted .
This has also simplified the interface controls.
![]() |
| A simpler interface |
![]() |
| A collection of shapes and sizes |
The interface now has dynamic level setting so that only high level events are triggered. Found a better way to detect events, helps with event visualisation . Added a Pingo Meter to the right hand side. ...
Update -20121112- New Features in my software.
The interface now has dynamic level setting so that only high level events are triggered. Also added a colourgram counter. A few more features in mind....
Update -20121025 - Mostly small and a few large meteors
![]() |
| Add caption |
Update 20121023 - If you don't want to do the programing with my little script.
Then this is the software that the big guys are using , ARGO can be found here
http://www.sdradio.eu/weaksignals/argo/index.html
It works with the Spaceweather audio stream just as mine does. Have fun......
UPDate: 20121022 - While watching the Audio trace from Space Weather Radio, this large burst happened by. Not sure what type of burst this was, usually they are spikes.
![]() |
| Large dispersed meteor trail.... |
![]() |
| What Roswell Radio saw......Image rights to Roswell Meteor |
Remember my post "Fun With Geminids 2010". It was a bit complicated to get working.
I have been using Processing.org software lately, lots of fun. There are a few examples of how to show an audio spectra gram in real time.
![]() |
| Me whistling... |
I wrapped a simple interface and triggered image saver around it.
I have it running now, but the big meteors just don't want to happen, so here's a few I have captured during testing.
![]() |
| Seven meteors across 1 hour... |
![]() |
| A medium size Meteor .. |
turn on the http://topaz.streamguys.tv/~spaceweather/ audio stream.
It will trigger on the larger Meteors, if you adjust the sensitivity correctly.
![]() |
| A larger longer Meteor trace.... |
And if you make a lot of noise, whistle or speak, it will trigger also.You will need to be vwary vwary quite when hunting radio meteors....
How to make for your self - Download Processing.org software, load the following text into a new sketch and save. Make a data directory under this sketch folder. Then run it.
The images appear in the sketch folder.
Have fun...... try to catch a few ..... Orionids Oct 21-22, 2012
.......hat tip to http://astroblogger.blogspot.com.au/...
---------------------\/\/\/\/\/\/ cut here \/\/\/\/\/\/ -----------------
/**
*** Radio Meteor Counter **** agw 20121021 ****
*** modified from
*
* LiveSpectrogram
* Takes successive FFTs and renders them onto the screen as grayscale, scrolling left.
*
* Dan Ellis dpwe@ee.columbia.edu 2010-01-15
*/
import ddf.minim.analysis.*;
import ddf.minim.*;
Minim minim;
AudioInput in;
FFT fft;
int colmax = 400;
int rowmax = 256;
int[][] sgram = new int[rowmax][colmax];
int[][] sgramb = new int[rowmax][colmax+1];
int col;
int leftedge;
int count = 0 ;
int sgwt = 0 ;
int sens = 155;
int lvl = 0;
int lvlm = 0;
int trig = 99 ;
String mode = "AUTO";
//
PFont fontA = createFont("Verdana", 12);
//
void setup()
{
size(colmax+ 150, rowmax+125, P3D);
// colorMode(RGB, 255);
background(0);
noStroke();
frameRate(100);
//
fontA = loadFont("CourierNew-12.vlw");
//
textFont(fontA, 12);
textMode(SCREEN);
textAlign(LEFT);
//
//
minim = new Minim(this);
in = minim.getLineIn(Minim.STEREO,2048);
fft = new FFT(in.bufferSize(), in.sampleRate());
fft.window(FFT.HAMMING);
}
void draw()
{
// background(0);
// perform a forward FFT on the samples in the input buffer
fft.forward(in.mix);
for(int i = 0; i < rowmax /* fft.specSize() */; i++)
{
// fill in the new column of spectral values
sgram[i][col] = (int)Math.round(Math.max(0,2*20*Math.log10(1000*fft.getBand(i))));
sgramb[i][col+1] = (int)Math.round(Math.max(0,2*20*Math.log10(1000*fft.getBand(i))));
//
if ( sgramb[i][col] >= sens) {
//
count++;
trig = 10;
//
}
fill(0, sgram[i][col], 255);
rect(460, 381, 5, -((sgram[i][col]*2)));
//
//
}
// next time will be the next column
col = col + 1;
// wrap back to the first column when we get to the end
if (col == colmax) { col = 0; }
// Draw points.
// leftedge is the column in the ring-filled array that is drawn at the extreme left
// start from there, and draw to the end of the array
for (int i = 0; i < colmax-leftedge; i++) {
for (int j = 0; j < rowmax; j++) {
//
// sgwt =(sgram[j][i+leftedge])-(sgramb[j][255]);
//stroke(sgwt);
stroke(sgram[j][i+leftedge]);
//
point(i,height-j);
//rect(i,height-j,2,2);
}
}
// Draw the rest of the image as the beginning of the array (up to leftedge)
for (int i = 0; i < leftedge; i++) {
for (int j = 0; j < rowmax; j++) {
//
stroke(sgram[j][i]);
point(i+colmax-leftedge,height-j);
//rect(i+colmax-leftedge,height-j,2,2);
}
}
// Next time around, we move the left edge over by one, to have the whole thing
// scroll left
leftedge = leftedge + 1;
// Make sure it wraps around
if (leftedge == colmax) { leftedge = 0; }
// Add frequency axis labels
int x = colmax + 5; // to right of spectrogram display
stroke(255);
line(x,0,x,height); // vertical line
// Make text appear centered relative to specified x,y point
textAlign(LEFT,CENTER);
for (float freq = 0.0; freq < in.sampleRate()/2; freq += 500.0) {
int y = height - fft.freqToIndex(freq); // which bin holds this frequency?
line(x,y,x+3,y); // add tick mark
text(Math.round(freq)+" Hz", x+5, y); // add text label
}
// zoom around mouse
//copy((mouseX-25), (mouseY-25), 50, 50, 1, 1, 100, 100);
fill(0, 0, 0);
rect(1, 1, 465, 120);
fill(255, 255, 255);
//
lvlm = mouseY-125;
if(lvlm < 0){ lvlm = 0 ;}
if (lvlm > 255) { lvlm =255;}
lvl = sgramb[lvlm][395] ;
//
text(mouseX+" "+mouseY+" "+lvl, 10,10 );
userinfo() ;
//text( date()+" "+time());
// 395 338 main ffequency location.
if (trig <= 10) {
copy(335, 125, 75, 256, 475, 125, 75, 256);
fill(0, 0, 0);
rect(475, 1, 72, 120);
fill(255, 255, 255);
text(mode+"\nSENS\n"+sens+"\nMtR\n"+count+"\n"+trig, 480,50);
trig-=1 ;
//
if (trig == 1) {
saveimg();
trig = 0 ;
}
}
if (trig <= 0) {
trig = 99 ;
}
}
void keyPressed() {
//
if (key =='x') {
exit();
}
if (key =='c') {
count = 0;
background(0);
text(mode+"\nSENS\n"+sens+"\nMtR\n"+count, 480,50);
}
if (key =='-') {
sens-=1 ;
background(0);
text(mode+"\nSENS\n"+sens+"\nMtR\n"+count, 480,50);
}
if (key =='=') {
sens+=1 ;
background(0);
text(mode+"\nSENS\n"+sens+"\nMtR\n"+count, 480,50);
}
if (key =='m') {
if (mode == "AUTO") { mode = "MANU"; }
else { mode = "AUTO"; }
background(0);
text(mode+"\nSENS\n"+sens+"\nMtR\n"+count, 480,50);
}
}
void saveimg() {
saveFrame("Meteor-####.jpg");
}
void userinfo() {
//
textFont(fontA, 14);
text("RADIO METEOR COUNTER\n"+"[E(x)it][Sensitivity -(-)/+(=)]\n"+"[(C)lear counter][(M)ode AUTO/MANU ]", 50, 50 );
textFont(fontA, 12);
//
}
void stop()
{
// always close Minim audio classes when you finish with them
in.close();
minim.stop();
super.stop();
}
---------------------/\/\/\/\/\ cut here /\/\/\/\/\-----------------
Subscribe to:
Posts (Atom)





















