Форум программистов, компьютерный форум CyberForum.ru

Исправить программу для спрямления (гиперболических) кривых - C++

Восстановить пароль Регистрация
 
Kerim_Geophysic
0 / 0 / 0
Регистрация: 08.01.2016
Сообщений: 108
13.04.2016, 03:17     Исправить программу для спрямления (гиперболических) кривых #1
впервые встречаюсь с необходимостью что-то изменить в программе вроде C++
немного знаком с Матлабом
значит есть программа, ее описание есть здесь http://ahay.org/blog/2013/04/08/prog...e-month-sfnmo/
суть в том, что есть три формулы, которые можно найти по ссылке выше. По умолчанию используется первая формула.
в этой формуле to- значение по вертикальной оси на рисунке (время), х- значение по горизонтали (расстояние).
входным параметрами являются to(время) и V(скорость).
программа написана так, что to(время) и V(скорость) надо задавать в файле, который непонятным для меня образом формируется. Я хочу задавать эти параметры двумя столбцами: первый- to(время), второй: V. То есть каждому времени соответствует скорость.
В приведенном ниже коде этот входной файл именуется как velocity
код программы:
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
#include <math.h>
 
#include <rsf.h>
 
#include "fint1.h"
 
static float *vel, *par, h;
static bool slow, squared;
 
static float nmo_map(float t, int it) {
    float v;
    
    v = vel[it];
    if (!squared) v *= v;
    v = slow ? h*v : h/v;
    t = t*t + v;
    if (t > 0.) t=sqrtf(t);
    return t;
}
 
static float shifted_nmo_map(float t, int it) {
    float v, q, s;
    
    v = vel[it];
    s = par[it];
 
    if (!squared) v *= v;
    v = slow ? h*v : h/v;
    q = t*t + s*v;
    if (q > 0.) {
    t=t*(1.-1./s)+sqrtf(q)/s;
    return t;
    } else {
    return q;
    }
}
 
static float taner_nmo_map(float t, int it) {
    float v, a;
    
    v = vel[it];
    a = par[it];
 
    if (!squared) v *= v;
    v = slow ? h*v : h/v;
    t = t*t + v/(1+a*v);
    if (t > 0.) t=sqrtf(t);
    return t;
}
 
