ocmath_find_markers_of_peak_centers_by_height_and_width

 

Description

Find feets of the peaks by user settings

Syntax

int ocmath_find_markers_of_peak_centers_by_height_and_width( UINT lSize, const double * px, const double * py, int nPeakNum, const double * pPeakX, int nMarkersNum, int * pMarkersIndices, double dMarkerHeight, double dMaxWidth, double dMinWidth )

Parameters

lSize
[input] size of px, py
px
[input] it contains curve's X coordinate's datas
py
[input] it contains curve's Y coordinate's datas
nPeakNum
[input] size of pPeakIndices
pPeakX
[input] X coordinates of peak centers
nMarkersNum
[input] size of pMarkersIndices, its initial value should larger than 2*nPeakNum
pMarkersIndices
[output] the indices of peaks' feets, its initial size is nMarkersNum, if nMarkersNum
is more than 2*nPeakNum, fill the excrescent elements with NANUM after pMarkersIndices[2*nPeakNum-1]
dMarkerHeight
[input] the estimated feet height of peaks, which will be marked
dMaxWidth
[input] the maximum width of peaks, which will be marked
dMinWidth
[input] the minimum width of peaks, which will be marked

Return

Return OE_NOERROR if succeed, otherwise, non-zero error code is returned.

Examples

EX1

//Assume in the current graph, curve's XY data is in the first data plot. This piece
//of code get the XY data of the curve from the first data plot and find the feets of
//the peaks. The result is output in a new worksheet and the feets will plot in the 
//original data plot.
void ocmath_find_markers_of_peak_centers_by_height_and_width_ex1( )
{
    GraphLayer gl = Project.ActiveLayer();
    if (!gl)
    {
        return;
    }
 
    //get plot data of curve from the first dataplot
    DataPlot dp = gl.DataPlots(0);        
    DataRange dr;
    vector vxData, vyData;
    if(dp.GetDataRange(dr))
    {
        DWORD dwPlotID;
        if(dr.GetData(DRR_GET_DEPENDENT | DRR_NO_FACTORS, 0, &dwPlotID, NULL, &vyData, &vxData) < 0)
        {
            printf("get_plot_data failed GetData");
            return;
        }
    }
 
    uint nDataSize = vxData.GetSize();
 
    vector vxPeaks, vyPeaks;
    vector<int> vnIndices;
 
    vxPeaks.SetSize(nDataSize);
    vyPeaks.SetSize(nDataSize);
    vnIndices.SetSize(nDataSize);
 
    ocmath_find_peaks_1st_derivative( &nDataSize, vxData, vyData, vxPeaks, vyPeaks, vnIndices, POSITIVE_DIRECTION | NEGATIVE_DIRECTION,3);
    vxPeaks.SetSize(nDataSize);
    vyPeaks.SetSize(nDataSize);
    vnIndices.SetSize(nDataSize);
    int lSize = vxData.GetSize();
    int peakNum=vnIndices.GetSize();
    double dMarkerHeight = 30;
    double dMaxWidth=30;
    double dMinWidth=1;
    vector<int> markeIndices;
    markeIndices.SetSize(2*peakNum);
    int nRet = ocmath_find_markers_of_peak_centers_by_height_and_width(vxData.GetSize(), vxData, vyData, peakNum, 
                                                        vxPeaks, 2*peakNum, markeIndices,  
                                                        dMarkerHeight, dMaxWidth, dMinWidth);
    if( nRet < OE_NOERROR )
    {
        printf("error code: %d\n", nRet);
        return;
    }
    
    vector vxMarker, vyMarker;
    int nSize = markeIndices.GetSize();
    vxMarker.SetSize(nSize);
    vyMarker.SetSize(nSize);
    for(int ii = 0; ii < markeIndices.GetSize(); ii++)
    {
        if(markeIndices[ii] >= 0)
        {
            vxMarker[ii] = vxData[markeIndices[ii]];
            vyMarker[ii] = vyData[markeIndices[ii]];
        }
        else
        {
            vxMarker[ii] = NANUM;
            vyMarker[ii] = NANUM;
        }        
    }
    vxMarker.Trim();
    vyMarker.Trim();
    
    //new a worksheet to show the result
    WorksheetPage wksPage;
    wksPage.Create();
    Worksheet wksResult = wksPage.Layers(0);
    //wksResult.AddCol("Indices");
    int nXCol, nYCol;
    nXCol = wksResult.AddCol("Marker X");
    nYCol = wksResult.AddCol("Marker Y");
    //wksResult.Columns(0).SetType(OKDATAOBJ_DESIGNATION_X);
    wksResult.Columns(1).SetType(OKDATAOBJ_DESIGNATION_X);
    wksResult.Columns(2).SetType(OKDATAOBJ_DESIGNATION_Y);
    
    //plot the feets in the original data plot
    DataRange drOut;
    drOut.Add("X", wksResult, 0, nXCol, -1, nXCol);
    drOut.Add("Y", wksResult, 0, nYCol, -1, nYCol);
    drOut.SetData(&vyMarker, &vxMarker);
    gl.AddPlot(drOut, IDM_PLOT_SCATTER);
}

Remark

Find feets of the peaks, whose height is greater than dMarkerHeight(means abs of the points' y values are greater than dMarkerHeight) and width is between dMinWidth and dMaxWidth, then record peaks' left and right feets in pMarkersIndices.

pMarkersIndices contains the indices of peaks' feets.

See Also

ocmath_find_peaks_1st_derivative, ocmath_find_peaks_2nd_derivative

header to Include

origin.h

Reference