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

C++

Войти
Регистрация
Восстановить пароль
 
Kerim_Geophysic
0 / 0 / 0
Регистрация: 08.01.2016
Сообщений: 134
#1

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

13.04.2016, 03:17. Просмотров 210. Ответов 0
Метки нет (Все метки)

впервые встречаюсь с необходимостью что-то изменить в программе вроде 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-как есть на самом деле)
0
Вложения
Тип файла: rar demo.rar (5.7 Кб, 2 просмотров)
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
13.04.2016, 03:17
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Исправить программу для спрямления (гиперболических) кривых (C++):

Необходимо разработать 3d рисовалку кривых Базье и прямых векторов и Инкрементировать в программу. - C++ Builder
Необходимо разработать 3d рисовалку кривых Базье и прямых векторов и Инкрементировать в программу. Суть программы: существует 3d...

Исправить программу для умножения матриц - C++
Привет. Создал программу для умножения матриц, используя метод Виноград. Вот что получилось #include &lt;iostream&gt; #include &lt;math.h&gt; using...

Исправить программу для определения чётных слов - C++
Было задание: Составить программу(без типа данных string), которая должна принимать слова, выводить только чётные из них, считать...

Исправить программу для определения площади заштрихованой фигуры - C++
Добрый день! Помогите исправить программу в C++. Вот такие условия: Составить программу для определения методом Монте-Карло площади...

Для графа определить его двудольность и вывести обе доли (исправить программу) - C++
помогите исправить программу! запускается, на файл вывода пуст... а когда раскоментирываю - не компилруется... вот код: #include...

Построить семейство кривых D(b0) для разных значений f0 (программа высчитывает D только при фиксированном f0) - C++
Здравствуйте!:) Прошу подсказать мне как построить семейство кривых D(b0) для разных значений f0(здесь программа высчитывает D только при...

Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
13.04.2016, 03:17
Привет! Вот еще темы с ответами:

Написать программу для построения графиков кривых второго порядка (гипербола) - Turbo Pascal
Написать программу для построения графиков кривых второго порядка(гипербола). При этом пользователь выбирает,для какой именно функции...

Составить программу построения известных математических кривых - Visual Basic
заданных уравнениями в полярных координатах.

Алгоритм обратных гиперболических функций - Алгоритмы
Здравствуйте! Нужна помощь в составлении алгоритмов для arsh(x), arch(х), arth(х), arcth(x). Хотя бы составленный алгоритм любого из них...

Составьте программу построения кривых, заданных в параметрической форме - PascalABC.NET
Составьте программу построения кривых, заданных в параметрической форме: 1. x=4cos^(3)(t), y=4sin^(3)(t) ; 2. x=4t, y=t^(3)+t^(2)+1


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

Или воспользуйтесь поиском по форуму:
Ответ Создать тему
Опции темы

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