// TODO: font and style
void AspUtil::writeFields(char *str, int x, int y, int w, int h) {
char fontName[64];
getOptName(TEXT,fontName);
set_graphics_font(fontName);
aip_drawString(str, x, y, 0, PARAM_COLOR);
set_graphics_font("Default");
}
int dpir(int argc, char *argv[], int retc, char *retv[])
/*************/
{ int groups;
int resets;
int i, min, ysize;
char name[64];
char tmpstr[MAXSTR];
(void) retc;
(void) retv;
if (argc > 1 && strcmp(argv[1],"off") == 0) {
set_dpir_flag(0,"");
redo_dpf_dpir();
RETURN;
} else if (argc > 1 && strcmp(argv[1],"turnoff") == 0) {
set_dpir_flag(0,"");
RETURN;
}
vertflag=0;
if ( ! P_getstring(GLOBAL, "integralOrient", tmpstr, 1, MAXSTR) )
{
vertflag = (strcmp(tmpstr,"vert") == 0);
}
if (argc > 1 && strcmp(argv[1],"vert") == 0) {
vertflag=1;
}
if(argc == 1 || strcmp(argv[1],"noli") != 0) {
execString("nli\n");
}
if (check_int(&resets))
{
Werrprintf("intmod must not be set to 'full'");
ABORT;
}
if (resets <= 1)
{
Werrprintf("no integral resets are defined");
ABORT;
}
plot = (argv[0][0] == 'p');
if (select_init(
0,
plot+1,
NO_FREQ_DIM,
NO_HEADERS,
DO_CHECK2D,
DO_SPECPARS,
NO_BLOCKPARS,
NO_PHASEFILE
)) {
ABORT;
}
if(strlen(argv[0]) > 4 && argv[0][4] == 'N') {
normInt = 1;
} else if(strlen(argv[0]) > 3 && argv[0][3] == 'N') {
normInt = 1;
}
if ((groups = getregions(resets)) == 0) {
ABORT;
}
min=mnumxpnts;
for (i = 1; i <= groups; i++) {
if(gb[i].lowlimit <= 0 || gb[i].highlimit <= 0
|| gb[i].lowlimit >= (mnumxpnts-xcharpixels)
|| gb[i].highlimit >= (mnumxpnts-xcharpixels)) continue;
if((gb[i].highlimit-gb[i].lowlimit) < min)
min = gb[i].highlimit-gb[i].lowlimit;
}
if(min != 0 && min < 3*xcharpixels) ysize = 3; else ysize=2;
if (!plot)
{
setwindows();
dispcalib = (float) (mnumypnts-ymin) / (float) wc2max;
CharSize = ycharpixels + ycharpixels / 5;
getOptName(INTEG_MARK,name);
set_line_thickness(name);
getOptName(INTEG_NUM,name);
set_graphics_font(name);
}
else
{
double size = G_getCharSize("IntegralNum");
charsize(size);
//charsize((double)0.7);
dispcalib = ppmm / ymultiplier;
CharSize = ycharpixels;
}
if (showPlotBox()) dpir_off=0.0;
else dpir_off=9.0;
if(dispcalib>0) {
int size= 12;
G_getSize("IntegralNum", &size);
i = (int)(vp*dispcalib);
if (vertflag && i < 4*size)
{
char cmd[64];
i = (int)(4*size/dispcalib)+1;
if(vertflag) sprintf(cmd,"vp=%d ds %s('vert')\n",i,argv[0]);
else sprintf(cmd,"vp=%d ds %s\n",i,argv[0]);
if(!plot) execString(cmd);
//Werrprintf("pir requires minimum vp of %d",i);
RETURN;
} else if (i < 3*size)
{
char cmd[64];
i = (int)(3*size/dispcalib)+1;
if(vertflag) sprintf(cmd,"vp=%d ds %s('vert')\n",i,argv[0]);
else sprintf(cmd,"vp=%d ds %s\n",i,argv[0]);
if(!plot) execString(cmd);
//Werrprintf("pir requires minimum vp of %d",i);
RETURN;
}
} else {
if (vp < 30)
{
char cmd[64];
i=30;
if(vertflag) sprintf(cmd,"vp=%d ds %s('vert')\n",i,argv[0]);
else sprintf(cmd,"vp=%d ds %s\n",i,argv[0]);
if(!plot) execString(cmd);
//Werrprintf("pir requires minimum vp of 30");
RETURN;
}
}
//color(PARAM_COLOR);
color(INTEG_MARK_COLOR);
if (debug1)
{
Wscrprintf("dispcalib= %g, dfpnt2= %d, CharSize= %d, bottom= %d\n",
dispcalib,dfpnt2,CharSize,
dfpnt2 - CharSize + (int) (dispcalib * (vp - dpir_off - 5.0)));
Wscrprintf("index highlimit lowlimit below\n");
for (i = 1; i <= groups; i++)
Wscrprintf("%d %d %d %s\n",i,gb[i].highlimit,gb[i].lowlimit,
(gb[i].below) ? "true" : "false");
}
if(vertflag) {
for (i = 1; i <= groups; i++) label_pirv(i);
} else {
for (i = 1; i <= groups; i++) label_pir(i, ysize);
}
if (!plot) {
char cmd[64];
if(argc > 1) {
if (isReal(argv[1]) ) sprintf(cmd,"%s(%s",argv[0],argv[1]);
else sprintf(cmd,"%s('%s'",argv[0],argv[1]);
for(i=2;i<argc;i++)
{
if (isReal(argv[i]) )
{
strcat(cmd,",");
strcat(cmd,argv[i]);
}
else
{
strcat(cmd,",'");
strcat(cmd,argv[i]);
strcat(cmd,"'");
}
}
strcat(cmd,")\n");
} else sprintf(cmd,"%s\n",argv[0]);
Wsetgraphicsdisplay("ds");
set_dpir_flag(1,cmd);
}
releaseAllWithId("dpir");
endgraphics();
disp_status(" ");
set_graphics_font("Default");
RETURN;
}
void AspInteg::display(spAspCell_t cell, spAspTrace_t trace, spAspDataInfo_t dataInfo,
int integFlag, double is, double off) {
set_transparency_level(transparency);
disFlag = integFlag;
float *traceData = trace->getFirstDataPtr();
if(!traceData) return;
int maxpts = trace->getTotalpts();
if(maxpts<1) return;
absValue = trace->getInteg(sCoord[0].x,sCoord[1].x);
double x1=pCoord[0].x=cell->val2pix(HORIZ,sCoord[0].x);
double x2=pCoord[1].x=cell->val2pix(HORIZ,sCoord[1].x);
// get integ data
if(m_data) delete[] m_data;
int p1 = trace->val2dpt(sCoord[0].x);
int p2 = trace->val2dpt(sCoord[1].x);
m_datapts = p2-p1+1;
m_data = new float[m_datapts];
float sum = 0;
int k=0;
for(int i=p1;i<=p2;i++) {
//sum += (*(traceData+i) - *(traceData+p1));
sum += (*(traceData+i));
(*(m_data+k)) = sum;
k++;
}
double vcali = cell->getCali(VERT);
m_scale = is/(double)maxpts;
m_yoff = off*vcali;
double vscale = dataInfo->getVscale();
if(vscale == 0) vscale = 1.0;
double scale = m_scale*vscale*vcali;
// add vp
double yoff = m_yoff + (dataInfo->getVpos())*vcali;
double y1=pCoord[0].y=cell->val2pix(VERT,(*(m_data))*m_scale) - m_yoff;
double y2=pCoord[1].y=cell->val2pix(VERT,(*(m_data+m_datapts-1))*m_scale) - m_yoff;
//AspUtil::drawMark((int)x1,(int)y1,ACTIVE_COLOR);
//AspUtil::drawMark((int)x2,(int)y2,ACTIVE_COLOR);
// check x1, x2 boundary
double px,pw,py,ph,px2;
cell->getPixCell(px,py,pw,ph);
px2=px+pw;
if(x1 < px && x2 < px) return;
if(x1 > px2 && x2 > px2) return;
if(x1 < px) x1=px;
if(x2 > px2) x2=px2;
int roiColor;
if(selected == ROI_SELECTED) roiColor=ACTIVE_COLOR;
else roiColor = INT_COLOR;
char thickName[64];
int thickness; // used by functions
int thick=0; // use thickness set by functions
string value;
if(integFlag & SHOW_INTEG) {
getOptName(INTEG_LINE,thickName);
//AspUtil::getDisplayOption(string(thickName)+"Thickness",value);
//thickness = atoi(value.c_str());
//set_spectrum_width(thickness);
//set_line_width(thickness);
set_spectrum_thickness(thickName,thickName,1.0);
set_line_thickness(thickName);
cell->drawPolyline(m_data,m_datapts,1,x1,py,x2-x1,ph,roiColor,scale,yoff);
int i = selectedHandle-1;
if(i>=0 && i < npts) {
AspUtil::drawMark((int)pCoord[i].x,(int)pCoord[i].y,ACTIVE_COLOR,thick);
}
}
labelX=labelY=labelW=labelH=0;
if((integFlag & SHOW_LABEL)) {
if(integFlag & SHOW_VERT_LABEL) rotate = 1;
else rotate=0;
string labelStr;
getLabel(dataInfo,labelStr,labelW,labelH);
labelX=(int)(labelLoc.x+(x1+x2)/2) - labelW/2;
labelY=(int)(labelLoc.y+(y1+y2)/2) - labelH;
if(labelY<(py+labelH)) labelY = (int)py+labelH;
char fontName[64];
if(LABEL_SELECTED) {
int labelColor;
setFont(labelColor);
AspUtil::drawString((char *)labelStr.c_str(), labelX,labelY, labelColor, "", rotate);
} else {
getOptName(INTEG_LABEL,fontName);
AspUtil::drawString((char *)labelStr.c_str(), labelX, labelY, -1, fontName,rotate);
}
labelY -= labelH;
}
if(selected == ROI_SELECTED) roiColor=ACTIVE_COLOR;
else roiColor = INTEG_MARK_COLOR;
getOptName(INTEG_MARK,thickName);
AspUtil::getDisplayOption(string(thickName)+"Thickness",value);
thickness = atoi(value.c_str());
if(thickness < 1) thickness=1;
//set_spectrum_width(thickness);
//set_line_width(thickness);
set_spectrum_thickness(thickName,thickName,1.0);
set_line_thickness(thickName);
if((integFlag & SHOW_VALUE) && (integFlag & SHOW_VERT_VALUE)) {
string labelStr;
int cwd, cht, ascent, descent;
getValue(dataInfo,labelStr,cht,cwd);
// overwrite cht, cwd with fixed length
char str[MAXSTR];
sprintf(str,"%.2f",100.00);
GraphicsWin::getTextExtents(str, 14, &ascent, &descent, &cht);
cwd = ascent + descent;
double y = py + ph - cht - 5*thickness;
Dpoint_t p1,p2;
p1.x=x1;
p2.x=x2;
p1.y=p2.y=y;
AspUtil::drawLine(p1,p2,roiColor,thick);
p1.x=p2.x=x1;
p2.y=y-5*thickness;
AspUtil::drawLine(p1,p2,roiColor,thick);
p1.x=p2.x=x2;
p2.y=y-5*thickness;
AspUtil::drawLine(p1,p2,roiColor,thick);
p1.x=p2.x=0.5*(x1+x2);
p2.y=y+5*thickness;
AspUtil::drawLine(p1,p2,roiColor,thick);
char fontName[64];
getOptName(INTEG_NUM,fontName);
y = py + ph - cwd/2;
AspUtil::drawString((char *)labelStr.c_str(), (int)p1.x-cwd/2, (int)y, -1, fontName,1);
} else if(integFlag & SHOW_VALUE) {
string labelStr;
int cwd, cht;
getValue(dataInfo,labelStr,cwd,cht);
double y = py + ph - cht - 15*thickness;
Dpoint_t p1,p2;
p1.x=x1;
p2.x=x2;
p1.y=p2.y=y;
AspUtil::drawLine(p1,p2,roiColor,thick);
p1.x=p2.x=x1;
p2.y=y-5*thickness;
AspUtil::drawLine(p1,p2,roiColor,thick);
p1.x=p2.x=x2;
p2.y=y-5*thickness;
AspUtil::drawLine(p1,p2,roiColor,thick);
p1.x=p2.x=0.5*(x1+x2);
p2.y=y+5*thickness;
AspUtil::drawLine(p1,p2,roiColor,thick);
char fontName[64];
getOptName(INTEG_NUM,fontName);
y = py + ph - cht/2;
AspUtil::drawString((char *)labelStr.c_str(), (int)p1.x-cwd/2, (int)y, -1, fontName,0);
}
getOptName(SPEC_LINE_MIN,thickName);
set_spectrum_width(1);
set_line_width(1);
set_transparency_level(0);
}
int dpf(int argc, char *argv[], int retc, char *retv[])
/*************/
{ int groups;
int noisemult_p = 0;
int numlines,maxlines;
int i;
int donll = TRUE;
int noneg = FALSE;
int top = FALSE;
char command[128];
char name[64];
(void) retc;
(void) retv;
Wturnoff_buttons();
leader_len = 20.0;
axish = FALSE;
axisp = FALSE;
axisp_freq = 1.0;
if (argc > 1)
{
for (i=1; i<argc; i++)
{
if (strcmp(argv[i],"off") == 0) {
set_dpf_flag(0,"");
redo_dpf_dpir();
RETURN;
} else if (strcmp(argv[i],"turnoff") == 0) {
set_dpf_flag(0,"");
RETURN;
} else if (strcmp(argv[i],"noll") == 0) donll = FALSE;
else if (strcmp(argv[i],"pos") == 0) noneg = TRUE;
else if (strcmp(argv[i],"top") == 0) top = TRUE;
else if (strcmp(argv[i],"axish") == 0) axish = TRUE;
else if (strcmp(argv[i],"axisp") == 0)
{
axisp = TRUE;
P_getreal(PROCESSED,"sfrq",&axisp_freq,1);
}
else if (strcmp(argv[i],"leader") == 0)
{
if (((i+1) < argc) && isReal(argv[i+1]) )
{
i++;
leader_len = stringReal(argv[i]);
}
}
else if (isReal(argv[i])) noisemult_p = i;
}
}
if (donll) /* if not donll, dpf uses last previous line listing */
{
if (noisemult_p != 0)
{
if (noneg)
strcpy(command,"nll('dpf','pos',");
else
strcpy(command,"nll('dpf',");
strcat(command,argv[noisemult_p]);
strcat(command,")\n");
}
else
{
if (noneg)
strcpy(command,"nll('dpf','pos')\n");
else
strcpy(command,"nll('dpf')\n");
}
execString(command);
}
if(P_getreal(CURRENT, "dpf_sc2",&dpf_sc2,1)) dpf_sc2=sc2;
if(P_getreal(CURRENT, "dpf_wc2",&dpf_wc2,1)) dpf_wc2=wc2;
if(dpf_wc2>wc2) dpf_wc2=wc2;
if(dpf_sc2<sc2) dpf_sc2=sc2;
/* if (init2d(1,1)) return(ERROR); */
scale = vs;
if (normflag)
scale *= normalize;
/* Wscrprintf("normflag=%d normalize=%g\n",normflag,normalize); */
plot = (argv[0][0] == 'p');
/*select_init(get_rev, dis_setup, fdimname, doheaders, docheck2d, dospecpars,
doblockpars, dophasefile)*/
/*if (init2d(0,plot + 1)) return(ERROR); */
if (select_init(
0,
plot+1,
NO_FREQ_DIM,
NO_HEADERS,
DO_CHECK2D,
DO_SPECPARS,
NO_BLOCKPARS,
NO_PHASEFILE
))
return(ERROR);
if ((numlines = getlines(noneg)) == 0) RETURN;
if (!plot)
{
setwindows();
dispcalib = (float) (mnumypnts-ymin) / (float) wc2max;
getOptName(PEAK_MARK,name);
set_line_thickness(name);
getOptName(PEAK_NUM,name);
set_graphics_font(name);
}
else
{
double size = G_getCharSize("PeakNum");
charsize(size);
//charsize((double)0.7);
dispcalib = ppmm / ymultiplier;
}
CharSize = ycharpixels + ycharpixels / 5;
maxlines = dnpnt / CharSize;
if (maxlines < numlines)
remove_lines(maxlines,&numlines);
if (numlines > 0)
{
//color(PARAM_COLOR);
color(PEAK_MARK_COLOR);
groups = groupcheck(numlines);
if (debug1)
{
Wscrprintf("index high low newhigh newlow highlimit lowlimit\n");
for (i = 1; i <= groups; i++)
Wscrprintf("%d %d %d %d %d %d %d\n",i,ga[i].high,ga[i].low,
ga[i].newhigh,ga[i].newlow,ga[i].highlimit,ga[i].lowlimit);
}
label_proc = (top) ? (PFV) label_top : (PFV) label_bot;
for (i = 1; i <= groups; i++)
label_group(i);
if (!plot) { // construct the command, and set the flag and command
char cmd[64];
if(argc > 1) {
if (isReal(argv[1]) ) sprintf(cmd,"%s(%s",argv[0],argv[1]);
else sprintf(cmd,"%s('%s'",argv[0],argv[1]);
for(i=2;i<argc;i++)
{
if (isReal(argv[i]) )
{
strcat(cmd,",");
strcat(cmd,argv[i]);
}
else
{
strcat(cmd,",'");
strcat(cmd,argv[i]);
strcat(cmd,"'");
}
}
strcat(cmd,")\n");
} else sprintf(cmd,"%s\n",argv[0]);
Wsetgraphicsdisplay("ds");
set_dpf_flag(1, cmd);
}
}
releaseAllWithId("dpf");
endgraphics();
disp_status(" ");
set_graphics_font("Default");
RETURN;
}
