I have made a Python script to plot the Doppler and geometry for a DSLWP-B observation. This can be useful for several things: predicting the downlink and uplink Doppler corrections, checking for the Doppler of the Moonbounce signal, checking the position of DSLWP-B in the sky, and studying Moonbounce geometry.
Author: destevez
Geometry for DSLWP-B Moonbounce
I have already spoken about the Moonbounce signal from DSLWP-B in several posts. To sum up, DSLWP-B is a Chinese satellite that is orbiting the Moon since May 25. The satellite has an Amateur payload that transmits GMSK and JT4G telemetry in the 70cm Amateur satellite band. This signal can be received by well equipped groundstations on Earth, including the 25m radiotelescope at Dwingeloo, in the Netherlands (and also by much smaller stations).
The people at Dwingeloo publish the recordings that they make of the RF signal. In two of these recordings, the signal from DSLWP-B is received not only via the direct path, but also through a reflection off the Moon’s surface. The reflected signal is around 25dB weaker, usually has a different Doppler shift, and has a Doppler spread of around 50 to 200Hz.
What I find most interesting about this is that of all the days that Dwingeloo has observed DSLWP-B, in only two of them (on 2018-10-07 and 2018-10-19) the Moonbounce signal has been visible. Mathematically, using a specular reflection on a sphere model, whenever the satellite is visible directly, there is also a ray from the spacecraft that reflects off the lunar surface and arrives at the groundstation (see the proof here). Therefore, I think that there must be something about the particular geometry of the days 7th and 19th that made the Moon reflections relatively stronger and hence detectable. Here I use GMAT to study the orbital geometry when the reflections were received.
DSLWP-B JT4G decoded via Moonbounce
I have already spoken about the Moonbounce signal from DSLWP-B that was received in Dwingeloo on 2018-10-07. I have matched it against the Doppler predictions and cross-correlated it against the direct signal. Since the reflected signal presented a high Doppler spread, decoding the GMSK data from the reflected signal would be very difficult or impossible.
On the other hand, JT4G is a digital mode designed for Earth-Moon-Earth microwave communications, so it is tolerant to high Doppler spreads. However, the reflections of the B0 transmitter at 435.4MHz, which contained the JT4G transmissions, were very weak, so I had not attempted to decode the JT4G Moonbounce signal.
On 2018-10-19, the Moonbounce signal from DSLWP-B was again visible in Dwingeloo’s recordings. I have used the 2018-10-19T17:53:35 435.4MHz recording and managed to decode the Moonbounce signal of one out of the five JT4G transmissions that appear in the recording.
To extract the data from the recording to WAV files that can be read by WSJT-X, I have used the following Jupyter notebook. Then I have used WSJT-X version 2.0.0-rc3 to try to decode the Moonbounce signal. Since the JT4 decoder only decodes a single signal at the selected frequency, it is enough to select the frequency of the Moonbounce signal in WSJT-X. The direct signal will not be decoded, even though it is also present in the WAV file.
The only transmission that I have managed to decode was made at 18:11 UTC. The two screenshots below show WSJT-X decoding the WAV file extracted from the recording.
Note the direct signal with a lowest tone at 1800Hz. The reflected signal is very faint, with a lowest tone at 700Hz. The Doppler spread of the reflected signal is large, approximately 200Hz, although it is difficult to judge from the spectrum.
When the WAV file is created, I also compensate for a linear frequency drift of 25Hz per minute due to Doppler, but this is not essential to obtain a valid decode.
The WAV file that produces a decode can be downloaded here. This file can be opened directly by WSJT-X.
Keeping track of DSLWP-B images
While writing my previous report, I had to keep track of when each of the images had been taken and when it was downloaded. This can be a bit tricky, since the images do not carry a timestamp. Here I describe the techniques that can bee used to identify and keep track of the images.
Download of DSLWP-B Moon and Earth pictures continues
In a previous post I spoke about the images of the Moon and Earth taken by the Chinese lunar-orbiting satellite DSLWP-B between October 6 and 10. Some images taken during these days hadn’t been downloaded yet. The activities have continued during this week by downloading the remaining images and taking and downloading new images. This is a report of the images downloaded between October 14 and 19.
DSLWP-B Moonbounce cross-correlation analysis
In one of my latest posts I commented on the Moonbounce signal of the Chinese lunar satellite DSLWP-B, as received in Dwingeloo. In the observation made in 2018-10-07 Cees Bassa discovered a signal in the waterfall of the Dwingeloo recordings that seemed to be a reflection off the Moon of DSLWP-B’s 70cm signal. My analysis showed that the Doppler of this signal was compatible with a specular reflection on the lunar surface.
In this post I study the cross-correlation of the Moonbounce signal against the direct signal. This gives some information about how the radio signals behave when reflecting off the Moon. Essentially, we compute the Doppler spread and time delay of the Moonbounce channel.
Report of the DSLWP-B Amateur observations of the Moon and Earth
The Chinese microsatellite DSLWP-B has been in lunar orbit since 25 May 2018. This satellite carries an Amateur radio payload which includes a small 640×480 CCD camera. The JPEG images taken by the camera can be transmitted using the SSDV protocol at 125 bits per second in the 70cm Amateur satellite band.
Update 17:00 UTC: Wei comments that the camera sensor is CMOS, not CCD, and it has 2592×1944 pixels. The image is resampled to 640×480 to save memory and bandwidth.
The orientation of the camera is fixed: the camera is mounted looking in the opposite direction of the solar panel, which is usually kept pointing directly to the Sun. Therefore, the camera is usually looking directly away from the Sun. The possibility of imaging celestial bodies such as the Moon and the Earth depends on the relative positions of these and the Sun.
During the first week of October there was a new Moon, which implied that it was possible to take images of the Moon and the Earth, as I have described in this post and this other post. This is a report of all the images taken and downloaded during the observation window.
DSLWP-B Moonbounce
If you have been following my latest posts, you will know that a series of observations with the DSLWP-B Inory eye camera have been scheduled over the last few days to try to take and download images of the Moon and Earth (see my last post). In a future post I will do a chronicle of these observations.
On October 6 an image of the Moon was taken to calibrate the exposure of the camera. This image was downlinked on the UTC morning of October 7. The download was commanded by Reinhard Kuehn DK5LA and received by the Dwingeloo radiotelescope.
Cees Bassa observed that in the waterfalls of the recordings made in Dwingeloo a weak Doppler-shifted signal of the DSLWP-B GMSK signal could be seen. This signal was a reflection off the Moon.
As far as I know, this is the first reported case of satellite-Moon-Earth (or SME) propagation, at least in Amateur radio. Here I do a Doppler analysis confirming that the signal is indeed reflected on the Moon surface and do some general remarks about the possibility of receiving the SME signal from DSLWP-B. Further analysis will be done in future posts.
DSLWP-B upcoming observations of Moon and Earth
Yesterday I looked at the photo planning for DSLWP-B, studying the appropriate times to take images with the DSLWP-B Inory eye camera so that there is a chance of getting images of the Moon or Earth. As I remarked, the Earth will be in view of the camera over the next few days, so this is a good time to plan and take images.
I asked Wei Mingchuan BG2BHC to compare his calculations with mine and shortly after this he emailed me his planning for observations between October 8 and 10. After measuring the field of view of the camera as 37×28 degrees, we can plot the angular distances between the Moon or Earth and the centre of the camera to check if the celestial body will be in the field of view of the camera.
The image below shows the angular distance between these celestial bodies and the centre of the camera. As we did in the photo planning post, we assume that the camera points precisely away from the Sun. Since the Moon and Earth (especially the Moon) have an angular size of several degrees, we plot the centre of these objects with a dashed line and the edges which are nearest and furthest from the camera centre with a solid line.
The field of view of the camera is represented with dotted red lines. Since the field of view is a rectangle, we have one mark for the minimum field of view, which is attained between the centre of the image and the centre of the top or bottom edge, and another mark for the maximum field of view, which is attained between the centre of the image and each of the corners.
The rotation of the spacecraft around the camera axis is not controlled precisely, so objects between the two red lines may or may not appear in the image depending on the rotation. Objects below the lower red line are guaranteed to appear if the pointing of the camera is correct and objects above the upper red line will not appear in the image, regardless of the rotation around the camera axis.
To calibrate the exposure of the camera, an image was taken yesterday on 2018-10-06 13:55 UTC. This time is marked in the figure above with an orange line. The image was downloaded this UTC morning. The download was commanded by Reinhard Kuehn DK5LA and received by Cees Bassa and the rest of the PI9CAM team in Dwingeloo. This image is shown below.
The image shows an over exposed Moon. Here we are interested in using the image to confirm the orientation of the camera. The distance between the centre of the image and the edge of the Moon is 240 pixels, which amounts to 14 degrees. The plot above gives a distance of 11 degrees between the edge of the Moon and the camera centre.
Thus, it seems that the camera is pointed off-axis by 3 degrees. This error is not important for scheduling camera photos, since an offset of a few degrees represents a small fraction of the total field of view and the largest error in predicting what will appear in the image is due to rotation of the spacecraft around the camera axis.
The observations planned by Wei for the upcoming days are shown in the plot above by green lines. The start of the observation is marked with a dashed line and the end of the observation (which is 2 hours later) is marked with a dotted line. The camera should take an image at the beginning of the observation and then we have 2 hours to download the image during the rest of the observation.
We see that Wei has taken care to schedule observations exactly on the next three times that the Moon will be closest to the camera centre. This gives the best chance of getting good images of the lunar surface (but the Moon will only fill the image partially, as in the picture shown above).
There are also two additional observations planned when the Moon is not in view. The first, on October 8 is guaranteed to give a good image of the Earth. The second, on October 10 will only give an image of the Earth if the rotation of the spacecraft is right.
The orbital calculations for the plot shown above have been done in GMAT. I have modified the photo_planning.script script to output a report with the coordinates of the Earth and the Moon in the Sun-pointing frame of reference (see the photo planning post).
The angle between the centre of the camera and the centre of the Earth of the Moon can be calculated as\[\arccos\left(\frac{x}{\sqrt{x^2+y^2+z^2}}\right),\]where \((x,y,z)\) are the coordinates of the celestial body in the Sun-pointing frame of reference. The apparent angular radius of the celestial body can be computed as\[\arcsin\left(\frac{r}{\sqrt{x^2+y^2+z^2}}\right),\]where \(r\) is the mean radius of the body.
These calculations and the plot have been made in this Jupyter notebook.
DSLWP-B camera field of view
In my previous post, I wondered about what was the field of view of the Inory eye camera in DSLWP-B. Wei Mingchuan BG2BHC has answered me on Twitter that the field of view of the camera is 14×18.5 degrees. However, he wasn’t clear about whether these figures are measured from the centre of the image to one side or between two opposite sides of the image. I guess that these values are measured from the centre to one side, since otherwise the total field of view of the camera seems too small.
Here I measure the field of view of the camera using the image of Mars and Capricornus taken on August 4, confirming that these numbers are measured from the centre of the image to one side, so the total field of view is 28×37 degrees.