Hfrhyu

Can smartphones with the same camera sensor have different image quality?

Limit for e and 1/e

Antler Helmet: Can it work?

The following signatures were invalid: EXPKEYSIG 1397BC53640DB551

Are my PIs rude or am I just being too sensitive?

Communication vs. Technical skills ,which is more relevant for today's QA engineer positions?

Aligning matrix of nodes with grid

When communicating altitude with a '9' in it, should it be pronounced "nine hundred" or "niner hundred"?

Determine whether f is a function, an injection, a surjection

Does a C shift expression have unsigned type? Why would Splint warn about a right-shift?

How to say that you spent the night with someone, you were only sleeping and nothing else?

Writing Thesis: Copying from published papers

Is there folklore associating late breastfeeding with low intelligence and/or gullibility?

How does modal jazz use chord progressions?

Windows 10: How to Lock (not sleep) laptop on lid close?

Estimated State payment too big --> money back; + 2018 Tax Reform

Interesting examples of non-locally compact topological groups

When is phishing education going too far?

Why does tar appear to skip file contents when output file is /dev/null?

New Order #5: where Fibonacci and Beatty meet at Wythoff

Can a zero nonce be safely used with AES-GCM if the key is random and never used again?

How to colour the US map with Yellow, Green, Red and Blue to minimize the number of states with the colour of Green

What are the performance impacts of 'functional' Rust?

Need a suitable toxic chemical for a murder plot in my novel

Updating geometry of point using OGR in Python?

.everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty margin-bottom:0;

2

I need to change the coordinates (geometry) of some points in a shapefile. I can't find the good method.

if ( zone[best_y, best_x] >= ( nb_pixel_area - toler) ) and ( zone[best_y, best_x] <= ( nb_pixel_area + toler) ):
 print('======= station deplacee')
 new_x = (( pix_x + best_x ) * hres ) + georef[0]
 new_y = (( pix_y + best_y ) * abs(hres) ) + georef[3]
 wkt = "POINT(%f %f)" % (new_x , new_y)
 point = ogr.CreateGeometryFromWkt(wkt)
 gauge.SetGeometryDirectly(point)

I was also not able to use the setX() and setY().

Has someone already done that action?

python ogr

share|improve this question

edited Apr 8 at 13:21

PolyGeo♦

54k1782246

asked Jan 27 '15 at 12:34

ChristineChristine

1709

• 
  
  
  

  
  1
  

  
  
  
  
  
  

  
  
  I'm not working under Qgis. Anyway, I've found the problem. This statement was missing : reloc_lyr.SetFeature(gauge)
  
  – Christine
  Jan 28 '15 at 7:28
  
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  If you are able to answer your own question then please expand your comment into a self-answer.
  
  – PolyGeo♦
  Apr 8 at 13:22
  

  
  
  
  

add a comment | 

2

I need to change the coordinates (geometry) of some points in a shapefile. I can't find the good method.

if ( zone[best_y, best_x] >= ( nb_pixel_area - toler) ) and ( zone[best_y, best_x] <= ( nb_pixel_area + toler) ):
 print('======= station deplacee')
 new_x = (( pix_x + best_x ) * hres ) + georef[0]
 new_y = (( pix_y + best_y ) * abs(hres) ) + georef[3]
 wkt = "POINT(%f %f)" % (new_x , new_y)
 point = ogr.CreateGeometryFromWkt(wkt)
 gauge.SetGeometryDirectly(point)

I was also not able to use the setX() and setY().

Has someone already done that action?

python ogr

share|improve this question

edited Apr 8 at 13:21

PolyGeo♦

54k1782246

asked Jan 27 '15 at 12:34

ChristineChristine

1709

• 
  
  
  

  
  1
  

  
  
  
  
  
  

  
  
  I'm not working under Qgis. Anyway, I've found the problem. This statement was missing : reloc_lyr.SetFeature(gauge)
  
  – Christine
  Jan 28 '15 at 7:28
  
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  If you are able to answer your own question then please expand your comment into a self-answer.
  
  – PolyGeo♦
  Apr 8 at 13:22
  

  
  
  
  

add a comment | 

I need to change the coordinates (geometry) of some points in a shapefile. I can't find the good method.

if ( zone[best_y, best_x] >= ( nb_pixel_area - toler) ) and ( zone[best_y, best_x] <= ( nb_pixel_area + toler) ):
 print('======= station deplacee')
 new_x = (( pix_x + best_x ) * hres ) + georef[0]
 new_y = (( pix_y + best_y ) * abs(hres) ) + georef[3]
 wkt = "POINT(%f %f)" % (new_x , new_y)
 point = ogr.CreateGeometryFromWkt(wkt)
 gauge.SetGeometryDirectly(point)

I was also not able to use the setX() and setY().

Has someone already done that action?

python ogr

share|improve this question

edited Apr 8 at 13:21

PolyGeo♦

54k1782246

asked Jan 27 '15 at 12:34

ChristineChristine

1709

