comparefits.c

图像处理的压缩算法

	wks.SetCell(iRow, 1, sStr1);
	for(iLoop = 2; iLoop <= 5; iLoop++)
	{
		wks.SetCell(iRow, iLoop, " ");
	}
	iRow++;

// Ruler
	wks.SetCell(iRow, 0, " ");
	for(iCol = 1;iCol <= 5;iCol++)
		wks.SetCell(iRow, iCol, "---");
	iRow++;

// Blank row
	for(iCol = 0;iCol <= 5;iCol++)
		wks.SetCell(iRow, iCol, " ");
	iRow++;

// Second Fit
// Indent and type a Fit results message
	sStr1.Format(FC_HEADER4, strCurve2);
	wks.SetCell(iRow, 0, " ");
	wks.SetCell(iRow, 1, sStr1);
	for(iCol=2 ; iCol <= 5; iCol++)
		wks.SetCell(iRow,iCol," ");
	iRow++;

// Ruler
	wks.SetCell(iRow, 0, " ");
	for(iCol = 1;iCol <= 5;iCol++)
		wks.SetCell(iRow, iCol, "---");
	iRow++;

// Headers
	wks.SetCell(iRow,0," ");
	sStr1 = FC_HEADER5;
	for(iCol = 1; iCol <= 5; iCol++)
	{
		sStr2 = sStr1.GetToken(iCol - 1, ' ');
		wks.SetCell(iRow, iCol, sStr2);
	}
	iRow++;
// Parameters
	for(iLoop = 1,iIdx = 4;iLoop <= (int) BData[3];iLoop++,iRow++)
	{
		wks.SetCell(iRow,0," ");
		bRes = BData.GetText(iIdx, sStr1);
		iIdx++;
		wks.SetCell(iRow, 1, sStr1);
		dTemp1 = BData[iIdx];
		iIdx++;
		wks.SetCell(iRow, 2, dTemp1);
		dTemp2 = BData[iIdx];
		iIdx++;
		wks.SetCell(iRow, 3, dTemp2);
		dTemp = dTemp2 * tTable(dFitCmpCL, BData[0]);
		wks.SetCell(iRow, 4, dTemp1 - dTemp);
		wks.SetCell(iRow, 5, dTemp1 + dTemp);
	}

// Blank row
	for(iCol = 0;iCol <= 5;iCol++)
		wks.SetCell(iRow, iCol, " ");
	iRow++;

// Special quantities
	wks.SetCell(iRow, 0, " ");
	sStr1 = "SSR = ";
	if(BData[1] == NANUM) sStr2 = "--"; else sStr2 = LocalizeDouble(BData[1], "%G");
	/// ML 2/6/2002 PROPER_CONCATENATION
	//lstrcat(sStr1,sStr2);
	sStr1 += sStr2;
	/// end PROPER_CONCATENATION
//	sStr1.Format("SSR = %f", BData[1]);
	wks.SetCell(iRow, 1, sStr1);
	for(iLoop = 2; iLoop <= 5; iLoop++)
	{
		wks.SetCell(iRow, iLoop, " ");
	}
	iRow++;
	wks.SetCell(iRow, 0, " ");
	sStr1 = "DOF = ";
	if(BData[0] == NANUM) sStr2 = "--"; else sStr2 = LocalizeDouble(BData[0], "%G");
	/// ML 2/6/2002 PROPER_CONCATENATION
	//lstrcat(sStr1,sStr2);
	sStr1 += sStr2;
	/// end PROPER_CONCATENATION
//	sStr1.Format("DOF = %f", BData[0]);
	wks.SetCell(iRow, 1, sStr1);
	for(iLoop = 2; iLoop <= 5; iLoop++)
	{
		wks.SetCell(iRow, iLoop, " ");
	}
	iRow++;
	wks.SetCell(iRow, 0, " ");
	sStr1 = "ChiSqr / DOF = ";
	if(BData[2] == NANUM) sStr2 = "--"; else sStr2 = LocalizeDouble(BData[2], "%G");
	/// ML 2/6/2002 PROPER_CONCATENATION
	//lstrcat(sStr1,sStr2);
	sStr1 += sStr2;
	/// end PROPER_CONCATENATION
