A veces veo srsName="EPSG:4326"en documentos GML (por ejemplo , https://postgis.net/docs/ST_GeomFromGML.html , aunque mi fuente particular no es pública). La intención (excepto quizás para el orden de los ejes que no es importante aquí) es clara y probablemente interoperable en varios sistemas,...
Estoy trabajando en un proyecto que requiere un búfer de 30 m alrededor de una característica lineal. La característica representa la costa y necesito el búfer para mostrar 30m tierra adentro desde la línea que representa la costa. ¿Hay alguna manera de hacer esto usando
¿Cómo mostrar etiquetas para puntos geojson en un mapa de folleto? Hay Leaflet.label que ahora está en desuso a favor de L.Tooltip , pero que solo muestra texto al pasar el mouse. Quiero mostrar las etiquetas de texto directamente en el mapa sin necesidad de la interacción del usuario. Entrada de...
Estoy intentando escribir un script que guardará una representación de varias capas usando el compositor de mapas. El problema con el que me encuentro es que el script se guarda antes de que qgis haya terminado de procesar todas las capas. Basado en varias otras respuestas ( 1 , 2 , 3 ), he...
Necesito realizar una nitidez panorámica de Gram-Schmidt , no tengo acceso a ArcGIS o ENVI, por lo que he estado usando QGIS . ¿Existe algún complemento o capacidad similar que esté disponible para realizar la nitidez panorámica de Gram-Schmidt en QGIS? Si no, ¿es esto posible en cualquier otro...
Estoy tratando de encontrar dónde varias líneas se cruzan con un polígono para dos geodataframes diferentes: from shapely.geometry import Polygon, LineString import geopandas as gpd polygon = Polygon([(0, 0), (1, 0), (1, 1), (0, 1), (0, 0)]) line1 = LineString([(0.5, 0.5), (0.7, 0.7)]) line2 =...
Tengo un conjunto de datos con más de 2000 funciones basadas en direcciones. Necesito acercarme a cada punto individualmente y moverlos un poco (sin desplazamiento constante, por lo que no hay otra forma de hacerlo). Seleccionar la siguiente característica en la tabla de atributos y presionar el...
Estoy escribiendo una función que tiene una geometría como parámetro que puede estar en diferentes proyecciones: CREATE OR REPLACE FUNCTION foo(in geometry) RETURNS boolean AS $$ DECLARE transformed_geom geometry := in; BEGIN IF ST_SRID(transformed_geom) != 32737 THEN transformed_geom :=...
Tengo una función que crea paneles fotovoltaicos solares representados como polígonos. Básicamente, crea una cuadrícula rectangular donde el usuario puede especificar los siguientes parámetros: Longitud Anchura Distancia horizontal Distancia vertical El código se basa en el complemento...
Necesito crear una cuadrícula de sector que se irradie desde un punto de interés. Aquí hay un ejemplo que laboriosamente puse a mano: ¿Encontré un par de soluciones similares en Cómo generar archivos de forma de tablero de dardos usando ArcGIS for Desktop? y Dibuje rosa de los vientos con QGIS...
Tengo dos archivos vectoriales en el panel de capas en QGIS (V2.8.2, Wien). Estoy tratando de seleccionar puntos de ambos archivos vectoriales simultáneamente para verlos en ambas tablas de atributos: ambas tablas están configuradas para mostrar solo las características seleccionadas. La única...
Estoy buscando una manera de fusionar los archivos de banda Sentinel 2 jp2 ( B02, B03, B04 ) y corregir los colores RGB. Todo debe hacerse con script bash o python. Para mi ejemplo, trabajo en estas imágenes . Idealmente, la solución estará cerca de este tutorial. Puedo fusionar las bandas con...
Estoy trabajando con archivos de forma censal (multipolígonos) y tengo otra tabla con puntos. Registramos en qué bloque censal se encuentra cada punto. Tengo un punto que devuelve falso para todas las consultas ST_Within y ST_Touches. Sin embargo, devuelve verdadero para ST_Intersects en un bloque....
Me gustaría saber cómo obtener todos los valores ráster dentro de un polígono usando gdal o pygeoprocessing, sin leer toda la cuadrícula como una matriz. pygeoprocessing y gdal pueden hacer estadísticas zonales, pero solo la función min, max, mean, stdev o count están disponibles desde dicha...
Quiero crear una capa temporal a partir de una capa existente y solo obtener esas características de la capa existente, que están seleccionadas. Sé que ArcGIS tiene esta funcionalidad que puede crear una capa a partir de entidades seleccionadas. Sin embargo, no pude encontrar ninguna solución de...
Estoy tratando de abrir un archivo shape con el cuaderno Jupyter, pero no puedo abrirlo. El archivo de forma se encuentra en la misma carpeta que mi libro de trabajo, y Python no puede identificarlo incluso si uso la ruta completa. import shapefile as shp import matplotlib.pyplot as plt sf =...
Tengo un marco de datos geopandas en EPSG: 4326 y crearía un nuevo marco de datos que consta de todas las filas que se encuentran dentro de un cuadro delimitador determinado. Primero obtengo el cuadro delimitador que me interesa (que en realidad es el cuadro delimitador de otro marco de...
Hace algún tiempo escuché sobre un nuevo DEM global llamado NASADEM que se lanzaría en 2017, y prometía ser el mejor DEM global hasta la fecha. ¿Qué es exactamente NASADEM y cuándo estará
¿En qué (decimal) lugar decimal de los valores de coordenadas QGIS define que cualquier conjunto de puntos está duplicado? Había estado pensando que QGIS maneja ~ 15 decimales; pero era la limitación solo porque había estado trabajando principalmente en Shapefiles. Sin embargo, el límite no se...
Estoy tratando de cambiar el tamaño de página en QGIS 3.0 Composer pero no puedo ver la pestaña Composición. He hecho clic derecho en un espacio vacío en la barra de herramientas pero no vi nada sobre la composición. Estoy usando Windows 7. Intenté mirar los documentos de QGIS pero todavía están en...