if ( zone[best_y, best_x] >= ( nb_pixel_area - toler) ) and ( zone[best_y, best_x] <= ( nb_pixel_area + toler) ):
 print('======= station deplacee')
 new_x = (( pix_x + best_x ) * hres ) + georef[0]
 new_y = (( pix_y + best_y ) * abs(hres) ) + georef[3]
 wkt = "POINT(%f %f)" % (new_x , new_y)
 point = ogr.CreateGeometryFromWkt(wkt)
 gauge.SetGeometryDirectly(point)

share|improve this question

edited Apr 8 at 13:21

PolyGeo♦

54k1782246

asked Jan 27 '15 at 12:34

ChristineChristine

1709

share|improve this question

edited Apr 8 at 13:21

PolyGeo♦

54k1782246

asked Jan 27 '15 at 12:34

ChristineChristine

1709

edited Apr 8 at 13:21

PolyGeo♦

54k1782246

edited Apr 8 at 13:21

PolyGeo♦

54k1782246

asked Jan 27 '15 at 12:34

ChristineChristine

1709

asked Jan 27 '15 at 12:34

ChristineChristine

1709

1 Answer
1

active

oldest

votes

0

Have you seen this snippet:

Modify Features

It is possible to either change feature’s geometry or to change some
attributes. The following example first changes values of attributes
with index 0 and 1, then it changes the feature’s geometry

fid = 100 # ID of the feature we will modify

if caps & QgsVectorDataProvider.ChangeAttributeValues:
 attrs = 0 : "hello", 1 : 123 
 layer.dataProvider().changeAttributeValues( fid : attrs )

if caps & QgsVectorDataProvider.ChangeGeometries:
 geom = QgsGeometry.fromPoint(QgsPoint(111,222))
 layer.dataProvider().changeGeometryValues( fid : geom )

share|improve this answer

answered Jan 27 '15 at 12:59

artwork21artwork21

31.1k555120

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



);

draft saved

draft discarded

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

Name

Email

Required, but never shown

StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fgis.stackexchange.com%2fquestions%2f131930%2fupdating-geometry-of-point-using-ogr-in-python%23new-answer', 'question_page');

);

Post as a guest

Name

Email

Required, but never shown

0

Have you seen this snippet:

Modify Features

It is possible to either change feature’s geometry or to change some
attributes. The following example first changes values of attributes
with index 0 and 1, then it changes the feature’s geometry

fid = 100 # ID of the feature we will modify

if caps & QgsVectorDataProvider.ChangeAttributeValues:
 attrs = 0 : "hello", 1 : 123 
 layer.dataProvider().changeAttributeValues( fid : attrs )

if caps & QgsVectorDataProvider.ChangeGeometries:
 geom = QgsGeometry.fromPoint(QgsPoint(111,222))
 layer.dataProvider().changeGeometryValues( fid : geom )

share|improve this answer

answered Jan 27 '15 at 12:59

artwork21artwork21

31.1k555120

add a comment | 

0

Have you seen this snippet:

Modify Features

It is possible to either change feature’s geometry or to change some
attributes. The following example first changes values of attributes
with index 0 and 1, then it changes the feature’s geometry

fid = 100 # ID of the feature we will modify

if caps & QgsVectorDataProvider.ChangeAttributeValues:
 attrs = 0 : "hello", 1 : 123 
 layer.dataProvider().changeAttributeValues( fid : attrs )

if caps & QgsVectorDataProvider.ChangeGeometries:
 geom = QgsGeometry.fromPoint(QgsPoint(111,222))
 layer.dataProvider().changeGeometryValues( fid : geom )

share|improve this answer

answered Jan 27 '15 at 12:59

artwork21artwork21

31.1k555120

add a comment | 

Have you seen this snippet:

Modify Features

It is possible to either change feature’s geometry or to change some
attributes. The following example first changes values of attributes
with index 0 and 1, then it changes the feature’s geometry

fid = 100 # ID of the feature we will modify

if caps & QgsVectorDataProvider.ChangeAttributeValues:
 attrs = 0 : "hello", 1 : 123 
 layer.dataProvider().changeAttributeValues( fid : attrs )

if caps & QgsVectorDataProvider.ChangeGeometries:
 geom = QgsGeometry.fromPoint(QgsPoint(111,222))
 layer.dataProvider().changeGeometryValues( fid : geom )

share|improve this answer

answered Jan 27 '15 at 12:59

artwork21artwork21

31.1k555120

fid = 100 # ID of the feature we will modify

if caps & QgsVectorDataProvider.ChangeAttributeValues:
 attrs = 0 : "hello", 1 : 123 
 layer.dataProvider().changeAttributeValues( fid : attrs )

if caps & QgsVectorDataProvider.ChangeGeometries:
 geom = QgsGeometry.fromPoint(QgsPoint(111,222))
 layer.dataProvider().changeGeometryValues( fid : geom )

share|improve this answer

answered Jan 27 '15 at 12:59

artwork21artwork21

31.1k555120

answered Jan 27 '15 at 12:59

artwork21artwork21

31.1k555120

answered Jan 27 '15 at 12:59

artwork21artwork21

31.1k555120