casa::ImageRegion* ImageRegionGenerator::_makeRegionPolygon( casa::ImageInterface<casa::Float> * image,
Carta::Lib::Regions::Polygon* polygon ){
casa::ImageRegion* imageRegion = nullptr;
if ( image && polygon ){
QPolygonF poly = polygon->qpolyf();
int cornerPointCount = poly.count();
casa::Vector<casa::Double> x( cornerPointCount );
casa::Vector<casa::Double> y( cornerPointCount );
const casa::CoordinateSystem &cs = image->coordinates( );
bool successful = true;
for ( int i = 0; i < cornerPointCount; ++i ) {
QPointF polyPt = poly.value( i );
casa::Vector<casa::Double> worldPt = _toWorld( cs, polyPt.x(), polyPt.y(), &successful );
if ( successful ){
x[i] = worldPt[0];
y[i] = worldPt[1];
}
else {
qDebug() << "i="<<i<<" could not convert x="<<polyPt.x()<<" y="<<polyPt.y()<<" to world.";
break;
}
}
if ( successful ){
casa::Vector<casa::Int> dispAxes(2);
const casa::CoordinateSystem &cs = image->coordinates( );
int directionIndex = cs.findCoordinate( casa::Coordinate::DIRECTION );
if ( directionIndex >= 0 ){
casa::Vector<casa::Int> dirPixelAxis = cs.pixelAxes(directionIndex);
dispAxes(0) = dirPixelAxis[0];
dispAxes(1) = dirPixelAxis[1];
const casa::String units( RAD_UNITS.toStdString().c_str() );
try {
casa::Quantum<casa::Vector<casa::Double> > qx( x, units );
casa::Quantum<casa::Vector<casa::Double> > qy( y, units );
casa::WCPolygon poly(qx, qy, casa::IPosition(dispAxes), cs);
imageRegion = new casa::ImageRegion(poly);
}
catch( casa::AipsError& error ) {
qDebug() << "Error making image region: "<<error.getMesg().c_str();
}
}
}
}
return imageRegion;
}
