Tengo una tabla postgis que contiene polígonos. Quiero transformar mis polígonos en polígonos de 4 puntos: recuperar el MBR más pequeño. Aquí viene una ilustración de lo que quiero: Mi polígono original es de color rosa. Supongo que será bastante fácil obtener el MBR (con la función...
Publiqué ¿Hay algún problema de proyección con los datos de hidrología del Atlas de Canadá de Recursos Naturales de Canadá (NRC)? que describe un cambio en algunos de los datos de NRC (no todos). Aparentemente, este cambio tiene que ver con la conversión entre coordenadas geográficas, latitud y...
Según la documentación, las opciones son: Sintaxis ExportToTIFF (map_document, out_tiff, {data_frame}, {df_export_width}, {df_export_height}, {resolution}, {world_file}, {color_mode}, {tiff_compression}, {geoTIFF_tags}) Sin embargo, esto funciona: import arcpy.mapping mxd =...
¿Cuál es la forma más sencilla de rotar y renderizar datos de Openstreetmap y aún así obtener un buen renderizado? Por ejemplo , Toronto , encajaría mejor girado ligeramente en una página en modo horizontal, por lo que el límite norte de la ciudad es paralelo al borde superior de la página...
Estoy tratando de usar la proyección cartográfica Bertin de 1953, como se describe en el libro "Semiology of graphics" de Jacques Bertin. Pero no puedo encontrar la configuración en ninguno de mis software SIG. ¿Alguien tiene una
Cuando busco imágenes de la tierra se divide en longitud y latitudes, parece que hay dos tipos principales: es un cuadrado perfecto y las latitudes en él crecen a una velocidad diferente desde el ecuador hasta arriba (y abajo) [media, parece que la longitud de 0 a 66 es más corta que el camino...
Estoy tratando de que OpenLayers muestre el cuadro delimitador del área del mapa. Intenté usar map.getExtent () porque debería ser exactamente lo que necesito, pero los valores devueltos son "extraños": Por ejemplo, si centro el mapa alrededor de Grenoble (Francia, 45.19 5.71) con nivel de zoom...
Estoy buscando asesoramiento para trabajar con conjuntos de datos mundiales donde el GCS tiene un rango positivo (0 grados a 360 grados). Los datos con los que estoy trabajando son datos oceanográficos de NetCDF y tienen valores de coordenadas positivos como se mencionó. Cuando se muestra en un GCS...
Veo aquí que hay un código fuente de Fortran para pasar de USNG a lon / lat: http://www.ngs.noaa.gov/PC_PROD/USNG/ ¿Hay algún código de Python o C / C ++ que haga lo mismo? Entiendo que esta no es una función proporcionada por Proj4 / pyproj / GDAL,
Hasta donde entiendo la literatura, hay varias maneras de convertir un conjunto de coordenadas cartesianas geocéntricas en coordenadas geodésicas Bowring, 1976 Borkowski, 1989 Lin y Wang, 1995 Toms, 1996 ... ¿Cuál sería la forma más eficiente de convertir las coordenadas cartesianas en...
Actualmente tengo un proyecto en el que rasterizo archivos de formas en rásteres de presencia / ausencia simples (1 y 0) y luego agrego cada uno a un total acumulado para crear un ráster de riqueza total al final. Los archivos de forma están en proyección WGS84, pero necesito que el ráster final...
Tengo una base de datos con muchos puntos en WGS84. Ahora estoy construyendo un caché que realiza consultas NN y puntos en rango usando un KDtree. El punto [sic] es que el radio de búsqueda se proporcionará en metros y que lat / lon no es un buen SRS para estas consultas geométricas. Estoy...
Mientras trazo el gadmmapa de Rusia R, tengo un problema en el meridiano de 180 grados, donde parte de la región rusa se traza por separado. Puede verlo en el diagrama usando el siguiente código (es el borde izquierdo y la parte está a la derecha). ¿Alguien puede sugerir cómo evitar este...
Tengo alrededor de 100-200 fotos tomadas en el terreno que me gustaría superponer en ArcGIS 10 que no tienen información espacial que necesito desarrollar de alguna manera usando una fórmula fotogramétrica. La cámara se conectó a 2 metros sobre el suelo en un poste GPS RTK. Relación uno a uno entre...
Me gustaría crear mapas de una "Tierra inclinada", como el siguiente ejemplo donde Australia está en el Polo Norte: http://i1048.photobucket.com/albums/s379/laskaris_mods/tilted_earth.jpg Puedo reproyectar mapas en formatos de imagen como .jpeg como este con el software que estoy usando, Geocart...
Tengo un archivo de forma de países del mundo que se proyecta como equidistante azimutal con Chicago, EE. UU. Como centro. Me gustaría cambiar la proyección de modo que Edimburgo, Reino Unido, sea el centro. Uno podría pensar que esto debería ser tan simple como cambiar 2 números en los parámetros...
Soy un poco nuevo en todo el mundo SIG, y aunque estoy haciendo grandes avances dibujando polígonos y cosas así en Bing Maps, actualmente estoy estancado. En nuestros datos espaciales de SQL, me encontré con datos MULTIPOLYGON, donde todos los valores de coordenadas son números muy grandes (6...
He leído cosas en Internet que afirman que NAD83 y WGS84 son efectivamente equivalentes (con pequeñas diferencias de elipsoide). Sin embargo, si hago conversiones geodésicas (WGS84) -> ECEF -> geodésicas (NAD83), usando solo los parámetros elipsoides WGS84 y GRS80 y ninguna transformación...
Estoy usando TileStache para servir mis propios mosaicos vectoriales usando una fuente de datos postgis. Primero, volví a proyectar los archivos de forma que tenía usando ArcMap (reproyectado en Web Mercator (Esfera auxiliar)) y los cargué en Postgres usando shp2pgsql (SRID: 900913). Puedo generar...
1 pregunta He encontrado una advertencia usando la función projectRaster () en el paquete ráster en R. A continuación se pega un ejemplo reproducible completo. Warning message: In rgdal::rawTransform(projto_int, projfrom, nrow(xy), xy[, 1], : 33940 projected point(s) not finite Mi pregunta...