7 votes

Remove clumps of pixels in R

Thanks to @gene and https://geoscripting-wur.github.io/AdvancedRasterAnalysis/ I can now answer my question (copied and modified): library(raster) # create some raster data r <- raster(ncols=12, ...
user avatar
  • 3,276
6 votes

Confidence in pan-sharpened classification?

In general, there are two approaches to classification: pixel-based and object-based: Pixel-based: Each spatial pixel is evaluated by itself against a set classification parameters. In this case, ...
user avatar
5 votes
Accepted

Erdas python scripting help

A good documentation can be found here (.pdf). Also look at the Erdas forum here. Another good resource is ERDAS IMAGINE 2014 Python Examples
user avatar
5 votes

Matching colours in already mosaiced single image?

The image needs to be cut apart and re-mosaiced after adjusting the colors in each piece. This can be done. As an example, I extracted the green band of the image. To make my work simple (the ...
user avatar
  • 67.8k
5 votes
Accepted

How to easily identify/classify coniferous forest from one date Landsat image?

NDVI is for vegetation/non vegetation discrimination. So if your vegetation is always coniferous forest, then it should be the most efficient method in your case. Otherwise you will have confusions ...
user avatar
  • 48.7k
5 votes
Accepted

Remove clumps of pixels in R

#reproducible example r <- raster(ncols=12, nrows=12) set.seed(0) r[] <- round(runif(ncell(r))*0.7 ) rc <- clump(r) #extract IDs of clumps according to some criteria clump9 = data.frame(...
user avatar
  • 9,921
4 votes

Matching colours in already mosaiced single image?

Histogram matching works by forcing the histogram of one image to match as closely as possible the histogram of a second target image. I'm afraid that it won't work on a single image. (There is the ...
user avatar
  • 10.8k
4 votes

1m cell size equivalent value in decimal degree

Assuming a spherical Earth of radius, R = 6,371,000 m, and a latitudinal (N-S) arc length, a = 1m, then the latitudinal arc angle, dφ = a / R = 0.000000157 rad = 0.0000089932 deg, but note what ...
user avatar
  • 8,604
4 votes

Using Erdas for land cover land use change detection

ATCOR for Erdas Imagine will convert DN to true reflectance--this step is critical for change detection analyses. Then you can use DeltaCue add-on in Erdas Imagine to detect land cover change. More ...
user avatar
  • 50.4k
4 votes
Accepted

Is possible to merge landsat 8 with High Resolution Satellite Imagery?

It is technically possible to use the pansharpening algorithm with different sensors, and all your tagged software have pansharpening tools (sometime . However, the quality of the outputs will depend: ...
user avatar
  • 48.7k
4 votes
Accepted

Existence of fiducials marks in a digital camera file

fiducial mark are used to define the coordinate system of the photograph. With film photograph, the paper moves under the objective and can get further distorted during storage and development, so you ...
user avatar
  • 48.7k
3 votes

Floating point to 8-bit unsigned raster conversion using ArcGIS for Desktop or ERDAS IMAGINE?

You will want to stretch your floating point data from (-1 to 1) to (0 - 255). You can accomplish this very easily in Erdas Imagine: Raster tab > Resolution Group > Radiometric > Rescale
user avatar
  • 50.4k
3 votes

How to rescale data ranges in ERDAS?

Rescaling an NDVI or EVI from -1 to 1, to 0 to 1, uses the Rescale function (under Raster, Radiometric). Clipping the top and bottom 0.5% is a percentage linear contrast stretch. To do this in ERDAS ...
user avatar
  • 8,926
3 votes
Accepted

Stereo Analyst & ArcGIS 10.1 - volume calculation

To calculate a volume, you need the 3D Analyst extension. 3D Analyst has several ways to calculate volumes. All methods require the data to be in a projected coordinate system. You can use TIN or ...
user avatar
  • 1,351
3 votes

Shifting cloud mask in Erdas/ArcMap/Py Script

There is a simple way you could do this using a python script, however it does not truly create a cloud shadow. For more information on creating a cloud shadow see this paper by Zhu and Woodcock (2014)...
user avatar
3 votes

Shapefile Conversion to ERDAS Signature editor

You can use the Signature Editor in ERDAS to use Shapefiles to create signature files. The Signature Editor is found under the supervised classification in the raster ribbon of ERDAS. You must have ...
user avatar
  • 708
3 votes
Accepted

Deriving slope angle value from DEM using ERDAS Imagine

Have you though of just using the free GDAL tool? http://www.gdal.org/gdaldem.html It can run in windows using OSGEO4W and can ouput .img files. For Imagine 8 it goes something like... From ...
user avatar
3 votes

Pro vs Con on using export to JPEG or JP2000 files for large imagery files

JPEG2000 Pros: Better compression than JPEG Supports both lossy and lossless compression, JPEG is lossy only Supports any number of bands, JPEG only supports 3 bands Supports more datatypes (...
user avatar
  • 53.4k
3 votes
Accepted

How to remove band from stack in ERDAS 2013?

Use the Layer Stack tool in reverse by using your 11 band stack as the input and selecting only the bands you want in the output. For example, in the attached screenshot I added bands 1, 2, 3 to the ...
user avatar
  • 50.4k
3 votes

What exactly does hill shade or Shaded relief , if performed on SRTM DEM will output?

Hillshade computes the local illumination from a light source located at infinity (like the sun). Basically, it yields the cosinus of the normal to the face of the terrain and the light ray. This can ...
user avatar
  • 48.7k
3 votes

Calculating NDWI in Erdas IMAGINE 2014

the original (GAO 1996) NDWI is based reflectance for the wavelength of 0.86 (NIR) and 1.24 µm (SWIR). This corresponds to bands 2 and 5 if you use MODIS. The equation for MODIS is thus (b2-b5)/(b2+b5)...
user avatar
  • 48.7k
3 votes

what is the basic programming language or code language for raster swipe tool?

In OpenLayers Swipe Control is easy enough to use, the following code should set you up: var map; var swipe; function init() { map = new OpenLayers.Map({ div : "map", ...
user avatar
3 votes

Knowing the pixel size for a raster before it was resampled

Upsample the image to the suspected resolution and Downsample it to the original. If the resulting and original are the same (or close enough), it is very likely that you have a lower resolution image....
user avatar
3 votes
Accepted

Smooth mosaicking of smaller orthomosaics to form a homogenous mosaic

That is noise / offset pixels or unexpected dark or light strips fringing an object. Data from airborne multispectral or hyperspectral imaging scanners shows a pronounced striping pattern due to ...
user avatar
  • 1,675
2 votes
Accepted

Stacking overlapping rasters, keeping cell values on top?

There are several methods to combine rasters in ArcGis; generally ERDAS is faster at raster functions than ArcGis but I'm so out-of-touch with ERDAS that I can't even remember how to mosaic. The ...
user avatar
2 votes
Accepted

Determining which "color" each image band represents?

As others have mentioned, the best practical way to determine which bands are which is to look at the source metadata, which is widely available for common products like ASTER and Rapideye. You can ...
user avatar
  • 50.4k
2 votes

ERDAS Imagine ecw compression settings

ECW V3 version (since the 2015 version at least) preserves the 16-bitness too, you used V2 and it is more compatible but cannot preserve it. Raster --> Radiometric --> LUT Stretch is how to do that 16 ...
user avatar
2 votes
Accepted

How to create an Erdas batch list file (.bls) using Python?

.bls files are relatively easy to create. It's the .bcf files that are very difficult to create in python because ERDAS has horrible batch creation and debugging. I've done both in Python before, and ...
user avatar
  • 11.5k
2 votes
Accepted

Creating training areas for parametric classifcation

One of the most common reason for the covariance matrix that is not invertible from samples is the presence of null variance for one band (all your pixel in one band have the same value). In this case,...
user avatar
  • 48.7k
2 votes

MODIS MOD13Q1 ndvi calculation

In ENVI go to "Band math" and tape: float(b1)/0.0001 look at http://www.exelisvis.fr/docs/BandMath.html , Actually there is no need for this, the reason why values ​​are multiplied by 0.0001 is only ...
user avatar
  • 21

Only top scored, non community-wiki answers of a minimum length are eligible