¿Cómo convierto un polígono en un solo punto o en todos los puntos?

Question

Quiero convertir los polígonos en puntos. Sé cómo se hace en ArcGis y creo que esto puede sonar trivial pero es un asunto serio.

Apreciaré cualquier indicio :)

Answer 1

Usted puede hacer polígono a centroide en QGIS usando:

Vector - Geometry Tools - Polygon Centroids

Answer 2

Usted puede hacer polígono a los nodos de la frontera en QGIS utilizando:

Vector - Geometry Tools - Nodes extraction

Answer 3

Ha determinado correctamente que $Y \sim \rm Exponential(\mathrm{e}^\theta)$ y la distribución determinada de $T$ Así que me concentraré en la computación

$$ \mathbb{E}( \log x, x \sim \rm Gamma(n, \frac{1}{\beta})) $$

La forma más sencilla es calcular $\mathbb{E}( x^s )$ y luego usar $\log(x) = \left. \frac{\mathrm{d}}{\mathrm{d} s} x^s \right\vert_{s=0}$ y $\log^2(x) = \left. \frac{\mathrm{d}^2}{\mathrm{d} s^2} x^s \right\vert_{s=0}$ .

$$ \mathbb{E}( x^s ) = \int_0^\infty x^s \cdot \frac{\beta^{n}}{\Gamma(n)} x^{n-1} \exp(-x \beta) \mathrm{d} x = \int_0^\infty \frac{\beta^{n}}{\Gamma(n)} x^{n + s-1} \exp(-x \beta) \mathrm{d} x = \beta^{-s} \frac{\Gamma(n+s)}{\Gamma(n)} $$

Ahora es fácil de entender:

$$ \mathbb{E}(\log T) = -\theta + \psi(n) \qquad \mathbb{V}(\log T) = \psi^{(1)}(n) $$ donde $\psi(n)$ es función digamma y $\psi^{(1)}(n)$ es su primera derivada.

Añadido : Para facilitar el cálculo de la varianza, hay que tener en cuenta que $\mathbb{V}(\log T) = \left. \frac{ \mathrm{d}^2}{ \mathrm{d}s^2 } \mu^{-s} \mathbb{E}(x^s)\right\vert_{s=0}$ , donde $\log \mu$ es la media de $\log T$ pero $\mu^{-s} \mathbb{E}(x^s) = \frac{\Gamma(n+s)}{\Gamma(n)} \exp(-s \psi(n)) = \exp( 1 + \psi(n) s + \frac{1}{2} \psi^\prime(n) s^2 + o(s^3)) \exp(-s \psi(n))$ . Por lo tanto, $ \mu^{-s} \mathbb{E}(x^s) = \exp(1 + \frac{1}{2} \psi^\prime(n) s^2 + o(s^3)) = 1 + \frac{1}{2} \psi^\prime(n) s^2 + o(s^3)$ .

Answer 4

Arcobjects:

IPolygon polygon = (IPolygon)geometry;

IPointCollection pointCollection = (IPointCollection)polygon;

Answer 5

Podrías empezar por aquí:

Característica al punto

import arcpy
from arcpy import env
env.workspace = "C:/data"
arcpy.FeatureToPoint_management("parcels.shp", "c:/data/output/parcels_center.shp", 
                            "CENTROID")

---OR---

Vértices de la característica al punto

aquí:

import arcpy
from arcpy import env
env.workspace = "C:/data"
arcpy.FeatureVerticesToPoints_management("parcels.shp",
                                     "c:/output/output.gdb/parcels_corner", 
                                     "ALL")