satellite-image

Problem importing USDA Crop Data Layer (GTIFF) in R using raster() package

由 岁酱吖の 提交于 2019-12-23 15:37:19
问题 I've downloaded the USDA Crop Data Layer (data available here, note large file: http://www.nass.usda.gov/research/Cropland/SARS1a.htm) and want to import it into R for analysis. I've imported other large GTIFF files (e.g., LANDSAT data), but get the following error when I call the function: r <- raster(file) Error in .rasterObjectFromFile(x, band = band, objecttype = "RasterLayer", : Cannot create a RasterLayer object from this file. where "file" is any of the .tif files provided by USDA (I'm

Big Satellite Image Processing

由 我只是一个虾纸丫 提交于 2019-12-19 10:11:35
问题 Im trying to run Mort Canty's http://mcanty.homepage.t-online.de/ Python iMAD implementation on bitemporal RapidEye Multispectral images. Which basically calculates the canonical correlation for the two images and then substracts them. The problem I'm having is that the images are of 5000 x 5000 x 5 (bands) pixels. If I try to run this on the whole image I get a memory error. Would the use of something like pyTables help me with this? What Mort Canty's code tries to do is that it loads the

Assign different color to NaN values in Matlab images

由 谁都会走 提交于 2019-12-13 05:23:59
问题 I am working with satellite images on Matlab and am writing a code for calculating Net Radiation. The output contains the image in form of a matrix with NaN values. When I try to plot it using: figure imshow(X); it assigns the NaN values the same color as 0. Is there any way I can assign those values to be a different color? Also, I want to save my matrix as an image, but in a format that its values don't get changed ... possibly ASCII. Is there a tool for converting matrix to ASCII? 回答1: If

Merging multiple rasters in R

由 試著忘記壹切 提交于 2019-12-04 09:21:11
问题 I've been trying to find a time-efficient way to merge multiple raster images in R. These are adjacent ASTER scenes from the southern Kilimanjaro region, and my target is to put them together to obtain one large image. This is what I got so far (object 'ast14dmo' representing a list of RasterLayer objects): # Loop through single ASTER scenes for (i in seq(ast14dmo.sd)) { if (i == 1) { # Merge current with subsequent scene ast14dmo.sd.mrg <- merge(ast14dmo.sd[[i]], ast14dmo.sd[[i+1]],

Merging multiple rasters in R

由 隐身守侯 提交于 2019-12-03 02:01:19
I've been trying to find a time-efficient way to merge multiple raster images in R. These are adjacent ASTER scenes from the southern Kilimanjaro region, and my target is to put them together to obtain one large image. This is what I got so far (object 'ast14dmo' representing a list of RasterLayer objects): # Loop through single ASTER scenes for (i in seq(ast14dmo.sd)) { if (i == 1) { # Merge current with subsequent scene ast14dmo.sd.mrg <- merge(ast14dmo.sd[[i]], ast14dmo.sd[[i+1]], tolerance = 1) } else if (i > 1 && i < length(ast14dmo.sd)) { tmp.mrg <- merge(ast14dmo.sd[[i]], ast14dmo.sd[[i

Big Satellite Image Processing

由 耗尽温柔 提交于 2019-12-01 11:06:08
Im trying to run Mort Canty's http://mcanty.homepage.t-online.de/ Python iMAD implementation on bitemporal RapidEye Multispectral images. Which basically calculates the canonical correlation for the two images and then substracts them. The problem I'm having is that the images are of 5000 x 5000 x 5 (bands) pixels. If I try to run this on the whole image I get a memory error. Would the use of something like pyTables help me with this? What Mort Canty's code tries to do is that it loads the images using gdal and then stores them in an 10 x 25,000,000 array. # initial weights wt = ones(cols*rows

Using pcolormesh for plotting an orbit data

由 痴心易碎 提交于 2019-11-29 05:14:37
I am trying to map a dataset with associated latitude and longitude. The details of the data I am using are given below: Variable Type Data/Info ------------------------------- lat ndarray 1826x960, type `float64` lon ndarray 1826x960, type `float64` data ndarray 1826x960, type `float64` I have created then a basemap: m = Basemap(projection='cyl', llcrnrlon=-180, urcrnrlon=180, llcrnrlat=-40, urcrnrlat=40, resolution='c') Now, on the basemap created, I'd plot the above mentioned dataset using pcolormesh: m.drawcoastlines() m.drawcountries x,y = m(lon,lat) m.pcolormesh(x,y,data) m.colorbar()

Using pcolormesh for plotting an orbit data

由 大憨熊 提交于 2019-11-27 18:57:32
问题 I am trying to map a dataset with associated latitude and longitude. The details of the data I am using are given below: Variable Type Data/Info ------------------------------- lat ndarray 1826x960, type `float64` lon ndarray 1826x960, type `float64` data ndarray 1826x960, type `float64` I have created then a basemap: m = Basemap(projection='cyl', llcrnrlon=-180, urcrnrlon=180, llcrnrlat=-40, urcrnrlat=40, resolution='c') Now, on the basemap created, I'd plot the above mentioned dataset using
订阅 satellite-image