Krioma.net Blogtag:www.krioma.net,2012:/blog//12012-02-04T14:38:27ZA site based on the thoughts of Samuel George - the guy behind Krioma.net. The thoughts mostly stem from his life experiences which mostly include astrophysics, astronomy, beer, football, music and computers.Movable Type 3.38FAS Newsletter #1tag:www.krioma.net,2012:/blog//1.17032012-02-04T14:21:24Z2012-02-04T14:38:27ZI feel like someone having in their hands issue number one of a comic... well maybe not like that but its still quite cool. I love looking back at how things started. The Federation of Astronomical Societies newsletter has been...
The Federation of Astronomical Societies newsletter has been running since June 1984 and over the last few months I've been sorting and scanning in a bunch of older records. This will eventually enhance the full online archive (before we only had back to 1999 online) of our newsletters.
The highlights of this issue are a write up of a night out at the Temple Observatory by a member of Coventry and Warwickshire Astronomical Society. Oh and if you were wondering, as I was, the telescope is still used. There is the first ever news round up, a stable of the FAS newsletter over the years and an article on how to make a star clock. I quite enjoyed seeing that the FAS subscription for the year was £7.50 - which would not be for a small society. Lets compare that to our medium society now of £44.50 (not including early payment) - its not that bad given we are looking back close to 3 decades.
Anyway, if you wish to read FAS Newsletter #1 you can get a pdf of it from the FAS website.]]>
Integration in Rtag:www.krioma.net,2011:/blog//1.16652012-01-26T08:45:01Z2012-01-26T08:50:15ZA reminder to myself of something I did quite a while go, how to do integration with [R]. dj integrand integrate(Vectorize(integrand), lower = 0, upper = 10) Of course the [R documentation] is of great use too....
[R].
Of course the [R documentation] is of great use too.
]]>
M42: The Orion Nebulatag:www.krioma.net,2012:/blog//1.16932012-01-15T00:11:43Z2012-01-15T00:13:36ZI love looking at the Orion nebula. Its always a lovely object to look at - its even better to take photos of it. I've not managed a properly deep image yet and the one that I'll show you below...
The same photo with different contrast levels chosen to highlight different parts of the Nebula.
]]>
Moon and plane over Birminghamtag:www.krioma.net,2012:/blog//1.16922012-01-08T20:45:20Z2012-01-08T21:08:39ZNot that long ago I managed to capture a lovely full Moon with a plane in the same field:...]]>
Convert 3gp to a better formattag:www.krioma.net,2012:/blog//1.17022012-01-02T11:21:24Z2012-01-02T11:32:51ZBy default my phone camera takes video in 3gp format - which isn't the best for sharing around... so here is how to convert it using ffmpeg on linux: Firstly you might need to install ffmpeg and some libraries: sudo...
By default my phone camera takes video in 3gp format - which isn't the best for sharing around... so here is how to convert it using ffmpeg on linux:
Firstly you might need to install ffmpeg and some libraries:
sudo apt-get install ffmpeg libavcodec-extra-52
You can just do:
ffmpeg -r 1 -i VIDEO0032.3gp -f avi -r 24 -vcodec libxvid -acodec libmp3lame movie.avi
but this gives some crapy quality - a better way is:
ffmpeg -r 1 -i VIDEO0034.3gp -f mp4 -acodec libmp3lame -ar 44100 -ab 128 -vcodec mpeg4 -maxrate 2000 -b 1500 -qmin 3 -qmax 5 -bufsize 4096 -g 300 -r 30000/1001 outfile2.mp4
AMI Small Array Movietag:www.krioma.net,2011:/blog//1.17012011-12-27T17:19:39Z2011-12-27T17:23:40ZA time lapse movie of the Arcminute Microkelvin Image (AMI) small array, I made by grabbing the webcam feed: Watch on youtube....
Arcminute Microkelvin Image (AMI) small array, I made by grabbing the webcam feed:
Watch on youtube.
]]>
Recover data with ddrescuetag:www.krioma.net,2011:/blog//1.16882011-12-24T12:13:44Z2011-12-24T12:11:32ZBeen handed a dead hard drive and want to try and recover some of the data? That's the problem I had the other day. So I turned to some good old linux tools to try and recover the data. One...
ddrescue. It works really well and given enough time can do a pretty solid job of data recovery - of course a better solution is proper backups but that's not always the way people go. Anyway, this is as much a note for myself but here is how to do some basic recovery using it:
sudo ddrescue -r 3 /dev/sde2 imaging loging
This images /dev/sed2 and produces the image file "imaging" with the log "logging".
To try and extract files you could use something like foremost
sudo foremost -w -i imaging -o /recovery/foremost
To make an audit of the files that can be recovered, recovering them with:
sudo foremost -i imaging -o /recovery/foremost2
]]>
Music of 2011tag:www.krioma.net,2011:/blog//1.16912011-12-22T11:44:48Z2011-12-23T08:52:48ZAt the end of each year, yes I know I'm a few days early but bah I have to move house at the start of the new year, I like to summarise the music I've been listening to. The last...
last.fm keeps hold of that info. I suspect that its missing about 20% of the music I listen to, which isn't scrobbled for one reason and another. It might be a little higher than that but I don't like to include the car journeys when we are using the CD player.
I was a little shocked to see the number 1 artist of the year, and indeed a few different names in the top 10. I'm putting this down to a Calgary scene influence... (number, artist, plays)
1 Daft Punk 1,163
2 Linkin Park 664
3 Avril Lavigne 552
4 My Chemical Romance 526
5 Tiësto 522
6 Rise Against 285
7 Angels & Airwaves 271
8 Dead by April 259
9 Lostprophets 239
10 Disturbed 238
I have to say I did love the Daft Punk Tron album so I'm guessing this is what dominates that. I also find Tiesto quite good for background music whilst I'm reading (in particular I really think Kaleidoscope fits in really well with the Garth Nix books I've been reading).
Not that much new for 2011 in this list though. Wonder if that will change in 2012. ]]>
A review of my 2011tag:www.krioma.net,2011:/blog//1.16902011-12-21T18:14:03Z2011-12-22T09:34:58Z2011 will always be one of the most memorable years in my life - it better be as I did get married :-) The year started off with me being back in the UK for my birthday which was great...
I went snowboarding for the first time:
Was painful and so far has turned out to be the last time. I just can't believe how poor my balance is.
In March my best-man Steve came over to visit me in Calgary and that was a lot of fun. Amongst the highlights of his trip over where seeing ice sculptures (and standing on the lake) at Lake Louise:
It wouldn't have been a proper trip over for him without taking in some hockey, so we went and
watched the outdoor WHL heritage game (that was damn cold!) and a couple of Flames games at the Saddledome:
We also went out and saw some world class snowboarding:
At the end of April I found it hard to get home, as stuff was falling off the building across the road from me!
In May, I said goodbye to Calgary as I had been offered a position at the University of Cambridge, back in the UK. I have to say I'm very glad to have come home and can spend a lot more time with my beautiful wife but I'll always have very fond memories of Calgary and miss lots about it - in particular the people. I wrote a detailed review of my time in Calgary. Though before I left Calgary I got to go to the wonderfully named small town of Vulcan:
Before I left Calgary I had a bit of a shock with a health issue, all was fine but lets just say I've since changed my work all hours, cut out go out drinking lots and living off caffeine.
In June I started at the Cavendish. I have to say I'm finding Cambridge a little different to the big
cities I'm used to. Oh and I got probably my best scientific result of the year accepted in June.
July was a very busy time. We were heavily involved in Wedding preparation though I did manage to stop for a few bits of fun. This included going to the British Grand Prix for the first time. Yes, it rained - it was a great time though.
July also involved my stag do, but we don't talk about that. It was great fun... we went go karting, had burgers and drank. Epic times.I even won a laser quest:
August... well this was a hectic month but all got more relaxed after the 19th. On the 19th we got married:
most of the better pictures aren't mine and well I feel bad about linking them so won't
(I wrote lots about the preparations before)
My wedding speech even made it onto youtube:
We took a short break to Bath as a honeymoon part 1.
In September I got to see my first supernova! (that's a star blowing up by the way)
and I took my first photo of a comet, Garrard:
October saw me start to settle into my new surroundings and I took my first trip out to the MRAO. I also had a paper on Ultra Steep Spectrum sources published.
At the of October we took our proper honeymoon out to Croatia, in particular Rovinj. It was absolutely beautiful.
The whole Istrian peninsula was full of interesting things, including lots of Roman sites. We also went inland to the Plitvice Lakes National Park - somewhere I'd love to spend a bit more time walking around.
In November I got to go to the University of Birmingham Vale Fireworks for the first time in 2 years so I was quite pleased - always
feels good to celebrate the saving of parliament by blowing up part of the country:
I got to take a trip over to Paris for ADASS XXI and I presented a poster on magnetic fields in GALFACTS. Elizabeth joined me after
the conference and we had a great weekend exploring Paris. Lots of fun.
At the end of the month we went out to Wast Hills and did some more of Messier catalogue but I went outside and tried out my camera at some night time constellation shots (Orion here above the dome):
and now its December. December has been a hectic month, and we are moving down to Cambridge in the new year so busy busy but we did manage to go out observing the other week and capture the Christmas tree cluser. Merry Christmas!
I'm looking forward to 2012... should be lots of fun. So far I have lots of observing planned, trip to India in Feb and tickets for the Olympics (alas just football). As a friend of mine said to me recently einen guten Rutsch ins neue Jahr
]]>
Observing last night - Jupiter + Dumbbelltag:www.krioma.net,2011:/blog//1.16892011-12-19T14:29:02Z2011-12-19T14:36:28ZIt wasn't particularly clear during our observing session, which was mostly aimed at taking a look at the ISS and Jupiter with our own eyes, but we focused the Meade out at Wast Hills and took a couple of images......
Jupiter with our own eyes, but we focused the Meade out at Wast Hills and took a couple of images...
If you change the contrast a little you can clearly see the Moons:
It was rather cloudy but we still managed a short exposure of the Dumbbell Nebula:
]]>
Creating FITS files in c++tag:www.krioma.net,2011:/blog//1.16812011-12-18T10:41:22Z2011-12-18T16:47:11ZA simple example, very similar to that given in the CFITSIO guidebook, on how to create a FITS file using CFITSIO. In this case I'm also building against some casacore libraries, but these aren't going to be used in this...
CFITSIO guidebook, on how to create a FITS file using CFITSIO. In this case I'm also building against some casacore libraries, but these aren't going to be used in this little code snippet but the idea is to use casacore todo further analysis. I'm hoping to post more here over time. Anyway the code (this can also be found as a filebuild_fits.cpp:
/* Create a FITS file, using cfitsio and some casacore libraries
by Samuel George
21-11-2011
Compile: g++ build_fits.cpp -o build_fits -lcasa_casa -lcfitsio
*/
#include // STL iostream
#include
#include
extern "C"{
#include
}
using std::cerr;
using std::cout;
using std::endl;
using std::string;
int main()
{
cout << "Create a FITS file" << endl;
int lenTime(10), status (0), lenFreq(20);
long naxis(2), naxes[2] = {lenTime,lenFreq};
long nelements (lenTime*lenFreq);
long fpixel (1), exposure (1500);
char comment[] ="Total Exposure Time";
cout << comment << endl;
float pixels[lenFreq][lenTime];
// create an array of pixels
try {
for (int ii(0); ii < lenTime; ii++){
for (int jj(0); jj
Save the code and compile like so:
g++ build_fits.cpp -o build_fits -lcasa_casa -lcfitsio
and then run with ./build_fits]]>
CASA leap seconds issuetag:www.krioma.net,2011:/blog//1.16822011-12-14T11:39:33Z2011-12-14T15:33:36ZI had an odd CASA issue yesterday. It was complaining about missing leap second information for TAI_UTC! Seemed really quite odd to me as I hadn't done anything to my installation... anyway simple fix is to update the latest leap...
cd $CASA_ROOT/data
rsync -avz rsync.aoc.nrao.edu::casadata .]]>
Spectral Extraction in Pythontag:www.krioma.net,2011:/blog//1.16872011-12-10T16:24:26Z2011-12-10T16:34:06ZI've wrote about how to read in FITS files in Python before, but I thought I'd readdress as I've been writing lots about fitting and wanted to build up to fitting properly calibrated data. So in this example I will...
how to read in FITS files in Python before, but I thought I'd readdress as I've been writing lots about fitting and wanted to build up to fitting properly calibrated data. So in this example I will add how to extract a spectra too. For this example I'm going to use a calibration lamp taken at the University of Birmingham Observatory. Both the script described here and the data will be given below. The output looks like:
So let's start with reading in a FITS file. Firstly ensure you have the most useful per-requisites ([Pyfits], [numpy], [scipy]) and import them:
We know that the spectra, by looking at the FITS file with DS9 has data in a small aperture. For now we will ignore any averaging. Lets take line 144 as we know there is data.
We have now extracted the data for a spectral column. I will follow this post up in the not to distant future with how to now identify the lines and run spectral calibration. Once you have this you can apply to a target spectra. Of course I've so far also ignored dark/bias counts but these can be easily dealt with.
The above can be downloaded as extract_spectra.py and the data Neon.fits (this is zipped, use gunzip).]]>
Christmas Tree Clustertag:www.krioma.net,2011:/blog//1.16862011-12-10T12:49:10Z2011-12-10T12:52:37ZProbably not the best picture I'll ever take of this source, I blame the bright moon, but here is a picture of the Christmas Tree Cluster that we took last night.... ... so "ho ho ho, merry Christmas"....
... so "ho ho ho, merry Christmas". ]]>
Multiple Gaussian Fitting in Pythontag:www.krioma.net,2011:/blog//1.16852011-12-09T15:04:53Z2011-12-09T15:14:23ZYesterday I showed you [how to fit a single Gaussian in some data]. Today lets deal with the case of two Gaussians. This came about due to some students trying to fit two Gaussian's to a shell star as the...
[how to fit a single Gaussian in some data]. Today lets deal with the case of two Gaussians. This came about due to some students trying to fit two Gaussian's to a shell star as the spectral line was altered from a simple Gaussian, actually there is a nice P-Cygni dip in there data so you should be able to recover the absorption line by this kind of fitting. Anyway, fitting 2 Gaussian's is basically the same thing as fitting one in python but with the added function. I'm thinking I'll work up how to deal with a whole spectrum including source finding. In this case what you have to deal with is that there are two sources and so a rough estimation of the peak position of both is crucial to the fit (well in the way it is implemented).
You basically want to end up with something like this:
Where the green points are the data, the blue dashed line the fit and the red line where the maximum is in the array. The fitting is done via a least-squares fitting routine.
In this example we will first start by generating some data, skipping the input from a file (but you can use the code from the 1 Gaussian example).
import scipy, numpy, pylab, asciidata
from numpy import *
from pylab import *
from scipy import *
from scipy.optimize import leastsq
#generate some data
gaussian = lambda x: 3*exp(-(10-x)**2/10.) + 1*exp(-(30-x)**2/10.)#change the parameters as you see fit
y_power = gaussian(arange(100))
x_pos = arange(100)
gauss_fit = lambda p, x: p[0]*(1/sqrt(2*pi*(p[2]**2)))*exp(-(x-p[1])**2/(2*p[2]**2))+p[3]*(1/sqrt(2*pi*(p[5]**2)))*exp(-(x-p[4])**2/(2*p[5]**2)) #1d Gaussian func
e_gauss_fit = lambda p, x, y: (gauss_fit(p,x) -y) #1d Gaussian fit
v0= [1,10,1,1,30,1] #inital guesses for Gaussian Fit. - just do it around the peaks
out = leastsq(e_gauss_fit, v0[:], args=(x_pos, y_power), maxfev=100000, full_output=1) #Gauss Fit
v = out[0] #fit parameters out
covar = out[1] #covariance matrix output
xxx = arange(min(x_pos),max(x_pos),x_pos[1]-x_pos[0])
ccc = gauss_fit(v,xxx) # this will only work if the units are pixel and not wavelength
fig = figure(figsize=(9, 9)) #make a plot
ax1 = fig.add_subplot(111)
ax1.plot(x_pos,y_power,'gs') #spectrum
ax1.plot(xxx,ccc,'b--') #fitted spectrum
ax1.axvline(x=xxx[where(ccc == max(ccc))[0]][0],color='r') #max position in data
setp(gca(), ylabel="power", xlabel="pixel position")
pylab.savefig("plotfitting.png")
print "p[0], a1: ", v[0]
print "p[1], mu1: ", v[1]
print "p[2], sigma1: ", v[2]
print "p[3], a2: ", v[0]
print "p[4], mu2: ", v[1]
print "p[5], sigma2: ", v[2]
The highlights of this issue are a write up of a night out at the Temple Observatory by a member of Coventry and Warwickshire Astronomical Society. Oh and if you were wondering, as I was, the telescope is still used. There is the first ever news round up, a stable of the FAS newsletter over the years and an article on how to make a star clock. I quite enjoyed seeing that the FAS subscription for the year was £7.50 - which would not be for a small society. Lets compare that to our medium society now of £44.50 (not including early payment) - its not that bad given we are looking back close to 3 decades.
Anyway, if you wish to read FAS Newsletter #1 you can get a pdf of it from the FAS website.]]>
]]>
