void Raster::refineRaster(std::size_t scaling)
{
double *new_raster_data(new double[_n_rows*_n_cols*scaling*scaling]);
for (std::size_t row(0); row<_n_rows; row++) {
for (std::size_t col(0); col<_n_cols; col++) {
const std::size_t idx(row*_n_cols+col);
for (std::size_t new_row(row*scaling); new_row<(row+1)*scaling; new_row++) {
const std::size_t idx0(new_row*_n_cols*scaling);
for (std::size_t new_col(col*scaling); new_col<(col+1)*scaling; new_col++) {
new_raster_data[idx0+new_col] = _raster_data[idx];
}
}
}
}
std::swap(_raster_data, new_raster_data);
_cell_size /= scaling;
_n_cols *= scaling;
_n_rows *= scaling;
delete [] new_raster_data;
}
void
TransfMat3x4::normalize_column (Proj::Axis axis) {
Proj::Pt2 new_col(column(axis));
new_col.normalize();
set_image_pt(axis, new_col);
}
