@@ -257,15 +257,16 @@ HDUList is the value.
257
257
.. doctest-remote-data ::
258
258
259
259
>>> from astroquery.esasky import ESASky
260
- >>> images = ESASky.get_images(position = " m51" , radius = " 20 arcmin" , missions = [' Herschel' , ' ISO-IR' ]) # doctest: +IGNORE_WARNINGS +IGNORE_OUTPUT
261
- Starting download of HERSCHEL data. (25 files)
262
- Downloading Observation ID: 1342188589 from http://archives.esac.esa.int/hsa/whsa-tap-server/data?RETRIEVAL_TYPE=STANDALONE&observation_oid=8618001&DATA_RETRIEVAL_ORIGIN=UI [Done]
263
- Downloading Observation ID: 1342188328 from http://archives.esac.esa.int/hsa/whsa-tap-server/data?RETRIEVAL_TYPE=STANDALONE&observation_oid=8637833&DATA_RETRIEVAL_ORIGIN=UI
260
+ >>> images = ESASky.get_images(position = " V* HT Aqr" , radius = " 15 arcmin" , missions = [' Herschel' , ' ISO-IR' ]) # doctest: +IGNORE_WARNINGS
261
+
262
+ Starting download of HERSCHEL data. (6 files)
263
+ Downloading Observation ID: 1342220557 from http://archives.esac.esa.int/hsa/whsa-tap-server/data?RETRIEVAL_TYPE=STANDALONE&observation_oid=8628906&DATA_RETRIEVAL_ORIGIN=UI [Done]
264
+ Downloading Observation ID: 1342221178 from http://archives.esac.esa.int/hsa/whsa-tap-server/data?RETRIEVAL_TYPE=STANDALONE&observation_oid=8628962&DATA_RETRIEVAL_ORIGIN=UI
264
265
...
265
266
>>> print (images) # doctest: +IGNORE_OUTPUT
266
267
{
267
268
'HERSCHEL': [{'70': [HDUList], '160': HDUList}, {'70': HDUList, '160': HDUList}, ...],
268
- 'ISO' : [HDUList, HDUList, HDUList, HDUList, ...]
269
+ 'ISO-IR ' : [HDUList, HDUList, HDUList, HDUList, ...]
269
270
...
270
271
}
271
272
@@ -274,12 +275,12 @@ parameter observation_id instead of target and position.
274
275
275
276
.. doctest-remote-data ::
276
277
277
- >>> images = ESASky.get_images(observation_ids = [" 100001010" , " 01500403" ], missions = [" SUZAKU" , " ISO-IR" ]) # doctest: +IGNORE_WARNINGS +IGNORE_OUTPUT
278
+ >>> images = ESASky.get_images(observation_ids = " 100001010" , missions = " SUZAKU" ) # doctest: +IGNORE_WARNINGS
279
+ >>> images = ESASky.get_images(observation_ids = [" 100001010" , " 01500403" ], missions = [" SUZAKU" , " ISO-IR" ])
278
280
279
281
Note that the fits files also are stored to disk. By default they are saved to the working directory but the location
280
282
can be chosen by the download_dir parameter.
281
283
282
-
283
284
Get maps
284
285
--------
285
286
@@ -291,15 +292,15 @@ position, radius and missions.
291
292
.. Skip testing examples with with hard-wired download dir values
292
293
.. doctest-remote-data ::
293
294
294
- >>> table_list = ESASky.query_region_maps(position = " m51 " , radius = " 20 arcmin" , missions = [' Herschel' , ' ISO-IR' ])
295
- >>> images = ESASky.get_maps(query_table_list = table_list, download_dir = " /home/user/esasky" ) # doctest: +SKIP
295
+ >>> table_list = ESASky.query_region_maps(position = " V* HT Aqr " , radius = " 15 arcmin" , missions = [' Herschel' , ' ISO-IR' ])
296
+ >>> images = ESASky.get_maps(query_table_list = table_list, download_dir = " /home/user/esasky" ) # doctest: +SKIP
296
297
297
298
This example is equivalent to:
298
299
299
300
.. doctest-remote-data ::
300
301
301
- >>> images = ESASky.get_images(position = " m51 " , radius = " 20 arcmin" , missions = [' Herschel' , ' ISO-IR' ],
302
- ... download_dir= " /home/user/esasky" ) # doctest: +SKIP
302
+ >>> images = ESASky.get_images(position = " V* HT Aqr " , radius = " 15 arcmin" , missions = [' Herschel' , ' ISO-IR' ],
303
+ download_dir="/home/user/esasky") # doctest: +SKIP
303
304
304
305
305
306
Get spectra
@@ -316,21 +317,21 @@ The methods returns a `dict` to separate the different missions. All mission exc
316
317
.. doctest-remote-data ::
317
318
318
319
>>> from astroquery.esasky import ESASky
319
- >>> spectra = ESASky.get_spectra(position = " m51" , radius = " 20 arcmin" , missions = [' Herschel' , ' XMM-NEWTON' ]) # doctest: +SKIP
320
- >>> spectra = ESASky.get_spectra(observation_ids = [" 02101201" , " z1ax0102t" ], missions = [" ISO-IR" , " HST-UV" ]) # doctest: +IGNORE_OUTPUT
320
+ >>> spectra = ESASky.get_spectra(position = " Gaia DR3 4512810408088819712" , radius = " 6.52 arcmin" ,
321
+ missions=['Herschel', 'XMM-NEWTON'])
322
+ >>> spectra = ESASky.get_spectra(observation_ids = [" 02101201" , " z1ax0102t" ], missions = [" ISO-IR" , " HST-UV" ])
321
323
322
324
or
323
325
324
326
.. doctest-remote-data ::
325
327
326
- >>> table_list = ESASky.query_region_spectra(position = " m51 " , radius = " 20 arcmin" ,
327
- ... missions= [' Herschel' , ' XMM-NEWTON' ])
328
- >>> spectra = ESASky.get_spectra_from_table(query_table_list = table_list)
328
+ >>> table_list = ESASky.query_region_spectra(position = " Gaia DR3 4512810408088819712 " , radius = " 6.52 arcmin" ,
329
+ missions=['Herschel', 'XMM-NEWTON'])
330
+ >>> spectra = ESASky.get_spectra_from_table(query_table_list = table_list, download_dir = " /home/user/esasky " )
329
331
dict: {
330
- 'HERSCHEL': {'1342211195 ': {'red' : {'HPSTBRRS' : HDUList}, 'blue' : {'HPSTBRBS': HDUList},
331
- '1342180796 ': {'WBS ' : {'WBS-H_LSB_5a' : HDUList}, 'HRS' : {'HRS-H_LSB_5a' : HDUList},
332
+ 'HERSCHEL': {'1342244919 ': {'red' : {'HPSTBRRS' : HDUList}, 'blue' : {'HPSTBRBS': HDUList},
333
+ '1342243607 ': {'SSW+SLW ' : {'SPSS' : HDUList},
332
334
...},
333
- 'HST-IR':[HDUList, HDUList, HDUList, HDUList, HDUList, ...],
334
335
'XMM-NEWTON' : [HDUList, HDUList, HDUList, HDUList, ...]
335
336
...
336
337
}
@@ -340,20 +341,12 @@ Here is another example for Herschel, since it is a bit special:
340
341
.. doctest-remote-data ::
341
342
342
343
>>> from astroquery.esasky import ESASky
343
- >>> result = ESASky.query_region_spectra(position = ' M51 ' , radius = ' 1arcmin ' , missions = [' HERSCHEL' ])
344
+ >>> result = ESASky.query_region_spectra(position = ' [SMB88] 6327 ' , radius = ' 1 arcmin ' , missions = [' HERSCHEL' ])
344
345
>>> herschel_result = result[' HERSCHEL' ]
345
346
>>> herschel_result[' observation_id' , ' target_name' , ' instrument' , ' observing_mode_name' , ' band' , ' duration' ].pprint()
346
347
>>> spectra = ESASky.get_spectra_from_table(query_table_list = [(' HERSCHEL' , herschel_result)], download_dir = ' Spectra_new' )
347
- >>> spectra[' HERSCHEL' ][' 1342211195' ][' red' ].keys()
348
- >>> spectra[' HERSCHEL' ][' 1342211195' ][' red' ][' HPSTBRRS' ].info()
349
- observation_id target_name instrument observing_mode_name band duration
350
- seconds
351
- -------------- ------------ ---------- ------------------------ ------------------- ---------
352
- 1342180479 M51 HIFI HifiMappingModeDBSRaster 555.45-571.56 GHz 8014.085
353
- 1342211195 m51 strip2-1 PACS PacsLineSpec 62.96-190.85 micron 1254.0
354
- 1342211188 m51 PACS PacsLineSpec 52.15-159.67 micron 5597.0
355
- 1342180478 M51 HIFI HifiMappingModeOTF 562.54-578.65 GHz 4012.741
356
- ...
348
+ >>> spectra[' HERSCHEL' ][' 1342249066' ][' SSW+SLW' ].keys()
349
+ >>> spectra[' HERSCHEL' ][' 1342249066' ][' SSW+SLW' ][' SPSS' ].info()
357
350
358
351
359
352
Solar System Object Crossmatch
@@ -424,7 +417,7 @@ Or download everything on an SSO by something like this:
424
417
.. doctest-remote-data ::
425
418
426
419
>>> from astroquery.esasky import ESASky
427
- >>> images= ESASky.get_images_sso(sso_name = " ganymede " )
420
+ >>> images= ESASky.get_images_sso(sso_name = " 2017 RN65 " )
428
421
429
422
430
423
This module also offers access to IMCCE's SsODNet resolver, which allows you to find solar and extra solar system
0 commit comments