Geting Scalars to DiscreteMarchingCubes

Development
6

Hello,

I was able to acess only point data with discretemarchingcubes and sufracenets3d through this code. The labels are correct.

CELL DATA CODE:

vtkCellData* cellData = polydataCopy->GetCellData();

if (!cellData) {std::cerr << "Error: Cell data is nullptr." << std::endl;}
 else{   std::cout << "Number of cell arrays: " << cellData->GetNumberOfArrays() << std::endl;}


vtkInformation* outputInfo = discrete_smooth_mesher->GetOutputPortInformation(0);
if (outputInfo->Has(vtkDataObject::CELL_DATA_VECTOR())) {
    vtkCellData* cellData = polydataCopy->GetCellData();} 
    else {std::cerr << "Error: No cell data available in the output." << std::endl;}

	

    vtkDataArray* scalarArray = cellData->GetArray(0);
    if(!scalarArray)vtkDataArray* scalarArray = cellData->GetArray("ImageScalars");
     

    if (scalarArray) {
    		    std::unordered_set<int> uniqueValues;

		    for (vtkIdType i = 0; i < scalarArray->GetNumberOfTuples(); ++i) {
     			double scalarValue = scalarArray->GetTuple1(i);
			int scalarValue_2 = static_cast<int>(scalarValue);
			
			uniqueValues.insert(scalarValue_2);}
		    

	    
	                // Print unique values
		    for (int value : uniqueValues) {
			std::cout << value << " ";
		    }


	}else {    std::cerr << "Error: Scalar array is nullptr." << std::endl;}

POINT DATA CODE:

	vtkPointData* pointData = polydataCopy->GetPointData();

	if (!pointData) {std::cerr << "Error: Point data is nullptr." << std::endl;
	} else {std::cout << "Number of point arrays: " << pointData->GetNumberOfArrays() << std::endl;

	    for (vtkIdType i = 0; i < pointData->GetNumberOfArrays(); ++i) {
		vtkDataArray* array = pointData->GetArray(i);
		const char* arrayName = array->GetName();
		std::cout << "Array " << i << ": " << (arrayName ? arrayName : "NoName") << std::endl;
	    }

	    vtkDataArray* scalarArray = pointData->GetScalars();

	    if (scalarArray) {
		// Assuming the scalar array is of type vtkDoubleArray
		std::unordered_set<int> uniqueValues;

		for (vtkIdType i = 0; i < scalarArray->GetNumberOfTuples(); ++i) {
		    int scalarValue = static_cast<int>(scalarArray->GetTuple1(i));
		    uniqueValues.insert(scalarValue);}

		// Print unique values
		for (int value : uniqueValues) {
		    std::cout << value << " ";}
	    } else {std::cerr << "Error: Scalar array is nullptr." << std::endl;}}

TERMINAL:

	Number of cell arrays: 0
	Error: No cell data available in the output.
	Error: Scalar array is nullptr.
	Number of point arrays: 2
	Array 0: scalars
	Array 1: Normals
	46 45 44 43 42 41 40 39 38 37 20 19 18 17 16 15 14 13 12 11 10 9 8 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36