Attention heat map

21 views
Skip to first unread message

C Oda

unread,
Nov 15, 2021, 6:58:44 AMNov 15
to astrometry
I'd love to use the Calibration data (Center, Size, Pixel Scale & Orientation) of publicly uploaded images to create a map of which areas of the sky have received the most attention, and with which accuracy :)

Has this ever been done?

Dustin Lang

unread,
Nov 15, 2021, 8:31:10 AMNov 15
to C Oda, astrometry
I made that map once, years ago.  I'd be happy to do a database dump for you to explore!

The database structure is managed by Django, via the models here, eg, here's Calibration

and the WCS values are all stored in the TanWCS table:

For making this heat-map, a rough approach would be to use the "SkyLocation" data,
which give the nearest power-of-two "HEALpix" pixel for each image.  (Here's where that gets populated:

If you let me know which tables you want, I can do the db dump.

cheers,
--dustin




On Mon, Nov 15, 2021 at 6:58 AM C Oda <ArndtS...@freenet.de> wrote:
I'd love to use the Calibration data (Center, Size, Pixel Scale & Orientation) of publicly uploaded images to create a map of which areas of the sky have received the most attention, and with which accuracy :)

Has this ever been done?

--
You received this message because you are subscribed to the Google Groups "astrometry" group.
To unsubscribe from this group and stop receiving emails from it, send an email to astrometry+...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/astrometry/5160e3f8-8f47-4e56-9361-a39869051b17n%40googlegroups.com.

C Oda

unread,
Nov 15, 2021, 8:46:26 AMNov 15
to astrometry
The TanWCS table seems to be the most informative, so I'll try to use that

In the meantime, I'll try to find out in which ways to map this :)

Thank you very much!

Dustin Lang

unread,
Nov 15, 2021, 2:09:14 PMNov 15
to C Oda, astrometry
I think you'll want to map it out on many scales.  Wide-field images are really a different kind of observation than single-deep-sky-object observations.

One thing you could do would be to put, say, a 100-pixel grid on the image and push each of those points through the WCS transformation to get RA,Dec.  Then accumulate those into, say, a Hammer-Aitoff or other all-sky projection.  That would give you basically a pixel count histogram.

Here's a database dump of the public solved image WCS values:

The columns are
        ra,dec,rad = tan.get_center_radecradius()
        print(imid, ra, dec, rad, tan.get_pixscale(),
              tan.get_orientation(), tan.get_parity(),
              tan.crval1, tan.crval2, tan.crpix1, tan.crpix2,
              tan.cd11, tan.cd12, tan.cd21, tan.cd22, tan.imagew, tan.imageh)

where 'tan' is a TanWCS,

(the one file is called "csv" but it's actually in space-separated format -- sorry)

cheers,
--dustin




C Oda

unread,
Nov 15, 2021, 6:18:48 PMNov 15
to astrometry
Here's my first attempt https://i.imgur.com/tTTjyQN.png !

I'll publish the code here: https://github.com/void4/astroheat

I'm pretty new to all this astro stuff, so more complex projections will take some time, I need to learn more about projections/mappings etc.
(I'm not even sure if this code is accurate but looking at the map the results make sense so far)

Clear Skies!
- coda

C Oda

unread,
Nov 15, 2021, 7:21:52 PMNov 15
to astrometry
Here's an even cooler looking heatmap using matplotlibs contourf() function!

https://raw.githubusercontent.com/void4/astroheat/main/contour_1637021811.png

Also using only 25000 of the 1899073 rows so far

Updated the repo with it. First, run main.py to get the old visualization but also export the grid data to grid.pickle. Then run display.py to produce this.
Reply all
Reply to author
Forward
0 new messages