landsat

3热度

2回答

我遇到了一个与R中的覆盖函数不同寻常的问题。我们试图用另一个图层中的值填充阴影像素。我可以让它工作得很好用栈如下 - library(raster) r1 <- raster(ncols=36, nrows=18) r1[] <- 1:ncell(r1) r1b <- r1a <- r1 r1_stack <- stack(r1, r1a, r1b) r2 <- setValues(r

0热度

2回答

当使用perl和earthexplorer时格式错误的JSON字符串

我目前正在与USGS的EarthExplorer一起根据空间坐标设置一些Landsat场景的批量下载。他们提供了一个非常好用的脚本，在https://earthexplorer.usgs.gov/inventory/example/json-download_data-pl这很好。我工作的集群上，尽管在正确安装的所有Perl模块，当我运行该脚本，我得到以下的输出： Running Script..

2热度

1回答

从AWS加载单波段Landsat 8栅格并将它们组合为单个多波段RDD的最简单方法是什么？

我使用geotrellis从Landsat 8中加载坐在S3上的geotiff栅格。但是，它们存储在每个频段的基础上。我可以使用S3GeoTiff类加载单独的频段，例如： val options = S3GeoTiffRDD.Options(getS3Client =() => S3Client.ANONYMOUS) val rddBlue = S3GeoTiffRDD.spatial("la

0热度

1回答

RGB复合Landsat-8图像

我试图从Landsat-8图像的各个波段创建RGB复合图像。尝试选项，如np.dstack，也是下面的格式 img[:,:,0] = b4 img[:,:,1] = b3 img[:,:,2] = b2 但是要得到一个错误 Error Creating RGB band

3热度

2回答

的Python：脚本导入多个陆地卫星图像，以草GIS

我想写一个Python代码，将通过调整这个代码导入LANDSAT卫星图像转换成草GIS：http://grass.osgeo.org/wiki/LANDSAT LANDSAT瓦片被下载的文件夹，每个包含7个tiff图像（波段1-7）。因此，我有一个包含几个子目录的目录（每个LANDSAT磁贴一个）。我目前的代码如下： #!/usr/bin/python import os import s

0热度

1回答

乐队如何读取数组只是返回我0，我该怎么办？

我是GDAL的新手。我使用的分类TM遥感影像，但我将这些代码行带数据只是给我0，有什么我不明白： raster_dataset = gdal.Open(‘LC81850552017065LGN00_B3.tif’, gdal.GA_ReadOnly) geo_transform = raster_dataset.GetGeoTransform() proj = raster_dataset.G

0热度

2回答

在栅格数据上实现回归树的方法 - python

我想在python中构建和实现一些栅格数据的回归树算法，并且似乎无法找到这样做的最佳方法。我将尝试解释我想要做的事情：我想要的输出是一个栅格图像，其值代表湖深，称之为depth.tif。我有一系列光栅图像，每一个都代表不同Landsat波段的反射率值，比如[B1.tif，B2.tif，...，B7.tif]，我想用它作为自变量来预测湖泊深度。对于我的训练数据，我有一个约6000点已知湖深的sh

1热度

1回答

获取值等于0的像素的坐标

我有一个202行和201列的栅格地图这个地图中有一些栅格，其像素值为0我想写一个函数，返回所有的像素值0格坐标我怎么能做到这一点我尝试使用，如果循环和while循环但它总是说TRUE/FALSE需要这里是我的示例代码 library(raster) library(rgdal) library(maptools) library(sp) setwd("E:/Landsat-d

0热度

1回答

在landsat上结合RGB通道并转换为png给出了一个可怕的结果

我正在使用Scala和Geotrellis库来执行Landsat-8图像的处理和分析，但是当我尝试将band 2,3和4（Red ，绿和蓝），并在这里是所得到的多频带TIFF（右）和PNG从它产生（左）： http://imgur.com/a/aA4XW 用于组合单频带的TIFF的代码是： class Tiler (outputPath : String, bandIds : Array[Stri

1热度

2回答

使用Python GDAL模块打开在同一个投影图像，但大小不同

我不知道这是正确的网站这一点，但，我使用gdal.Open（）模块Python 2.7版打开Landsat GeoTIFF图像。它们具有相同的UTM地图投影，但图像大小不同。如何修复图像大小相同？我想直接进行逐像素比较。