4contour画图

最新推荐文章于 2024-08-09 15:34:32 发布

糖糖不要甜_a

最新推荐文章于 2024-08-09 15:34:32 发布

阅读量502

点赞数 1

文章标签： python

本文链接：https://blog.csdn.net/weixin_53421104/article/details/121472471

版权

本文介绍如何使用Python的matplotlib库来绘制等高线图，通过contour函数展示数据的分布和变化，适用于地理、物理、工程等多个领域的数据可视化。

摘要由CSDN通过智能技术生成

from os import sep
from astropy.io import fits
from astropy.wcs import WCS
from reproject import reproject_interp
import matplotlib.pyplot as plt
import numpy as np
from astropy.io import fits as pf
from copy import deepcopy
from regions import PixCoord
from regions import EllipseSkyRegion, EllipsePixelRegion
from astropy import units as u
from astropy.coordinates import SkyCoord
import pandas as pd
import sys
from matplotlib import colors
import matplotlib
import os
import matplotlib.cm as cm
from matplotlib.ticker import MultipleLocator
from scipy.interpolate import griddata

def draw_contour(source_name, source1,source2,levels,vm):#定义contour函数

  hdu = fits.open('/home/wy/Desktop/show/4c/'+source_name+'/''v.fits')[0]#打开fits文件
  hdr = (hdu.header)


  hdu1 = fits.open('/home/wy/Desktop/show/4c/'+source_name+'/'+source1)[0]
  hdr0 = deepcopy(hdu1.header)#深度拷贝头文件
  hdr1 = pf.Header({ 'SIMPLE':   True,
                    'BITPIX':    -32,
                     'NAXIS':       2,#改成两个轴，只保留需要的前两个轴
                    'NAXIS1':  hdr0['NAXIS1'],
                    'NAXIS2':  hdr0['NAXIS2'],
                    'CDELT1':  hdr0['CDELT1'],