获得两类相关点之间的线性关系

2024-08-31 08:44

本文主要是介绍获得两类相关点之间的线性关系,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

例如我是有两幅图,但是他们的单位不一样(一个温度,一个央斯基)原则上他们只是相差一个因子。可以利用一些线性拟合的关系来获得那个因子。其实代码和那个计算谱指数差不多,只是做了一些容易犯错误的修改(斜率的计算上)

from astropy.io import fits as pf
from astropy.wcs import WCS
import sys
from copy import deepcopy
from astropy.coordinates import SkyCoord
from reproject import reproject_interp
import matplotlib.pyplot as plt
import numpy as np
from astropy import units as u
import math
from statistics import mode
from astropy.convolution import convolve, Gaussian2DKernel
from astropy.modeling import Fittable1DModel, Parameter, models, fitting#例如data和data_cite 记得先平滑到一个分辨率上
cal_spectral(data[mark1==1].flatten()[::1], data_cite[mark1==1].flatten()[::1], 1.248, 1.248)def cal_alpha(data1,data2):n = data1.sizex = np.zeros(n) y = np.zeros(n)kk = 0for i in range(n):if np.isnan(data1[i]) or np.isnan(data2[i]): continuex[kk] = data1[i]y[kk] = data2[i]kk += 1if kk<10: return np.nan,np.nan,np.nan,np.nanp_init = models.Polynomial1D(1)fit_p = fitting.LinearLSQFitter()p = fit_p(p_init, x[:kk], y[:kk])a, b = p.c0.value, p.c1.valuereturn a, bdef cal_spectral(data1, data2, freq1, freq2):a1, b1 = cal_alpha(data1, data2)a2, b2 = cal_alpha(data2, data1)print(b1, 1/b2)alpha_1 = (b1 + 1/b2) / 2.d_alpha_1 = abs(1/b2 - b1)a1_p = a1 * 1.b1_p = b1 * 1.a2_p = -a2/b2b2_p = 1./b2#画图plt.rcParams.update({'font.size': 25})fig = plt.figure(figsize=(15,11))ax = fig.add_subplot(1, 1, 1)ax.plot(data1, data2, 'o')ax.plot(data1, a1_p + b1_p * data1)ax.plot(data1, a2_p + b2_p * data1)ss = r'$-$' + '%4.2f' % abs(alpha_1) + r'$\pm$' + '%4.2f' % d_alpha_1ax.text(0.1, 0.9, ss, horizontalalignment='left', verticalalignment='center', transform=ax.transAxes)ax.set_xlabel('T (K, at %5.3f GHz)' % freq1)ax.set_ylabel('T(K, Extrapolating E. M. Berkhuijsen observ-\n-ational data to a frequency of 1.248 GHz)')ax.xaxis.labelpad = 15ax.yaxis.labelpad = 8plt.savefig('Test_Temperature calibration.png', bbox_inches = 'tight')plt.show()

这篇关于获得两类相关点之间的线性关系的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!


原文地址:
本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.chinasem.cn/article/1123424

相关文章

CSS3中的字体及相关属性详解

《CSS3中的字体及相关属性详解》:本文主要介绍了CSS3中的字体及相关属性,详细内容请阅读本文,希望能对你有所帮助... 字体网页字体的三个来源:用户机器上安装的字体,放心使用。保存在第三方网站上的字体,例如Typekit和Google,可以link标签链接到你的页面上。保存在你自己Web服务器上的字

java Long 与long之间的转换流程

《javaLong与long之间的转换流程》Long类提供了一些方法,用于在long和其他数据类型(如String)之间进行转换,本文将详细介绍如何在Java中实现Long和long之间的转换,感... 目录概述流程步骤1:将long转换为Long对象步骤2:将Longhttp://www.cppcns.c

解决tomcat启动时报Junit相关错误java.lang.ClassNotFoundException: org.junit.Test问题

《解决tomcat启动时报Junit相关错误java.lang.ClassNotFoundException:org.junit.Test问题》:本文主要介绍解决tomcat启动时报Junit相... 目录tomcat启动时报Junit相关错误Java.lang.ClassNotFoundException

Maven中引入 springboot 相关依赖的方式(最新推荐)

《Maven中引入springboot相关依赖的方式(最新推荐)》:本文主要介绍Maven中引入springboot相关依赖的方式(最新推荐),本文给大家介绍的非常详细,对大家的学习或工作具有... 目录Maven中引入 springboot 相关依赖的方式1. 不使用版本管理(不推荐)2、使用版本管理(推

Python的time模块一些常用功能(各种与时间相关的函数)

《Python的time模块一些常用功能(各种与时间相关的函数)》Python的time模块提供了各种与时间相关的函数,包括获取当前时间、处理时间间隔、执行时间测量等,:本文主要介绍Python的... 目录1. 获取当前时间2. 时间格式化3. 延时执行4. 时间戳运算5. 计算代码执行时间6. 转换为指

Java中Date、LocalDate、LocalDateTime、LocalTime、时间戳之间的相互转换代码

《Java中Date、LocalDate、LocalDateTime、LocalTime、时间戳之间的相互转换代码》:本文主要介绍Java中日期时间转换的多种方法,包括将Date转换为LocalD... 目录一、Date转LocalDateTime二、Date转LocalDate三、LocalDateTim

golang获取当前时间、时间戳和时间字符串及它们之间的相互转换方法

《golang获取当前时间、时间戳和时间字符串及它们之间的相互转换方法》:本文主要介绍golang获取当前时间、时间戳和时间字符串及它们之间的相互转换,本文通过实例代码给大家介绍的非常详细,感兴趣... 目录1、获取当前时间2、获取当前时间戳3、获取当前时间的字符串格式4、它们之间的相互转化上篇文章给大家介

Vue中组件之间传值的六种方式(完整版)

《Vue中组件之间传值的六种方式(完整版)》组件是vue.js最强大的功能之一,而组件实例的作用域是相互独立的,这就意味着不同组件之间的数据无法相互引用,针对不同的使用场景,如何选择行之有效的通信方式... 目录前言方法一、props/$emit1.父组件向子组件传值2.子组件向父组件传值(通过事件形式)方

JavaScript Array.from及其相关用法详解(示例演示)

《JavaScriptArray.from及其相关用法详解(示例演示)》Array.from方法是ES6引入的一个静态方法,用于从类数组对象或可迭代对象创建一个新的数组实例,本文将详细介绍Array... 目录一、Array.from 方法概述1. 方法介绍2. 示例演示二、结合实际场景的使用1. 初始化二

Python实现PDF与多种图片格式之间互转(PNG, JPG, BMP, EMF, SVG)

《Python实现PDF与多种图片格式之间互转(PNG,JPG,BMP,EMF,SVG)》PDF和图片是我们日常生活和工作中常用的文件格式,有时候,我们可能需要将PDF和图片进行格式互转来满足... 目录一、介绍二、安装python库三、Python实现多种图片格式转PDF1、单张图片转换为PDF2、多张图