//	sStr1.Format("ChiSqr / DOF = %f", BData[2] / BData[0]);
	wks.SetCell(iRow, 1, sStr1);
	for(iLoop = 2; iLoop <= 5; iLoop++)
	{
		wks.SetCell(iRow, iLoop, " ");
	}
	iRow++;

// Ruler
	wks.SetCell(iRow, 0, " ");
	for(iCol = 1;iCol <= 5;iCol++)
		wks.SetCell(iRow, iCol, "---");
	iRow++;

// Blank row
	for(iCol = 0;iCol <= 5;iCol++)
		wks.SetCell(iRow, iCol, " ");
	iRow++;

// F-Test Results
	wks.SetCell(iRow, 0, " ");
	wks.SetCell(iRow, 1, FC_HEADER6);
	for(iLoop = 2; iLoop <= 5; iLoop++)
	{
		wks.SetCell(iRow, iLoop, " ");
	}
	iRow++;

// Ruler
	wks.SetCell(iRow, 0, " ");
	for(iCol = 1;iCol <= 5;iCol++)
		wks.SetCell(iRow, iCol, "---");
	iRow++;

// F and P Values
	wks.SetCell(iRow, 0, " ");
	sStr1 = FC_HEADER5B;
	if(dfval == NANUM) sStr2 = "NAN"; else sStr2 = LocalizeDouble(dfval, "%G");
	/// ML 2/6/2002 PROPER_CONCATENATION
	//lstrcat(sStr1,sStr2);
	sStr1 += sStr2;
	/// end PROPER_CONCATENATION
	sStr2 = ", P = ";
	/// ML 2/6/2002 PROPER_CONCATENATION
	//lstrcat(sStr1,sStr2);
	sStr1 += sStr2;
	/// end PROPER_CONCATENATION
	if(dpval == NANUM) sStr2 = "NAN"; else sStr2 = LocalizeDouble(dpval, "%G");
	/// ML 2/6/2002 PROPER_CONCATENATION
	//lstrcat(sStr1,sStr2);
	sStr1 += sStr2;
	/// end PROPER_CONCATENATION
//	sStr1.Format("F-Value = %f, P = %f", dfval, dpval);
	wks.SetCell(iRow, 1, sStr1);
	for(iLoop = 2; iLoop <= 5; iLoop++)
	{
		wks.SetCell(iRow, iLoop, " ");
	}
	iRow++;

// Ruler
	wks.SetCell(iRow, 0, " ");
	for(iCol = 1;iCol <= 5;iCol++)
		wks.SetCell(iRow, iCol, "---");
	iRow++;

// Conclusion
	LT_get_var("OFitCmpCL", &dFitCmpCL);
	dTemp = 1.0 - dFitCmpCL / 100.0;
	if(dpval > dTemp)
	{
		sStr1.Format(FC_FOOTER1, dTemp);
	}
	else
	{
		sStr1.Format(FC_FOOTER2, dTemp);
	}
	wks.SetCell(iRow, 0, " ");
	wks.SetCell(iRow, 1, sStr1);
	for(iLoop = 2; iLoop <= 5; iLoop++)
	{
		wks.SetCell(iRow, iLoop, " ");
	}
	iRow++;

// Blank row
	for(iCol = 0;iCol <= 5;iCol++)
		wks.SetCell(iRow, iCol, " ");
	iRow++;

