Eu tenho uma tabela postgis contendo polígonos. Quero transformar meus polígonos em polígonos de 4 pontos: recuperar o menor MBR. Aqui vem uma ilustração do que eu quero: Meu polígono original é na cor rosa. Eu acho que será muito fácil obter o MBR (com as funções ST_Envelope e BBox). Mas...
Publiquei Existe um problema de projeção nos dados de hidrologia do Atlas de Recursos Naturais do Canadá (NRC) no Canadá? descrevendo uma mudança em alguns dos dados do NRC (não todos). Aparentemente, essa mudança tem a ver com a conversão entre coordenadas geográficas, latitude e longitude e...
De acordo com a documentação, as opções são: Sintaxe ExportToTIFF (map_document, out_tiff, {data_frame}, {df_export_width}, {df_export_height}, {resolution}, {world_file}, {color_mode}, {tiff_compression}, {geoTIFF_tags}) No entanto, isso funciona: import arcpy.mapping mxd =...
Qual é a maneira mais simples de girar e renderizar dados do Openstreetmap e ainda obter uma boa renderização? Por exemplo , Toronto , caberia girar um pouco melhor em uma página no modo paisagem, para que o limite norte da cidade seja paralelo à borda superior da página impressa. Exemplo: Toronto...
Estou tentando usar a projeção do mapa Bertin de 1953, conforme descrito no livro "Semiologia dos gráficos", de Jacques Bertin. Mas não consigo encontrar as configurações em nenhum dos meus softwares GIS. Alguém tem uma solução
Quando procuro imagens da Terra, as longitudes e latitudes parecem ser os dois principais tipos: é um quadrado perfeito e as latitudes nele crescem uma taxa diferente do equador para cima (para baixo) [média, parece que o comprimento de 0 a 66 é mais curto que o caminho de 60 para o polo (90)]...
Estou tentando obter a caixa delimitadora da área do mapa exibida pelo OpenLayers. Tentei usar o map.getExtent () porque deveria ser exatamente o que eu precisava, mas os valores retornados são "estranhos": Por exemplo, se eu centralizar o mapa em torno de Grenoble (França, 45,19 5,71) com o...
Estou procurando conselhos sobre como trabalhar com conjuntos de dados mundiais nos quais o GCS tem um intervalo positivo (0 a 360 graus). Os dados com os quais estou trabalhando são dados oceanográficos do NetCDF e possuem valores de coordenadas positivos, conforme mencionado. Quando exibido em um...
Vejo aqui que existe o código-fonte Fortran para ir de USNG para lon / lat: http://www.ngs.noaa.gov/PC_PROD/USNG/ Existe algum código python ou C / C ++ que faz o mesmo? Entendo que isso não é uma função fornecida pelo Proj4 / pyproj / GDAL,
Até onde eu entendo a literatura, existem várias maneiras de converter um conjunto de coordenadas cartesianas geoecêntricas em coordenadas geodésicas Bowring, 1976 Borkowski, 1989 Lin e Wang, 1995 Toms, 1996 ... Qual seria a maneira mais eficiente de converter coordenadas cartesianas em...
Atualmente, tenho um projeto em que rasterizo os shapefiles em rasters de presença / ausência simples (1 e 0) e, em seguida, adiciono cada um a um total em execução para criar um raster de riqueza total no final. Os shapefiles estão na projeção WGS84, mas eu preciso que a varredura final esteja na...
Eu tenho um banco de dados com muitos pontos no WGS84. Agora estou construindo um cache que executa NN e aponta consultas de intervalo usando um KDtree. O ponto [sic] é que o raio de pesquisa será fornecido em metros e que lat / lon não é um bom SRS para essas consultas geométricas. Estou...
Ao traçar o gadmmapa da Rússia R, tenho um problema no meridiano de 180 graus, onde parte da região russa é plotada separadamente. Você pode vê-lo no gráfico usando o código a seguir (é a borda esquerda e a parte está à direita). Alguém pode sugerir como evitar esse problema? Aqui está uma versão...
Eu tenho cerca de 100-200 fotos de campo que gostaria de sobrepor no ArcGIS 10 que não possuem informações espaciais que eu preciso desenvolver de alguma forma usando uma fórmula fotogramétrica. A câmera foi acoplada 2 metros acima do solo em um poste de GPS RTK. Relação individual entre pontos de...
Eu gostaria de criar mapas de uma "Terra inclinada", como o exemplo a seguir, onde a Austrália está no Polo Norte: http://i1048.photobucket.com/albums/s379/laskaris_mods/tilted_earth.jpg Posso reprojetar mapas em formatos de imagem como .jpeg como este com o software que estou usando, Geocart 3....
Eu tenho um shapefile dos países do mundo que é projetado como Equidistante Azimutal com Chicago, EUA como o centro. Gostaria de mudar a projeção para que Edimburgo, Reino Unido, seja o centro. Alguém poderia pensar que isso deve ser tão simples quanto alterar 2 números nos parâmetros do sistema de...
Sou um pouco novo em todo o mundo GIS e, embora esteja fazendo grandes progressos desenhando polígonos e outros no Bing Maps, atualmente estou preso. Em nossos dados espaciais SQL, me deparei com dados MULTIPOLYGON, onde todos os valores de coordenadas são números muito grandes (6 e 7 dígitos para...
Li em toda a Internet coisas que afirmam que o NAD83 e o WGS84 são efetivamente equivalentes (com pequenas diferenças elipsóides). No entanto, se eu fizer conversões geodésicas (WGS84) -> ECEF -> geodésicas (NAD83), usando apenas os parâmetros elipsóides WGS84 e GRS80 e nenhuma transformação...
Estou usando o TileStache para servir meus próprios blocos de vetor usando uma fonte de dados postgis. Primeiro, eu re-projetei os shapefiles que eu tinha usando o ArcMap (reprojetado no Web Mercator (Auxilary Sphere)) e os carreguei no Postgres usando shp2pgsql (SRID: 900913). Eu sou capaz de...
1. Pergunta Encontrei um aviso usando a função projectRaster () no pacote raster em R. Um exemplo completo e reproduzível é colado abaixo. Warning message: In rgdal::rawTransform(projto_int, projfrom, nrow(xy), xy[, 1], : 33940 projected point(s) not finite Minha pergunta é: Esse aviso é um...