Coğrafi Bilgi Sistemleri

Kriging yapmak istediğimde, datatable'ımda kullandığım değerlere bağlı olarak sadece bazen çalışıyor. Krige işlevinin bir sonucu olarak alıyorum var1.pred: NA NA NA ...ve var1.var: NA NA NA ...(ancak yalnızca datatable'ımda "yanlış" değerleri kullandığımda). Örneğin: o çalışır her zaman (şimdiye kadar) Ben sadece 10 değerler kullandığınızda o çalışır Ben 50 değerler kullandığınızda …

10 r  interpolation  kriging  gstat 

Windows 7'de R'ye okumak istediğim çok büyük bir ESRI şekil dosyası var (> 8000 çokgen, .shp dosyası> 32MB). Rgdal ve readOGR () yöntemine aşinayım. Bununla birlikte, bu şekil dosyası boyutuyla, tüm şekil dosyasını belleğe okumak çok yavaş bir işlemdir. Ayrıca, şekil dosyasında bulunan 8.000'in yalnızca birkaç kümelenmiş çokgeniyle (<100) ilgileniyorum. …

10 shapefile  r  rgdal 

QGIS ile pek tecrübem yok ve mmqgis eklentisi ile oluşturulan yerleri birbirine bağlayan çizgilerle bir harita üzerinde çalışıyorum. Bu direkt çizgileri estetik amaçlar için biraz kıvrımlı yapmanın kolay bir yolu var mı? Aşağıdaki harita gibi bir şey elde etmeye çalışıyorum ama birkaç hub ile. PostGIS ve diğer gelişmiş araçlarla bunun …

10 qgis  qgis-plugins 

Ölümcül trafik kazalarını temsil eden noktaları olan bir katman var ve alan "dates_14_D" tarihlerini içerir, ancak türü dize. Bu dize alanını tarih alanına dönüştürmek istiyorum ancak Tamam düğmesi etkin değil. Sorun ne? QGis kullanıyorum. Aşağıdaki resme bakın

10 qgis  field-calculator  convert  fields-attributes  date 

Yeni bir alan eklemek ve her özellik için değerleri hesaplamak için PyQGIS kullanmak istiyorum. Alan Hesaplayıcı seçeneğine benzer. 'Alan Hesaplayıcı' ifadem örneğin: y(start_point($geometry)) from PyQt4.QtCore import QVariant from qgis.core import QgsField, QgsExpression, QgsFeature vl = iface.activeLayer() vl.startEditing() #step 1 myField = QgsField( 'myNewColumn', QVariant.Float ) vl.addAttribute( myField ) idx = …

10 pyqgis  fields-attributes 

Sahip olduğumuz yüksek çözünürlüklü bir görüntü için fayans üretmek için biraz uğraşıyorum. Şu anki görüntümüz GeoTiff dosyası olarak kaydedilen çok büyük (+ 20 GB) bir görüntüdür. Gdal2tiles komut satırı yardımcı programını kullanarak fayans oluşturmak ve daha sonra fayans sağlamak için TMS görüntü sağlayıcı kullanarak Cesium, açmak ve görüntülemek istiyorum. Gdalinfo …

10 gdal  gdalwarp  cesium  gdal2tiles 

Haritaların sayfanın kenarına kadar ilerlediği bir kitap için haritalar yapıyorum. Yazıcı, görüntülerin çevresinde 3 mm kanama alanı önerdi, bu nedenle sayfam 127x177 mm ise, görüntüler 133x183 mm olmalıdır. Bu taşma payı alanı, kitap sayfaları kesildiğinde hatalara izin vermektir ve haritalarda önemli bir bilgi olmamalıdır. Bu yüzden her tarafta bu ekstra …

10 qgis  labeling  print-composer 

ArcMap 10.2.2'yi Gelişmiş lisansla kullanıyorum. ArcMap'ı diğer temel analizler için kullandım, ancak program için oldukça yeniyim ve sadece temel kodlama becerilerine sahibim, ArcMap'te SQL kod oluşturucularını kullanan çok temel sorguların dışında SQL deneyimi yok . Türlerin coğrafi aralıklarının çokgenleri var ve bu türlerin taşıdığı patojenlerin / hastalıkların bir listesini içeren …

10 arcgis-desktop  arcgis-10.2  attribute-joins  spatial-join 

Kapalı . Bu sorunun daha fazla odaklanması gerekiyor . Şu anda cevapları kabul etmiyor. Bu soruyu geliştirmek ister misiniz? Soruyu, yalnızca bu yayını düzenleyerek tek bir soruna odaklanacak şekilde güncelleyin . 3 yıl önce kapalı . 1 hektarlık bir bloktaki ananas bitkilerinin sayısını nasıl hesaplayabilirim? Bu algoritmaları duydum ama nasıl …

10 remote-sensing  imagery  density  feature-extraction  unmanned-aerial-vehicle 

Uydu yörünge biçimleri için çokgenler üretmeye çalışıyorum. Şimdiye kadar [lat, long] her bir alanın kenarını temsil eden iki çizgi oluşturmak için bir yöntemim var. Bazı dallar uluslararası datelini geçiyor ve etrafı sarıyorlar: Bunu şu şekilde çözebildim ogr2ogr -wrapdateline: ogr2ogr -wrapdateline -f "ESRI Shapefile" test.shp orbits.shp Muhtemelen çizgileri bölüyor Şimdi her …

10 python  ogr  wgs84  shapely  antimeridian 

Bir MSSQL DB SRID 4269 olarak depolanan devlet sınırlarına karşı ABD nüfus sayımı verileri var . Bu veri kümesine karşı hesaplamalar için kullanılan giriş verileri SRID 4326 olarak saklanır . Bildiğim kadarıyla anlayabilirsiniz olarak SpatialReference.org , SRID 4269 SRID 4326 sadece bir alt kümesidir ve hiçbir yeniden projeksiyon gereklidir koordinat. …

10 coordinate-system  wgs84  srid  nad83 

Ggplot2 kullanarak uzamsal verilerimi çizerken sorun yaşıyorum. Harita spplot kullanılarak çizildiğinde iyi görünüyor, bu yüzden yırtılmanın takviye aşamasında gerçekleştiğini varsayıyorum. Kod aşağıdaki gibidir: #install the packages library(rgdal) library(mapproj) library(raster) library(rgeos) library(ggplot2) library(plyr) if (!require(gpclib)) install.packages("gpclib", type="source") gpclibPermit() setwd("C:/Users/My Documents") #read in laa to regional mapping #must aggregate to higher level …

10 r  ggplot2  rgeos 

ArcGIS pythonscript ile bir dizi polyilnes içeren bir şekil dosyasının tüm uç noktalarını birbirine bağlayan bir çokgen oluşturmaya çalışıyorum, çokgendeki düğümlerin sırası önemli olduğu için bunu yaparken sorun yaşıyorum. Resimdeki gri poligonu yeşil çizgilerden elde etmek istiyorum

10 arcpy  polygon  line 

Bir özellik sınıfındaki tüm satırları silmek istiyorum, ancak 1 - 5 arasındaki satırlar olarak kalıyorum. Diğer yandan ilk beş satır olarak kalmak ve diğerlerini silmek istiyorum. Arama imlecini ve güncelleme imlecini kullanmam gerektiğini biliyorum ama kullanamadım. ArcPy kullanarak satırları nasıl silebilirim?

10 arcpy  cursor 

Vektör geometrilerini basitleştirmek için QGIS 2.8 kullanmaya çalışıyorum. Ancak, Geometrileri basitleştir aracında Toleransı basitleştir değerinin önemini anlamıyorum. Dokümantasyon bunu açıklamıyor.

10 qgis  geoprocessing  simplify