search

radio-astro-tools / spectral-cube

Library for reading and analyzing astrophysical spectral data cubes
http://spectral-cube.rtfd.org
BSD 3-Clause "New" or "Revised" License
98 stars 65 forks source link

find_line problem #908

Closed sylviankahane closed 6 months ago

sylviankahane commented 6 months ago

filename_1 = download_file("https://almascience.eso.org/dataPortal/member.uid___A001_X1465_X3a33.BrickMaser_sci.spw71.cube.I.manual.image.pbcor.fits", cache=True)

cube1 = SpectralCube.read(filename_1)

cube1.find_lines(chemical_name=' H2CS ').show_in_notebook() . . . File ~\AppData\Local\Programs\Python\Python312\Lib\site-packages\astropy\io\ascii\core.py:1457, in BaseReader.read(self, table) 1449 if len(str_vals) != n_cols: 1450 errmsg = ( 1451 f"Number of header columns ({n_cols}) inconsistent with " 1452 f"data columns ({len(str_vals)}) at data line {i}\n" 1453 f"Header values: {[x.name for x in cols]}\n" 1454 f"Data values: {str_vals}" 1455 ) -> 1457 raise InconsistentTableError(errmsg) 1459 for j, col in enumerate(cols): 1460 col.str_vals.append(str_vals[j])

InconsistentTableError: Number of header columns (1) inconsistent with data columns (2) at data line 13 Header values: ['<-- start of nrao-header.html -->'] Data values: ['<a href= http', '//www.nrao.edu/ ">']

keflavich commented 6 months ago

That error indicates that your installation of astroquery is out of date.

sylviankahane commented 6 months ago

Thank you

Upgrading astroquery worked. I have another problem now with cube_reproj from the same tutorial:

cube1vel

(SpectralCube with shape=(75, 250, 250) and unit=Jy / beam:

n_x: 250 type_x: RA---SIN unit_x: deg range: 266.534072 deg: 266.554577 deg n_y: 250 type_y: DEC--SIN unit_y: deg range: -28.713958 deg: -28.695975 deg n_s: 75 type_s: VRAD unit_s: km / s range: -43.509 km / s: 104.685 km / s

cube2vel_spatialspectralsmooth

SpectralCube with shape=(75, 420, 420) and unit=Jy / beam:

n_x: 420 type_x: RA---SIN unit_x: deg range: 266.537002 deg: 266.551600 deg n_y: 420 type_y: DEC--SIN unit_y: deg range: -28.711371 deg: -28.698569 deg n_s: 75 type_s: VRAD unit_s: km / s range: -43.509 km / s: 104.685 km / s

The following raises an error:

cube2vel_reproj = cube2vel_spatialspectralsmooth.reproject(cube1vel.header) cube2vel_reproj

.

.

.

ValueError: All values in reprojected cube are nan. This can be caused by an error in which coordinates do not 'round-trip'. Try setting roundtrip_coords=False. You might also check whether the WCS transformation produces valid pixel->world and world->pixel coordinates in each axis.

Sincerely

Syvlvian Kahane

On Fri, May 24, 2024 at 3:43 PM Adam Ginsburg @.***> wrote:

That error indicates that your installation of astroquery is out of date.

— Reply to this email directly, view it on GitHub https://github.com/radio-astro-tools/spectral-cube/issues/908#issuecomment-2129445159, or unsubscribe https://github.com/notifications/unsubscribe-auth/ADAD74A4CD3VVTEMJL6DQLLZD4YV5AVCNFSM6AAAAABIHJ4H6SVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDCMRZGQ2DKMJVHE . You are receiving this because you authored the thread.Message ID: @.***>

-- Sylvian (gmail: @.***)