7
$\begingroup$

I was experimenting with the code from this question when I ran into another problem with regions.

Ω = RegionDifference[Rectangle[{0, 0}, {10, 10}], Rectangle[{4, 4}, {8, 8}]];
RegionPlot[Ω]

plot

Ω1 = TransformedRegion[Ω, RotationTransform[45 °, {5, 5}]];
RegionQ[Ω1]

True

RegionPlot[Ω1]

RegionPlot::invplotreg: TransformedRegion[RegionDifference[Rectangle[{0, 0}, {10, 10}], Rectangle[{4, 4}, {8, 8}]], TransformationFunction[...]] is not a valid region to plot. >>

What is the difference between a "valid region to plot" and a region that satisfies RegionQ? Or, perhaps, to put it better, am I seeing a bug in RegionPlot or just an incomplete implementation?

I note that 

RegionPlot[TransformedRegion[Rectangle[], RotationTransform[45 °, {.5, .5}]], 
  PlotRange -> All]

works as expected, so it would seem RegionPlot can handle rotations for some class of inputs.

Improve this question
$\endgroup$
7
  • $\begingroup$ I can reproduce this. Might try tracing it some time later (for testing my tracer :) $\endgroup$
    – Silvia
    Commented Jul 19, 2014 at 12:27
  • 3
    $\begingroup$ For me even the first basic example in the TransformedRegion documentation doesn't work properly. And the second example hangs the kernel... $\endgroup$
    – Simon Woods
    Commented Jul 19, 2014 at 13:24
  • $\begingroup$ @SimonWoods For the first example, looks like it can use a PlotRange -> All. (+1) $\endgroup$
    – Silvia
    Commented Jul 19, 2014 at 16:47
  • 1
    $\begingroup$ @SimonWoods, if you could send this and your OS version to the support that would be useful to track it down. Thanks. I have filed fact that an expr with Head TransformedRegion does not work as a suggestion for improvement. $\endgroup$
    – user21
    Commented Jul 21, 2014 at 8:51
  • 1
    $\begingroup$ @user21. re: "an expr with Head TransformedRegion does not work." I don't think that is the issue since TransformedRegion[Rectangle[], RotationTransform[45 °, {.5, .5}]] has head TransformedRegion, but can be plotted by RegionPlot. $\endgroup$
    – m_goldberg
    Commented Jul 21, 2014 at 9:07

1 Answer 1

Reset to default
6
$\begingroup$

Here is a workaround:

r1 = RegionDifference[Rectangle[{0, 0}, {10, 10}], 
   Rectangle[{4, 4}, {8, 8}]];
r2 = TransformedRegion[r1, RotationTransform[45 \[Degree], {5, 5}]];
mr = DiscretizeRegion[r2]

enter image description here

And then:

RegionPlot[mr]

enter image description here

Improve this answer
$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.