«gdalwarp» 태그된 질문

경도가 0에서 360 사이 인 좌표계를 가진 지오메트리 래스터 이미지가 있습니다. 이미지의 가로 가운데는 180 도입니다. 아래 이미지를보십시오 : -180180 경도 범위의 EPSG : 4326 SRS로 변환하고 싶습니다. 그리고 이미지의 중심이 그리니치 자오선 (0)에 있기를 원합니다. 나는이 srs가 매우 널리 사용되는 것 같아요. 결과는 다음과 같습니다. 그래서 gdalwarp 명령을 사용하여 …

31 coordinate-system  gdal  geotiff-tiff  proj  gdalwarp 

shapefile 형상에 따라 Raster 파일을 자르려고합니다. 다음 코드를 사용하고 있습니다 gdalwarp -cutline INPUT.shp INPUT.tif OUTPUT.tif 그러나 쉐이프 파일 지오메트리 범위를 벗어난 검은 색이됩니다. 나는 여기에 예를주고있다. 클립하고 싶은 첫 번째 이미지에서. 두 번째 이미지는 결과 래스터이지만 지오메트리의 외부가 null 인 세 번째 래스터를 원합니다. 첫 번째 이미지 : 입력 래스터 …

26 gdal  gdalwarp 

GDAL python 바인딩을 사용하여 다시 프로젝트 / 다시 샘플링하려고하지만 명령 줄 유틸리티의 결과와 약간 다른 결과가 나타납니다 gdalwarp. 더 짧은 예는 아래 업데이트를 참조하십시오. 이 스크립트는 Python 접근 방식을 보여줍니다. from osgeo import osr, gdal import numpy def reproject_point(point, srs, target_srs): ''' Reproject a pair of coordinates from one spatial …

20 python  gdal  gdalwarp  numpy 

여러 래스터 gdalwarp를 통해 그리드에 정렬하고 정렬 한 후 -tap출력 래스터가 원래 래스터보다 상당히 큰 것을 알았습니다. 상당히 철저한 웹 검색 으로이 Trac 문제 가 나타났습니다 . 프랭크 워머 담은 그 이유를 설명했다. "주의 깊게 검토 할 때 gdal_translate는 TIFFWriteScanline () 인터페이스를 사용하여 GTiffDataset :: CreateCopy? () 내에서 출력 파일을 …

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

나는 약 100MB마다 약 20 개의 dem tiff를 병합하고 싶었고 gdalwarp (가상 파일 옵션에 대해 알고 있음)를 사용했습니다. gdalwarp $(list_of_tiffs) merged.tiff 그러나이 명령은 매우 오랜 시간이 걸렸으며 15 분 후에 세트에서 5 번째 이미지를 처리하고있었습니다. 종료합니다. 그런 다음 gdal_merge.py 스크립트를 모자이크 처리에도 사용할 수 있다는 것을 알았습니다. gdal_merge.py $(list_of_tiffs) 3 …

14 gdal  gdalwarp  gdal-merge 

다중 대역 지오메트리 래스터 파일을 동일한 (더 작은) 영역으로 자르는 간단한 유틸리티를 작성 중입니다. gdalwarp를 사용하여 단일 다각형 클리핑 shapefile을 사용하여 파일을 쉽게자를 수 있습니다. gdalwarp -cutline clipper.shp -crop_to_cutline input.tif output.tif 그러나 클립하려는 실제 영역은 항상 shapefile이 아닌 다른 geotiff 래스터 파일로 정의됩니다. 래스터의 범위를 클리핑 파일로 사용할 수 있다면 …

14 python  raster  gdal  clip  gdalwarp 

QGIS 2.XX에서 드롭 다운 메뉴 표시 줄 도구는 하단에 편집 가능한 전화 상자가 있으며 GUI에서 사용할 수없는 명령은 연필 버튼을 클릭하여 수동으로 입력 할 수 있습니다. 3.0에서 메뉴 표시 줄 도구는 처리 도구 상자에있는 GDAL 도구와 동일하지만 더 의미가 있지만 도구 작업에 간단한 명령을 추가 할 수있는 방법은 없습니다. 내가 …

13 gdalwarp  qgis-3 

QGis Save Raster 의 코드를 렌더링 된 이미지로 사용하려고했습니다 . i래스터 레이어입니다. pipelayer = i pipeextent = pipelayer.extent() pipewidth, pipeheight = (pipelayer.width(), pipelayer.height()) piperenderer = pipelayer.renderer() pipeprovider = pipelayer.dataProvider() crs = pipelayer.crs().toWkt() pipe = QgsRasterPipe() pipe.set(pipeprovider.clone()) pipe.set(piperenderer.clone()) pipedFile = os.path.join(tempfile.gettempdir(), safeLayerName + '_pipe.tif') print pipedFile file_writer = QgsRasterFileWriter(pipedFile) file_writer.writeRaster(pipe, pipewidth, …

12 qgis  gdal  pyqgis  gdalwarp  gdal-translate 

나는 gdal로 일부 이미지를 모자이크 처리하고 모자이크 중간의 날카로운 가장자리를 제거하기 위해 각 이미지의 가장자리쪽으로 페이딩 / 점진적 알파 레이어를 사용하여 최종 결과를 개선하고 싶습니다. 내가 겪고있는 문제는 점진적 알파 레이어가있는 각 개별 이미지의 일부가 아래 표시된 것처럼 반투명하지 않고 최종 모자이크에서 이미지 아래에 이미지를 숨기고 있다는 것입니다. 이상적 으로이 …

12 gdal  geotiff-tiff  gdalwarp  gdal-translate  gdalbuildvrt 

gdalwarp를 사용하여 날짜 표시 줄 근처에서 SRTM 타일을 조작하고 있습니다 (예 : 반 자오선). SRTM 타일은 자오선과 매우 약간 (1/2 픽셀) 겹칩니다. gdalinfo를 사용하여 이것을 볼 수 있습니다 : gdalinfo S16W180.hgt Driver: SRTMHGT/SRTMHGT File Format Files: S16W180.hgt Size is 1201, 1201 [...] Lower Left (-180.0004167, -16.0004167) (180d 0' 1.50"W, 16d …

11 srtm  gdalwarp 

gdalwarp를 사용하여 Geostationary에서 Lambert Conformal로 다시 매핑하려고합니다. 내 입력 데이터는 netcdf에 있으며 지리적 좌표 (도)이며 다시 매핑 된 데이터를 netcdf에 출력하고 싶습니다. 입력 netcdf 데이터에 해당하는 vrt 파일을 작성했습니다. Gdalwarp는 netcdf 파일을 출력하지만 출력 데이터는 모두 0이며 다음과 같은 오류가 발생합니다. Creating output file that is 5120P x 5120L. Processing …

11 gdal  coordinate-system  proj  gdalwarp  netcdf 

pythonand를 사용하여 래스터를 지리 참조하고 싶습니다 GDAL. 저의 현재 접근 방식은 전화 gdal_translate및 gdalwarp사용 os.system및 추악한 지상 통제 지점 목록입니다. 나는 이것을 기본적으로 할 수있는 방법을 정말로 원합니다 python. 이것은 현재 사용중인 프로세스입니다. import os os.system('gdal_translate -of GTiff -gcp 1251.92 414.538 -7.9164e+06 5.21094e+06 -gcp 865.827 107.699 -7.91651e+06 5.21104e+06 "inraster.tif" "outraster1.tif"') …

10 python  raster  gdal  coordinate-system  gdalwarp 

우리가 가지고있는 고해상도 이미지의 타일을 생성하기 위해 약간의 노력을 기울였습니다. 현재 가지고있는 이미지는 GeoTiff 파일로 저장된 매우 큰 (+ 20GB) 이미지입니다. gdal2tiles 명령 행 유틸리티를 사용하여 타일을 생성 한 다음 TMS 이미지 제공자를 사용하여 타일을 제공하기 위해 Cesium에서 열어서보고 싶습니다. gdalinfo를 사용하여 이미지에 대한 세부 정보는 다음과 같습니다. Driver: GTiff/GeoTIFF …

10 gdal  gdalwarp  cesium  gdal2tiles 

래스터를 투영하려고합니다. R에는 이것에 대한 projectRaster()기능이 있습니다 (완전히 재현 가능한 예 아래). # example Raster require(raster) r <- raster(xmn=-110, xmx=-90, ymn=40, ymx=60, ncols=40, nrows=40) r <- setValues(r, 1:ncell(r)) projection(r) # project to newproj <- "+init=epsg:4714" # using raster package to reproject pr1 <- projectRaster(r, crs = CRS(newproj), method = 'bilinear') …

9 raster  gdal  r  coordinate-system  gdalwarp 

파이썬과 QGIS 2.0을 사용하고 있습니다. 하나의 다각형 기능으로 폴더의 래스터를 자르려고합니다. "PyQGIS"를 사용하는 것은 처음입니다. 어쨌든, 나는 간단한 스크립트가 작동하지 않으면 어떤 제안이라도 대단히 감사하겠습니다! import qgis.core, qgis,utils QgsApplication.setPrefixPath("C:/OSGeo4W64/apps/qgis", True) QgsApplication.initQgis() CLIP= "C:/Users/unim/Documents/Umberto/Universita/PhD/Guglielmin/Permafrost/Alta_Valtellina/Landsat_ita/study_area_foscagno.shp" INPUT_FOLDER="C:/Users/unimi/Documents/Umberto/Universita/PhD/Guglielmin/Permafrost/Alta_Valtellina/Landsat_ita/LE71930282000259EDC00" OUTPUT= "C:/Users/unim/Documents/Umberto/Universita/PhD/Guglielmin/Permafrost/Alta_Valtellina/Landsat_ita/foscagno_pyqgis/" for RASTER in INPUT_FOLDER.tif do echo "Processing $RASTER" gdalwarp -q -cutline CLIP -crop_to_cutline -of …

9 gdal  pyqgis  gdalwarp