Python 和 Origin 实验室之间的交互

ra 格式翻译：hh mm ss；十进制： 度 毫米 秒：

import pandas as pd
import astropy as ap
import numpy as np
from astropy import units as u
from astropy.coordinates import SkyCoord

df = pd.read_csv('WD.csv')
ra1 = df['icrsra'] #прямое восхождение
dec1 = df['icrsdec'] #склонение
ra2 = ['']*len(ra1)
dec2 = ['']*len(ra1)

for i in range(len(ra1)):
  if not pd.isnull(ra1[i]) and not pd.isnull(dec1[i]): # это убирает пустые значения, по необходимости можно изменить/убрать
    k = ra1[i]

    k = k.split()
    ra2[i] = float(k[0])*15+float(k[1])/4+float(k[2])/4/60
    k = dec1[i].split(' ')
    dec2[i] = float(k[0])+float(k[1])/60+float(k[2])/3600
  else: # ненужные, при необходимости изменить
    ra2[i] = 0
    dec2[i] = 0
c = SkyCoord(ra2*u.degree,dec2*u.degree, frame = 'icrs')

结果，我们得到 c 作为天空坐标。最好将其转换为正常格式：

sky_coords = []
for i in range(len(c)):
    sky_coords.append([c[i].l.degree,c[i].b.degree])

sky_coords 最终成为一个常规列表，如 [l, b]

my_df = pd.DataFrame(mylist, columns=['col1'])
my_df.to_csv('out33.csv')

这是将列表输出到单独的文件中

现在的起源。基本定义 f(x)：

function double Sos(double m, double q)
{
double sum;
double fx;
sum = 0; 
for(i = 1; i<45; i=i+1)
{
  fx = 2.7255*(1+wcol(1)[i]);
  sum = sum + (wcol(2)[i]-fx)^2/(wcol(3)[i])^2;
}

return sum;
}

适用的：

wcol(7)[7] = Sos(0.03,0.2) 

对于循环（嵌套在此处），如果循环不起作用，则检查 for(mi=4;mi<42;mi=mi+1) 后是否没有 SEMOLON：

for(mi = 4; mi<42; mi=mi+1)
{ m = (mi-1)/100;
  for(qi = 20; qi<41; qi=qi+1)
        {
         q = qi/100;
         wcol(mi)[qi] = Sos(m,q);
        }
}

示例如果：

for(i = 1; i<250; i=i+1)
{
    if (wcol(8)[i]>365)
    {wcol(8)[i] = wcol(8)[i] - 365;};
    else
    { };
}

线路

变量输入：符号“NOT SUITABLE, ONLY”

wcol(7)[7]$ = "privet";

string y$ = "poka";

使用字符串：

wcol(7)[7]$ = y$ + "privet";

这将使 wcol[1][899] 正常：（否则它将像文本一样，它将是 wcol(z)[k]$*6）

wcol(1)[899] = wcol(z)[k]$*6;

在以下条件下一切正常：

if (wcol(1)[1]$ == "2023-May-08")
{
    wcol(7)[9] = 1;
}