One thought was that these were farm fields in satellite images. I just had no
clue why one would subdivide them. Maybe a training and testing portion of
these polygons? At a few dozen meters the curvature might be negligible on a
sufficiently large spheroid. If 45-55% split is approximate then go through the
entire polygon area pixel by pixel and use rnorm() to generate a random value
that puts the pixel in A if <0.4 and in B it >0.4. This has a better chance of
working if the polygon is centered in the image so that curvature weighs
equally.
For this post my thought was that it was either homework because the question
was so simple, or it was not providing enough information to answer the real
question. The generalized question is to divide an n-sided polygon mapped to
the surface of a spheroid with semi-axes A, B, C into 45-55% split into two
parts 45 and 55%. If multiple polygons were spread over the image and the image
curvature was large relative to the image size I don't see how you could
recover polygons close to the edge of the spheroid unless you knew beforehand
the curvature at that point. Possibly in the most general sense the problem is
unsolvable because angles near the horizon will be nearly straight lines.
Another part of asking questions is to provide enough detail that others may
arrive at a creative answer. How about a mechanical solution. Print the polygon
onto paper. Cut the polygon out. Cut the polygon in half. Weigh the halves to
get actual split. Scan each piece in. You can print several copies and repeat
until you get a 45-55 split as close as your equipment will measure. I have
solved the problem as asked, though it did not involve R. I split the polygon
problem solved or do I need to do something with the pieces?
I had posted a simple R solution, but was then asked if the solution would work
for a trapezoid with 21 bottom and 18 top. My solution was not generalizable
but the question was for a rectangle (four 90 degree angles and parallel
opposite sides) with sides 18 and 200 meters.
You could draw polygons on a balloon using a sharpie marker (or similar).
Photograph the balloon at different distances, and process the images. At least
this way you have a system to test that part of the program. Correct pixel area
for curvature, estimate area of polygon, and then figure out what it means to
split the polygon. If the non-regular polygon is not evenly divisible 45-55
what happens to the remainder?
Is this problem better handled using GIS techniques?
Tim
-----Original Message-----
From: R-help <r-help-boun...@r-project.org> On Behalf Of Avi Gross via R-help
Sent: Wednesday, April 27, 2022 11:27 AM
To: r-help@r-project.org
Subject: Re: [R] Convert one rectangle into two rectangles
[External Email]
Just FYI, Jim, I was sent private mail by Javad that turns the problem aroundso
not only are there no rectangles, but the problem is not 2-D. He is working
with Spatial Polygons representing areas on the surface of a sphere (presumably
Earth) and wants to subdivide them. This is a VERY different topic and there
are packages and other references that might apply to his needs. I have no idea
why the 60/40 split by area.
Yes, as a very simplified idea, I understand why he proposed dividing a
rectangle proportionately but realistically what he is working with is a bit
more like a trapezoid which is also bent into a third dimension.
So those of us wanting to help on the original problem can stop and, speaking
for myself, I am going to approach many questions posed here more carefully to
see if they are well-thought-out or are some kind of fishing expedition. And,
since it has been made very clear multiple times that the scope of this forum
is narrow and not meant to help with HW, and I have no idea how to verify it is
not, ...
-----Original Message-----
From: Jim Lemon <drjimle...@gmail.com>
To: javad bayat <j.bayat...@gmail.com>; r-help mailing list
<r-help@r-project.org>
Sent: Wed, Apr 27, 2022 5:59 am
Subject: Re: [R] Convert one rectangle into two rectangles
Hi javad,
Let's think about it a little. There are only two ways you can divide a
rectangle into two other rectangles.
The dividing line must be straight and parallel to one side or the other or you
won't get two rectangles.
If there are no other constraints, you can take your pick.
A E B
-------------------------------------
| | |
-------------------------------------
C F D
AE = CF = 0.45 * AB
Alternatively the same equation can be used on the shorter sides.
If you want to cut it across the length, you can do it with an abacus.
While the IrregLong package can produce an abacus plot, I was unable to find an
R-based abacus.
For that I suggest:
https://urldefense.proofpoint.com/v2/url?u=https-3A__toytheater.com_abacus_&d=DwIFaQ&c=sJ6xIWYx-zLMB3EPkvcnVg&r=9PEhQh2kVeAsRzsn7AkP-g&m=Tbmiv3MGuN8f7gYSHUOKBKkBAx5goWhCbuyiJMDVP2c3HV-_FfjZ8__7NHiO0K9V&s=Fz2RhXtevzROx9IOzX2dn3FYg972J9cnuWoWW7MqZBo&e=
Jim
On Wed, Apr 27, 2022 at 1:42 AM javad bayat <j.bayat...@gmail.com> wrote:
>
> Dear R users;
> I want to convert a rectangle into 2 new rectangles. I have a
> rectangle (18 m * 200 m), and I want to subtract 45 percent of it's
> area as a new rectangle. Is it possible to do this in R?
> Sincerely.
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://urldefense.proofpoint.com/v2/url?u=https-3A__stat.ethz.ch_mail
> man_listinfo_r-2Dhelp&d=DwIFaQ&c=sJ6xIWYx-zLMB3EPkvcnVg&r=9PEhQh2kVeAs
> Rzsn7AkP-g&m=Tbmiv3MGuN8f7gYSHUOKBKkBAx5goWhCbuyiJMDVP2c3HV-_FfjZ8__7N
> HiO0K9V&s=AZ_JuNx5fiOamGA0k4TSWI5sC8DZZdomay2LjrgyMtM&e=
> PLEASE do read the posting guide
> https://urldefense.proofpoint.com/v2/url?u=http-3A__www.R-2Dproject.or
> g_posting-2Dguide.html&d=DwIFaQ&c=sJ6xIWYx-zLMB3EPkvcnVg&r=9PEhQh2kVeA
> sRzsn7AkP-g&m=Tbmiv3MGuN8f7gYSHUOKBKkBAx5goWhCbuyiJMDVP2c3HV-_FfjZ8__7
> NHiO0K9V&s=VPNX6IvYr7F4Vti-KkGsoS-sxBZHfX4T6JTrzo7_aOw&e=
> and provide commented, minimal, self-contained, reproducible code.
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://urldefense.proofpoint.com/v2/url?u=https-3A__stat.ethz.ch_mailman_listinfo_r-2Dhelp&d=DwIFaQ&c=sJ6xIWYx-zLMB3EPkvcnVg&r=9PEhQh2kVeAsRzsn7AkP-g&m=Tbmiv3MGuN8f7gYSHUOKBKkBAx5goWhCbuyiJMDVP2c3HV-_FfjZ8__7NHiO0K9V&s=AZ_JuNx5fiOamGA0k4TSWI5sC8DZZdomay2LjrgyMtM&e=
PLEASE do read the posting guide
https://urldefense.proofpoint.com/v2/url?u=http-3A__www.R-2Dproject.org_posting-2Dguide.html&d=DwIFaQ&c=sJ6xIWYx-zLMB3EPkvcnVg&r=9PEhQh2kVeAsRzsn7AkP-g&m=Tbmiv3MGuN8f7gYSHUOKBKkBAx5goWhCbuyiJMDVP2c3HV-_FfjZ8__7NHiO0K9V&s=VPNX6IvYr7F4Vti-KkGsoS-sxBZHfX4T6JTrzo7_aOw&e=
and provide commented, minimal, self-contained, reproducible code.
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://urldefense.proofpoint.com/v2/url?u=https-3A__stat.ethz.ch_mailman_listinfo_r-2Dhelp&d=DwIFaQ&c=sJ6xIWYx-zLMB3EPkvcnVg&r=9PEhQh2kVeAsRzsn7AkP-g&m=Tbmiv3MGuN8f7gYSHUOKBKkBAx5goWhCbuyiJMDVP2c3HV-_FfjZ8__7NHiO0K9V&s=AZ_JuNx5fiOamGA0k4TSWI5sC8DZZdomay2LjrgyMtM&e=
PLEASE do read the posting guide
https://urldefense.proofpoint.com/v2/url?u=http-3A__www.R-2Dproject.org_posting-2Dguide.html&d=DwIFaQ&c=sJ6xIWYx-zLMB3EPkvcnVg&r=9PEhQh2kVeAsRzsn7AkP-g&m=Tbmiv3MGuN8f7gYSHUOKBKkBAx5goWhCbuyiJMDVP2c3HV-_FfjZ8__7NHiO0K9V&s=VPNX6IvYr7F4Vti-KkGsoS-sxBZHfX4T6JTrzo7_aOw&e=
and provide commented, minimal, self-contained, reproducible code.
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
