Adding elevation to point in PyQGIS? Planned maintenance scheduled April 17/18, 2019 at 00:00UTC (8:00pm US/Eastern) Announcing the arrival of Valued Associate #679: Cesar Manara Unicorn Meta Zoo #1: Why another podcast?How to get the Z values of existing lineStringZ with pyqgis?How to get an elevation profile for a gps track?Elevation and point dataGPS point elevation vs Contours layer elevationAdd height/elevation in generated sql files from osm2poNeed to cluster the plots points into groups based on their elevationPostgis Self-intersection on geometry<->raster comparisonFinding elevation of points using PyQGIS?Extracting raster values to point feature class failing in ArcMap?Method asPoint().x() doesn't work in QGIS 3.4.3 but did in 3.4.1Adding tile layer via PyQGIS
What would be the ideal power source for a cybernetic eye?
How to draw this diagram using TikZ package?
3 doors, three guards, one stone
Is there a Spanish version of "dot your i's and cross your t's" that includes the letter 'ñ'?
If 'B is more likely given A', then 'A is more likely given B'
Storing hydrofluoric acid before the invention of plastics
Is high blood pressure ever a symptom attributable solely to dehydration?
When to stop saving and start investing?
Did Kevin spill real chili?
How can I fade player when goes inside or outside of the area?
How does a Death Domain cleric's Touch of Death feature work with Touch-range spells delivered by familiars?
What is this single-engine low-wing propeller plane?
Antler Helmet: Can it work?
How much radiation do nuclear physics experiments expose researchers to nowadays?
When is phishing education going too far?
Is it true to say that an hosting provider's DNS server is what links the entire hosting environment to ICANN?
What do you call a plan that's an alternative plan in case your initial plan fails?
How to find all the available tools in macOS terminal?
List *all* the tuples!
The logistics of corpse disposal
Why is "Captain Marvel" translated as male in Portugal?
How can players work together to take actions that are otherwise impossible?
Does accepting a pardon have any bearing on trying that person for the same crime in a sovereign jurisdiction?
How do I keep my slimes from escaping their pens?
Adding elevation to point in PyQGIS?
Planned maintenance scheduled April 17/18, 2019 at 00:00UTC (8:00pm US/Eastern)
Announcing the arrival of Valued Associate #679: Cesar Manara
Unicorn Meta Zoo #1: Why another podcast?How to get the Z values of existing lineStringZ with pyqgis?How to get an elevation profile for a gps track?Elevation and point dataGPS point elevation vs Contours layer elevationAdd height/elevation in generated sql files from osm2poNeed to cluster the plots points into groups based on their elevationPostgis Self-intersection on geometry<->raster comparisonFinding elevation of points using PyQGIS?Extracting raster values to point feature class failing in ArcMap?Method asPoint().x() doesn't work in QGIS 3.4.3 but did in 3.4.1Adding tile layer via PyQGIS
.everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty margin-bottom:0;
I am trying to add elevation to a 2D point in a QGIS 3 plugin, preferably in a concise and straightforward way. Not finding this easy.
I think part of the problem is that in QGIS 3, fromPoint()
has been renamed to fromPointXY()
. I am not very clear on how to create 3D point geometry.
I have tried a number of variations without finding one which works and have not been able to find a suitable example online.
Basically what I want to achieve is PointXY + elevation -> Point
, which I would expect to be something like this:
z = 1.0 # or whatever
for l in selectedLayers: # for each selected layer
if l.geometryType() == QgsWkbTypes.PointGeometry: # point layer
l.startEditing()
feat = l.getFeatures()
for f in feat: # each feature
# Get existing geometry...
geom = f.geometry()
pt = geom.asPoint() # returns PointXY
# Replace with 3D geometry...
pt3D = QgsPoint(pt.x(), pt.y(), z)
geom = QgsGeometry.fromPointXY(pt3D) # error, because pt is 3D
f.setGeometry(geom)
l.updateFeature(f)
l.commitChanges()
How do I do this?
For the benefit of others this is now resolved using the answer(s) below.
I am now adding the elevation like this:
f.setGeometry( QgsPoint(pt.x(), pt.y(), z) ) # x, y, z
It turns out that my problem was not so much adding the elevation, as knowing when I had successfully added it.
elevation pyqgis-3
add a comment |
I am trying to add elevation to a 2D point in a QGIS 3 plugin, preferably in a concise and straightforward way. Not finding this easy.
I think part of the problem is that in QGIS 3, fromPoint()
has been renamed to fromPointXY()
. I am not very clear on how to create 3D point geometry.
I have tried a number of variations without finding one which works and have not been able to find a suitable example online.
Basically what I want to achieve is PointXY + elevation -> Point
, which I would expect to be something like this:
z = 1.0 # or whatever
for l in selectedLayers: # for each selected layer
if l.geometryType() == QgsWkbTypes.PointGeometry: # point layer
l.startEditing()
feat = l.getFeatures()
for f in feat: # each feature
# Get existing geometry...
geom = f.geometry()
pt = geom.asPoint() # returns PointXY
# Replace with 3D geometry...
pt3D = QgsPoint(pt.x(), pt.y(), z)
geom = QgsGeometry.fromPointXY(pt3D) # error, because pt is 3D
f.setGeometry(geom)
l.updateFeature(f)
l.commitChanges()
How do I do this?
For the benefit of others this is now resolved using the answer(s) below.
I am now adding the elevation like this:
f.setGeometry( QgsPoint(pt.x(), pt.y(), z) ) # x, y, z
It turns out that my problem was not so much adding the elevation, as knowing when I had successfully added it.
elevation pyqgis-3
1
That contains 2 x ')' and 1 x '(' so syntactically I wouldn't expect that to work as written. It seems clear from the documentation that QgSPointXY is specifically a 2D point and does not support a z dimension. If you mean replace 'geom = QgsGeometry.fromPointXY(pt3D)' with something like 'geom = QgsGeometry.fromPointXY( QgsPointXY(1, 2), 3)' then no, that doesn't work.
– wotnot
Apr 9 at 7:31
1
QgsPointXY is specifically a 2D point and does not support a z dimension. This seems clear from the QGIS API documentation and I have tried it. 'pt3D = QgsPointXY((1,2), 3)' results in a TypeError.
– wotnot
Apr 9 at 7:44
Please write a self-answer in the area reserved for answers rather than in the area reserved for questions.
– PolyGeo♦
Apr 11 at 11:34
add a comment |
I am trying to add elevation to a 2D point in a QGIS 3 plugin, preferably in a concise and straightforward way. Not finding this easy.
I think part of the problem is that in QGIS 3, fromPoint()
has been renamed to fromPointXY()
. I am not very clear on how to create 3D point geometry.
I have tried a number of variations without finding one which works and have not been able to find a suitable example online.
Basically what I want to achieve is PointXY + elevation -> Point
, which I would expect to be something like this:
z = 1.0 # or whatever
for l in selectedLayers: # for each selected layer
if l.geometryType() == QgsWkbTypes.PointGeometry: # point layer
l.startEditing()
feat = l.getFeatures()
for f in feat: # each feature
# Get existing geometry...
geom = f.geometry()
pt = geom.asPoint() # returns PointXY
# Replace with 3D geometry...
pt3D = QgsPoint(pt.x(), pt.y(), z)
geom = QgsGeometry.fromPointXY(pt3D) # error, because pt is 3D
f.setGeometry(geom)
l.updateFeature(f)
l.commitChanges()
How do I do this?
For the benefit of others this is now resolved using the answer(s) below.
I am now adding the elevation like this:
f.setGeometry( QgsPoint(pt.x(), pt.y(), z) ) # x, y, z
It turns out that my problem was not so much adding the elevation, as knowing when I had successfully added it.
elevation pyqgis-3
I am trying to add elevation to a 2D point in a QGIS 3 plugin, preferably in a concise and straightforward way. Not finding this easy.
I think part of the problem is that in QGIS 3, fromPoint()
has been renamed to fromPointXY()
. I am not very clear on how to create 3D point geometry.
I have tried a number of variations without finding one which works and have not been able to find a suitable example online.
Basically what I want to achieve is PointXY + elevation -> Point
, which I would expect to be something like this:
z = 1.0 # or whatever
for l in selectedLayers: # for each selected layer
if l.geometryType() == QgsWkbTypes.PointGeometry: # point layer
l.startEditing()
feat = l.getFeatures()
for f in feat: # each feature
# Get existing geometry...
geom = f.geometry()
pt = geom.asPoint() # returns PointXY
# Replace with 3D geometry...
pt3D = QgsPoint(pt.x(), pt.y(), z)
geom = QgsGeometry.fromPointXY(pt3D) # error, because pt is 3D
f.setGeometry(geom)
l.updateFeature(f)
l.commitChanges()
How do I do this?
For the benefit of others this is now resolved using the answer(s) below.
I am now adding the elevation like this:
f.setGeometry( QgsPoint(pt.x(), pt.y(), z) ) # x, y, z
It turns out that my problem was not so much adding the elevation, as knowing when I had successfully added it.
elevation pyqgis-3
elevation pyqgis-3
edited Apr 11 at 11:34
PolyGeo♦
54k1782246
54k1782246
asked Apr 9 at 6:09
wotnotwotnot
999
999
1
That contains 2 x ')' and 1 x '(' so syntactically I wouldn't expect that to work as written. It seems clear from the documentation that QgSPointXY is specifically a 2D point and does not support a z dimension. If you mean replace 'geom = QgsGeometry.fromPointXY(pt3D)' with something like 'geom = QgsGeometry.fromPointXY( QgsPointXY(1, 2), 3)' then no, that doesn't work.
– wotnot
Apr 9 at 7:31
1
QgsPointXY is specifically a 2D point and does not support a z dimension. This seems clear from the QGIS API documentation and I have tried it. 'pt3D = QgsPointXY((1,2), 3)' results in a TypeError.
– wotnot
Apr 9 at 7:44
Please write a self-answer in the area reserved for answers rather than in the area reserved for questions.
– PolyGeo♦
Apr 11 at 11:34
add a comment |
1
That contains 2 x ')' and 1 x '(' so syntactically I wouldn't expect that to work as written. It seems clear from the documentation that QgSPointXY is specifically a 2D point and does not support a z dimension. If you mean replace 'geom = QgsGeometry.fromPointXY(pt3D)' with something like 'geom = QgsGeometry.fromPointXY( QgsPointXY(1, 2), 3)' then no, that doesn't work.
– wotnot
Apr 9 at 7:31
1
QgsPointXY is specifically a 2D point and does not support a z dimension. This seems clear from the QGIS API documentation and I have tried it. 'pt3D = QgsPointXY((1,2), 3)' results in a TypeError.
– wotnot
Apr 9 at 7:44
Please write a self-answer in the area reserved for answers rather than in the area reserved for questions.
– PolyGeo♦
Apr 11 at 11:34
1
1
That contains 2 x ')' and 1 x '(' so syntactically I wouldn't expect that to work as written. It seems clear from the documentation that QgSPointXY is specifically a 2D point and does not support a z dimension. If you mean replace 'geom = QgsGeometry.fromPointXY(pt3D)' with something like 'geom = QgsGeometry.fromPointXY( QgsPointXY(1, 2), 3)' then no, that doesn't work.
– wotnot
Apr 9 at 7:31
That contains 2 x ')' and 1 x '(' so syntactically I wouldn't expect that to work as written. It seems clear from the documentation that QgSPointXY is specifically a 2D point and does not support a z dimension. If you mean replace 'geom = QgsGeometry.fromPointXY(pt3D)' with something like 'geom = QgsGeometry.fromPointXY( QgsPointXY(1, 2), 3)' then no, that doesn't work.
– wotnot
Apr 9 at 7:31
1
1
QgsPointXY is specifically a 2D point and does not support a z dimension. This seems clear from the QGIS API documentation and I have tried it. 'pt3D = QgsPointXY((1,2), 3)' results in a TypeError.
– wotnot
Apr 9 at 7:44
QgsPointXY is specifically a 2D point and does not support a z dimension. This seems clear from the QGIS API documentation and I have tried it. 'pt3D = QgsPointXY((1,2), 3)' results in a TypeError.
– wotnot
Apr 9 at 7:44
Please write a self-answer in the area reserved for answers rather than in the area reserved for questions.
– PolyGeo♦
Apr 11 at 11:34
Please write a self-answer in the area reserved for answers rather than in the area reserved for questions.
– PolyGeo♦
Apr 11 at 11:34
add a comment |
2 Answers
2
active
oldest
votes
You don't need create a QgsGeometry object.
Add minimal sample for setGeometry 3D
feature = QgsFeature()
feature.setAttributes([lon, lat, alt])
feature.setGeometry(QgsPoint(lon, lat, alt))
layer.addFeatures([feature])
The original code is in my github,where a lot of layers are in 3d
https://github.com/All4Gis/QGISFMV/blob/35eb3f59b4d6423a375d7a04d54b6765dcc5c986/code/utils/QgsFmvLayers.py#L515
And for get to z value from feature,you can use something like that:
feature.geometry().vertexAt(0).z()
Thanks. Understood on the QgsGeometry object. It may be that I am successfully adding the elevation but I just cannot access the z value to check it. If I try to get the geometry using asPoint() it returns a 2D point and any z value is discarded. The API documentation is clear on this. So... how do I get the full coordinates of a 3D point?? I want to do something like this: 'QgsMessageLog.logMessage(str(f.geometry().z()))' but that doesn't work. I am getting 'AttributeError: 'QgsGeometry' object has no attribute 'z''.
– wotnot
Apr 9 at 14:36
update my answer
– Fran Raga
Apr 9 at 14:46
add a comment |
If you are looking at the Geometry of a feature you will have to do something like:
val_x = geom.get().x()
val_y = geom.get().y()
val_z = geom.get().z()
How to get the Z values of existing lineStringZ with pyqgis?
add a 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%2f318193%2fadding-elevation-to-point-in-pyqgis%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
You don't need create a QgsGeometry object.
Add minimal sample for setGeometry 3D
feature = QgsFeature()
feature.setAttributes([lon, lat, alt])
feature.setGeometry(QgsPoint(lon, lat, alt))
layer.addFeatures([feature])
The original code is in my github,where a lot of layers are in 3d
https://github.com/All4Gis/QGISFMV/blob/35eb3f59b4d6423a375d7a04d54b6765dcc5c986/code/utils/QgsFmvLayers.py#L515
And for get to z value from feature,you can use something like that:
feature.geometry().vertexAt(0).z()
Thanks. Understood on the QgsGeometry object. It may be that I am successfully adding the elevation but I just cannot access the z value to check it. If I try to get the geometry using asPoint() it returns a 2D point and any z value is discarded. The API documentation is clear on this. So... how do I get the full coordinates of a 3D point?? I want to do something like this: 'QgsMessageLog.logMessage(str(f.geometry().z()))' but that doesn't work. I am getting 'AttributeError: 'QgsGeometry' object has no attribute 'z''.
– wotnot
Apr 9 at 14:36
update my answer
– Fran Raga
Apr 9 at 14:46
add a comment |
You don't need create a QgsGeometry object.
Add minimal sample for setGeometry 3D
feature = QgsFeature()
feature.setAttributes([lon, lat, alt])
feature.setGeometry(QgsPoint(lon, lat, alt))
layer.addFeatures([feature])
The original code is in my github,where a lot of layers are in 3d
https://github.com/All4Gis/QGISFMV/blob/35eb3f59b4d6423a375d7a04d54b6765dcc5c986/code/utils/QgsFmvLayers.py#L515
And for get to z value from feature,you can use something like that:
feature.geometry().vertexAt(0).z()
Thanks. Understood on the QgsGeometry object. It may be that I am successfully adding the elevation but I just cannot access the z value to check it. If I try to get the geometry using asPoint() it returns a 2D point and any z value is discarded. The API documentation is clear on this. So... how do I get the full coordinates of a 3D point?? I want to do something like this: 'QgsMessageLog.logMessage(str(f.geometry().z()))' but that doesn't work. I am getting 'AttributeError: 'QgsGeometry' object has no attribute 'z''.
– wotnot
Apr 9 at 14:36
update my answer
– Fran Raga
Apr 9 at 14:46
add a comment |
You don't need create a QgsGeometry object.
Add minimal sample for setGeometry 3D
feature = QgsFeature()
feature.setAttributes([lon, lat, alt])
feature.setGeometry(QgsPoint(lon, lat, alt))
layer.addFeatures([feature])
The original code is in my github,where a lot of layers are in 3d
https://github.com/All4Gis/QGISFMV/blob/35eb3f59b4d6423a375d7a04d54b6765dcc5c986/code/utils/QgsFmvLayers.py#L515
And for get to z value from feature,you can use something like that:
feature.geometry().vertexAt(0).z()
You don't need create a QgsGeometry object.
Add minimal sample for setGeometry 3D
feature = QgsFeature()
feature.setAttributes([lon, lat, alt])
feature.setGeometry(QgsPoint(lon, lat, alt))
layer.addFeatures([feature])
The original code is in my github,where a lot of layers are in 3d
https://github.com/All4Gis/QGISFMV/blob/35eb3f59b4d6423a375d7a04d54b6765dcc5c986/code/utils/QgsFmvLayers.py#L515
And for get to z value from feature,you can use something like that:
feature.geometry().vertexAt(0).z()
edited Apr 9 at 14:46
answered Apr 9 at 13:45
Fran RagaFran Raga
3,35911123
3,35911123
Thanks. Understood on the QgsGeometry object. It may be that I am successfully adding the elevation but I just cannot access the z value to check it. If I try to get the geometry using asPoint() it returns a 2D point and any z value is discarded. The API documentation is clear on this. So... how do I get the full coordinates of a 3D point?? I want to do something like this: 'QgsMessageLog.logMessage(str(f.geometry().z()))' but that doesn't work. I am getting 'AttributeError: 'QgsGeometry' object has no attribute 'z''.
– wotnot
Apr 9 at 14:36
update my answer
– Fran Raga
Apr 9 at 14:46
add a comment |
Thanks. Understood on the QgsGeometry object. It may be that I am successfully adding the elevation but I just cannot access the z value to check it. If I try to get the geometry using asPoint() it returns a 2D point and any z value is discarded. The API documentation is clear on this. So... how do I get the full coordinates of a 3D point?? I want to do something like this: 'QgsMessageLog.logMessage(str(f.geometry().z()))' but that doesn't work. I am getting 'AttributeError: 'QgsGeometry' object has no attribute 'z''.
– wotnot
Apr 9 at 14:36
update my answer
– Fran Raga
Apr 9 at 14:46
Thanks. Understood on the QgsGeometry object. It may be that I am successfully adding the elevation but I just cannot access the z value to check it. If I try to get the geometry using asPoint() it returns a 2D point and any z value is discarded. The API documentation is clear on this. So... how do I get the full coordinates of a 3D point?? I want to do something like this: 'QgsMessageLog.logMessage(str(f.geometry().z()))' but that doesn't work. I am getting 'AttributeError: 'QgsGeometry' object has no attribute 'z''.
– wotnot
Apr 9 at 14:36
Thanks. Understood on the QgsGeometry object. It may be that I am successfully adding the elevation but I just cannot access the z value to check it. If I try to get the geometry using asPoint() it returns a 2D point and any z value is discarded. The API documentation is clear on this. So... how do I get the full coordinates of a 3D point?? I want to do something like this: 'QgsMessageLog.logMessage(str(f.geometry().z()))' but that doesn't work. I am getting 'AttributeError: 'QgsGeometry' object has no attribute 'z''.
– wotnot
Apr 9 at 14:36
update my answer
– Fran Raga
Apr 9 at 14:46
update my answer
– Fran Raga
Apr 9 at 14:46
add a comment |
If you are looking at the Geometry of a feature you will have to do something like:
val_x = geom.get().x()
val_y = geom.get().y()
val_z = geom.get().z()
How to get the Z values of existing lineStringZ with pyqgis?
add a comment |
If you are looking at the Geometry of a feature you will have to do something like:
val_x = geom.get().x()
val_y = geom.get().y()
val_z = geom.get().z()
How to get the Z values of existing lineStringZ with pyqgis?
add a comment |
If you are looking at the Geometry of a feature you will have to do something like:
val_x = geom.get().x()
val_y = geom.get().y()
val_z = geom.get().z()
How to get the Z values of existing lineStringZ with pyqgis?
If you are looking at the Geometry of a feature you will have to do something like:
val_x = geom.get().x()
val_y = geom.get().y()
val_z = geom.get().z()
How to get the Z values of existing lineStringZ with pyqgis?
answered Apr 9 at 15:12
Cary HCary H
1559
1559
add a comment |
add a 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%2f318193%2fadding-elevation-to-point-in-pyqgis%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
1
That contains 2 x ')' and 1 x '(' so syntactically I wouldn't expect that to work as written. It seems clear from the documentation that QgSPointXY is specifically a 2D point and does not support a z dimension. If you mean replace 'geom = QgsGeometry.fromPointXY(pt3D)' with something like 'geom = QgsGeometry.fromPointXY( QgsPointXY(1, 2), 3)' then no, that doesn't work.
– wotnot
Apr 9 at 7:31
1
QgsPointXY is specifically a 2D point and does not support a z dimension. This seems clear from the QGIS API documentation and I have tried it. 'pt3D = QgsPointXY((1,2), 3)' results in a TypeError.
– wotnot
Apr 9 at 7:44
Please write a self-answer in the area reserved for answers rather than in the area reserved for questions.
– PolyGeo♦
Apr 11 at 11:34