// Line Fit comparison
	if(strFunc.CompareNoCase("Line") == 0)
	{
		wks.SetCell(iRow, 0, " ");
		wks.SetCell(iRow, 1, FC_LINEARHEADER1);
		for(iCol = 2;iCol <= 5;iCol++)
			wks.SetCell(iRow, iCol, " ");
		iRow++;
		// Ruler
		wks.SetCell(iRow, 0, " ");
		for(iCol = 1;iCol <= 5;iCol++)
			wks.SetCell(iRow, iCol, "---");
		iRow++;
		// Intercepts
		dTemp1 = AData[6] * AData[6] / AData[2];
		dTemp2 = BData[6] * BData[6] / BData[2];
		dTemp3 = (AData[0] * AData[2] + BData[0] * BData[2]) / (AData[0] + BData[0]);
		LT_get_var("FitCmpCL", &dFitCmpCL);
		dTemp = tTable(dFitCmpCL, AData[0] + BData[0]) * sqrt(dTemp3) * sqrt(dTemp1 + dTemp2);
		dTemp1 = AData[5] - BData[5] - dTemp;
		dTemp2 = AData[5] - BData[5] + dTemp;
		dTemp3 = AData[8] - BData[8] - dTemp;
		dTemp4 = AData[8] - BData[8] + dTemp;
		sStr1.Format(FC_LINEARHEADER2, dTemp1, dTemp2);
		sStr2.Format(" %f", NANUM); // GRD - Replace missing value with dashes
		sStr1.Replace(sStr2," --");
		sStr1.Replace(sStr2," --");
		wks.SetCell(iRow, 0, " ");
		wks.SetCell(iRow, 1, sStr1);
		for(iCol = 2;iCol <= 5;iCol++)
			wks.SetCell(iRow, iCol, " ");
		iRow++;
		wks.SetCell(iRow, 0, " ");
		if(dTemp1 > dTemp2)
		{
			dTemp = dTemp1;
			dTemp1 = dTemp2;
			dTemp2 = dTemp;
		}
		LT_get_var("OFitCmpCL", &dFitCmpCL);
		dFitCmpCL = 1.0 - dFitCmpCL / 100.0;
		// GRD - Handle missing values
		if((dTemp1 <= 0 && dTemp2 >= 0) || (dTemp1 == NANUM || dTemp2 == NANUM))
		{
			sStr1.Format(FC_LINEARFOOTER1A, dFitCmpCL);
		}
		else
		{
			sStr1.Format(FC_LINEARFOOTER1B, dFitCmpCL);
		}
		wks.SetCell(iRow, 1, sStr1);
		for(iCol = 2;iCol <= 5;iCol++)
			wks.SetCell(iRow, iCol, " ");
		iRow++;

		// Blank row
		for(iCol = 0;iCol <= 5;iCol++)
			wks.SetCell(iRow, iCol, " ");
		iRow++;

		// Slopes
		sStr1.Format(FC_LINEARHEADER3, dTemp3, dTemp4);
		sStr2.Format(" %f", NANUM); // GRD - Replace missing value with dashes
		sStr1.Replace(sStr2," --");
		sStr1.Replace(sStr2," --");
		wks.SetCell(iRow, 0, " ");
		wks.SetCell(iRow, 1, sStr1);
		for(iCol = 2;iCol <= 5;iCol++)
			wks.SetCell(iRow, iCol, " ");
		iRow++;
		wks.SetCell(iRow, 0, " ");
		if(dTemp3 > dTemp4)
		{
			dTemp = dTemp2;
			dTemp3 = dTemp4;
			dTemp4 = dTemp;
		}
		// GRD - Handle missing values
		if((dTemp3 <= 0 && dTemp4 >= 0) || (dTemp3 == NANUM || dTemp4 == NANUM))
		{
			sStr1.Format(FC_LINEARFOOTER2A, dFitCmpCL);
		}
		else
		{
			sStr1.Format(FC_LINEARFOOTER2B, dFitCmpCL);
		}
		wks.SetCell(iRow, 1, sStr1);
		for(iCol = 2;iCol <= 5;iCol++)
			wks.SetCell(iRow, iCol, " ");
		iRow++;

	}

// Blank row
	for(iCol = 0;iCol <= 5;iCol++)
		wks.SetCell(iRow, iCol, " ");
	iRow++;

	LT_execute("type.wks(%M);type;type.EndResults();");
}