int main (int argc, char* argv[])
{
    fint1 nmo;
    bool half;
    int ix,ih, nt,nx,nw, nh, noff, nmask, CDPtype, mute, *mask;
    float dt, t0, h0, dh, dy, str;
    float *trace, *off;
    mapfunc nmofunc;
    sf_file cmp, nmod, velocity, offset, msk, het;
 
    sf_init (argc,argv);
    cmp      = sf_input("in");
    velocity = sf_input("velocity");
    nmod     = sf_output("out");
 
    if (NULL != sf_getstring("s")) {
        het = sf_input("s");
    nmofunc = shifted_nmo_map;  
    } else if (NULL != sf_getstring("a")) {
        het = sf_input("a");
        nmofunc = taner_nmo_map;
    } else {
    het = NULL;
    nmofunc =  nmo_map;
    } 
 
    if (SF_FLOAT != sf_gettype(cmp)) sf_error("Need float input");
    if (!sf_histint  (cmp,"n1",&nt)) sf_error("No n1= in input");
    if (!sf_histfloat(cmp,"d1",&dt)) sf_error("No d1= in input");
    if (!sf_histfloat(cmp,"o1",&t0)) sf_error("No o1= in input");
    if (!sf_histint  (cmp,"n2",&nh)) sf_error("No n2= in input");
    nx = sf_leftsize(cmp,2);
 
    if (!sf_getbool("half",&half)) half=true;
    /* if y, the second axis is half-offset instead of full offset */
    if (!sf_getfloat("str",&str)) str=0.5;
    /* maximum stretch allowed */
 
    if (!sf_getint("mute",&mute)) mute=12;
    /* mute zone */
 
    CDPtype=1;
    off = sf_floatalloc(nh);
 
    if (NULL != sf_getstring("offset")) {
    offset = sf_input("offset");
    if (SF_FLOAT != sf_gettype(offset)) sf_error("Need float offset");
    noff = sf_filesize(offset);
    if (noff == nh) {
        sf_floatread (off,nh,offset);
    } else if (noff != nh*nx) {
        sf_error("Wrong dimensions in offset");
    }
    } else {
    if (!sf_histfloat(cmp,"d2",&dh)) sf_error("No d2= in input");
    if (!sf_histfloat(cmp,"o2",&h0)) sf_error("No o2= in input");
 
    if (sf_histfloat(cmp,"d3",&dy) && !sf_getint("CDPtype",&CDPtype)) {
        CDPtype=half? 0.5+dh/dy : 0.5+0.5*dh/dy;
        if (CDPtype < 1) {
        CDPtype=1;
        } else if (1 != CDPtype) {
        sf_histint(cmp,"CDPtype",&CDPtype);
            sf_warning("CDPtype=%d",CDPtype);
        }
    }       
 
    for (ih = 0; ih < nh; ih++) {
        off[ih] = h0 + ih*dh; 
    }
 
    noff = nh;
    offset = NULL;
    }
    
    if (NULL != sf_getstring("mask")) {
    msk = sf_input("mask");
 
    if (SF_INT != sf_gettype(msk)) sf_error("Need integer mask");
    nmask = sf_filesize(msk);
    mask = sf_intalloc(nh);
 
    if (nmask == nh) {
        sf_intread (mask,nh,msk);
    } else if (nmask != nh*nx) {
        sf_error("Wrong dimensions in mask");
    }
    } else {
    nmask = nh;
 
    msk = NULL;
    mask = NULL;
    }
 
    if (!sf_getbool("slowness",&slow)) slow=false;
    /* if y, use slowness instead of velocity */
 
    if (!sf_getbool("squared",&squared)) squared=false;
    /* if y, the slowness or velocity is squared */
 
    if (!sf_getfloat ("h0",&h0)) h0=0.;
    /* reference offset */
    if (half) h0 *= 2.;
    if (!sf_getint("extend",&nw)) nw=4;
    /* trace extension */
 
    trace = sf_floatalloc(nt);
    vel   = sf_floatalloc(nt);
    if (NULL != het) {
    par = sf_floatalloc(nt);
    } else {
    par = NULL;
    }
 
    nmo = fint1_init (nw, nt, mute);
    
    for (ix = 0; ix < nx; ix++) {
    sf_warning("CMP %d of %d;",ix+1,nx);
 
    sf_floatread (vel,nt,velocity); 
    if (NULL != het) sf_floatread(par,nt,het);
    if (NULL != offset && noff != nh) sf_floatread (off,nh,offset);
    if (NULL != msk && nmask != nh) sf_intread (mask,nh,msk);
 
    for (ih = 0; ih < nh; ih++) {
        sf_floatread (trace,nt,cmp);
        
        /* skip dead traces */
        if (NULL != msk && 0==mask[ih]) {
        sf_floatwrite (trace,nt,nmod);
        continue;
        }
        
        fint1_set(nmo,trace);
 
        h = off[ih];
        if (NULL == offset) h += (dh/CDPtype)*(ix%CDPtype); 
        if (half) h *= 2;
        h = h*h - h0*h0;
        
        stretch(nmo,nmofunc,nt,dt,t0,nt,dt,t0,trace,str);
        sf_floatwrite (trace,nt,nmod);
    }
    }
    sf_warning(".");
 
    exit (0);
}
 
/*  $Id$     */
как мне изменить этот файл, чтобы velocity задавался в файле в котором я хочу? (см прикрепленные файлы, pick.rsf- как я хочу, picks.rsf-как есть на самом деле)
Вложения
Тип файла: rar demo.rar (5.7 Кб, 2 просмотров)
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
13.04.2016, 03:17     Исправить программу для спрямления (гиперболических) кривых
Посмотрите здесь:

C++ исправить программу
C++ Исправить программу
Необходимо разработать 3d рисовалку кривых Базье и прямых векторов и Инкрементировать в программу. C++ Builder
Исправить программу C++
C++ Исправить программу
Исправить программу C++
Исправить программу С++ C++
Для графа определить его двудольность и вывести обе доли (исправить программу) C++
Построить семейство кривых D(b0) для разных значений f0 (программа высчитывает D только при фиксированном f0) C++
C++ Исправить программу для умножения матриц
Исправить программу для определения площади заштрихованой фигуры C++
Исправить программу для определения чётных слов C++

Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Ответ Создать тему
Опции темы

Текущее время: 23:12. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru