Запитання з тегом «gdal»

Після використання gdalwarpдля проектування та вирівнювання до сітки (через -tap) декілька растрових я помітив, що вихідний растр був значно більшим, ніж оригінальний растр. Досить ретельний веб-пошук виявив цю проблему Trac : Френк Вармердам пояснив причину: "При уважному огляді різниця у розглянутому файлі полягає в тому, що gdal_translate використовує інтерфейс TIFFWriteScanline …

19 gdal  coordinate-system  geotiff-tiff  gdalwarp  file-size 

Так, GDAL нещодавно додав нову функцію, яка дозволяє випадково читати файли з відра S3. Я хочу обрізати зображення GDAL з декількох плиток зображення, не завантажуючи весь файл. Я бачив дуже рідкісну документацію про те, як налаштувати та отримати доступ до відра S3, хоча GDAL, і я трохи не збентежений, як …

19 raster  gdal  gdal-merge  amazon-ec2  amazon-s3 

Гаразд, трохи нерозумного питання, але: Я чув, як люди вимовляють бібліотеку GDAL (Бібліотека абстракції геопросторових даних), як goooo-dal. Моя природна схильність полягала в тому, щоб вимовляти його як джи-дал, як у букві 'g'-dal. Будь-яке пояснення того, як goo-dal став популярним? Чи був залучений goo-gle?

18 gdal 

У Python, використовуючи GDAL, я витягнув проекцію растру як рядок WKT таким чином: wkt = dataset.GetProjection() # wkt is 'PROJCS["GDA_1994_Transverse_Mercator",GEOGCS["GDA_1994",DATUM["GDA_1994",SPHEROID["GRS_1980",6378137,298.2572221010002],TOWGS84[0,0,0,0,0,0,0]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",117],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",10000000],UNIT["Meter",1]]' Використовуючи рядок WKT, я можу створити екземпляр SpatialReference наступним чином: src = osr.SpatialReference() src.ImportFromWkt(wkt) Це легко-горох. Я можу досить легко витягти кілька параметрів проекції з srcтаких, як зона UTM тощо. …

18 python  gdal  coordinate-system  python-2.7 

У мене є растр (насправді USGS DEM), і мені потрібно розділити його на менші шматки, як показано на малюнку нижче. Це було досягнуто в ArcGIS 10.0 за допомогою інструмента Split Raster. Я хотів би це зробити методом FOSS. Я дивився на GDAL, думаючи, що це зроблять це (якось із gdal_translate), …

18 raster  gdal  grass  dem  open-source-gis 

Ця публікація є спільнотою Wiki . Відредагуйте наявні відповіді, щоб покращити цю публікацію. Наразі не приймає нових відповідей. Я розробляю деякий код у Python, і хочу використовувати прив'язки python GDAL / OGR. Який мінімум мені потрібно встановити для використання прив'язок пітонів GDAL та OGR? Який найпростіший спосіб встановити GDAL / …

18 python  gdal  ogr 

Я намагаюся змінити значення поля файлу форми. Однак це значення не змінює. Що не так з моїм кодом? import ogr driver = ogr.GetDriverByName('ESRI Shapefile') fn = 'dist.shp' dataSource = driver.Open(fn, 0) layer = dataSource.GetLayer() feature = layer.GetNextFeature() dist = 233 while feature: feature.SetField("dist", dist) layer.SetFeature(feature) feature = layer.GetNextFeature() dataSource.Destroy()

18 gdal  ogr  fields-attributes 

У мене формуляр з багатокутниками. І у мене є глобальний растровий файл. Я хочу перекласти багатокутники формифільму на растрову сітку та обчислити середнє значення растру для кожного багатокутника. Як я можу це зробити за допомогою GDAL, записуючи результати у файл форми?

18 python  raster  shapefile  gdal  spatial-statistics 

Наразі я використовую такий сценарій, щоб додати деякі дані атрибутів із таблиці до безлічі індивідуальних файлів: library(rgdal) specieslist <- read.csv("SpeciesList1.txt", header=F) attdata <- read.table("TestAtt.csv", sep = ",", header=T) for (n in 1:dim(specieslist)[1]) { speciesname <- specieslist[n,1] shp <- readOGR("Mesoamerica_modified_polygons", speciesname) shp$ENGL_NAME<-attdata[n,2] writeOGR(shp, "PolygonsV2", speciesname, driver="ESRI Shapefile") } В кінці я …

18 shapefile  gdal  r  ogr  attribute-table 

Я досить часто в коді Python GDAL бачу, що люди закривають набори даних наприкінці сценарію. Чому має сенс закрити набір даних у Python GDAL? Чи є наслідки, якщо я цього не роблю? import gdal # open dataset ds = gdal.Open('test.tif') # close dataset ds = None

17 python  gdal 

Я перепроектую растри в python за допомогою GDAL. Мені потрібно спроектувати декілька тифів від географічних координат WGS 84 до Web Mercator WGS 1984 (допоміжна сфера), щоб згодом використовувати їх у Openlayers разом із OpenStreetMap та, можливо, картами Google. Я використовую Python 2.7.5 та GDAL 1.10.1 звідси і перетворюю координати, використовуючи …

17 python  gdal  coordinate-system  epsg  web-mercator 

Я готую дані для включення до вузла геомережі, де карти відображаються через геосервер. Проблема полягає в тому, що клітини, у яких "немає даних", відображаються як геотифи як чорні, і я хотів би, щоб колір був встановлений на білий або прозорий. Я використовую ArcView і не маю доступу ні до просторового …

17 arcgis-desktop  geoserver  gdal  geotiff-tiff  geonetwork 

Я намагаюся обробити деякі біокліматичні растрові файли, такі як можна завантажити з http://www.worldclim.org/current (набір bioclim). Вони, схоже, мають значення -3.4e+38нодату, встановлені відповідно до QGIS (дивлячись на вихід gdalinfo, це -3.39999999999999996e+38). Здається, що інструменти gdal не в змозі впоратися з цим значенням nodata, а qgis, схоже, не може розпізнати його. У …

17 qgis  raster  gdal 

Хтось знає, який найпростіший спосіб встановити GDAL 2.0 на MacOSX? Звичайні підозрювані бінарні розповсюдження Kingchaos та домашній мов, схоже, застрягли в GDAL 1.11 Мені потрібна підтримка властивостей функції TopoJSON, які, здається, були включені в 2.0 EDIT 2017/01/07: двійковий розподіл kingchaos тепер забезпечує розподіл GDAL 2.1

16 gdal  installation  ogr  mac  homebrew 

Зараз я працюю над імпортом кліматичних даних CANGRID (надаються як Surfer Grid ascii, ".grd" файли) в ArcGIS. Сітка має 95 рядів розміром 125 стовпців. Метадані містять ширину / тривалість походження (лівий нижній кут), розмір комірки (50 км), а також відзначає проекцію як полярну стереографічну з центральним меридіаном (110 градусів Шт) …

16 arcgis-desktop  python  gdal  coordinate-system