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; }