Points within polygons in different projectionsOverlaying layers with different projections on AndroidCorrect coordinate systems, but trustworthy distance between points?How do you clip a gnomonic or stereographic azimuthal projection?Overlay polygons on raster in different projectionsIssue rendering polygons at global scale when using the Robinson coordinate reference system in QGIS 2.10Change coordinate system to add buffer in metersHow to determine if lat/long is within SABS school boundary?Project raster tool stays white at iteration in ModelBuilder?Vizualise different projections of a vector using QGISConversion between unknown map projections?
Can an open source licence be revoked if it violates employer's IP?
Which are the methodologies for interpreting Vedas?
Can an escape pod land on Earth from orbit and not be immediately detected?
Can you open the door or die? v2
Am I being scammed by a sugar daddy?
Do Veracrypt encrypted volumes have any kind of brute force protection?
How to represent jealousy in a cute way?
Why would a car salesman tell me not to get my credit pulled again?
Why are backslashes included in this shell script?
Purpose of cylindrical attachments on Power Transmission towers
Does WiFi affect the quality of images downloaded from the internet?
What do I need to do, tax-wise, for a sudden windfall?
Tiffeneau–Demjanov rearrangement products
Idiom for 'person who gets violent when drunk"
Is it possible to have battery technology that can't be duplicated?
Are athlete's college degrees discounted by employers and graduate school admissions?
Is it good practice to create tables dynamically?
Why did the Death Eaters wait to reopen the Chamber of Secrets?
Someone who is granted access to information but not expected to read it
Why does there seem to be an extreme lack of public trashcans in Taiwan?
Am I allowed to determine tenets of my contract as a warlock?
In Pandemic, why take the extra step of eradicating a disease after you've cured it?
Print "N NE E SE S SW W NW"
Simple log rotation script
Points within polygons in different projections
Overlaying layers with different projections on AndroidCorrect coordinate systems, but trustworthy distance between points?How do you clip a gnomonic or stereographic azimuthal projection?Overlay polygons on raster in different projectionsIssue rendering polygons at global scale when using the Robinson coordinate reference system in QGIS 2.10Change coordinate system to add buffer in metersHow to determine if lat/long is within SABS school boundary?Project raster tool stays white at iteration in ModelBuilder?Vizualise different projections of a vector using QGISConversion between unknown map projections?
.everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty margin-bottom:0;
I have a polygon and a point feature class in some projection. The points are contained within the polygons
The points are always inside the polygons no matter the projection destination I'm projecting points and polygons to (including 'unprojecting' to geographic).
Putting visually JGH example I have a big polygon in WKID =4326 with a point outside of it...
...when I project them (say in 3347 to follow JGH example) I get this...
where the point now is inside the polygon. I used the ArcGIS Pro 'project' function.
When I use project setting the optional parameter 'Preserve Shape' as in the following picture.
...I get something definitely better like...
where (hopefully) the point stays outside the polygon.
'Preserve shape' parameter did what outlined by J.R. Extra vertices that provide a better discretization of the polygon in the destination projection....otherwise everything is reduced to projecting 4 vertices and drawing straight lines in the destination space to connect them (producing the indesiderable output of the second picture above)
Am I getting this right ?
coordinate-system arcgis-pro reprojection-mathematics
add a comment |
I have a polygon and a point feature class in some projection. The points are contained within the polygons
The points are always inside the polygons no matter the projection destination I'm projecting points and polygons to (including 'unprojecting' to geographic).
Putting visually JGH example I have a big polygon in WKID =4326 with a point outside of it...
...when I project them (say in 3347 to follow JGH example) I get this...
where the point now is inside the polygon. I used the ArcGIS Pro 'project' function.
When I use project setting the optional parameter 'Preserve Shape' as in the following picture.
...I get something definitely better like...
where (hopefully) the point stays outside the polygon.
'Preserve shape' parameter did what outlined by J.R. Extra vertices that provide a better discretization of the polygon in the destination projection....otherwise everything is reduced to projecting 4 vertices and drawing straight lines in the destination space to connect them (producing the indesiderable output of the second picture above)
Am I getting this right ?
coordinate-system arcgis-pro reprojection-mathematics
add a comment |
I have a polygon and a point feature class in some projection. The points are contained within the polygons
The points are always inside the polygons no matter the projection destination I'm projecting points and polygons to (including 'unprojecting' to geographic).
Putting visually JGH example I have a big polygon in WKID =4326 with a point outside of it...
...when I project them (say in 3347 to follow JGH example) I get this...
where the point now is inside the polygon. I used the ArcGIS Pro 'project' function.
When I use project setting the optional parameter 'Preserve Shape' as in the following picture.
...I get something definitely better like...
where (hopefully) the point stays outside the polygon.
'Preserve shape' parameter did what outlined by J.R. Extra vertices that provide a better discretization of the polygon in the destination projection....otherwise everything is reduced to projecting 4 vertices and drawing straight lines in the destination space to connect them (producing the indesiderable output of the second picture above)
Am I getting this right ?
coordinate-system arcgis-pro reprojection-mathematics
I have a polygon and a point feature class in some projection. The points are contained within the polygons
The points are always inside the polygons no matter the projection destination I'm projecting points and polygons to (including 'unprojecting' to geographic).
Putting visually JGH example I have a big polygon in WKID =4326 with a point outside of it...
...when I project them (say in 3347 to follow JGH example) I get this...
where the point now is inside the polygon. I used the ArcGIS Pro 'project' function.
When I use project setting the optional parameter 'Preserve Shape' as in the following picture.
...I get something definitely better like...
where (hopefully) the point stays outside the polygon.
'Preserve shape' parameter did what outlined by J.R. Extra vertices that provide a better discretization of the polygon in the destination projection....otherwise everything is reduced to projecting 4 vertices and drawing straight lines in the destination space to connect them (producing the indesiderable output of the second picture above)
Am I getting this right ?
coordinate-system arcgis-pro reprojection-mathematics
coordinate-system arcgis-pro reprojection-mathematics
edited May 30 at 10:03
PolyGeo♦
54.3k1784253
54.3k1784253
asked May 28 at 16:26
SupereshekSupereshek
262
262
add a comment |
add a comment |
2 Answers
2
active
oldest
votes
Some software will use great circle arcs to connect unprojected vertices (sometimes when using a special data type, like PostGIS geography
), while projected (or unprojected using the geometry
datatype) vertices are connected using straight lines.
This can result in a point being inside a polygon expressed as geography
but outside of it if expressed as geometry
The following example uses PostGIS. The polygon goes up to latitude 50, the point is at latitude 51.
WITH poly AS (select ST_GeomFromText('polygon((0 0, 50 0, 50 50, 0 50, 0 0))',4326) as geom),
pnt AS (select ST_GeomFromText('point(25 51)',4326) as geom)
SELECT ST_INTERSECTS(poly.geom,pnt.geom) intersect_geometry,
ST_INTERSECTS(poly.geom::geography,pnt.geom::geography) intersect_geography
FROM poly, pnt;
intersect_geometry | intersect_geography
--------------------+---------------------
f | t
Edit
Extending on @JR comment, here is an example when projecting a huge polygon to a Lambert Conformal Conic projection (3347) and checking the point intersection near the central meridian: we can see a difference of 4 degrees of latitude between the two!
WITH poly AS (select ST_GeomFromText('polygon((-70 45, -130 45, -130 50, -70 50, -70 45))',4326) as geom),
pnt AS (select ST_GeomFromText('point(-100 54)',4326) as geom)
SELECT ST_INTERSECTS(poly.geom,pnt.geom) intersect_geometry,
ST_INTERSECTS(st_transform(poly.geom,3347),st_transform(pnt.geom,3347)) intersect_reproject
FROM poly, pnt;
intersect_geometry | intersect_reproject
--------------------+---------------------
f | t
Same in SQL server...
– Supereshek
May 29 at 9:55
add a comment |
Yes. projections will never reproject a point that was inside a polygon in one projection to be outside it in another--unless there's some sort of precision error. I'm not sure what this property is called in geography, but I just realized it's essentially relativistic invariance, which basically says that as time dilates for us, and our coordinate systems are compressed or stretched, no observer in any given frame will disagree on causal ordering of events. Likewise, no matter in what projection an "observer lives", no one will disagree on what points are in what polygons, even if they disagree on exactly how far the points are from the edges of the polygons, and stuff like that.
Well, except for horizon clipping, and discontinuous projections, which might throw a monkey wench in topology.
– Vince
May 28 at 16:53
3
Not sure but I think the question come from the way polygon may be incorrectly reprojected as in a big polygon defined by only corner vertex. When reprojecting you will project the corner then reconstruct the polygon by linking the corner with a straight line. In this case a point close to the border could be seen on the other side after reprojecting but that's because the projection is wrong along the length of the side (to prevent that you need to add vertex on the polygon side before reprojecting)
– J.R
May 28 at 16:54
I hadn't considered horizons or discontinous projections. I guess like in physics geography does have singularities where certain space-time coordinate systems break down; poles are another example.
– 0mn1
May 28 at 19:37
Yes...as J.R. is saying my question comes from the way a polygon may be incorrectly re-projected. I think ArcGIS 'Project' accounts for this if you use "preserve_shape" optional parameter. It adds extra vertices in order to get a better projected shape.
– Supereshek
May 29 at 9:54
@J.R I see. Is it correct to say that this happens because the PostGIS geometry is a different projection than the geography? If I connect to points on unprojected data with a straight line, then do the same thing with correct geodesic calculations and add that (apparently) curved line to the map, then there will be points between them. But the straight line was never the correct line for unprojected lon-lat data--it's not like the line is shifting from one side of a point to the other between two projections, rather the geometry projection is wrong.
– 0mn1
May 29 at 14:00
|
show 1 more comment
Your Answer
StackExchange.ready(function()
var channelOptions =
tags: "".split(" "),
id: "79"
;
initTagRenderer("".split(" "), "".split(" "), channelOptions);
StackExchange.using("externalEditor", function()
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled)
StackExchange.using("snippets", function()
createEditor();
);
else
createEditor();
);
function createEditor()
StackExchange.prepareEditor(
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: false,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: null,
bindNavPrevention: true,
postfix: "",
imageUploader:
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
,
onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
);
);
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fgis.stackexchange.com%2fquestions%2f324171%2fpoints-within-polygons-in-different-projections%23new-answer', 'question_page');
);
Post as a guest
Required, but never shown
2 Answers
2
active
oldest
votes
2 Answers
2
active
oldest
votes
active
oldest
votes
active
oldest
votes
Some software will use great circle arcs to connect unprojected vertices (sometimes when using a special data type, like PostGIS geography
), while projected (or unprojected using the geometry
datatype) vertices are connected using straight lines.
This can result in a point being inside a polygon expressed as geography
but outside of it if expressed as geometry
The following example uses PostGIS. The polygon goes up to latitude 50, the point is at latitude 51.
WITH poly AS (select ST_GeomFromText('polygon((0 0, 50 0, 50 50, 0 50, 0 0))',4326) as geom),
pnt AS (select ST_GeomFromText('point(25 51)',4326) as geom)
SELECT ST_INTERSECTS(poly.geom,pnt.geom) intersect_geometry,
ST_INTERSECTS(poly.geom::geography,pnt.geom::geography) intersect_geography
FROM poly, pnt;
intersect_geometry | intersect_geography
--------------------+---------------------
f | t
Edit
Extending on @JR comment, here is an example when projecting a huge polygon to a Lambert Conformal Conic projection (3347) and checking the point intersection near the central meridian: we can see a difference of 4 degrees of latitude between the two!
WITH poly AS (select ST_GeomFromText('polygon((-70 45, -130 45, -130 50, -70 50, -70 45))',4326) as geom),
pnt AS (select ST_GeomFromText('point(-100 54)',4326) as geom)
SELECT ST_INTERSECTS(poly.geom,pnt.geom) intersect_geometry,
ST_INTERSECTS(st_transform(poly.geom,3347),st_transform(pnt.geom,3347)) intersect_reproject
FROM poly, pnt;
intersect_geometry | intersect_reproject
--------------------+---------------------
f | t
Same in SQL server...
– Supereshek
May 29 at 9:55
add a comment |
Some software will use great circle arcs to connect unprojected vertices (sometimes when using a special data type, like PostGIS geography
), while projected (or unprojected using the geometry
datatype) vertices are connected using straight lines.
This can result in a point being inside a polygon expressed as geography
but outside of it if expressed as geometry
The following example uses PostGIS. The polygon goes up to latitude 50, the point is at latitude 51.
WITH poly AS (select ST_GeomFromText('polygon((0 0, 50 0, 50 50, 0 50, 0 0))',4326) as geom),
pnt AS (select ST_GeomFromText('point(25 51)',4326) as geom)
SELECT ST_INTERSECTS(poly.geom,pnt.geom) intersect_geometry,
ST_INTERSECTS(poly.geom::geography,pnt.geom::geography) intersect_geography
FROM poly, pnt;
intersect_geometry | intersect_geography
--------------------+---------------------
f | t
Edit
Extending on @JR comment, here is an example when projecting a huge polygon to a Lambert Conformal Conic projection (3347) and checking the point intersection near the central meridian: we can see a difference of 4 degrees of latitude between the two!
WITH poly AS (select ST_GeomFromText('polygon((-70 45, -130 45, -130 50, -70 50, -70 45))',4326) as geom),
pnt AS (select ST_GeomFromText('point(-100 54)',4326) as geom)
SELECT ST_INTERSECTS(poly.geom,pnt.geom) intersect_geometry,
ST_INTERSECTS(st_transform(poly.geom,3347),st_transform(pnt.geom,3347)) intersect_reproject
FROM poly, pnt;
intersect_geometry | intersect_reproject
--------------------+---------------------
f | t
Same in SQL server...
– Supereshek
May 29 at 9:55
add a comment |
Some software will use great circle arcs to connect unprojected vertices (sometimes when using a special data type, like PostGIS geography
), while projected (or unprojected using the geometry
datatype) vertices are connected using straight lines.
This can result in a point being inside a polygon expressed as geography
but outside of it if expressed as geometry
The following example uses PostGIS. The polygon goes up to latitude 50, the point is at latitude 51.
WITH poly AS (select ST_GeomFromText('polygon((0 0, 50 0, 50 50, 0 50, 0 0))',4326) as geom),
pnt AS (select ST_GeomFromText('point(25 51)',4326) as geom)
SELECT ST_INTERSECTS(poly.geom,pnt.geom) intersect_geometry,
ST_INTERSECTS(poly.geom::geography,pnt.geom::geography) intersect_geography
FROM poly, pnt;
intersect_geometry | intersect_geography
--------------------+---------------------
f | t
Edit
Extending on @JR comment, here is an example when projecting a huge polygon to a Lambert Conformal Conic projection (3347) and checking the point intersection near the central meridian: we can see a difference of 4 degrees of latitude between the two!
WITH poly AS (select ST_GeomFromText('polygon((-70 45, -130 45, -130 50, -70 50, -70 45))',4326) as geom),
pnt AS (select ST_GeomFromText('point(-100 54)',4326) as geom)
SELECT ST_INTERSECTS(poly.geom,pnt.geom) intersect_geometry,
ST_INTERSECTS(st_transform(poly.geom,3347),st_transform(pnt.geom,3347)) intersect_reproject
FROM poly, pnt;
intersect_geometry | intersect_reproject
--------------------+---------------------
f | t
Some software will use great circle arcs to connect unprojected vertices (sometimes when using a special data type, like PostGIS geography
), while projected (or unprojected using the geometry
datatype) vertices are connected using straight lines.
This can result in a point being inside a polygon expressed as geography
but outside of it if expressed as geometry
The following example uses PostGIS. The polygon goes up to latitude 50, the point is at latitude 51.
WITH poly AS (select ST_GeomFromText('polygon((0 0, 50 0, 50 50, 0 50, 0 0))',4326) as geom),
pnt AS (select ST_GeomFromText('point(25 51)',4326) as geom)
SELECT ST_INTERSECTS(poly.geom,pnt.geom) intersect_geometry,
ST_INTERSECTS(poly.geom::geography,pnt.geom::geography) intersect_geography
FROM poly, pnt;
intersect_geometry | intersect_geography
--------------------+---------------------
f | t
Edit
Extending on @JR comment, here is an example when projecting a huge polygon to a Lambert Conformal Conic projection (3347) and checking the point intersection near the central meridian: we can see a difference of 4 degrees of latitude between the two!
WITH poly AS (select ST_GeomFromText('polygon((-70 45, -130 45, -130 50, -70 50, -70 45))',4326) as geom),
pnt AS (select ST_GeomFromText('point(-100 54)',4326) as geom)
SELECT ST_INTERSECTS(poly.geom,pnt.geom) intersect_geometry,
ST_INTERSECTS(st_transform(poly.geom,3347),st_transform(pnt.geom,3347)) intersect_reproject
FROM poly, pnt;
intersect_geometry | intersect_reproject
--------------------+---------------------
f | t
edited May 29 at 11:21
answered May 28 at 17:33
JGHJGH
14.4k21439
14.4k21439
Same in SQL server...
– Supereshek
May 29 at 9:55
add a comment |
Same in SQL server...
– Supereshek
May 29 at 9:55
Same in SQL server...
– Supereshek
May 29 at 9:55
Same in SQL server...
– Supereshek
May 29 at 9:55
add a comment |
Yes. projections will never reproject a point that was inside a polygon in one projection to be outside it in another--unless there's some sort of precision error. I'm not sure what this property is called in geography, but I just realized it's essentially relativistic invariance, which basically says that as time dilates for us, and our coordinate systems are compressed or stretched, no observer in any given frame will disagree on causal ordering of events. Likewise, no matter in what projection an "observer lives", no one will disagree on what points are in what polygons, even if they disagree on exactly how far the points are from the edges of the polygons, and stuff like that.
Well, except for horizon clipping, and discontinuous projections, which might throw a monkey wench in topology.
– Vince
May 28 at 16:53
3
Not sure but I think the question come from the way polygon may be incorrectly reprojected as in a big polygon defined by only corner vertex. When reprojecting you will project the corner then reconstruct the polygon by linking the corner with a straight line. In this case a point close to the border could be seen on the other side after reprojecting but that's because the projection is wrong along the length of the side (to prevent that you need to add vertex on the polygon side before reprojecting)
– J.R
May 28 at 16:54
I hadn't considered horizons or discontinous projections. I guess like in physics geography does have singularities where certain space-time coordinate systems break down; poles are another example.
– 0mn1
May 28 at 19:37
Yes...as J.R. is saying my question comes from the way a polygon may be incorrectly re-projected. I think ArcGIS 'Project' accounts for this if you use "preserve_shape" optional parameter. It adds extra vertices in order to get a better projected shape.
– Supereshek
May 29 at 9:54
@J.R I see. Is it correct to say that this happens because the PostGIS geometry is a different projection than the geography? If I connect to points on unprojected data with a straight line, then do the same thing with correct geodesic calculations and add that (apparently) curved line to the map, then there will be points between them. But the straight line was never the correct line for unprojected lon-lat data--it's not like the line is shifting from one side of a point to the other between two projections, rather the geometry projection is wrong.
– 0mn1
May 29 at 14:00
|
show 1 more comment
Yes. projections will never reproject a point that was inside a polygon in one projection to be outside it in another--unless there's some sort of precision error. I'm not sure what this property is called in geography, but I just realized it's essentially relativistic invariance, which basically says that as time dilates for us, and our coordinate systems are compressed or stretched, no observer in any given frame will disagree on causal ordering of events. Likewise, no matter in what projection an "observer lives", no one will disagree on what points are in what polygons, even if they disagree on exactly how far the points are from the edges of the polygons, and stuff like that.
Well, except for horizon clipping, and discontinuous projections, which might throw a monkey wench in topology.
– Vince
May 28 at 16:53
3
Not sure but I think the question come from the way polygon may be incorrectly reprojected as in a big polygon defined by only corner vertex. When reprojecting you will project the corner then reconstruct the polygon by linking the corner with a straight line. In this case a point close to the border could be seen on the other side after reprojecting but that's because the projection is wrong along the length of the side (to prevent that you need to add vertex on the polygon side before reprojecting)
– J.R
May 28 at 16:54
I hadn't considered horizons or discontinous projections. I guess like in physics geography does have singularities where certain space-time coordinate systems break down; poles are another example.
– 0mn1
May 28 at 19:37
Yes...as J.R. is saying my question comes from the way a polygon may be incorrectly re-projected. I think ArcGIS 'Project' accounts for this if you use "preserve_shape" optional parameter. It adds extra vertices in order to get a better projected shape.
– Supereshek
May 29 at 9:54
@J.R I see. Is it correct to say that this happens because the PostGIS geometry is a different projection than the geography? If I connect to points on unprojected data with a straight line, then do the same thing with correct geodesic calculations and add that (apparently) curved line to the map, then there will be points between them. But the straight line was never the correct line for unprojected lon-lat data--it's not like the line is shifting from one side of a point to the other between two projections, rather the geometry projection is wrong.
– 0mn1
May 29 at 14:00
|
show 1 more comment
Yes. projections will never reproject a point that was inside a polygon in one projection to be outside it in another--unless there's some sort of precision error. I'm not sure what this property is called in geography, but I just realized it's essentially relativistic invariance, which basically says that as time dilates for us, and our coordinate systems are compressed or stretched, no observer in any given frame will disagree on causal ordering of events. Likewise, no matter in what projection an "observer lives", no one will disagree on what points are in what polygons, even if they disagree on exactly how far the points are from the edges of the polygons, and stuff like that.
Yes. projections will never reproject a point that was inside a polygon in one projection to be outside it in another--unless there's some sort of precision error. I'm not sure what this property is called in geography, but I just realized it's essentially relativistic invariance, which basically says that as time dilates for us, and our coordinate systems are compressed or stretched, no observer in any given frame will disagree on causal ordering of events. Likewise, no matter in what projection an "observer lives", no one will disagree on what points are in what polygons, even if they disagree on exactly how far the points are from the edges of the polygons, and stuff like that.
answered May 28 at 16:37
0mn10mn1
784
784
Well, except for horizon clipping, and discontinuous projections, which might throw a monkey wench in topology.
– Vince
May 28 at 16:53
3
Not sure but I think the question come from the way polygon may be incorrectly reprojected as in a big polygon defined by only corner vertex. When reprojecting you will project the corner then reconstruct the polygon by linking the corner with a straight line. In this case a point close to the border could be seen on the other side after reprojecting but that's because the projection is wrong along the length of the side (to prevent that you need to add vertex on the polygon side before reprojecting)
– J.R
May 28 at 16:54
I hadn't considered horizons or discontinous projections. I guess like in physics geography does have singularities where certain space-time coordinate systems break down; poles are another example.
– 0mn1
May 28 at 19:37
Yes...as J.R. is saying my question comes from the way a polygon may be incorrectly re-projected. I think ArcGIS 'Project' accounts for this if you use "preserve_shape" optional parameter. It adds extra vertices in order to get a better projected shape.
– Supereshek
May 29 at 9:54
@J.R I see. Is it correct to say that this happens because the PostGIS geometry is a different projection than the geography? If I connect to points on unprojected data with a straight line, then do the same thing with correct geodesic calculations and add that (apparently) curved line to the map, then there will be points between them. But the straight line was never the correct line for unprojected lon-lat data--it's not like the line is shifting from one side of a point to the other between two projections, rather the geometry projection is wrong.
– 0mn1
May 29 at 14:00
|
show 1 more comment
Well, except for horizon clipping, and discontinuous projections, which might throw a monkey wench in topology.
– Vince
May 28 at 16:53
3
Not sure but I think the question come from the way polygon may be incorrectly reprojected as in a big polygon defined by only corner vertex. When reprojecting you will project the corner then reconstruct the polygon by linking the corner with a straight line. In this case a point close to the border could be seen on the other side after reprojecting but that's because the projection is wrong along the length of the side (to prevent that you need to add vertex on the polygon side before reprojecting)
– J.R
May 28 at 16:54
I hadn't considered horizons or discontinous projections. I guess like in physics geography does have singularities where certain space-time coordinate systems break down; poles are another example.
– 0mn1
May 28 at 19:37
Yes...as J.R. is saying my question comes from the way a polygon may be incorrectly re-projected. I think ArcGIS 'Project' accounts for this if you use "preserve_shape" optional parameter. It adds extra vertices in order to get a better projected shape.
– Supereshek
May 29 at 9:54
@J.R I see. Is it correct to say that this happens because the PostGIS geometry is a different projection than the geography? If I connect to points on unprojected data with a straight line, then do the same thing with correct geodesic calculations and add that (apparently) curved line to the map, then there will be points between them. But the straight line was never the correct line for unprojected lon-lat data--it's not like the line is shifting from one side of a point to the other between two projections, rather the geometry projection is wrong.
– 0mn1
May 29 at 14:00
Well, except for horizon clipping, and discontinuous projections, which might throw a monkey wench in topology.
– Vince
May 28 at 16:53
Well, except for horizon clipping, and discontinuous projections, which might throw a monkey wench in topology.
– Vince
May 28 at 16:53
3
3
Not sure but I think the question come from the way polygon may be incorrectly reprojected as in a big polygon defined by only corner vertex. When reprojecting you will project the corner then reconstruct the polygon by linking the corner with a straight line. In this case a point close to the border could be seen on the other side after reprojecting but that's because the projection is wrong along the length of the side (to prevent that you need to add vertex on the polygon side before reprojecting)
– J.R
May 28 at 16:54
Not sure but I think the question come from the way polygon may be incorrectly reprojected as in a big polygon defined by only corner vertex. When reprojecting you will project the corner then reconstruct the polygon by linking the corner with a straight line. In this case a point close to the border could be seen on the other side after reprojecting but that's because the projection is wrong along the length of the side (to prevent that you need to add vertex on the polygon side before reprojecting)
– J.R
May 28 at 16:54
I hadn't considered horizons or discontinous projections. I guess like in physics geography does have singularities where certain space-time coordinate systems break down; poles are another example.
– 0mn1
May 28 at 19:37
I hadn't considered horizons or discontinous projections. I guess like in physics geography does have singularities where certain space-time coordinate systems break down; poles are another example.
– 0mn1
May 28 at 19:37
Yes...as J.R. is saying my question comes from the way a polygon may be incorrectly re-projected. I think ArcGIS 'Project' accounts for this if you use "preserve_shape" optional parameter. It adds extra vertices in order to get a better projected shape.
– Supereshek
May 29 at 9:54
Yes...as J.R. is saying my question comes from the way a polygon may be incorrectly re-projected. I think ArcGIS 'Project' accounts for this if you use "preserve_shape" optional parameter. It adds extra vertices in order to get a better projected shape.
– Supereshek
May 29 at 9:54
@J.R I see. Is it correct to say that this happens because the PostGIS geometry is a different projection than the geography? If I connect to points on unprojected data with a straight line, then do the same thing with correct geodesic calculations and add that (apparently) curved line to the map, then there will be points between them. But the straight line was never the correct line for unprojected lon-lat data--it's not like the line is shifting from one side of a point to the other between two projections, rather the geometry projection is wrong.
– 0mn1
May 29 at 14:00
@J.R I see. Is it correct to say that this happens because the PostGIS geometry is a different projection than the geography? If I connect to points on unprojected data with a straight line, then do the same thing with correct geodesic calculations and add that (apparently) curved line to the map, then there will be points between them. But the straight line was never the correct line for unprojected lon-lat data--it's not like the line is shifting from one side of a point to the other between two projections, rather the geometry projection is wrong.
– 0mn1
May 29 at 14:00
|
show 1 more comment
Thanks for contributing an answer to Geographic Information Systems Stack Exchange!
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
To learn more, see our tips on writing great answers.
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fgis.stackexchange.com%2fquestions%2f324171%2fpoints-within-polygons-in-different-projections%23new-answer', 'question_page');
);
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown